Formation of G‐wires,bimolecular and tetramolecular quadruplex: Cation‐induced structural polymorphs of G‐rich DNA sequence of human SYTX gene

An exceptional property of auto‐folding into a range of intra‐ as well as intermolecular quadruplexes by guanine‐rich oligomers (GROs) of promoters, telomeres and various other genomic locations is still one of the most attractive areas of research at present times. The main reason for this attention is due to their established in vivo existence and biological relevance. Herein, the structural status of a 20‐nt long G‐rich sequence with two G5 stretches (SG20) is investigated using various biophysical and biochemical techniques. Bioinformatics analysis suggested the presence of a 17‐nt stretch of this SG20 sequence in the intronic region of human SYTX (Synaptotagmin 10) gene. The SYTX gene helps in sensing out the Ca²⁺ ion, causing its intake in the pre‐synaptic neuron. A range of various topologies like bimolecular, tetramolecular and guanine‐wires (nano‐wires) was exhibited by the studied sequence, as a function of cations (Na⁺/K⁺) concentration. UV‐thermal denaturation, gel electrophoresis, and circular dichroism (CD) spectroscopy showed correlations and established a cation‐dependent structural switch. The G‐wire formation, in the presence of K⁺, may further be explored for its possible relevance in nano‐biotechnological applications.     

Biophysical probing of the binding properties of a Cu(II) complex with G‐quadruplex DNA: an experimental and computational study

Telomerase inhibition through G‐quadruplex stabilization by small molecules is of great interest as a novel anticancer therapeutic strategy. Here, we show that newly synthesized Cu‐complex binds to G‐quadruplex DNA and induces changes in its stability. This biophysical interaction was investigated in vitro using spectroscopic, voltammetric and computational techniques. The binding constant for this complex to G‐quadruplex using spectroscopic and electrochemical methods is in the order of 10⁵. The binding stoichiometry was investigated using spectroscopic techniques and corresponded to a ratio of 1: 1. Fluorescence titration results reveal that Cu‐complex is quenched in the presence of G‐quadruplex DNA. Analysis of the fluorescence emission at different temperatures shows that ΔH° > 0, ΔS° > 0 and ΔG° < 0, and indicates that hydrophobic interactions played a major role in the binding processes. MD simulation results suggested that this ligand could stabilize the G‐quadruplex structure. An optimized docked model of the G‐quadruplex–ligand mixture confirmed the experimental results. Based on the results, we conclude that Cu‐complex as an anticancer candidate can bind and stabilize the G‐quadruplex DNA structure. Copyright © 2015 John Wiley & Sons, Ltd.     

Characterization of the G‐quadruplex structure of a catalytic DNA with peroxidase activity

It has been reported that the complexes formed by hemin and some G‐quadruplexes can be developed as a new class of DNAzyme with peroxidase activity. This kind of DNAzyme has received a great deal of attention. But to date, the actual G‐quadruplex structure that can provide hemin with enhanced peroxidase activity is in doubt. Herein, the G‐quadruplex structure of CatG4, a 21‐nucleotide DNA oligomer which was previously reported to bind hemin and the resulting complex exhibiting enhanced peroxidase activity, was characterized by fluorescence and circular dichroism measurements. The results suggest that the catalytically active form of CatG4 may be a unimolecular parallel quadruplex rather than a unimolecular chair‐type antiparallel quadruplex or a multistranded parallel quadruplex. In addition, the fluorescence analysis of labeled oligonucleotides may be developed as a supplementary tool for the study of DNA conformations. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 331–339, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Recognition of human telomeric G‐quadruplex DNA by berberine analogs: effect of substitution at the 9 and 13 positions of the isoquinoline moiety

