Validation of a highly sensitive HaloTag-based assay to evaluate the potency of a novel class of allosteric β-Galactosidase correctors

Site-directed Enzyme Enhancement Therapy (SEE-Tx®) technology is a disease-agnostic drug discovery tool that can be applied to any protein target of interest with a known three-dimensional structure. We used this proprietary technology to identify and characterize the therapeutic potential of structurally targeted allosteric regulators (STARs) of the lysosomal hydrolase β-galactosidase (β-Gal), which is deficient due to gene mutations in galactosidase beta 1 (GLB1)-related lysosomal storage disorders (LSDs). The biochemical HaloTag cleavage assay was used to monitor the delivery of wildtype (WT) β-Gal and four disease-related β-Gal variants (p.Ile51Thr, p.Arg59His, p.Arg201Cys and p.Trp273Leu) in the presence and absence of two identified STAR compounds. In addition, the ability of STARs to reduce toxic substrate was assessed in a canine fibroblast cell model. In contrast to the competitive pharmacological chaperone N-nonyl-deoxygalactonojirimycin (NN-DGJ), the two identified STAR compounds stabilized and substantially enhanced the lysosomal transport of wildtype enzyme and disease-causing β-Gal variants. In addition, the two STAR compounds reduced the intracellular accumulation of exogenous GM1 ganglioside, an effect not observed with the competitive chaperone NN-DGJ. This proof-of-concept study demonstrates that the SEE-Tx® platform is a rapid and cost-effective drug discovery tool for identifying STARs for the treatment of LSDs. In addition, the HaloTag assay developed in our lab has proved valuable in investigating the effect of STARs in promoting enzyme transport and lysosomal delivery. Automatization and upscaling of this assay would be beneficial for screening STARs as part of the drug discovery process.

Extracting Atomic Contributions to Binding Free Energy Using Molecular Dynamics Simulations with Mixed Solvents (MDmix)

BACKGROUND

Mixed solvents MD (MDmix) simulations have proved to be a useful and increasingly accepted technique with several applications in structure-based drug discovery. One of the assumptions behind the methodology is the transferability of free energy values from the simulated cosolvent molecules to larger drug-like molecules. However, the binding free energy maps (ΔGbind) calculated for the different moieties of the cosolvent molecules (e.g. a hydroxyl map for the ethanol) are largely influenced by the rest of the solvent molecule and do not reflect the intrinsic affinity of the moiety in question. As such, they are hardly transferable to different molecules.

METHOD

To achieve transferable energies, we present here a method for decomposing the molecular binding free energy into accurate atomic contributions.

RESULT

We demonstrate with two qualitative visual examples how the corrected energy maps better match known binding hotspots and how they can reveal hidden hotspots with actual drug design potential.

CONCLUSION

Atomic decomposition of binding free energies derived from MDmix simulations provides transferable and quantitative binding free energy maps.

Seasonal trends and spatial relations between environmental/meteorological factors and leishmaniosis sand fly vector abundances in Central Spain

This paper reports on an entomological survey performed over the period 2006-2008 in Central Spain (mainly in the Madrid province) where canine leishmaniosis (CanL) is endemic. The study area was selected on the grounds of its wide altitude range, which determines both broad climate and vegetation ranges that could affect sand fly distributions. This area was surveyed from NE to SW across its mountain range (Sistema Central) and plateau area using sticky traps mainly on embankments. In 2006 and 2007, 123 sites were sampled (9557 sand flies captured) to establish possible relations between environmental or meteorological factors and vector densities (Phlebotomus perniciosus and Phlebotomus ariasi). The factors correlated with higher vector densities were: a sample site between villages or at the edge of a village, the lack of a paved road, a rural habitat, an east or south-facing wall or wall sheltered from the wind, the presence of livestock or birds, a holm-oak wood vegetation, a lower summer mean temperature and lower annual mean precipitation. This study was followed by a seasonal survey conducted at 16 selected sites (14,353 sand flies) sampled them monthly from May to November 2008. P. perniciosus showed a diphasic seasonal trend with two abundance peaks in July and September whereas P. ariasi showed a monophasic trend with one peak in August. Comparing with data from studies performed in 1991 in the same area, vector densities are significantly higher. A possible explanation for this is that the vectors (mainly P. ariasi) are moving towards higher altitudes perhaps because of global change. This increasing trend could have an impact on CanL and its geographical distribution.

