PubMed Health⌕ Search

Biomedical subjects

M S Searle

Publications and source records attributed to M S Searle.

At least 19 recordsLinked to original sources

Loss of ubiquitin binding is a unifying mechanism by which mutations of SQSTM1 cause Paget's disease of bone.

Ubiquitin-associated (UBA) domain mutations of SQSTM1 are an important cause of Paget's disease of bone (PDB), which is a human skeletal disorder characterized by abnormal bone turnover. We previously showed that, when introduced into the full-length SQSTM1 protein, the disease-causing P392L, M404V, G411S, and G425R missense mutations and the E396X truncating mutation (representative of all of the SQSTM1 truncating mutations) cause a generalized loss of monoubiquitin binding and impaired K48-linked polyubiquitin binding at physiological temperature. Here, we show that the remaining three known PDB missense mutations, P387L, S399P, and M404T, have similar deleterious effects on monoubiquitin binding and K48-linked polyubiquitin binding by SQSTM1. The P387L mutation affects an apparently unstructured region at the N terminus of the UBA domain, some five residues from the start of the first helix, which is dispensable for polyubiquitin binding by the isolated UBA domain. Our findings support the proposal that the disease mechanism in PDB with SQSTM1 mutations involves a common loss of ubiquitin binding function of SQSTM1 and implicate a sequence extrinsic to the compact globular region of the UBA domain as a critical determinant of ubiquitin recognition by the full-length SQSTM1 protein.

Adaptor Proteins, Signal Transducing↗

p62 mutations, ubiquitin recognition and Paget's disease of bone.