G‐quadruplex forming sequences are widely distributed in human genome and serve as novel targets for regulating gene expression and chromosomal maintenance. They offer unique targets for anticancer drug development. Here, the interaction of berberine (BC) and two of its analogs bearing substitution at 9 and 13‐position with human telomeric G‐quadruplex DNA sequence has been investigated by biophysical techniques. Both the analogs exhibited several‐fold higher binding affinity than berberine. The Scatchard binding isotherms revealed non‐cooperative binding. 9‐ω‐amino hexyl ether analog (BC1) showed highest affinity (1.8 × 10⁶ M^?1) while the affinity of the 13‐phenylpropyl analog (BC2) was 1.09 × 10⁶ M^?1. Comparative fluorescence quenching and polarization anisotropy of the emission spectra gave evidence for a stronger stacking interaction of the analogs compared to berberine. The thiazole orange displacement assay has clearly established that the analogs were more effective in displacing the end stacked dye in comparison to berberine. However, the binding of the analogs did not induce any major structural perturbation in the G‐quadruplex structure, but led to higher thermal stability. Energetics of the binding indicated that the association of the analogs was exothermic and predominantly entropy driven phenomenon. Increasing the temperature resulted in weaker binding; the enthalpic contribution increased and the entropic contribution decreased. A small negative heat capacity change with significant enthalpy–entropy compensation established the involvement of multiple weak noncovalent interactions in the binding process. The 9‐ω‐amino hexyl ether analog stabilized the G‐quadruplex structure better than the 13‐phenyl alkyl analog. Copyright © 2015 John Wiley & Sons, Ltd.     

Targeting human telomeric G‐quadruplex DNA with antitumour natural alkaloid aristololactam‐β‐D‐glucoside and its comparison with daunomycin

Study on anticancer agents that act via stabilization of telomeric G‐quadruplex DNA has emerged as novel and exciting field for anticancer drug discovery. The interaction of carbohydrate containing anticancer alkaloid aristololactam‐β‐D‐glucoside (ADG) with human telomeric G‐quadruplex DNA sequence was characterized by different biophysical techniques. The binding parameters were compared with daunomycin (DAN), a well‐known chemotherapeutic drug. The Scatchard binding isotherms revealed noncooperative binding for both with the binding affinity values of (1.01 ± 0.05) × 10⁶ and (1.78 ± 0.18) × 10⁶ M⁻¹ for ADG and DAN, respectively. Circular dichroism, ferrocyanide quenching study, anisotropy study, thiazole orange displacement, optical melting, differential scanning calorimetry study, and molecular docking study suggest significant stacking and stabilizing efficiency of ADG with comparison to DAN. The energetics of the interaction for ADG and DAN revealed that both reactions were predominantly entropy driven. Negative heat capacity values were obtained from the temperature dependence of the enthalpy change. The standard molar Gibbs energy change exhibited only marginal alterations with temperature suggesting the occurrence of enthalpy‐entropy compensation. These findings indicate that ADG can act as a stabilizer of telomeric G‐quadruplex DNA and thereby can be considered as a potential telomerase inhibitor.     

Folding thermodynamics of the hybrid‐1 type intramolecular human telomeric G‐quadruplex

Guanine‐rich DNA sequences that may form G‐quadruplexes are located in strategic DNA loci with the ability to regulate biological events. G‐quadruplexes have been under intensive scrutiny owing to their potential to serve as novel drug targets in emerging anticancer strategies. Thermodynamic characterization of G‐quadruplexes is an important and necessary step in developing predictive algorithms for evaluating the conformational preferences of G‐rich sequences in the presence or the absence of their complementary C‐rich strands. We use a combination of spectroscopic, calorimetric, and volumetric techniques to characterize the folding/unfolding transitions of the 26‐meric human telomeric sequence d[A₃G₃(T₂AG₃)₃A₂]. In the presence of K⁺ ions, the latter adopts the hybrid‐1 G‐quadruplex conformation, a tightly packed structure with an unusually small number of solvent‐exposed atomic groups. The K⁺‐induced folding of the G‐quadruplex at room temperature is a slow process that involves significant accumulation of an intermediate at the early stages of the transition. The G‐quadruplex state of the oligomeric sequence is characterized by a larger volume and compressibility and a smaller expansibility than the coil state. These results are in qualitative agreement with each other all suggesting significant dehydration to accompany the G‐quadruplex formation. Based on our volume data, 432 ± 19 water molecules become released to the bulk upon the G‐quadruplex formation. This large number is consistent with a picture in which DNA dehydration is not limited to water molecules in direct contact with the regions that become buried but involves a general decrease in solute–solvent interactions all over the surface of the folded structure. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 216–227, 2014.     