Structural properties of g,t-parallel duplexes

The structure of G,T-parallel-stranded duplexes of DNA carrying similar amounts of adenine and guanine residues is studied by means of molecular dynamics (MD) simulations and UV- and CD spectroscopies. In addition the impact of the substitution of adenine by 8-aminoadenine and guanine by 8-aminoguanine is analyzed. The presence of 8-aminoadenine and 8-aminoguanine stabilizes the parallel duplex structure. Binding of these oligonucleotides to their target polypyrimidine sequences to form the corresponding G,T-parallel triplex was not observed. Instead, when unmodified parallel-stranded duplexes were mixed with their polypyrimidine target, an interstrand Watson-Crick duplex was formed. As predicted by theoretical calculations parallel-stranded duplexes carrying 8-aminopurines did not bind to their target. The preference for the parallel-duplex over the Watson-Crick antiparallel duplex is attributed to the strong stabilization of the parallel duplex produced by the 8-aminopurines. Theoretical studies show that the isomorphism of the triads is crucial for the stability of the parallel triplex.

Discovery of 5-HT6 receptor ligands based on virtual HTS

Based on a pharmacophore alignment on a 5-HT(6) ligand applying 4SCan technology, a new lead series was identified and further structurally investigated. K(i)s down to 8 nM were achieved.

Theoretical study of the Hoogsteen-Watson-Crick junctions in DNA

A series of d (AT)(n) oligonucleotides containing mixtures of normal B-type Watson-Crick and antiparallel Hoogsteen helices have been studied using molecular dynamics simulation techniques to analyze the structural and thermodynamic impact of the junction between Watson-Crick and antiparallel Hoogsteen structures. Analysis of molecular dynamics simulations strongly suggests that for all oligonucleotides studied the antiparallel Hoogsteen appears as a reasonable conformation, only slightly less stable than the canonical B-type Watson-Crick one. The junctions between the Watson-Crick and Hoogsteen structures introduces a priori a sharp discontinuity in the helix, because the properties of each type of conformation are very well preserved in the corresponding fragments. However, and quite counterintuitively, junctions do not largely distort the duplex in structural, dynamics or energetic terms. Our results strongly support the possibility that small fragments of antiparallel Hoogsteen duplex might be embedded into large fragments of B-type Watson-Crick helices, making possible protein-DNA interactions that are specific of the antiparallel Hoogsteen conformation.

Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations?

The counterion distribution around a DNA dodecamer (5'-CGCGAATTCGCG-3') is analyzed using both standard and novel techniques based on state of the art molecular dynamics simulations. Specifically, we have explored the population of Na(+) in the minor groove of DNA duplex, and whether or not a string of Na(+) can replace the spine of hydration in the narrow AATT minor groove. The results suggest that the insertion of Na(+) in the minor groove is a very rare event, but that when once the ion finds specific sites deep inside the groove it can reside there for very long periods of time. According to our simulation the presence of Na(+) inside the groove does not have a dramatic influence in the structure or dynamics of the duplex DNA. The ability of current MD simulations to obtain equilibrated pictures of the counterion atmosphere around DNA is critically discussed.

Theoretical Study of the Guanine → 6-Thioguanine Substitution in Duplexes, Triplexes, and Tetraplexes

Molecular dynamics and thermodynamic integration calculations have been carried out on a set of G-rich single-strand, duplex, triplex, and quadruplex DNAs to study the structural and stability changes connected with the guanine --> 6-thioguanine (G --> S) mutation. The presence of 6-thioguanine leads to a shift of the geometry from the B/A intermediate to the pure B-form in duplex DNA. The G --> S mutation does not largely affect the structure of the antiparallel triplex when it is located at the reverse-Hoogsteen position, but leads to a non-negligible local distortion in the structure when it is located at the Watson-Crick position. The G --> S mutation leads to destabilization of all studied structures: the lowest effect has been observed for the G --> S mutation in the reverse-Hoogsteen strand of the triplex, a medium effect has been observed in the Watson-Crick strand of the triplex and duplex, and the highest influence of the G -->S mutation has been found for the quadruplex structures.

Effect of bulky lesions on DNA: solution structure of a DNA duplex containing a cholesterol adduct

