The 8-17 DNAzyme can operate in a single active structure regardless of metal ion cofactor.

DNAzymes - synthetic enzymes made of DNA - have long attracted attention as RNA-targeting therapeutic agents. Yet, as of now, no DNAzyme-based drug has been approved, partially due to our lacking understanding of their molecular mode of action. In this work we report the solution structure of 8-17 DNAzyme bound to a Zn2+ ion solved through NMR spectroscopy. Surprisingly, it turned out to be very similar to the previously solved Pb2+-bound form (catalytic domain RMSD = 1.28 Å), despite a long-standing literature consensus that Pb2+ recruits a different DNAzyme fold than other metal ion cofactors. Our follow-up NMR investigations in the presence of other ions - Mg2+, Na+, and Pb2+ - suggest that at DNAzyme concentrations used in NMR all these ions induce a similar tertiary fold. Based on these findings, we propose a model for 8-17 DNAzyme interactions with metal ions postulating the existence of only a single catalytically-active structure, yet populated to a different extent depending on the metal ion cofactor. Our results provide structural information on the 8-17 DNAzyme in presence of non-Pb2+ cofactors, including the biologically relevant Mg2+ ion.

Comparison of the Backbone Dynamics of Dehaloperoxidase-Hemoglobin Isoenzymes.

Dehaloperoxidase (DHP) is a multifunctional hemeprotein with a functional switch generally regulated by the chemical class of the substrate. Its two isoforms, DHP-A and DHP-B, differ by only five amino acids and have an almost identical protein fold. However, the catalytic efficiency of DHP-B for oxidation by a peroxidase mechanism ranges from 2- to 6-fold greater than that of DHP-A depending on the conditions. X-ray crystallography has shown that many substrates and ligands have nearly identical binding in the two isoenzymes, suggesting that the difference in catalytic efficiency could be due to differences in the conformational dynamics. We compared the backbone dynamics of the DHP isoenzymes at pH 7 through heteronuclear relaxation dynamics at 11.75, 16.45, and 19.97 T in combination with four 300 ns MD simulations. While the overall dynamics of the isoenzymes are similar, there are specific local differences in functional regions of each protein. In DHP-A, Phe35 undergoes a slow chemical exchange between two conformational states likely coupled to a swinging motion of Tyr34. Moreover, Asn37 undergoes fast chemical exchange in DHP-A. Given that Phe35 and Asn37 are adjacent to Tyr34 and Tyr38, it is possible that their dynamics modulate the formation and migration of the active tyrosyl radicals in DHP-A at pH 7. Another significant difference is that both distal and proximal histidines have a 15-18% smaller S2 value in DHP-B, thus their greater flexibility could account for the higher catalytic activity. The distal histidine grants substrate access to the distal pocket. The greater flexibility of the proximal histidine could also accelerate H2O2 activation at the heme Fe by increased coupling of an amino acid charge relay to stabilize the ferryl Fe(IV) oxidation state in a Poulos-Kraut "push-pull"-type peroxidase mechanism.

In Vitro Enzyme Kinetics and NMR-Based Product Elucidation for Glutathione S-Conjugation of the Anticancer Unsymmetrical Bisacridine C-2028 in Liver Microsomes and Cytosol: Major Role of Glutathione S-Transferase M1-1 Isoenzyme.

This work is the next step in studying the interplay between C-2028 (anticancer-active unsymmetrical bisacridine developed in our group) and the glutathione S-transferase/glutathione (GST/GSH) system. Here, we analyzed the concentration- and pH-dependent GSH conjugation of C-2028 in rat liver microsomes and cytosol. We also applied three recombinant human GST isoenzymes, which altered expression was found in various tumors. The formation of GSH S-conjugate of C-2028 in liver subfractions followed Michaelis-Menten kinetics. We found that C-2028 was conjugated with GSH preferentially by GSTM1-1, revealing a sigmoidal kinetic model. Using a colorimetric assay (MTT test), we initially assessed the cellular GST/GSH-dependent biotransformation of C-2028 in relation to cytotoxicity against Du-145 human prostate cancer cells in the presence or absence of the modulator of GSH biosynthesis. Pretreatment of cells with buthionine sulfoximine resulted in a cytotoxicity decrease, suggesting a possible GSH-mediated bioactivation process. Altogether, our results confirmed the importance of GSH conjugation in C-2028 metabolism, which humans must consider when planning a treatment strategy. Finally, nuclear magnetic resonance spectroscopy elucidated the structure of the GSH-derived product of C-2028. Hence, synthesizing the compound standard necessary for further advanced biological and bioanalytical investigations will be achievable.