G‐quadruplex DNA‐based asymmetric catalysis of michael addition: Effects of sonication,ligands, and co‐solvents

There is an escalating interest of using double stranded DNA molecules as a chiral scaffold to construct metal‐biomacromolecule hybrid catalysts for asymmetric synthesis. Several recent studies also evaluated the use of G‐quadruplex DNA‐based catalysts for asymmetric Diels‐Alder and Friedel‐Crafts reactions. However, there is still a lack of understanding of how different oligonucleotides, salts (such as NaCl and KCl), metal ligands and co‐solvents affect the catalytic performance of quadruplex DNA‐based hybrid catalysts. In this study, we aim to systematically evaluate these key factors in asymmetric Michael addition reactions, and to examine the conformational and molecular changes of DNA by circular dichroism (CD) spectroscopy and gel electrophoresis. We achieved up to 95% yield and 50% enantiomeric excess (ee) when the reaction of 2‐acylimidazole 1a and dimethylmalonate was catalyzed by 5′‐G₃(TTAG₃)₃?3′ (G4DNA1) in 20 mM MOPS (pH 6.5) containing 50 mM KCl and 40 µM [Cu(dmbipy)(NO₃)₂], and G4DNA1 was pre‐sonicated in ice bath for 10 min prior to the reaction. G‐quadruplex‐based hybrid catalysts provide a new tool for asymmetric catalysis, but future mechanistic studies should be sought to further improve the catalytic efficiency. The current work presents a systematic study of asymmetric Michael addition catalyzed by G‐quadruplex catalysts constructed via non‐covalent complexing, and an intriguing finding of the effect of pre‐sonication on catalytic efficiency. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:891–898, 2016     

pH‐dependent conformational dynamics of beta‐secretase 1: A molecular dynamics study

Beta‐secretase 1 (BACE‐1) is an aspartyl protease implicated in the overproduction of β‐amyloid fibrils responsible for Alzheimer disease. The process of β‐amyloid genesis is known to be pH dependent, with an activity peak between solution pH of 3.5 and 5.5. We have studied the pH‐dependent dynamics of BACE‐1 to better understand the pH dependent mechanism. We have implemented support for graphics processor unit (GPU) accelerated constant pH molecular dynamics within the AMBER molecular dynamics software package and employed this to determine the relative population of different aspartyl dyad protonation states in the pH range of greatest β‐amyloid production, followed by conventional molecular dynamics to explore the differences among the various aspartyl dyad protonation states. We observed a difference in dynamics between double‐protonated, mono‐protonated, and double‐deprotonated states over the known pH range of higher activity. These differences include Tyr 71‐aspartyl dyad proximity and active water lifetime. This work indicates that Tyr 71 stabilizes catalytic water in the aspartyl dyad active site, enabling BACE‐1 activity.     

Cholesterol‐β1AR interaction versus cholesterol‐β2AR interaction