The three-dimensional solution structure of two DNA decamers of sequence d(CCACXGGAAC)-(GTTCCGGTGG) with a modified nucleotide containing a cholesterol derivative (X) in its C1 '(chol)alpha or C1 '(chol)beta diastereoisomer form has been determined by using NMR and restrained molecular dynamics. This DNA derivative is recognized with high efficiency by the UvrB protein, which is part of the bacterial nucleotide excision repair, and the alpha anomer is repaired more efficiently than the beta one. The structures of the two decamers have been determined from accurate distance constraints obtained from a complete relaxation matrix analysis of the NOE intensities and torsion angle constraints derived from J-coupling constants. The structures have been refined with molecular dynamics methods, including explicit solvent and applying the particle mesh Ewald method to properly evaluate the long range electrostatic interactions. These calculations converge to well defined structures whose conformation is intermediate between the A- and B-DNA families as judged by the root mean square deviation but with sugar puckerings and groove shapes corresponding to a distorted B-conformation. Both duplex adducts exhibit intercalation of the cholesterol group from the major groove of the helix and displacement of the guanine base opposite the modified nucleotide. Based on these structures and molecular dynamics calculations, we propose a tentative model for the recognition of damaged DNA substrates by the UvrB protein.

Antiparallel triple helices. Structural characteristics and stabilization by 8-amino derivatives