The interactions of monomeric acridines and unsymmetrical bisacridines (UAs) with DNA duplexes: an insight provided by NMR and MD studies.

Members of a novel class of anticancer compounds, exhibiting high antitumor activity, i.e. the unsymmetrical bisacridines (UAs), consist of two heteroaromatic ring systems. One of the ring systems is an imidazoacridinone moiety, with the skeleton identical to the structural base of Symadex. The second one is a 1-nitroacridine moiety, hence it may be regarded as Nitracrine's structural basis. These monoacridine units are connected by an aminoalkyl linker, which vary in structure. In theory, these unsymmetrical dimers should act as double-stranded DNA (dsDNA) bis-intercalators, since the monomeric units constituting the UAs were previously reported to exhibit an intercalating mode of binding into dsDNA. On the contrary, our earlier, preliminary studies have suggested that specific and/or structurally well-defined binding of UAs into DNA duplexes might not be the case. In this contribution, we have revisited and carefully examined the dsDNA-binding properties of monoacridines C-1305, C-1311 (Symadex), C-283 (Ledakrin/Nitracrine) and C-1748, as well as bisacridines C-2028, C-2041, C-2045 and C-2053 using advanced NMR techniques, aided by molecular modelling calculations and the analysis of UV-VIS spectra, decomposed by chemometric techniques. These studies allowed us to explain, why the properties of UAs are not a simple sum of the features exhibited by the acridine monomers.

Solution Structure of a Lanthanide-binding DNA Aptamer Determined Using High Quality pseudocontact shift restraints.

In this contribution we report the high-resolution NMR structure of a recently identified lanthanide-binding aptamer (LnA). We demonstrate that the rigid lanthanide binding by LnA allows for the measurement of anisotropic paramagnetic NMR restraints which to date remain largely inaccessible for nucleic acids. One type of such restraints - pseudocontact shifts (PCS) induced by four different paramagnetic lanthanides - was extensively used throughout the current structure determination study and the measured PCS turned out to be exceptionally well reproduced by the final aptamer structure. This finding opens the perspective for a broader application of paramagnetic effects in NMR studies of nucleic acids through the transplantation of the binding site found in LnA into other DNA/RNA systems.

Acid-Base Equilibrium and Self-Association in Relation to High Antitumor Activity of Selected Unsymmetrical Bisacridines Established by Extensive Chemometric Analysis.

Unsymmetrical bisacridines (UAs) represent a novel class of anticancer agents previously synthesized by our group. Our recent studies have demonstrated their high antitumor potential against multiple cancer cell lines and human tumor xenografts in nude mice. At the cellular level, these compounds affected 3D cancer spheroid growth and their cellular uptake was selectively modulated by quantum dots. UAs were shown to undergo metabolic transformations in vitro and in tumor cells. However, the physicochemical properties of UAs, which could possibly affect their interactions with molecular targets, remain unknown. Therefore, we selected four highly active UAs for the assessment of physicochemical parameters under various pH conditions. We determined the compounds' pKa dissociation constants as well as their potential to self-associate. Both parameters were determined by detailed and complex chemometric analysis of UV-Vis spectra supported by nuclear magnetic resonance (NMR) spectroscopy. The obtained results indicate that general molecular properties of UAs in aqueous media, including their protonation state, self-association ratio, and solubility, are strongly pH-dependent, particularly in the physiological pH range of 6 to 8. In conclusion, we describe the detailed physicochemical characteristics of UAs, which might contribute to their selectivity towards tumour cells as opposed to their effect on normal cells.

Iso-Partricin, an Aromatic Analogue of Amphotericin B: How Shining Light on Old Drugs Might Help Create New Ones.