Two 8‐µs all‐atom molecular dynamics simulations have been performed on the two highly homologous G protein‐coupled receptor (GPCR) subtypes, β₁‐ and β₂‐adrenergic receptors, which were embedded in a lipid bilayer with randomly dispersed cholesterol molecules. During the simulations, cholesterol molecules accumulate to different surface regions of the two receptors, suggesting the subtype specificity of cholesterol–β‐adrenergic receptor interaction and providing some clues to the physiological difference of the two subtypes. Meanwhile, comparison between the two receptors in interacting with cholesterols shed some new light on general determinants of cholesterol binding to GPCRs. Our results indicate that although the concave surface, charged residues and aromatic residues are important, neither of these stabilizing factors is indispensable for a cholesterol interaction site. Different combinations of these factors lead to the diversified binding modes of cholesterol binding to the receptors. Our long‐time simulations, for the first time, revealed the pathway of a cholesterol molecule entering the consensus cholesterol motif (CCM) site, and the binding process of cholesterol to CCM is accompanied by a side chain flipping of the conserved Trp4.50. Moreover, the simulation results suggest that the I‐/V‐/L‐rich region on the extracellular parts of helix 6 might be an alternatively conserved cholesterol‐binding site for the class‐A GPCRs. Proteins 2014; 82:760–770. © 2013 Wiley Periodicals, Inc.     

G‐LoSA: An efficient computational tool for local structure‐centric biological studies and drug design

Molecular recognition by protein mostly occurs in a local region on the protein surface. Thus, an efficient computational method for accurate characterization of protein local structural conservation is necessary to better understand biology and drug design. We present a novel local structure alignment tool, G‐LoSA. G‐LoSA aligns protein local structures in a sequence order independent way and provides a GA‐score, a chemical feature‐based and size‐independent structure similarity score. Our benchmark validation shows the robust performance of G‐LoSA to the local structures of diverse sizes and characteristics, demonstrating its universal applicability to local structure‐centric comparative biology studies. In particular, G‐LoSA is highly effective in detecting conserved local regions on the entire surface of a given protein. In addition, the applications of G‐LoSA to identifying template ligands and predicting ligand and protein binding sites illustrate its strong potential for computer‐aided drug design. We hope that G‐LoSA can be a useful computational method for exploring interesting biological problems through large‐scale comparison of protein local structures and facilitating drug discovery research and development. G‐LoSA is freely available to academic users at http://im.compbio.ku.edu/GLoSA/ .     

The stability of the G‐matrix: The role of spatial heterogeneity

The temporal stability of the genetic variance‐covariance matrix ( G ) has been discussed for a long time in the evolutionary literature. A common assumption in all studies, including empirical ones, is that spatial heterogeneity is minor such that the population can be represented by a single mean and variance. We use the well‐established allocation‐acquisition model to analyze the effect of relaxing of this assumption, simulating a case where the population is divided into patches with a variance in quality between patches. This variance can in turn differ between years. We found that changes in spatial variance in patch quality over years can make the G ‐matrix vary substantially over years and that the estimated genetic correlations, evolvability, and response to selection are different dependent on whether spatial heterogeneity is taken into account or not. This will have profound implications for our ability to predict evolutionary change and understanding of the evolutionary process.     

Introducing folding stability into the score function for computational design of RNA‐binding peptides boosts the probability of success

A computational strategy that integrates our peptide search algorithm with atomistic molecular dynamics simulation was used to design rational peptide drugs that recognize and bind to the anticodon stem and loop domain (ASL^Lys3) of human for the purpose of interrupting HIV replication. The score function of the search algorithm was improved by adding a peptide stability term weighted by an adjustable factor λ to the peptide binding free energy. The five best peptide sequences associated with five different values of λ were determined using the search algorithm and then input in atomistic simulations to examine the stability of the peptides' folded conformations and their ability to bind to ASL^Lys3. Simulation results demonstrated that setting an intermediate value of λ achieves a good balance between optimizing the peptide's binding ability and stabilizing its folded conformation during the sequence evolution process, and hence leads to optimal binding to the target ASL^Lys3. Thus, addition of a peptide stability term significantly improves the success rate for our peptide design search. Proteins 2016; 84:700–711. © 2016 Wiley Periodicals, Inc.     

Predicting binding modes of reversible peptide‐based inhibitors of falcipain‐2 consistent with structure–activity relationships