The structural, dynamical, and recognition properties of antiparallel DNA triplexes formed by the antiparallel d(G#G.C), d(A#A.T), and d(T#A.T) motifs (the pound sign and dot mean reverse-Hoogsteen and Watson-Crick hydrogen bonds, respectively) are studied by means of "state of the art" molecular dynamics simulations. Once the characteristics of the helix are defined, molecular dynamics and thermodynamic integration calculations are used to determine the expected stabilization of the antiparallel triplex caused by the introduction of 8-aminopurines. Finally, oligonucleotides containing 8-aminopurine derivatives are synthesized and tested experimentally using several approaches in a variety of systems. A very large stabilization of the triplex is found experimentally, as predicted by simulations. These results open the possibility for the use of oligonucleotides carrying 8-aminopurines to bind single-stranded nucleic acids by formation of antiparallel triplexes.

Theoretical study of a new DNA structure: the antiparallel Hoogsteen duplex

The structure of a new form of duplex DNA, the antiparallel Hoogsteen duplex, is studied in polyd(AT) sequences by means of state-of-the-art molecular dynamics simulations in aqueous solution. The structure, which was found to be stable in all of the simulations, has many similarities with the standard Watson-Crick duplex in terms of general structure, flexibility, and molecular recognition patterns. Accurate MM-PB/SA (and MM-GB/SA) analysis shows that the new structure has an effective energy similar to that of the B-type duplex, while it is slightly disfavored by intramolecular entropic considerations. Overall, MD simulations strongly suggest that the antiparallel Hoogsteen duplex is an accessible structure for a polyd(AT) sequence, which might compete under proper experimental conditions with normal B-DNA. MD simulations also suggest that chimeras containing Watson-Crick duplex and Hoogsteen antiparallel helices might coexist in a common structure, but with the differential characteristics of both type of structures preserved.

Theoretical study of alkyl-pi and aryl-pi interactions. Reconciling theory and experiment

Quantum mechanical and quantum mechanical/molecular mechanical calculations in conjunction with continuum solvation models have been used to analyze CH-pi interactions in model systems of aryl- and alkyl-aromatic interactions, as well as in a model folding system designed to study those interactions. High level calculations reproduced accurately the interaction of CH-pi interactions in both alkyl- and aryl-based model systems. Dispersion effects dominate the interaction, but the electrostatics term is also relevant for aryl CH-pi interactions. Theoretical calculations were also used to examine the influence of CH-pi interactions in determining the conformational flexibility of folding models. Finally, a critical comparison of the results obtained from high level calculations on model systems and the experimental data derived for folding models in apolar solvents was carried out, which allowed us to reconcile the apparent discrepancy found between both data.

Hoogsteen-based parallel-stranded duplexes of DNA. Effect of 8-amino-purine derivatives

The structure of parallel-stranded duplexes of DNA-containing a mixture of guanines (G) and adenines (A) is studied by means of molecular dynamics (MD) simulation, as well as NMR and circular dichroism (CD) spectroscopy. Results demonstrate that the structure is based on the Hoogsteen motif rather than on the reverse Watson-Crick one. Molecular dynamics coupled to thermodynamic integration (MD/TI) calculations and melting experiments allowed us to determine the effect of 8-amino derivatives of A and G and of 8-amino-2'-deoxyinosine on the stability of parallel-stranded duplexes. The large stabilization of the parallel-stranded helix upon 8-amino substitution agrees with a Hoogsteen pairing, confirming MD, NMR, and CD data, and suggests new methods to obtain DNA triplexes for antigene and antisense purposes.

Theoretical studies of d(A:T)-based parallel-stranded DNA duplexes

Poly d(A:T) parallel-stranded DNA duplexes based on the Hoogsteen and reverse Watson-Crick hydrogen bond pairing are studied by means of extensive molecular dynamics (MD) simulations and molecular mechanics coupled to Poisson-Boltzmann (MM-PB/SA) calculations. The structural, flexibility, and reactivity characteristics of Hoogsteen and reverse Watson-Crick parallel duplexes are described from the analysis of the trajectories. Theoretical calculations show that the two parallel duplexes are less stable than the antiparallel Watson-Crick duplex. The difference in stability between antiparallel and parallel duplexes increases steadily as the length of the duplex increases. The reverse Watson-Crick arrangement is slightly more stable than the Hoogsteen duplex, the difference being also increased linearly with the length of the duplex. A subtle balance of intramolecular and solvation terms is responsible for the preference of a given helical structure.

Parallel-stranded hairpins containing 8-aminopurines. Novel efficient probes for triple-helix formation

We describe novel oligomers with a greater propensity to form triplexes than oligomers containing only natural bases. They consist of a polypyrimidine sequence linked head-to-head with a polypurine sequence carrying one or several 8-aminoadenine or 8-aminoguanines. The presence of 8-aminopurines also stabilised the parallel-stranded duplex structure.

The effect of amino groups on the stability of DNA duplexes and triplexes based on purines derived from inosine

The effect of amino groups attached at positions 2 and 8 of the hypoxanthine moiety in the structure, reactivity and stability of DNA duplexes and triplexes is studied by means of quantum mechanical calculations, as well as extended molecular dynamics (MD) and thermodynamic integration (MD/TI) simulations. Theoretical estimates of the change in stability related to 2'-deoxyguanosine (G) --> 2'-deoxyinosine (I) --> 8-amino-2'-deoxyinosine (8AI) mutations have been experimentally verified, after synthesis of the corresponding compounds. An amino group placed at position 2 stabilizes the duplex, as expected, and surprisingly also the triplex. The presence of an amino group at position 8 of the hypoxanthine moiety stabilizes the triplex but, surprisingly, destabilizes the duplex. The subtle electronic redistribution occurring upon the introduction of an amino group on the purine seems to be responsible for this surprising behavior. Interesting 'universal base' properties are found for 8AI.

Is polarization important in cation-pi interactions?

The importance of cation->aromatic polarization effects on cation-pi interactions has been explored. Theoretical calculations demonstrate that polarization is a large contribution to cation-aromatic interactions, and particularly to cation-pi interactions. For a series of compounds with a similar aromatic core, polarization is constant and makes small influence in the relative cation-binding energies. However, when the aromatic core changes polarization contributions might be very different. We found that the generalized molecular interaction potential with polarization is a very fast and powerful tool for the prediction of cation binding of aromatic compounds.

Antiviral effects of xanthate D609 on the human respiratory syncytial virus growth cycle

The antiviral compound tricyclo-decan-9-yl-xanthogenate (D609) inhibits respiratory syncytial (RS) virus growth in human epithelial (Hep 2) cells. D609 treatment resulted in a decrease in the accumulation of viral proteins, in the phosphorylation of the viral phosphoprotein, and in the amount of extracellular antigens and infectious particles. The relative accumulation of viral proteins was also unbalanced, however no differences were found in the amount of viral RNA with plus or minus polarity. In addition nucleocapsids formation was not inhibited. These observations suggested that this antiviral compound affects the relative proportion of viral proteins and the phosphorylation of P protein. Both features appear to be important in RS virus morphogenesis.