Partricin is a heptaene macrolide antibiotic complex that exhibits exceptional antifungal activity, yet poor selective toxicity, in the pathogen/host system. It consists of two compounds, namely partricin A and B, and both of these molecules incorporate two cis-type bonds within their heptaenic chromophores: 28Z and 30Z. In this contribution, we have proven that partricins are susceptible to a chromophore-straightening photoisomerization process. The occurring 28Z→28E and 30Z→30E switches are irreversible in given conditions, and they are the only structural changes observed during the experiment. The obtained all-trans partricin's derivatives, namely iso-partricins A and B, exhibit very promising features, potentially resulting in the improvement of their selective toxicity.

Ipertrofan Revisited-The Proposal of the Complete Stereochemistry of Mepartricin A and B.

Being a methyl ester of partricin, the mepartricin complex is the active substance of a drug called Ipertrofan (Tricandil), which was proven to be useful in treatment of benign prostatic hyperplasia and chronic nonbacterial prostatitis/chronic pelvic pain syndrome. Nevertheless, no direct structural evidence on the stereochemistry of its components has been presented to date. In this contribution, we have conducted detailed, NMR-driven stereochemical studies on mepartricins A and B, aided by molecular dynamics simulations. The absolute configuration of all the stereogenic centers of mepartricin A and B was defined as 3R, 7R, 9R, 11S, 13S, 15R, 17S, 18R, 19S, 21R, 36S, 37R, and 38S, and proposed as 41R. The geometry of the heptaenic chromophore of both compounds has been established as 22E, 24E, 26E, 28Z, 30Z, 32E, and 34E. Our studies on mepartricin ultimately proved that partricins A and B are structurally identical to the previously described main components of the aureofacin complex: gedamycin and vacidin, respectively. The knowledge of the stereochemistry of this drug is a fundamental matter not only in terms of studies on its molecular mode of action, but also for potential derivatization, aiming at improvement of its pharmacological properties.

The origin of the high stability of 3'-terminal uridine tetrads: contributions of hydrogen bonding, stacking interactions, and steric factors evaluated using modified oligonucleotide analogs.

RNA G-quadruplexes fold almost exclusively into parallel-stranded structures and thus display much less structural diversity than their DNA counterparts. However, also among RNA G-quadruplexes peculiar structural elements can be found which are capable of reshaping the physico-chemical properties of the folded structure. A striking example is provided by a uridine tetrad (U-tetrad) placed on the 3'-terminus of the tetramolecular G-quadruplex. In this context, the U-tetrad adopts a unique conformation involving chain reversal and is responsible for a tremendous stabilization of the G-quadruplex (ΔTm up to 30°C). In this report, we attempt to rationalize the origin of this stabilizing effect by concurrent structural, thermal stability, and molecular dynamics studies of a series of G-quadruplexes with subtle chemical modifications at their 3'-termini. Our results provide detailed insights into the energetics of the "reversed" U-tetrad motif and the requirements for its formation. They point to the importance of the 2'OH to phosphate hydrogen bond and preferential stacking interactions for the formation propensity and stability of the motif.

A strong preference for the TA/TA dinucleotide step discovered for an acridine-based, potent antitumor dsDNA intercalator, C-1305: NMR-driven structural and sequence-specificity studies.

Triazoloacridinone C-1305, a potent antitumor agent recommended for Phase I clinical trials, exhibits high activity towards a wide range of experimental colon carcinomas, in many cases associated with complete tumor regression. C-1305 is a well-established dsDNA intercalator, yet no information on its mode of binding into DNA is available to date. Herein, we present the NMR-driven and MD-refined reconstruction of the 3D structures of the d(CGATATCG)2:C-1305 and d(CCCTAGGG)2:C-1305 non-covalent adducts. In both cases, the ligand intercalates at the TA/TA site, forming well-defined dsDNA:drug 1:1 mol/mol complexes. Orientation of the ligand within the binding site was unambiguously established by the DNA/ligand proton-proton NOE contacts. A subsequent, NMR-driven study of the sequence-specificity of C-1305 using a series of DNA duplexes, allowed us to confirm a strong preference towards TA/TA dinucleotide steps, followed by the TG/CA steps. Interestingly, no interaction at all was observed with duplexes containing exclusively the AT/AT, GG/CC and GA/TC steps.

