Ligand-induced unfolding mechanism of an RNA G-quadruplex

The cationic porphyrin, TMPyP4, is a well-established DNA G-quadruplex (G4) binding ligand that can stabilize different topologies via multiple binding modes. However, TMPyP4 has completely opposite destabilizing and unwinding effect on RNA G4 structures. The structural mechanisms that mediate RNA G4 unfolding remains unknown. Here, we report on the TMPyP4-induced RNA G4 unfolding mechanism studied by well-tempered metadynamics (WT-MetaD) with supporting biophysical experiments. The simulations predict a two-state mechanism of TMPyP4 interaction via a groove-bound and a top-face bound conformation. The dynamics of TMPyP4 stacking on the top tetrad disrupts Hoogsteen H-bonds between guanine bases resulting in the consecutive TMPyP4 intercalation from top-to-bottom G-tetrads. The results reveal a striking correlation between computational and experimental approaches and validate WT-MetaD simulations as a powerful tool for studying RNA G4-ligand interactions

Exploring the Dynamics of Propeller Loops in Human Telomeric DNA Quadruplexes Using Atomistic Simulations

We have carried out a series of extended unbiased molecular dynamics (MD) simulations (up to 10 \u3bcs long, 3c162 \u3bcs in total) complemented by replica-exchange with the collective variable tempering (RECT) approach for several human telomeric DNA G-quadruplex (GQ) topologies with TTA propeller loops. We used different AMBER DNA force-field variants and also processed simulations by Markov State Model (MSM) analysis. The slow conformational transitions in the propeller loops took place on a scale of a few \u3bcs, emphasizing the need for long simulations in studies of GQ dynamics. The propeller loops sampled similar ensembles for all GQ topologies and for all force-field dihedral-potential variants. The outcomes of standard and RECT simulations were consistent and captured similar spectrum of loop conformations. However, the most common crystallographic loop conformation was very unstable with all force-field versions. Although the loss of canonical \u3b3-trans state of the first propeller loop nucleotide could be related to the indispensable bsc0 \u3b1/\u3b3 dihedral potential, even supporting this particular dihedral by a bias was insufficient to populate the experimentally dominant loop conformation. In conclusion, while our simulations were capable of providing a reasonable albeit not converged sampling of the TTA propeller loop conformational space, the force-field description still remained far from satisfactory

Dynamic Profiling of β-Coronavirus 3CL M^proProtease Ligand-Binding Sites

Data availability statement: The trajectories of Mpro simulations and models of the metastable states can be downloaded from 10.5281/zenodo.4782284.β-coronavirus (CoVs) alone has been responsible for three major global outbreaks in the 21st century. The current crisis has led to an urgent requirement to develop therapeutics. Even though a number of vaccines are available, alternative strategies targeting essential viral components are required as a backup against the emergence of lethal viral variants. One such target is the main protease (Mpro) that plays an indispensable role in viral replication. The availability of over 270 Mpro X-ray structures in complex with inhibitors provides unique insights into ligand–protein interactions. Herein, we provide a comprehensive comparison of all nonredundant ligand-binding sites available for SARS-CoV2, SARS-CoV, and MERS-CoV Mpro. Extensive adaptive sampling has been used to investigate structural conservation of ligand-binding sites using Markov state models (MSMs) and compare conformational dynamics employing convolutional variational auto-encoder-based deep learning. Our results indicate that not all ligand-binding sites are dynamically conserved despite high sequence and structural conservation across β-CoV homologs. This highlights the complexity in targeting all three Mpro enzymes with a single pan inhibitor.There was no funding for this wor

Selectivity of major isoquinoline alkaloids from Chelidonium majus towards telomeric G-quadruplex: A study using a transition-FRET (t-FRET) assay