Falcipain‐2 (FP‐2) is a major hemoglobinase of Plasmodium falciparum, considered an important drug target for the development of antimalarials. A previous study reported a novel series of 20 reversible peptide‐based inhibitors of FP‐2. However, the lack of tridimensional structures of the complexes hinders further optimization strategies to enhance the inhibitory activity of the compounds. Here we report the prediction of the binding modes of the aforementioned inhibitors to FP‐2. A computational approach combining previous knowledge on the determinants of binding to the enzyme, docking, and postdocking refinement steps, is employed. The latter steps comprise molecular dynamics simulations and free energy calculations. Remarkably, this approach leads to the identification of near‐native ligand conformations when applied to a validation set of protein‐ligand structures. Overall, we proposed substrate‐like binding modes of the studied compounds fulfilling the structural requirements for FP‐2 binding and yielding free energy values that correlated well with the experimental data. Proteins 2017; 85:1666–1683. © 2017 Wiley Periodicals, Inc.     

A tale of two paralogs: human Transformer2 proteins with differential RNA-binding affinities

The Transformer2 (Tra2) proteins in humans are homologues of the Drosophila Tra2 protein. One of the two RNA-binding paralogs, Tra2β, has been very well-studied over the past decade, but not much is known about Tra2α. It was very recently shown that the two proteins demonstrate the phenomenon of paralog compensation. Here, we provide a structural basis for this genetic backup circuit, using molecular modelling and dynamics studies. We show that the two proteins display similar binding specificities, but differential affinities to a short GAA-rich RNA stretch. Starting from the 6-nucleotide RNA in the solution structure, close to 4000 virtual mutations were modelled on RNA and the domain–RNA interactions were studied after energy minimisation to convergence. Separately, another known 13-nucleotide stretch was docked and the domain–RNA interactions were observed through a 100-ns dynamics trajectory. We have also demonstrated the ‘compensatory’ mechanism at the level of domains in one of the domain repeat-containing RNA-binding proteins.     

Dynamic behavior of the post‐SET loop region of NSD1: Implications for histone binding and drug development

NSD1 is a SET‐domain histone methyltransferase that methylates lysine 36 of histone 3. In the crystal structure of NSD1, the post‐SET loop is in an autoinhibitory position that blocks binding of the histone peptide as well as the entrance to the lysine‐binding channel. The conformational dynamics preceding histone binding and the mechanism by which the post‐SET loop moves to accommodate the target lysine is currently unknown, although potential models have been proposed. Using molecular dynamics simulations, we have identified potential conformations of the post‐SET loop differing from those of previous studies, as well as proposed a model of peptide‐bound NSD1. Our simulations illustrate the dynamic behavior of the post‐SET loop and the presence of a few distinct conformations. In every case, the post‐SET loop remains in an autoinhibitory position blocking the peptide‐binding cleft, suggesting that another interaction is required to optimally position NSD1 in an active conformation. This finding provides initial evidence for a mechanism by which NSD1 preferentially binds nucleosomal substrates.     

Binding characteristics and interactive region of 2‐phenylpyrazolo[1,5‐c]quinazoline with DNA

The interaction between 2‐phenylpyrazolo[1,5‐c]quinazoline (PQ) and DNA under physiological conditions was investigated using multi‐spectroscopic techniques, atomic force microscopy and gel electrophoresis. The thermodynamic parameters were estimated and were discussed in detail. The results of fluorescence‐quenching experiments indicated that the main interactive force between PQ and DNA was a hydrophobic interaction and that it was a static quenching process. Potassium iodide and single‐strand (ss)DNA quenching studies, together with circular dichroism spectra implied groove binding of PQ with DNA. Atomic force microscopy and gel electrophoresis experiments suggested that there were no major conformational changes in DNA upon interaction with PQ. In addition, UV/vis absorption titration of DNA bases confirmed that PQ bound with DNA mainly through a minor groove interaction and preferentially interacted with adenine and thymine. We anticipate that this work will provide useful information for the application of quinazoline derivatives in the fields of medicinal and pharmaceutical chemistry. Copyright © 2014 John Wiley & Sons, Ltd.     