Thermodynamic and structural contributions of the 6-thioguanosine residue to helical properties of RNA.

Thionucleotides, especially 4-thiouridine and 6-thioguanosine, are photosensitive molecules that photocrosslink to both proteins and nucleic acids, and this feature is a major reason for their application in various investigations. To get insight into the thermodynamic and structural contributions of 6-thioguanosine to the properties of RNA duplexes a systematic study was performed. In a series of RNA duplexes, selected guanosine residues located in G-C base pairs, mismatches (G-G, G-U, and G-A), or 5' and 3'-dangling ends were replaced with 6-thioguanosine. Generally, the presence of 6-thioguanosine diminishes the thermodynamic stability of RNA duplexes. This effect depends on its position within duplexes and the sequence of adjacent base pairs. However, when placed at a dangling end a 6-thioguanosine residue actually exerts a weak stabilizing effect. Furthermore, the structural effect of 6-thioguanosine substitution appears to be minimal based on NMR and Circular Dichroism (CD) data.

Unraveling the structural basis for the exceptional stability of RNA G-quadruplexes capped by a uridine tetrad at the 3' terminus.

Uridine tetrads (U-tetrads) are a structural element encountered in RNA G-quadruplexes, for example, in the structures formed by the biologically relevant human telomeric repeat RNA. For these molecules, an unexpectedly strong stabilizing influence of a U-tetrad forming at the 3' terminus of a quadruplex was reported. Here we present the high-resolution solution NMR structure of the r(UGGUGGU)4 quadruplex which, in our opinion, provides an explanation for this stabilization. Our structure features a distinctive, abrupt chain reversal just prior to the 3' uridine tetrad. Similar "reversed U-tetrads" were already observed in the crystalline phase. However, our NMR structure coupled with extensive explicit solvent molecular dynamics (MD) simulations identifies some key features of this motif that up to now remained overlooked. These include the presence of an exceptionally stable 2'OH to phosphate hydrogen bond, as well as the formation of an additional K+ binding pocket in the quadruplex groove.

Assessing protein conformational landscapes: integration of DEER data in Maximum Occurrence analysis.

The properties of the conformational landscape of a biomolecule are of capital importance to understand its function. It is widely accepted that a statistical ensemble is far more representative than a single structure, especially for proteins with disordered regions. While experimental data provide the most important handle on the conformational variability that the system is experiencing, they usually report on either time or ensemble averages. Since the available conformations largely outnumber the (independent) available experimental data, the latter can be equally well reproduced by a variety of ensembles. We have proposed the Maximum Occurrence (MaxOcc) approach to provide an upper bound of the statistical weight of each conformation. This method is expected to converge towards the true statistical weights by increasing the number of independent experimental datasets. In this paper we explore the ability of DEER (Double Electron Electron Resonance) data, which report on the distance distribution between two spin labels attached to a biomolecule, to restrain the MaxOcc values and its complementarity to previously introduced experimental techniques such as NMR and Small-Angle X-ray Scattering. We here present the case of Ca2+ bound calmodulin (CaM) as a test case and show that DEER data impose a sizeable reduction of the conformational space described by high MaxOcc conformations.

Identification of productive and futile encounters in an electron transfer protein complex.

Well-defined, stereospecific states in protein complexes are often in exchange with an ensemble of more dynamic orientations: the encounter states. The structure of the stereospecific complex between cytochrome P450cam and putidaredoxin was solved recently by X-ray diffraction as well as paramagnetic NMR spectroscopy. Other than the stereospecific complex, the NMR data clearly show the presence of additional states in the complex in solution. In these encounter states, populated for a small percentage of the time, putidaredoxin assumes multiple orientations and samples a large part of the surface of cytochrome P450cam. To characterize the nature of the encounter states, an extensive paramagnetic NMR dataset has been analyzed using the Maximum Occurrence of Regions methodology. The analysis reveals the location and maximal spatial extent of the additional states needed to fully explain the NMR data. Under the assumption of sparsity of the size of the conformational ensemble, several minor states can be located quite precisely. The distribution of these minor states correlates with the electrostatic potential map around cytochrome P450cam. Whereas some minor states are on isolated positively charged patches, others are connected to the stereospecific site via positively charged paths. The existence of electrostatically favorable pathways between the stereospecific interaction site and the different minor states or lack thereof suggests a means to discriminate between productive and futile encounter states.

