CAG RNAs induce DNA damage and apoptosis by silencing NUDT16 expression in polyglutamine degeneration.

DNA damage plays a central role in the cellular pathogenesis of polyglutamine (polyQ) diseases, including Huntington's disease (HD). In this study, we showed that the expression of untranslatable expanded CAG RNA per se induced the cellular DNA damage response pathway. By means of RNA sequencing (RNA-seq), we found that expression of the Nudix hydrolase 16 (NUDT16) gene was down-regulated in mutant CAG RNA-expressing cells. The loss of NUDT16 function results in a misincorporation of damaging nucleotides into DNAs and leads to DNA damage. We showed that small CAG (sCAG) RNAs, species generated from expanded CAG transcripts, hybridize with CUG-containing NUDT16 mRNA and form a CAG-CUG RNA heteroduplex, resulting in gene silencing of NUDT16 and leading to the DNA damage and cellular apoptosis. These results were further validated using expanded CAG RNA-expressing mouse primary neurons and in vivo R6/2 HD transgenic mice. Moreover, we identified a bisamidinium compound, DB213, that interacts specifically with the major groove of the CAG RNA homoduplex and disfavors the CAG-CUG heteroduplex formation. This action subsequently mitigated RNA-induced silencing complex (RISC)-dependent NUDT16 silencing in both in vitro cell and in vivo mouse disease models. After DB213 treatment, DNA damage, apoptosis, and locomotor defects were rescued in HD mice. This work establishes NUDT16 deficiency by CAG repeat RNAs as a pathogenic mechanism of polyQ diseases and as a potential therapeutic direction for HD and other polyQ diseases.

Identification of 13 Guanidinobenzoyl- or Aminidinobenzoyl-Containing Drugs to Potentially Inhibit TMPRSS2 for COVID-19 Treatment.

Positively charged groups that mimic arginine or lysine in a natural substrate of trypsin are necessary for drugs to inhibit the trypsin-like serine protease TMPRSS2 that is involved in the viral entry and spread of coronaviruses, including SARS-CoV-2. Based on this assumption, we identified a set of 13 approved or clinically investigational drugs with positively charged guanidinobenzoyl and/or aminidinobenzoyl groups, including the experimentally verified TMPRSS2 inhibitors Camostat and Nafamostat. Molecular docking using the C-I-TASSER-predicted TMPRSS2 catalytic domain model suggested that the guanidinobenzoyl or aminidinobenzoyl group in all the drugs could form putative salt bridge interactions with the side-chain carboxyl group of Asp435 located in the S1 pocket of TMPRSS2. Molecular dynamics simulations further revealed the high stability of the putative salt bridge interactions over long-time (100 ns) simulations. The molecular mechanics/generalized Born surface area-binding free energy assessment and per-residue energy decomposition analysis also supported the strong binding interactions between TMPRSS2 and the proposed drugs. These results suggest that the proposed compounds, in addition to Camostat and Nafamostat, could be effective TMPRSS2 inhibitors for COVID-19 treatment by occupying the S1 pocket with the hallmark positively charged groups.

Detecting Protein-Ligand Interaction from Integrated Transient Induced Molecular Electronic Signal (i-TIMES).