Functional analyses of PDB (Paget's disease of bone)-associated mutants of the p62 [also known as SQSTM1 (sequestosome 1)] signalling adaptor protein represent an interesting paradigm for understanding not only the disease mechanism in this skeletal disorder, but also the critical determinants of ubiquitin recognition by an ubiquitin-binding protein. The 11 separate PDB mutations identified to date all affect the C-terminal region of p62 containing the UBA domain (ubiquitin-associated domain), a ubiquitin-binding element. All of these mutations have deleterious effects on ubiquitin binding by p62 in vitro, and there is evidence of an inverse relationship between ubiquitin-binding function and disease severity. The effects on ubiquitin-binding function of most of the mutations can be attributed to either reduced UBA domain stability, and/or the mutations affecting the presumed ubiquitin-binding interface of the UBA domain. However, a subset of the mutations are more difficult to rationalize; several of these affect sequences of p62 outside of the minimal ubiquitin-binding region, providing insights into non-UBA domain sequences within the host protein which mediate ubiquitin-binding affinity. The p62 mutations are presumed to result in activation of (osteoclast) NF-kappaB (nuclear factor kappaB) signalling. Understanding how loss of ubiquitin-binding function of p62 impacts on signal transduction events in osteoclasts will undoubtedly further our understanding of the disease mechanism in PDB at the molecular level.

Adaptor Proteins, Signal Transducing↗

Structural and functional studies of mutations affecting the UBA domain of SQSTM1 (p62) which cause Paget's disease of bone.

Mutations affecting the UBA (ubiquitin-associated) domain of SQSTM1 (Sequestosome 1) (p62) are a common cause of Paget's disease of bone. The missense mutations resolve into those which retain [P392L (Pro(392)-->Leu), G411S] or abolish (M404V, G425R) the ability of the isolated UBA domain to bind Lys-48-linked polyubiquitin. These effects can be rationalized with reference to the solution structure of the UBA domain, which we have determined by NMR spectroscopy. The UBA domain forms a characteristic compact three-helix bundle, with a hydrophobic patch equivalent to that previously implicated in ubiquitin binding by other UBA domains. None of the mutations affect overall folding of the UBA domain, but both M404V and G425R involve residues in the hydrophobic patch, whereas Pro-392 and Gly-411 are more remote. A simple model assuming the isolated UBA domain is functioning as a compact monomer can explain the effects of the mutations on polyubiquitin binding. The P392L and G411S mutations do however have subtle local effects on secondary structure, which may become more relevant in full-length SQSTM1. Identification of the in vivo ubiquitylated substrates of SQSTM1 will be most informative in determining the functional significance of the SQSTM1-ubiquitin interaction, and consequences of the disease-associated mutations.

Adaptor Proteins, Signal Transducing↗

Cooperativity in drug-DNA recognition: a molecular dynamics study.

NMR studies have shown that the minor groove-binding ligand Hoechst 33258 binds to the two T4/A4 tracts within the duplex d(CTTTTCGAAAAG)2 in a highly cooperative manner, such that in titration experiments no intermediate 1:1 complex can be detected. The NMR-derived structures of the free DNA and the 2:1 complex have been obtained, but can shed little light on what the origins of this cooperativity may be. Here we present the results of a series of molecular dynamics simulations on the free DNA, the 1:1 complex, and the 2:1 complex, which have been designed to enable us to calculate thermodynamic parameters associated with the molecular recognition events. The results of the molecular dynamics studies confirm that structural factors alone cannot explain the cooperativity observed, indeed when enthalpic and hydration factors are looked at in isolation, the recognition process is predicted to be slightly anticooperative. However, when changes in configurational entropy are taken into account as well, the overall free energy differences are such that the calculated cooperativity is in good agreement with that observed experimentally. The results indicate the power of molecular dynamics methods to provide reasonable explanations for phenomena that are difficult to explain on the basis of static models alone, and provide a nice example of the concept of "allostery without conformational change".

Bisbenzimidazole↗

Effects of amino acid phi,psi propensities and secondary structure interactions in modulating H alpha chemical shifts in peptide and protein beta-sheet.

H alpha chemical shifts are often used as indicators of secondary structure formation in protein structural analysis and peptide folding studies. On the basis of NMR analysis of model beta-sheet and alpha-helical peptides, together with a statistical analysis of protein structures for which NMR data are available, we show that although the gross pattern of H alpha chemical shifts reflects backbone torsion angles, longer range effects from distant amino acids are the dominant factor determining experimental chemical shifts in beta-sheets of peptides and proteins. These show context-dependent variations that aid structural assignment and highlight anomalous shifts that may be of structural significance and provide insights into beta-sheet stability.

Amino Acid Sequence↗

Insights into the stability of native and partially folded states of ubiquitin: effects of cosolvents and denaturants on the thermodynamics of protein folding.

The thermodynamics of the native<-->A state and native<-->unfolded transitions for ubiquitin have been characterized in detail using the denaturants methanol and guanidinium chloride (Gdn.HCl) both separately and in combination. Gdn.HCl destabilizes the partially folded alcohol-induced A state such that the effects of alcoholic solvents on the native<-->unfolded transition can be investigated directly via a two-state model. The combined denaturing effects of methanol and Gdn.HCl appear to conform to a simple additive model. We show that ubiquitin folds and unfolds cooperatively in all cases, forming the same "native" state; however, the thermodynamics of the N<-->U transition change dramatically between alcoholic and Gdn.HCl solutions, with folding in aqueous methanol associated with large negative enthalpy and entropy terms at 298 K with a gradual falloff in DeltaC(p) at higher methanol concentrations, as previously reported for the N<-->A transition (Woolfson, D. N., Cooper, A., Harding, M. M., Williams, D. H., and Evans, P. A. (1993) J. Mol. Biol. 229, 502-511.). Both the N<-->U and the N<-->A transitions are enthalpy driven to a similar extent. We conclude that under these conditions van der Waals interactions in the packing of the nonpolar protein core, which is common to both the N<-->U and the N<-->A transitions, appear to drive folding in the absence of entropic effects associated with release of ordered solvent (hydrophobic effect). Solvent transfer studies of hydrocarbons into alcoholic solvents, with and without Gdn.HCl, are consistent with a large enthalpic driving force for burial of a nonpolar surface, with a linear dependence of protein stability (DeltaG(N)(<-->)(U)) on cosolvent concentration reflected in a similar linear dependence of hydrocarbon solubility. The data demonstrate that the hydrophobic effect is not a prerequisite for specific stabilization of the native state or the A state and that van der Waals packing of the nonpolar core appears to be the dominant factor in stabilization of the native state.

Animals↗

Cooperative assembly of a nativelike ubiquitin structure through peptide fragment complexation: energetics of peptide association and folding.

Peptide fragments corresponding to the N- and C-terminal portions of bovine ubiquitin, U(1-35) and U(36-76), are shown by NMR to associate in solution to form a complex of modest stability (Kassn approximately 1.4 x 10(5) M(-1) at pH 7.0), with NMR features characteristic of a nativelike structure. The complex undergoes cold denaturation, with temperature-dependent estimates of stability from NMR indicating a DeltaC(p) degrees for fragment complexation in good agreement with that determined for native ubiquitin, suggesting that fragment association results in the burial of a similar hydrophobic surface area. The stability of the complex shows appreciable pH dependence, suggesting that ionic interactions on the surface of the protein contribute significantly. However, denaturation studies of native ubiquitin in the presence of guanidine hydrochloride (Gdn.HCl) show little pH dependence, suggesting that ionic interactions may be "screened" by the denaturant, as recently suggested. Examination of the conformation of the isolated peptide fragments has shown evidence for a low population of nativelike structure in the N-terminal beta-hairpin (residues 1-17) and weak nascent helical propensity in the helical fragment (residues 21-35). In contrast, the C-terminal peptide (36-76) shows evidence in aqueous solution, from some Halpha chemical shifts, for nonnative phi and psi angles; nonnative alpha-helical structure is readily induced in the presence of organic cosolvents, indicating that tertiary interactions in both native ubiquitin and the folded fragment complex strongly dictate its structural preference. The data suggest that the N-terminal fragment (1-35), where interaction between the helix and hairpin requires the minimum loss of conformational entropy, may provide the nucleation site for fragment complexation.

Animals↗

Templating peptide folding on the surface of a micelle: nucleating the formation of a beta-hairpin.

NMR spectroscopy is used to show that a 20-residue beta-hairpin peptide sequence derived from ferredoxin I, with a Pro-Asp two-residue type I turn which is uncommon in beta-hairpins, is unstructured in aqueous solution but shows NOE evidence for partial folding in the presence of sodium dodecylsulphate micelles. The peptide has a number of lysine residues in the N-terminal beta-strand capable of interacting with the micelle surface and templating the partial folding of the hairpin by reducing the entropic cost of ordering the peptide backbone.

Amino Acid Motifs↗

Sequence-dependent variation in DNA minor groove width dictates orientational preference of Hoechst 33258 in A-tract recognition: solution NMR structure of the 2:1 complex with d(CTTTTGCAAAAG)(2).

The solution structure of the dodecamer duplex d(CTTTTGCAAAAG)(2)and its 2:1 complex with the bis -benzimidazole Hoechst 33258 has been investigated by NMR and NOE-restrained molecular dynamics (rMD) simulations. Drug molecules are bound in each of the two A-tracts with the bulky N-methylpiperazine ring of each drug located close to the central TG (CA) step, binding essentially to the narrow minor groove of each A-tract. MD simulations over 1 ns, using an explicit solvation model, reveal time-averaged sequence-dependent narrowing of the minor groove from the 3'-end towards the 5'-end of each TTTT sequence. Distinct junctions at the TpG (CpA) steps, characterised by large positive roll, low helical and propeller twists and rapid AT base pair opening rates, add to the widening of the groove at these sites and appear to account for the bound orientation of the two drug molecules with the N-methylpiperazine ring binding in the wider part of the groove close to the junctions. Comparisons between the free DNA structure and the 2:1 complex (heavy atom RMSD 1.55 A) reveal that these sequence-dependent features persist in both structures. NMR studies of the sequence d(GAAAAGCTTTTC)(2), in which the A-tracts have been inverted with the elimination of the TpG junctions, results in loss of orientational specificity of Hoechst 33258 and formation of multiple bound species in solution, consistent with the drug binding in a number of different orientations.

Allosteric Regulation↗

Folding of a beta-hairpin peptide derived from the N-terminus of ubiquitin. Conformational preferences of beta-turn residues dictate non-native beta-strand interactions.

The role of the non-native beta-turn sequence (NPDG) in nucleating the folding of a beta-hairpin peptide derived from the N-terminus of ubiquitin, has been examined by NMR and CD spectroscopy. The NPDG sequence, while representing a common two-residue type I turn sequence in proteins, folds to give a G1-bulged type I turn in the context of a beta-hairpin peptide, to the exclusion of other possible conformations. The turn conformation results in misalignment of the two beta strands and a beta hairpin with non-native side chain interactions. A truncated 12-residue analogue of the hairpin, in which the majority of residues in the N-terminal beta strand have been deleted, shows some weak propensity to fold into a G-bulged type I turn conformation in the absence of interstrand stabilizing interactions. The NPDG turn sequence pays some of the entropic cost in initiating folding allowing interstrand interactions, which in this case arise from the non-native pairing of residue side chains, to stabilize a significant population of the folded state. Examination of the relative abundance of the Pro-Asp type I turn, with G in the +B1 position, vs. the type I G-bulged turn PXG, in a database of high resolution structures, reveals 48 instances of PXG bulged turns for which X = Asp is by far the most common residue with 20 occurrences. Strikingly, there are no examples of a type I PD turn with G at the +B1 position, in good agreement with our experimental observations that the PDG G-bulged turn is populated preferentially in solution.

Amino Acid Sequence↗

Dissecting the stability of a beta-hairpin peptide that folds in water: NMR and molecular dynamics analysis of the beta-turn and beta-strand contributions to folding.

NMR studies of the folding and conformational properties of a beta-hairpin peptide, several peptide fragments of the hairpin, and sequence-modified analogues, have enabled the various contributions to beta-hairpin stability in water to be dissected. Temperature and pH-induced unfolding studies indicate that the folding-unfolding equilibrium approximates to a two-state model. The hairpin is highly resistant to denaturation and is still significantly folded in 7 M urea at 298 K. Thermodynamic analysis shows the hairpin to fold in water with a significant change in heat capacity, however, DeltaCp degrees in 7 M urea is reduced. V/Y-->A mutations on one strand of the hairpin reduce folding to <10 %, consistent with a hydrophobic stabilisation model. We show that in a truncated peptide (residues 6-16) lacking the hydrophobic residues on one beta-strand, the type I' Asn-Gly turn in the sequence SINGKK is significantly populated in water in the absence of interstrand hydrophobic contacts. Unrestrained molecular dynamics simulations of unfolding, using an explicit solvation model, show that the conformation of the NG turn persists for longer than the AG analogue, which has a much lower propensity for type I' turn formation from a data base analysis of preferred turns. The origin of the high stability of the Asn-Gly turn is not entirely clear; data base analysis of 66 NG turns, together with molecular dynamics simulations, reveals no participation of the Asn side-chain in turn-stabilising interactions with the peptide backbone. However, hydration analysis of the molecular dynamics simulations reveals a pocket of "high density" water bridging between the Asn side-chain and peptide main-chain that suggests solvent-mediated interactions may play an important role in modulating phi,psi propensities in the NG turn region.

Amino Acid Sequence↗

Solution structure and dynamics of the A-T tract DNA decamer duplex d(GGTAATTACC)2: implications for recognition by minor groove binding drugs.

The structure of the DNA decamer duplex d(GGTAATTACC)(2) has been determined using NMR distance restraints and molecular dynamics simulations of 500 ps to 1 ns in aqueous solution at 300 K. Using both canonical A and canonical B starting structures [root-mean-square deviation (RMSD) 4.6 A; 1 A=10(-10) m], with and without experimental restraints, we show that all four simulations converge to a similar envelope of final conformations with B-like helical parameters (pairwise RMSD 1.27-2.03 A between time-averaged structures). While the two restrained simulations reach a stable trajectory after 300-400 ps, the unrestrained trajectories take longer to equilibrate. We have analysed the dynamic aspects of these structures (sugar pucker, helical twist, roll, propeller twist and groove width) and show that the minor groove width in the AATT core of the duplex fluctuates significantly, sampling both wide and narrow conformations. The structure does not have the highly pre-organized narrow minor groove generally regarded as essential for recognition and binding by small molecules, suggesting that ligand binding carries with it a significant component of 'induced-fit'. Our simulations show that there are significant differences in structure between the TpA step (where p=phosphate) and the ApA and ApT steps, where a large roll into the major groove at the TpA step appears to be an important factor in widening the minor groove at this position.

Allosteric Site↗

Structure, dynamics and hydration of the nogalamycin-d(ATGCAT)2Complex determined by NMR and molecular dynamics simulations in solution.

The structure of the 1:1 nogalamycin:d(ATGCAT)2 complex has been determined in solution from high-resolution NMR data and restrained molecular dynamics (rMD) simulations using an explicit solvation model. The antibiotic intercalates at the 5'-TpG step with the nogalose lying along the minor groove towards the centre of the duplex. Many drug-DNA nuclear Overhauser enhancements (NOEs) in the minor groove are indicative of hydrophobic interactions over the TGCA sequence. Steric occlusion prevents a second nogalamycin molecule from binding at the symmetry-related 5'-CpA site, leading to the conclusion that the observed binding orientation in this complex is the preferred orientation free of the complication of end-effects (drug molecules occupy terminal intercalation sites in all X-ray structures) or steric interactions between drug molecules (other NMR structures have two drug molecules bound in close proximity), as previously suggested. Fluctuations in key structural parameters such as rise, helical twist, slide, shift, buckle and sugar pucker have been examined from an analysis of the final 500 ps of a 1 ns rMD simulation, and reveal that many sequence-dependent structural features previously identified by comparison of different X-ray structures lie within the range of dynamic fluctuations observed in the MD simulations. Water density calculations on MD simulation data reveal a time-averaged pattern of hydration in both the major and minor groove, in good agreement with the extensive hydration observed in two related X-ray structures in which nogalamycin is bound at terminal 5'-TpG sites. However, the pattern of hydration determined from the sign and magnitude of NOE and ROE cross-peaks to water identified in 2D NOESY and ROESY experiments identifies only a few "bound" water molecules with long residence times. These solvate the charged bicycloaminoglucose sugar ring, suggesting an important role for water molecules in mediating drug-DNA electrostatic interactions within the major groove. The high density of water molecules found in the minor groove in X-ray structures and MD simulations is found to be associated with only weakly bound solvent in solution.

Base Sequence↗

Molecular recognition between a new pentacyclic acridinium salt and DNA sequences investigated by optical spectroscopic techniques, proton nuclear magnetic resonance spectroscopy, and molecular modeling.

A pentacyclic acridine, 1H-2,3-dihydroindolizino[7,6,5-kl]acridinium chloride (1), related in structure to tetra- and pentacyclic marine natural products, has previously been shown to induce apoptosis in breast and non-small-cell lung tumor cell lines and shows significant differences in biological potency and antitumor profile from other intercalating agents based on the acridine framework. We report on the molecular recognition of the acridinium salt with DNA, quantified by optical spectroscopic methods, and have compared these results with the clinical agent amsacrine (m-AMSA). The results point to an intercalative association between 1 and G-C-rich sequences of DNA. We have synthesized a hexamer duplex d(ACGCGT)2, presenting two potential 5'-CpG recipient sites, and have investigated in detail by NMR and molecular modeling methods the orientational preferences of 1, particularly with regard to the pyrrolidine ring system. On the basis of the intermolecular nuclear Overhauser effect (NOE) data, four possible intercalation models were considered; no single model produced a significantly better fit than any of the others. The best fit to the experimental data was obtained by considering a dynamic equilibrium between the different intercalated orientations with the drug maximizing pi-overlap with the G-C base pairs at the intercalation site. We found little evidence for any degree of groove specificity imparted by the pyrrolidine ring. If these simulations have biological relevance they suggest that, at most, the agent induces only a transitory hot spot in the DNA which, evidently, is sufficient to be sensed by damage-recognition mechanisms of the cell.

Acridines↗

DNA minor groove recognition by bis-benzimidazole analogues of Hoechst 33258: insights into structure-DNA affinity relationships assessed by fluorescence titration measurements.

Fluorescence titration measurements have been used to examine the binding interaction of a number of analogues of the bis -benzimidazole DNA minor groove binding agent Hoechst 33258 with the decamer duplex d(GCAAATTTGC)2. The method of continuous variation in ligand concentration (Job plot analysis) reveals a 1:1 binding stoichiometry for all four analogues; binding constants are independent of drug concentration (in the range [ligand] = 0.1-5 microM). The four analogues studied were chosen in order to gain some insight into the relative importance of a number of key structural features for minor groove recognition, namely (i) steric bulk of the N -methylpiperazine ring, (ii) ligand hydrophobicity, (iii) isohelicity with the DNA minor groove and (iv) net ligand charge. This was achieved, first, by replacing the bulky, non-planar N -methylpiperazine ring with a less bulky planar charged imidazole ring permitting binding to a narrower groove, secondly, by linking the N -methylpiperazine ring to the phenyl end of the molecule to give the molecule a more linear, less isohelical conformation and, finally, by introducing a charged imidazole ring in place of the phenolic OH making it dicationic, enabling the contribution of the additional electrostatic interaction and extended conformation to be assessed. Delta G values were measured at 20 degrees C in the range -47.6 to -37.5 kJ mol-1 and at a number of pH values between 5.0 and 7.2. We find a very poor correlation between Delta G values determined by fluorescence titration and effects of ligand binding on DNA melting temperatures, concluding that isothermal titration methods provide the most reliable method of determining binding affinities. Our results indicate that the bulky N -methylpiperazine ring imparts a large favourable binding interaction, despite its apparent requirement for a wider minor groove, which others have suggested arises in a large part from the hydrophobic effect. The binding constant appears to be insensitive to the isohelical arrangement of the constituent rings which in these analogues gives the same register of hydrogen bonding interactions with the floor of the groove.

Base Sequence↗

Structure/activity studies of the anti-MUC1 monoclonal antibody C595 and synthetic MUC1 mucin-core-related peptides and glycopeptides.

MUC1 mucin is a large complex glycoprotein expressed on normal epithelial cells in humans and overexpressed and under or aberrantly glycosylated on many malignant cancer cells which consequently allows recognition of the protein core by antibodies. In order to understand how glycosylation may modulate or regulate antibody binding of mucin protein core epitopes, we have analyzed the antibody C595 (epitope RPAP) for its structure, stability, and its binding to a series of synthetic peptides and glycopeptides by a number of spectroscopic methods. Thermal and pH denaturation studies followed by changes in the CD spectrum of the antibody indicate critical involvement of specific residues to the stability of the antibody. Fluorescence binding studies indicate that alpha-N-acetylgalactosamine (GalNAc) glycosylation of a MUC1 mucin synthetic peptide TAPPAHGVT9SAPDTRPAPGS20T21APPA at threonine residues 9 and 21 and serine residue 20 enhanced the binding of antibody. The structural effects of GalNAc glycosylation on the conformation of the MUC1 peptide were studied. CD of the peptides and glycopeptides in a cryogenic mixture cooled to approximately -97 degrees C revealed that a left-handed polyproline II helix (PPII) is adopted by the peptides in solution, which appears to be further stabilized by addition of the GalNAc residues. Consistent with the PPII helical structure, which has no intra-amide hydrogen bonds, high-field NMR spectroscopy of the glycopeptide revealed no sequential dNN, medium-range, or long-range nuclear Overhauser effect (NOE) connectivities. These studies indicate that stabilization of the PPII helix by GalNAc glycosylation present the epitope of C595 antibody with a favorable conformation for binding. Furthermore, they illustrate that glycosylation of the MUC1 tumor marker protein with a simple O-linked saccharide expressed in many cancers, can enhance the binding of the clinically relevant C595 antibody.

Amino Acid Sequence↗

Modulation of intrinsic phi,psi propensities of amino acids by neighbouring residues in the coil regions of protein structures: NMR analysis and dissection of a beta-hairpin peptide.

Analysis of residues in coil regions of protein structures presents a novel approach to deconvoluting the various competing factors which determine the intrinsic phi,psi propensities of amino acids free from the regular interactions associated with beta-strands and alpha-helices. We have considered the role of context on phi,psi preferences by examining the effects of neighbouring residues in modulating coil propensities within a data base of 512 high-resolution, low-homology structures. In the general case, when flanking residues are beta-branched or aromatic (Val, Ile, Tyr and Phe) the beta-propensity (Pbeta) increases significantly, largely due to steric effects between flanking residues. More subtle residue-specific effects are apparant when Pbeta values are examined in detail, showing "random coil" conformations to be highly sequence-dependent. The effects of flanking residues on phi distributions have been used to calculate context-dependent average 3JNH-Halpha coupling constants. We have examined these findings in the context of the folding of a model 16-residue beta-hairpin peptide, "mutant" hairpin (VSI-->KSK sequence change) and the isolated C-terminal beta-strand fragments of both hairpins. We find a better correlation between 3JNH-Halpha values derived from the data base model and those determined experimentally when context-dependent phi distributions are considered. The individual C-terminal beta-strand sequences (GKKITVSI versus GKKITKSK) of the two hairpins are predisposed to different extents to formation of an extended beta-like conformation. Conformational "predisposition" in this context may contribute significantly to beta-hairpin stability.

Amino Acid Sequence↗

DNA minor groove recognition by a tetrahydropyrimidinium analogue of hoechst 33258: NMR and molecular dynamics studies of the complex with d(GGTAATTACC)2.

Hoechst 43254 (H43254), a 2,3,4,5-tetrahydropyrimidin-1-ium analogue of the bis-benzimidazole minor groove binding agent Hoechst 33258 (H33258), has been studied by NMR and restrained molecular dynamics in its complex with d(GGTAATTACC)2. We investigate the origin of the enhanced complex stability afforded by the replacement of the N-methylpiperazine ring of H33258 with the tetrahydropyrimidinium ring of H43254, the latter presenting the opportunity for specific minor groove-directed recognition through a pyrimidinium NH. A set of 25 drug-DNA NOEs define the binding site with some precision and are used as part of the structural analysis using restrained molecular dynamics simulations considering explicit solvation and the treatment of electrostatic interactions using the particle mesh Ewald method within AMBER 4.1. Starting with three different initial structures with the drug located at different sites in the groove (pairwise RMSD 4.3-12.6 A) we arrive at three very similar structures (pairwise RMSD 0.80-1.34 A) representing one converged binding site at the centre of the AATT tract. Two of the three structures show the tetrahydropyrimidinium ring to be suitably positioned for an -NH to adenine N3 hydrogen bond suggesting that electrostatic interactions may play an important role in the enhanced affinity as well as imparting additional A-T specificity. The NMR data show that the pyrimidinium NH interaction is dynamic since signal averaging from the two sides of the ring indicate rapid rotations in the bound form.

Base Sequence↗