How to tackle protein structural data from solution and solid state: An integrated approach.

Long-range NMR restraints, such as diamagnetic residual dipolar couplings and paramagnetic data, can be used to determine 3D structures of macromolecules. They are also used to monitor, and potentially to improve, the accuracy of a macromolecular structure in solution by validating or "correcting" a crystal model. Since crystal structures suffer from crystal packing forces they may not be accurate models for the macromolecular structures in solution. However, the presence of real differences should be tested for by simultaneous refinement of the structure using both crystal and solution NMR data. To achieve this, the program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic and paramagnetic NMR data and/or diamagnetic residual dipolar couplings. Inconsistencies between crystal structures and solution NMR data, if any, may be due either to structural rearrangements occurring on passing from the solution to solid state, or to a greater degree of conformational heterogeneity in solution with respect to the crystal. In the case of multidomain proteins, paramagnetic restraints can provide the correct mutual orientations and positions of domains in solution, as well as information on the conformational variability experienced by the macromolecule.

Inter-helical conformational preferences of HIV-1 TAR-RNA from maximum occurrence analysis of NMR data and molecular dynamics simulations.

Detecting conformational heterogeneity in biological macromolecules is a key for the understanding of their biological function. We here provide a comparison between two independent approaches to assess conformational heterogeneity: molecular dynamics simulations, performed without inclusion of any experimental data, and maximum occurrence (MaxOcc) distribution over the topologically available conformational space. The latter only reflects the extent of the averaging and identifies regions which are most compliant with the experimentally measured NMR Residual Dipolar Couplings (RDCs). The analysis was performed for the HIV-1 TAR RNA, consisting of two helical domains connected by a flexible bulge junction, for which four sets of RDCs were available as well as an 8.2 µs all-atom molecular dynamics simulation. A sample and select approach was previously applied to extract from the molecular dynamics trajectory conformational ensembles in agreement with the four sets of RDCs. The MaxOcc analysis performed here identifies the most likely sampled region in the conformational space of the system which, strikingly, overlaps well with the structures independently sampled in the molecular dynamics calculations and even better with the RDC selected ensemble.

Information content of long-range NMR data for the characterization of conformational heterogeneity.

Long-range NMR data, namely residual dipolar couplings (RDCs) from external alignment and paramagnetic data, are becoming increasingly popular for the characterization of conformational heterogeneity of multidomain biomacromolecules and protein complexes. The question addressed here is how much information is contained in these averaged data. We have analyzed and compared the information content of conformationally averaged RDCs caused by steric alignment and of both RDCs and pseudocontact shifts caused by paramagnetic alignment, and found that, despite the substantial differences, they contain a similar amount of information. Furthermore, using several synthetic tests we find that both sets of data are equally good towards recovering the major state(s) in conformational distributions.

Exploring regions of conformational space occupied by two-domain proteins.

The presence of heterogeneity in the interdomain arrangement of several biomolecules is required for their function. Here we present a method to obtain crucial clues to distinguish between different kinds of protein conformational distributions based on experimental NMR data. The method explores subregions of the conformational space and provides both upper and lower bounds of probability for the system to be in each subregion.

Site-selective photoemission from delocalized valence shells induced by molecular rotation.

Due to the generally delocalized nature of molecular valence orbitals, valence-shell spectroscopies do not usually allow to specifically target a selected atom in a molecule. However, in X-ray electron spectroscopy, the photoelectron momentum is large and the recoil angular momentum transferred to the molecule is larger when the photoelectron is ejected from a light atom compared with a heavy one. This confers an extreme sensitivity of the rotational excitation to the ionization site. Here we show that, indeed, the use of high-energy photons to photoionize valence-shell electrons of hydrogen chloride offers an unexpected way to decrypt the atomic composition of the molecular orbitals due to the rotational dependence of the photoionization profiles. The analysis of the site-specific rotational envelopes allows us to disentangle the effects of the two main mechanisms of rotational excitation, based on angular momentum exchange between the molecule and either the incoming photon or the emitted electron.