Molecular simulations study of novel 1,4‐dihydropyridines derivatives with a high selectivity for Cav3.1 calcium channel

1,4‐Dihydropyridines (DHPs) have been developed to treat hypertension, angina, and nerve system disease. They are thought to mainly target the L‐type calcium channels, but low selectivity prompts them to block Cav1.2 and Cav3.1 channels simultaneously. Recently, some novel DHPs with different hydrophobic groups have been synthesized and among them M12 has a higher selectivity for Cav3.1. However, the structural information about Cav3.1‐DHPs complexes is not available in the experiment. Thus, we combined homology modeling, molecular docking, molecular dynamics simulations, and binding free energy calculations to quantitatively elucidate the inhibition mechanism of DHPs. The calculated results indicate that our model is in excellent agreement with experimental results. On the basis of conformational analysis, we identify the main interactions between DHPs and calcium channels and further elaborate on the different selectivity of ligands from the micro perspective. In conjunction with energy distribution, we propose that the binding sites of Cav3.1‐DHPs is characterized by several interspersed hydrophobic amino acid residues on the IIIS6 and IVS6 segments. We also speculate the favorable function groups on prospective DHPs. Besides, our model provides important information for further mutagenesis experiments.     

Flexible functional interactions between G‐protein subunits contribute to the specificity of plant responses

Plants being sessile integrate information from a variety of endogenous and external cues simultaneously to optimize growth and development. This necessitates the signaling networks in plants to be highly dynamic and flexible. One such network involves heterotrimeric G‐proteins comprised of Gα, Gβ, and Gγ subunits, which influence many aspects of growth, development, and stress response pathways. In plants such as Arabidopsis, a relatively simple repertoire of G‐proteins comprised of one canonical and three extra‐large Gα, one Gβ and three Gγ subunits exists. Because the Gβ and Gγ proteins form obligate dimers, the phenotypes of plants lacking the sole Gβ or all Gγ genes are similar, as expected. However, Gα proteins can exist either as monomers or in a complex with Gβγ, and the details of combinatorial genetic and physiological interactions of different Gα proteins with the sole Gβ remain unexplored. To evaluate such flexible, signal‐dependent interactions and their contribution toward eliciting a specific response, we have generated Arabidopsis mutants lacking specific combinations of Gα and Gβ genes, performed extensive phenotypic analysis, and evaluated the results in the context of subunit usage and interaction specificity. Our data show that multiple mechanistic modes, and in some cases complex epistatic relationships, exist depending on the signal‐dependent interactions between the Gα and Gβ proteins. This suggests that, despite their limited numbers, the inherent flexibility of plant G‐protein networks provides for the adaptability needed to survive under continuously changing environments.     

Binding hotspots on K‐ras: Consensus ligand binding sites and other reactive regions from probe‐based molecular dynamics analysis

We have used probe‐based molecular dynamics (pMD) simulations to search for interaction hotspots on the surface of the therapeutically highly relevant oncogenic K‐Ras G12D. Combining the probe‐based query with an ensemble‐based pocket identification scheme and an analysis of existing Ras‐ligand complexes, we show that (i) pMD is a robust and cost‐effective strategy for binding site identification, (ii) all four of the previously reported ligand binding sites are suitable for structure‐based ligand design, and (iii) in some cases probe binding and expanded sampling of configurational space enable pocket expansion and increase the likelihood of site identification. Furthermore, by comparing the distribution of hotspots in nonpocket‐like regions with known protein‐ and membrane‐interacting interfaces, we propose that pMD has the potential to predict surface patches responsible for protein‐biomolecule interactions. These observations have important implications for future drug design efforts and will facilitate the search for potential interfaces responsible for the proposed transient oligomerization or interaction of Ras with other biomolecules in the cellular milieu. Proteins 2015; 83:898–909. © 2015 Wiley Periodicals, Inc.