Quantitative information about protein-ligand interactions is central to drug discovery. To obtain the quintessential reaction dissociation constant, ideally measurements of reactions should be performed without perturbations by molecular labeling or immobilization. The technique of transient induced molecular electrical signal (TIMES) has provided a promising technique to meet such requirements, and its performance in a microfluidic environment further offers the potential for high throughput and reduced consumption of reagents. In this work, we further the development by using integrated TIMES signal (i-TIMES) to greatly enhance the accuracy and reproducibility of the measurement. While the transient response may be of interest, the integrated signal directly measures the total amount of surface charge density resulted from molecules near the surface of electrode. The signals enable quantitative characterization of protein-ligand interactions. We have demonstrated the feasibility of i-TIMES technique using different biomolecules including lysozyme, N,N',N″-triacetylchitotriose (TriNAG), aptamer, p-aminobenzamidine (pABA), bovine pancreatic ribonuclease A (RNaseA), and uridine-3'-phosphate (3'UMP). The results show i-TIMES is a simple and accurate technique that can bring tremendous value to drug discovery and research of intermolecular interactions.

Label-free analytical performances of a peptide-based QCM biosensor for trypsin.

A label-free biosensor was fabricated for the detection of trypsin by using a peptide-functionalized quartz crystal microbalance gold electrode. The synthetized peptide chains were immobilized tightly on the QCM electrode via a self-assembly method, which formed a thin and approximate rigid layer of peptides. The detection signal was achieved by calculating the mass changes on the QCM electrode because the peptide chains could be specifically cleaved in the carboxyl terminuses of arginine and lysine by trypsin. When gold nanoparticles were coupled to the peptide chains, the sensing signal would be amplified 10.9 times. Furthermore, the sensor interface shows a lower resonance resistance change when the peptide chain is immobilized horizontally. Independent detections in parallel on different electrodes have a wide linear range. Under the optimum conditions, the signal-amplified biosensor allowed the measurement of trypsin over the range of 0-750 ng mL-1 with a detection limit of 8.6 ng mL-1. Moreover, for screening the inhibitor of trypsin, the IC50 values were obtained to be 1.85 µg mL-1 for benzamidine hydrochloride and 20.5 ng mL-1 for the inhibitor from soybean.

Crystal structure of a human GABAA receptor.

Type-A γ-aminobutyric acid receptors (GABAARs) are the principal mediators of rapid inhibitory synaptic transmission in the human brain. A decline in GABAAR signalling triggers hyperactive neurological disorders such as insomnia, anxiety and epilepsy. Here we present the first three-dimensional structure of a GABAAR, the human ß3 homopentamer, at 3 Å resolution. This structure reveals architectural elements unique to eukaryotic Cys-loop receptors, explains the mechanistic consequences of multiple human disease mutations and shows an unexpected structural role for a conserved N-linked glycan. The receptor was crystallized bound to a previously unknown agonist, benzamidine, opening a new avenue for the rational design of GABAAR modulators. The channel region forms a closed gate at the base of the pore, representative of a desensitized state. These results offer new insights into the signalling mechanisms of pentameric ligand-gated ion channels and enhance current understanding of GABAergic neurotransmission.

Kallikrein 5 inhibitors identified through structure based drug design in search for a treatment for Netherton Syndrome.

Netherton syndrome (NS) is a rare and debilitating severe autosomal recessive genetic skin disease with high mortality rates particularly in neonates. NS is caused by loss-of-function SPINK5 mutations leading to unregulated kallikrein 5 (KLK5) and kallikrein 7 (KLK7) activity. Furthermore, KLK5 inhibition has been proposed as a potential therapeutic treatment for NS. Identification of potent and selective KLK5 inhibitors would enable further exploration of the disease biology and could ultimately lead to a treatment for NS. This publication describes how fragmentation of known trypsin-like serine protease (TLSP) inhibitors resulted in the identification of a series of phenolic amidine-based KLK5 inhibitors 1. X-ray crystallography was used to find alternatives to the phenol interaction leading to identification of carbonyl analogues such as lactam 13 and benzimidazole 15. These reversible inhibitors, with selectivity over KLK1 (10-100 fold), provided novel starting points for the guided growth towards suitable tool molecules for the exploration of KLK5 biology.

Small molecule binders recognize DNA microstructural variations via an induced fit mechanism.

Regulation of gene-expression by specific targeting of protein-nucleic acid interactions has been a long-standing goal in medicinal chemistry. Transcription factors are considered "undruggable" because they lack binding sites well suited for binding small-molecules. In order to overcome this obstacle, we are interested in designing small molecules that bind to the corresponding promoter sequences and either prevent or modulate transcription factor association via an allosteric mechanism. To achieve this, we must design small molecules that are both sequence-specific and able to target G/C base pair sites. A thorough understanding of the relationship between binding affinity and the structural aspects of the small molecule-DNA complex would greatly aid in rational design of such compounds. Here we present a comprehensive analysis of sequence-specific DNA association of a synthetic minor groove binder using long timescale molecular dynamics. We show how binding selectivity arises from a combination of structural factors. Our results provide a framework for the rational design and optimization of synthetic small molecules in order to improve site-specific targeting of DNA for therapeutic uses in the design of selective DNA binders targeting transcription regulation.

Interactive molecular dynamics in virtual reality from quantum chemistry to drug binding: An open-source multi-person framework.

As molecular scientists have made progress in their ability to engineer nanoscale molecular structure, we face new challenges in our ability to engineer molecular dynamics (MD) and flexibility. Dynamics at the molecular scale differs from the familiar mechanics of everyday objects because it involves a complicated, highly correlated, and three-dimensional many-body dynamical choreography which is often nonintuitive even for highly trained researchers. We recently described how interactive molecular dynamics in virtual reality (iMD-VR) can help to meet this challenge, enabling researchers to manipulate real-time MD simulations of flexible structures in 3D. In this article, we outline various efforts to extend immersive technologies to the molecular sciences, and we introduce "Narupa," a flexible, open-source, multiperson iMD-VR software framework which enables groups of researchers to simultaneously cohabit real-time simulation environments to interactively visualize and manipulate the dynamics of molecular structures with atomic-level precision. We outline several application domains where iMD-VR is facilitating research, communication, and creative approaches within the molecular sciences, including training machines to learn potential energy functions, biomolecular conformational sampling, protein-ligand binding, reaction discovery using "on-the-fly" quantum chemistry, and transport dynamics in materials. We touch on iMD-VR's various cognitive and perceptual affordances and outline how these provide research insight for molecular systems. By synergistically combining human spatial reasoning and design insight with computational automation, technologies such as iMD-VR have the potential to improve our ability to understand, engineer, and communicate microscopic dynamical behavior, offering the potential to usher in a new paradigm for engineering molecules and nano-architectures.

Imino proton NMR guides the reprogramming of A•T specific minor groove binders for mixed base pair recognition.

Sequence-specific binding to DNA is crucial for targeting transcription factor-DNA complexes to modulate gene expression. The heterocyclic diamidine, DB2277, specifically recognizes a single G•C base pair in the minor groove of mixed base pair sequences of the type AAAGTTT. NMR spectroscopy reveals the presence of major and minor species of the bound compound. To understand the principles that determine the binding affinity and orientation in mixed sequences of DNA, over thirty DNA hairpin substrates were examined by NMR and thermal melting. The NMR exchange dynamics between major and minor species shows that the exchange is much faster than compound dissociation determined from biosensor-surface plasmon resonance. Extensive modifications of DNA sequences resulted in a unique DNA sequence with binding site AAGATA that binds DB2277 in a single orientation. A molecular docking result agrees with the model representing rapid flipping of DB2277 between major and minor species. Imino spectral analysis of a (15)N-labeled central G clearly shows the crucial role of the exocyclic amino group of G in sequence-specific recognition. Our results suggest that this approach can be expanded to additional modules for recognition of more sequence-specific DNA complexes. This approach provides substantial information about the sequence-specific, highly efficient, dynamic nature of minor groove binding agents.

Multiple Ligand Unbinding Pathways and Ligand-Induced Destabilization Revealed by WExplore.

We report simulations of full ligand exit pathways for the trypsin-benzamidine system, generated using the sampling technique WExplore. WExplore is able to observe millisecond-scale unbinding events using many nanosecond-scale trajectories that are run without introducing biasing forces. The algorithm generates rare events by dividing the coordinate space into regions, on-the-fly, and balancing computational effort between regions through cloning and merging steps, as in the weighted ensemble method. The averaged exit flux yields a ligand exit rate of 180 µs, which is within an order of magnitude of the experimental value. We obtain broad sampling of ligand exit pathways, and visualize our findings using conformation space networks. The analysis shows three distinct exit channels, two of which are formed through large, rare motions of the loop regions in trypsin. This broad set of ligand-bound poses is then used to investigate general properties of ligand binding: we observe both a direct stabilizing effect of ligand-protein interactions and an indirect destabilizing effect on intraprotein interactions that is induced by the ligand. Significantly, the crystallographic binding poses are distinguished not only because their ligands induce large stabilizing effects, but also because they induce relatively low indirect destabilizations.

Funnel metadynamics as accurate binding free-energy method.

A detailed description of the events ruling ligand/protein interaction and an accurate estimation of the drug affinity to its target is of great help in speeding drug discovery strategies. We have developed a metadynamics-based approach, named funnel metadynamics, that allows the ligand to enhance the sampling of the target binding sites and its solvated states. This method leads to an efficient characterization of the binding free-energy surface and an accurate calculation of the absolute protein-ligand binding free energy. We illustrate our protocol in two systems, benzamidine/trypsin and SC-558/cyclooxygenase 2. In both cases, the X-ray conformation has been found as the lowest free-energy pose, and the computed protein-ligand binding free energy in good agreement with experiments. Furthermore, funnel metadynamics unveils important information about the binding process, such as the presence of alternative binding modes and the role of waters. The results achieved at an affordable computational cost make funnel metadynamics a valuable method for drug discovery and for dealing with a variety of problems in chemistry, physics, and material science.

Fluorescence-labeled bis-benzamidines as fluorogenic DNA minor-groove binders: photophysics and binding dynamics.

In recent decades there has been great interest in the design of highly sensitive sequence-specific DNA binders. The eligibility of the binder depends on the magnitude of the fluorescence increase upon binding, related to its photophysics, and on its affinity and specificity, which is, in turn, determined by the dynamics of the binding process. Therefore, progress in the design of DNA binders requires both thorough photophysical studies and precise determination of the association and dissociation rate constants involved. We have studied two bis-benzamidine (BBA) derivatives labeled by linkers of various lengths with the dye Oregon Green (OG). These fluorogenic binders show a dramatic fluorescence enhancement upon binding to the minor groove of double-stranded (ds) DNA, as well as significant improvement in their sequence specificity versus the parent BBA, although with decreased affinity constants. Detailed photophysical analysis shows that static and dynamic quenching of the OG fluorescence by BBA through photoinduced electron transfer is suppressed upon insertion of BBA into the minor groove of DNA. Fluorescence correlation spectroscopy yields precise dynamic rate constants that prove that the association process of these fluorogenic binders to dsDNA is very similar to that of BBA alone and that their lower affinity is mainly a consequence of their weaker attachment to the minor groove and the resultant faster dissociation process. The conclusions of this study will allow us to go one step further in the design of new DNA binders with tunable fluorescence and binding properties.

Affinity screening using competitive binding with fluorine-19 hyperpolarized ligands.

Fluorine-19 NMR and hyperpolarization form a powerful combination for drug screening. Under a competitive equilibrium with a selected fluorinated reporter ligand, the dissociation constant (K(D)) of other ligands of interest is measurable using a single-scan Carr-Purcell-Meiboom-Gill (CPMG) experiment, without the need for a titration. This method is demonstrated by characterizing the binding of three ligands with different affinities for the serine protease trypsin. Monte Carlo simulations show that the highest accuracy is obtained when about one-half of the bound reporter ligand is displaced in the binding competition. Such conditions can be achieved over a wide range of affinities, allowing for rapid screening of non-fluorinated compounds when a single fluorinated ligand for the binding pocket of interest is known.

Pharmacokinetic comparison to determine the mechanisms underlying the differential efficacies of cationic diamidines against first- and second-stage human African trypanosomiasis.

Human African trypanosomiasis (HAT), a neglected tropical disease, is fatal without treatment. Pentamidine, a cationic diamidine, has been used to treat first-stage (hemolymphatic) HAT since the 1940s, but it is ineffective against second-stage (meningoencephalitic, or central nervous system [CNS]) infection. Novel diamidines (DB75, DB820, and DB829) have shown promising efficacy in both mouse and monkey models of first-stage HAT. However, only DB829 cured animals with second-stage infection. In this study, we aimed to determine the mechanisms underlying the differential efficacies of these diamidines against HAT by conducting a comprehensive pharmacokinetic characterization. This included the determination of metabolic stability in liver microsomes, permeability across MDCK and MDR1-MDCK cell monolayers, interaction with the efflux transporter MDR1 (P-glycoprotein 1 or P-gp), drug binding in plasma and brain, and plasma and brain concentration-time profiles after a single dose in mice. The results showed that DB829, an azadiamidine, had the highest systemic exposure and brain-to-plasma ratio, whereas pentamidine and DB75 had the lowest. None of these diamidines was a P-gp substrate, and the binding of each to plasma proteins and brain differed greatly. The brain-to-plasma ratio best predicted the relative efficacies of these diamidines in mice with second-stage infection. In conclusion, pharmacokinetics and CNS penetration influenced the in vivo efficacies of cationic diamidines against first- and second-stage HAT and should be considered when developing CNS-active antitrypanosomal diamidines.

Functional characterization of protein variants encoded by nonsynonymous single nucleotide polymorphisms in MARC1 and MARC2 in healthy Caucasians.

Human molybdenum-containing enzyme mitochondrial amidoxime reducing component (mARC), cytochrome b5 type B, and NADH cytochrome b5 reductase form an N-reductive enzyme system that is capable of reducing N-hydroxylated compounds. Genetic variations are known, but their functional relevance is unclear. Our study aimed to investigate the incidence of nonsynonymous single nucleotide polymorphisms (SNPs) in the mARC genes in healthy Caucasian volunteers, to determine saturation of the protein variants with molybdenum cofactor (Moco), and to characterize the kinetic behavior of the protein variants by in vitro biotransformation studies. Genotype frequencies of six SNPs in the mARC genes (c.493A>G, c.560T>A, c.736T>A, and c.739G>C in MARC1; c.730G>A and c.735T>G in MARC2) were determined by pyrosequencing in a cohort of 340 healthy Caucasians. Protein variants were expressed in Escherichia coli. Saturation with Moco was determined by measurement of molybdenum by inductively coupled mass spectrometry. Steady state assays were performed with benzamidoxime. The six variants were of low frequency in this Caucasian population. Only one homozygous variant (c.493A; MARC1) was detected. All protein variants were able to bind Moco. Steady state assays showed statistically significant decreases of catalytic efficiency values for the mARC-2 wild type compared with the mARC-1 wild type (P < 0.05) and for two mARC-2 variants compared with the mARC-2 wild type (G244S, P < 0.05; C245W, P < 0.05). After simultaneous substitution of more than two amino acids in the mARC-1 protein, N-reductive activity was decreased 5-fold. One homozygous variant of MARC1 was detected in our sample. The encoded protein variant (A165T) showed no different kinetic parameters in the N-reduction of benzamidoxime.

Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations.

The understanding of protein-ligand binding is of critical importance for biomedical research, yet the process itself has been very difficult to study because of its intrinsically dynamic character. Here, we have been able to quantitatively reconstruct the complete binding process of the enzyme-inhibitor complex trypsin-benzamidine by performing 495 molecular dynamics simulations of free ligand binding of 100 ns each, 187 of which produced binding events with an rmsd less than 2 Å compared to the crystal structure. The binding paths obtained are able to capture the kinetic pathway of the inhibitor diffusing from solvent (S0) to the bound (S4) state passing through two metastable intermediate states S2 and S3. Rather than directly entering the binding pocket the inhibitor appears to roll on the surface of the protein in its transition between S3 and the final binding pocket, whereas the transition between S2 and the bound pose requires rediffusion to S3. An estimation of the standard free energy of binding gives ΔG° = -5.2 ± 0.4 kcal/mol (cf. the experimental value -6.2 kcal/mol), and a two-states kinetic model k(on) = (1.5 ± 0.2) × 10(8) M(-1) s(-1) and k(off) = (9.5 ± 3.3) × 10(4) s(-1) for unbound to bound transitions. The ability to reconstruct by simple diffusion the binding pathway of an enzyme-inhibitor binding process demonstrates the predictive power of unconventional high-throughput molecular simulations. Moreover, the methodology is directly applicable to other molecular systems and thus of general interest in biomedical and pharmaceutical research.

On the reproducibility of protein crystal structures: five atomic resolution structures of trypsin.

Structural studies of proteins usually rely on a model obtained from one crystal. By investigating the details of this model, crystallographers seek to obtain insight into the function of the macromolecule. It is therefore important to know which details of a protein structure are reproducible or to what extent they might differ. To address this question, the high-resolution structures of five crystals of bovine trypsin obtained under analogous conditions were compared. Global parameters and structural details were investigated. All of the models were of similar quality and the pairwise merged intensities had large correlation coefficients. The C(α) and backbone atoms of the structures superposed very well. The occupancy of ligands in regions of low thermal motion was reproducible, whereas solvent molecules containing heavier atoms (such as sulfur) or those located on the surface could differ significantly. The coordination lengths of the calcium ion were conserved. A large proportion of the multiple conformations refined to similar occupancies and the residues adopted similar orientations. More than three quarters of the water-molecule sites were conserved within 0.5 Šand more than one third were conserved within 0.1 Å. An investigation of the protonation states of histidine residues and carboxylate moieties was consistent for all of the models. Radiation-damage effects to disulfide bridges were observed for the same residues and to similar extents. Main-chain bond lengths and angles averaged to similar values and were in agreement with the Engh and Huber targets. Other features, such as peptide flips and the double conformation of the inhibitor molecule, were also reproducible in all of the trypsin structures. Therefore, many details are similar in models obtained from different crystals. However, several features of residues or ligands located in flexible parts of the macromolecule may vary significantly, such as side-chain orientations and the occupancies of certain fragments.

Investigation of the effect of structure modification of furamidine on the DNA minor groove binding and antiprotozoal activity.

New analogs of the antiprotozoal agent Furamidine were prepared utilizing Stille coupling reactions and amidation of the bisnitrile intermediate using lithium bis-trimethylsilylamide. Both the phenyl groups and the furan moiety of furamidine were replaced by heterocycles including thiophene, selenophene, indole or benzimidazole. Based upon the ΔTm and the CD results, the new compounds showed strong binding to the DNA minor groove. The new analogues are also more active both in vitro and in vivo than furamidine. Compounds 7a, 7b, and 7f showed the highest activity in vivo by curing 75% of animals, and this merits further evaluation.

Crystallographic and kinetic evidence of allostery in a trypsin-like protease.

Protein allostery is based on the existence of multiple conformations in equilibrium linked to distinct functional properties. Although evidence of allosteric transitions is relatively easy to identify by functional studies, structural detection of a pre-existing equilibrium between alternative conformations remains challenging even for textbook examples of allosteric proteins. Kinetic studies show that the trypsin-like protease thrombin exists in equilibrium between two conformations where the active site is either collapsed (E*) or accessible to substrate (E). However, structural demonstration that the two conformations exist in the same enzyme construct free of ligands has remained elusive. Here we report the crystal structure of the thrombin mutant N143P in the E form, which complements the recently reported structure in the E* form, and both the E and E* forms of the thrombin mutant Y225P. The side chain of W215 moves 10.9 Å between the two forms, causing a displacement of 6.6 Å of the entire 215-217 segment into the active site that in turn opens or closes access to the primary specificity pocket. Rapid kinetic measurements of p-aminobenzamidine binding to the active site confirm the existence of the E*-E equilibrium in solution for wild-type and the mutants N143P and Y225P. These findings provide unequivocal proof of the allosteric nature of thrombin and lend strong support to the recent proposal that the E*-E equilibrium is a key property of the trypsin fold.