Background Natural bioproducts are invaluable resources in drug discovery. Isoquinoline alkaloids of Chelidonium majus constitute a structurally diverse family of natural products that are of great interest, one of them being their selectivity for human telomeric G-quadruplex structure and telomerase inhibition. Methods The study focuses on the mechanism of telomerase inhibition by stabilization of telomeric G-quadruplex structures by berberine, chelerythrine, chelidonine, sanguinarine and papaverine. Telomerase activity and mRNA levels of hTERT were estimated using quantitative telomere repeat amplification protocol (q-TRAP) and qPCR, in MCF-7 cells treated with different groups of alkaloids. The selectivity of the main isoquinoline alkaloids of Chelidonium majus towards telomeric G-quadruplex forming sequences were explored using a sensitive modified thermal FRET-melting measurement in the presence of the complementary oligonucleotide CT22. We assessed and monitored G-quadruplex topologies using circular dichroism (CD) methods, and compared spectra to previously well-characterized motifs, either alone or in the presence of the alkaloids. Molecular modeling was performed to rationalize ligand binding to the G-quadruplex structure. Results The results highlight strong inhibitory effects of chelerythrine, sanguinarine and berberine on telomerase activity, most likely through substrate sequestration. These isoquinoline alkaloids interacted strongly with telomeric sequence G-quadruplex. In comparison, chelidonine and papaverine had no significant interaction with the telomeric quadruplex, while they strongly inhibited telomerase at transcription level of hTERT. Altogether, all of the studied alkaloids showed various levels and mechanisms of telomerase inhibition. Conclusions We report on a comparative study of anti-telomerase activity of the isoquinoline alkaloids of Chelidonium majus. Chelerythrine was most effective in inhibiting telomerase activity by substrate sequesteration through G-quadruplex stabilization. General significance Understanding structural and molecular mechanisms of anti-cancer agents can help in developing new and more potent drugs with fewer side effects. Isoquinolines are the most biologically active agents from Chelidonium majus, which have shown to be telomeric G-quadruplex stabilizers and potent telomerase inhibitors

Interaction of anesthetic supplement thiopental with human serum albumin

Thiopental (TPL) is a commonly used barbiturate anesthetic. Its binding with human serum albumin (HSA) was studied to explore the anesthetic-induced protein dysfunction. The basic binding interaction was studied by UV-absorption and fluorescence spectroscopy. An increase in the binding affinity (K) and in the number of binding sites (n) with the increasing albumin concentration was observed. The interaction was conformation-dependent and the highest for the F isomer of HSA, which implicates its slow elimination. The mode of binding was characterized using various thermodynamic parameters. Domain II of HSA was found to possess a high affinity binding site for TPL. The effect of micro-metal ions on the binding affinity was also investigated. The molecular distance, r, between donor (HSA) and acceptor (TPL) was estimated by fluorescence resonance energy transfer (FRET). Correlation between the stability of the TPL-N and TPL-F complexes and drug distribution is discussed. The structural changes in the protein investigated by circular dichroism (CD) and Fourier transform infrared (FT-IR) spectroscopy reflect perturbation of the albumin molecule and provide an explanation for the heterogeneity of action of this anesthetic

Phosphorylation of Histone Deacetylase 8: Structural and Mechanistic Analysis of the Phosphomimetic S39E Mutant

International audienc

Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G‑Quadruplexes

G-quadruplexes (GQs) are key noncanonical DNA and RNA architectures stabilized by desolvated monovalent cations present in their central channels. We analyze extended atomistic molecular dynamics simulations (∼580 μs in total) of GQs with 11 monovalent cation parametrizations, assessing GQ overall structural stability, dynamics of internal cations, and distortions of the G-tetrad geometries. Majority of simulations were executed with the SPC/E water model; however, test simulations with TIP3P and OPC water models are also reported. The identity and parametrization of ions strongly affect behavior of a tetramolecular d­[GGG]₄ GQ, which is unstable with several ion parametrizations. The remaining studied RNA and DNA GQs are structurally stable, though the G-tetrad geometries are always deformed by bifurcated H-bonding in a parametrization-specific manner. Thus, basic 10-μs-scale simulations of fully folded GQs can be safely done with a number of cation parametrizations. However, there are parametrization-specific differences and basic force-field errors affecting the quantitative description of ion-tetrad interactions, which may significantly affect studies of the ion-binding processes and description of the GQ folding landscape. Our d­[GGG]₄ simulations indirectly suggest that such studies will also be sensitive to the water models. During exchanges with bulk water, the Na⁺ ions move inside the GQs in a concerted manner, while larger relocations of the K⁺ ions are typically separated. We suggest that the Joung-Cheatham SPC/E K⁺ parameters represent a safe choice in simulation studies of GQs, though variation of ion parameters can be used for specific simulation goals