Molecular modelling of the interaction of carbocyclic analogues of netropsin and distamycin with d(CGCGAATTCGCG)2

A molecular mechanics and molecular dynamics approach was used to examine the structure of complexes formed between the d(CGCGAATTCGCG)2 duplex and netropsin, distamycin, and four carbocyclic analogues of netropsin and distamycin (1-4). The resulting structures of the ligand-DNA model complexes and their energetics were examined. It is predicted that the compounds 1-4 should have a decreased affinity for the minor groove of AT-rich regions in comparison to netropsin and distamycin. From the energetic analysis it appears that van der Waals and electrostatic interactions are more important than specific hydrogen bonds in stabilizing the ligand-duplex complexes. We predict that compounds 1 and 2 are effectively isohelical with the DNA minor groove. The superior DNA-binding afforded by 1 and 2 in comparison to 3 and 4 results from their more effective penetration into the minor groove and smaller perturbation of molecular structure upon complex formation.     

Structural investigation of the hedamycin:d(ACCGGT)2 complex by NMR and restrained molecular dynamics

Hedamycin, a member of the pluramycin family of drugs, displays a range of biological responses including antitumor and antimicrobial activity. The mechanism of action is via direct interaction with DNA through intercalation between the bases of the oligonucleotide and alkylation of a guanine residue at 5'-PyG-3' sites. There appears to be some minor structural differences between two earlier studies on the interaction of hedamycin with 5'-PyG-3' sites. In this study, a high-resolution NMR analysis of the hedamycin:d(ACCGGT)2 complex was undertaken in order to investigate the effect of replacing the thymine with a guanine at the preferred 5'-CGT-3' site. The resultant structure was compared with earlier work, with particular emphasis placed on the drug conformation. The structure of the hedamycin:d(ACCGGT)2 complex has many features in common with the two previous NMR structures of hedamycin:DNA complexes but differed in the conformation and orientation of the N,N-dimethylvancosamine saccharide of hedamycin in one of these structures. The preferential binding of hedamycin to 5'-CG-3' over 5'-TG-3' binding sites is explained in terms of the orientation and location of the N,N-dimethylvancosamine saccharide in the minor groove.     

Interactions of quinoxaline antibiotic and DNA: the molecular structure of a triostin A-d(GCGTACGC) complex

The crystal structure of a DNA octamer d(GCGTACGC) complexed to an antitumor antibiotic, triostin A, has been solved and refined to 2.2 A resolution by x-ray diffraction analysis. The antibiotic molecule acts as a true bis intercalator surrounding the d(CpG) sequence at either end of the unwound right-handed DNA double helix. As previously observed in the structure of triostin A-d(CGTACG) complex (A.H.-J. Wang, et. al., Science, 225, 1115-1121 (1984)), the alanine amino acid residues of the drug molecule form sequence-specific hydrogen bonds to guanines in the minor groove. The two central A.T base pairs are in Hoogsteen configuration with adenine in the syn conformation. In addition, the two terminal G.C base pairs flanking the quinoxaline rings are also held together by Hoogsteen base pairing. This is the first observation in an oligonucleotide of. Hoogsteen G.C base pairs where the cytosine is protonated. The principal functional components of a bis-intercalative compound are discussed.     

Structural analysis of the complex of a distamycin analogue with the Dickerson dodecamer 13C labeled at 5'-carbons using NMR spectroscopy.

Structural analysis of the complex of a distamycin analogue (Tallimustine) with the Dickerson dodecamer d(C*G*C*G*A*A*T*T*C*G*C*G) [N*:[5'-(13)C]nucleotide] was performed by NMR spectroscopy and the results will be described in detail.     

The solution structure of the Ga(III)-bleomycin A2 complex resolved by NMR and molecular modeling; interaction with d(CCAGGCCTGG)

The solution structure of the Ga(III)-bleomycin A2 complex (GaBLM) has been determined using 2D NMR methods in combination with molecular dynamics calculations. Complete assignment of the amide and amine protons, observation of 80 NOEs and measurement of 15 (3)JH(-H) coupling constants provided us with a well-defined structure using a restrained simulated annealing protocol. On the basis of distance and dihedral angle constraints agreement, along with potential energy considerations, the favored model is a five-coordinate complex with the primary amine of beta-aminoalanine holding the axial position of a distorted tetragonal pyramid. The disaccharide moiety of GaBLM is not a ligand, sharing the same side of the equatorial plane with the axial amine ligand. Titration of the self-complementary oligonucleotide d(CCAGGCCTGG) with GaBLM results in the formation of only one 1:1 complex in slow exchange on the NMR time scale. Our data indicate that the bithiazole moiety intercalates between the C6*G15 and C7*G14 base pairs, in a similar mode to that reported by earlier studies. Structural implications and comparisons to other metallo-bleomycins are discussed.     

Interaction of berenil with the EcoRI dodecamer d(CGCGAATTCGCG)2 in solution studied by NMR

The conformation of the EcoRI dodecamer d(CGCGAATTCGCG)2 has been examined in solution by 1H and 31P NMR. Spin-spin coupling constants and nuclear Overhauser (NOE) enhancement spectroscopy show that all deoxyriboses lie in the south domain, with a small admixture of the north conformation (0-20%). The time dependence of the nuclear Overhauser enhancements also reveals a relatively uniform conformation at the glycosidic bonds (average angle, chi = -114 degrees). The average helical twist is 36.5 degrees (9.8 base pairs per turn). Tilt angles are small (in the range 0 to -10 degrees), and roll angles are poorly determined. Unlike single-crystal X-ray studies of the same sequence, there is no evidence for asymmetry in the structure. Both the NOE intensities and 31P relaxation data imply conformational anomalies at the C3-G4/C9-G10 and the A5-A6/T7-T8 steps. Berenil binds in 1:1 stoichiometry to the dodecamer with high affinity (Kd = 1 microM at 298 K) and causes substantial changes in chemical shifts of the sugar protons of nucleotides Ado 5-Cyt 9 and of the H2 resonances of the two Ado residues. No significant asymmetry appears to be induced in the DNA conformation on binding, and there is no evidence for intercalation, although the binding site is not centrosymmetric. NOEs are observed between the aromatic protons of berenil and the H1' of both Thy 7 and Thy 8, as well as to Ado 5 and Ado 6 H2. These results firmly establish that berenil binds via the minor groove and closely approaches the nucleotides Ado 6, Thy 7, and Thy 8. On the basis of quantitative NOE spectroscopy and measurements of spin-spin coupling constants, changes in the conformations of the nucleotides are found to be small. Using the observed NOEs between the ligand and the DNA together with the derived glycosidic torsion angles, we have built models that satisfy all of the available solution data. The berenil molecule binds at the 5'-AAT (identical to 5'-ATT on the complementary strand) site such that (i) favorable hydrogen bonds are formed between the charged amidinium groups and the N3 atoms of Ado 6 and Ado 18 and (ii) the ligand is closely isohelical with the floor of the minor groove.     

Conformational analysis of the nonapeptide leuprorelin using NMR and molecular modeling

The nonapeptide Leuprorelin, one of the LHRH agonists, was studied by means of 2D nuclear magnetic resonance spectroscopy and molecular modeling. NOESY spectra in aqueous/deuterated methanol solution (50%H₂O/CD₃OD) at low temperature (268 K) revealed short-range nOe connectivities (i, i+1), characteristic of flexibility of the molecule. The H^N–H^N sequential connectivities observed provide evidence that the sequence has the propensity to form a bend involving residues 5 and 6 and the N-terminal segment. The -proton chemical shifts compared to random coil and additional data from the amide proton temperature coefficients support this assumption. One long-range nOe cross peak between H₂ –H^NEth is indicative of proximity between C- and N-termini.     

Conformational analysis of the nonapeptide leuprorelin using NMR and molecular modeling

Summary The nonapeptide Leuprorelin, one of the LHRH agonists, was studied by means of 2D nuclear magnetic resonance spectroscopy and molecular modeling. NOESY spectra in aqueous/deuterated methanol solution (50% H₂O/CD₃OD) at low temperature (268 K) revealed short-range nOe connectivities (i, i+1), characteristic of flexibility of the molecule. The H ^N -H ^N sequential connectivities observed provide evidence that the sequence has the propensity to form a bend involving residues 5 and 6 and the N-terminal segment. The α-proton chemical shifts compared to random coil and additional data from the amide proton temperature coefficients support this assumption. One long-range nOe cross peak between H ₂ ^α -H ^NEth is indicative of proximity between C- and N-termini.     

The molecular structure of the complex of Hoechst 33258 and the DNA dodecamer d(CGCGAATTCGCG).

The crystal structure of the complex between the dodecamer d(CGCGAATTCGCG) and a synthetic dye molecule Hoechst 33258 was solved by X-ray diffraction analysis and refined to an R-factor of 15.7% at 2.25 A resolution. The crescent-shaped Hoechst compound is found to bind to the central four AATT base pairs in the narrow minor groove of the B-DNA double helix. The piperazine ring of the drug has its flat face almost parallel to the aromatic bisbenzimidazole ring and lies sideways in the minor groove. No evidence of disordered structure of the drug is seen in the complex. The binding of Hoechst to DNA is stabilized by a combination of hydrogen bonding, van der Waals interaction and electrostatic interactions. The binding preference for AT base pairs by the drug is the result of the close contact between the Hoechst molecule and the C2 hydrogen atoms of adenine. The nature of these contacts precludes the binding of the drug to G-C base pairs due to the presence of N2 amino groups of guanines. The present crystal structural information agrees well with the data obtained from chemical footprinting experiments.     

Structural comparison of the B-DNA dodecamers d(CGCGTTAACGCG) and d(CGCGAATTCGCG) with T2A2 and A2T2 tracts.

The x-ray structure of the deoxy oligonucleotide dodecamer d(CGCGTTAACGCG) recently determined in our laboratory shows that the helical parameters of the central TTAA segment are significantly different compared to the central AATT in d(CGCGAATTCGCG). The roll in the central TA step of the T2A2 dodecamer opens towards the minor groove while the AT step of the A2T2 dodecamer opens towards the major groove. Also, the roll angles at the steps 4 and 8 (GT and AC in T2A2) and (GA and TC in A2T2) are in opposite directions. The high cup and helical twist angles at the central base-pair of T2A2 decreases the base stacking interactions compared to A2T2. Tilt angles within the tetranucleotide segments TTAA and AATT have opposite signs. In spite of the local differences caused by the sequence inversion (TTAA----AATT), the two dodecamers exhibit similar overall bending. The top third is more bent than the bottom third relative to the central segment. This asymmetric bending in the two dodecamers is mainly due to crystal packing interactions.     

Structural analysis of the complex of phenoxazone antibiotic actinocyl-bis-(2-deoxytetranucleotide 5'-d(TpGpCpA) by the methods of two dimensional 1H-NMR-spectroscopy and molecular mechanics

The spatial structures of intercalated complexes of synthetic phenoxazone antibiotic actinocyl-bis-(2-dimethylaminoethyl) amide with self-complementary deoxytetranucleotide 5'-d(TpGpCpA) have been investigated. Analysis has been made using two-dimensional NMR (2D-NOESY) data in aqueous solution and molecular mechanics simulation. Distinctive features of the conformation of drug-DNA complexes have been determined at two possible orientations of the chromophore of phenoxazone antibiotic at the intercalation site.     

Structural investigation of syringomycin-E using molecular dynamics simulation and NMR

Syringomycin-E (SR-E) is a cyclic lipodepsinonapeptide produced by certain strains of the bacterium Pseudomonas syringae pv. syringae. It shows inhibitory effects against many fungal species, including human pathogens. Its primary biological target is the plasma membrane, where it forms channels comprised of at least six SR-E molecules. The high-resolution structure of SR-E and the structure of the channels are currently not known. In this paper, we investigate in atomic detail the molecular features of SR-E in water by NMR and in water and octane by molecular dynamics simulation (MD). We built a model of the peptide and examined its structure in water and octane in 200 ns MD simulations both with and without distance restraints derived from NMR NOE data. The resulting trajectories show good agreement with the measured NOEs and circular dichroism data from the literature and provide atomistic models of SR-E that are an important step toward a better understanding of the antifungal and antibacterial activity of this peptide.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Stability and motion of a hairpin and the corresponding mismatched duplex: a theoretical exploration using molecular mechanics and normal mode analysis of 2D NMR results on d(GCCGCAGC)

The oligomer d(GCCGCAGC) can adopt two different conformations: i) a duplex with two mismatched A.C base pairs and ii) a hairpin with two C.G base pairs and a single stranded loop. We report molecular mechanics, normal mode analysis, and thermodynamic stability calculations for both structures. We show that the energy-minimized structure and harmonic-dynamics results are in complete agreement with the observed NOE spectrum and imino proton exchange data. We conclude that the high stability of the hairpin structure over the duplex at low salt concentration is due to the higher vibrational entropy contribution to the system free energy by the single stranded loop and to the lack of minor groove phosphate/phosphate electrostatic repulsions that tend to destabilize the duplex.     

Novel allosteric effectors of the tryptophan synthase alpha(2)beta(2) complex identified by computer-assisted molecular modeling

Tryptophan synthase is a pyridoxal 5'-phosphate-dependent alpha(2)beta(2) complex catalyzing the formation of L-tryptophan. The functional properties of one subunit are allosterically regulated by ligands of the other subunit. Molecules tailored for binding to the alpha-active site were designed using as a starting model the three-dimensional structure of the complex between the enzyme from Salmonella typhimurium and the substrate analog indole-3-propanol phosphate. On the basis of molecular dynamics simulations, indole-3-acetyl-X, where X is glycine, alanine, valine and aspartate, and a few other structurally related compounds were found to be good candidates for ligands of the alpha-subunit. The binding of the designed compounds to the alpha-active site was evaluated by measuring the inhibition of the alpha-reaction of the enzyme from Salmonella typhimurium. The inhibition constants were found to vary between 0.3 and 1.7 mM. These alpha-subunit ligands do not bind to the beta-subunit, as indicated by the absence of effects on the rate of the beta-reaction in the isolated beta(2) dimer. A small inhibitory effect on the activity of the alpha(2)beta(2) complex was caused by indole-3-acetyl-glycine and indole-3-acetyl-aspartate whereas a small stimulatory effect was caused by indole-3-acetamide. Furthermore, indole-3-acetyl-glycine, indole-3-acetyl-aspartate and indole-3-acetamide perturb the equilibrium of the catalytic intermediates formed at the beta-active site, stabilizing the alpha-aminoacrylate Schiff base. These results indicate that (i) indole-3-acetyl-glycine, indole-3-acetyl-aspartate and indole-3-acetamide bind to the alpha-subunit and act as allosteric effectors whereas indole-3-acetyl-valine and indole-3-acetyl-alanine only bind to the alpha-subunit, and (ii) the terminal phosphate present in the already known allosteric effectors of tryptophan synthase is not strictly required for the transmission of regulatory signals.     

2D NMR studies and 3D structure of the parallel-stranded duplex oligonucleotide Acrm5-alpha-d(TCTAAACTC)-beta-d(AGATTTGAG) via complete relaxation matrix analysis of the NOE effects and molecular mechanics calculations

The three-dimensional structure of the duplex formed by the association of the unnatural oligonucleotide alpha-d(TCTAAACTC) covalently linked to an acridine derivative (m5Acr) with its natural and parallel complementary sequence beta-d(AGATTTGAG) was investigated by nuclear magnetic resonance spectroscopy and constrained molecular mechanics calculations. All the nonexchangeable and exchangeable resonances were assigned in this duplex. The structure was refined by using interproton distances determined by NOE measurements. The NOE values were converted into distances by using the complete 190 x 190 relaxation matrix. The unnatural duplex Acrm5-alpha-d(TCTAAACTC)-beta-d(AGATTTGAG) forms a parallel right-handed helix with Watson-Crick base pairing; the alpha and beta deoxyriboses adopt a 3'-exo conformation. The acridine moiety was found stacked up the C9-G9 base pair. The structure of the first seven base pairs of this duplex was found similar to that of the duplex alpha-d(TCTAAAC)-beta-d(AGATTTG), which we had already investigated [Lancelot, G., et al. (1989) Biochemistry 28, 7871-7878]. Since these structures were generated by using experimental NOE values obtained independently on macromolecules whose global correlation time was different (3.8 and 2.2 ns), we conclude that this comparison is a good test of the viability of our method to generate three-dimensional structures of oligonucleotides in solution. Starting from different initial conformations, we show that the NOE constraints allow one to reach the same final restrained conformation, taking into account implicitly the solvent effect.     

Synthesis of 2-(2-adamantyl)piperidines and structure anti-influenza virus A activity relationship study using a combination of NMR spectroscopy and molecular modeling

The 2-(2-adamantyl)piperidines 13 and 15a-c were synthesized and evaluated for anti-influenza virus A and B activity. The parent N-H compound 13 was 3–4 times more active than amantadine and rimantadine against H₂N₂ influenza A. N-alkylation of 13 resulted in derivatives 15a-c that were devoid of biological activity. This dramatic reduction in biological activity may be attributed to the different conformational properties between N-H and N-alkyl piperidines, as deduced from the combination of computational chemistry and NMR spectroscopy.     

Deuterium NMR investigation of backbone dynamics in the synthetic oligonucleotide [d(CGCGAATTCGCG)]2

Backbone dynamics in the [5',5"-2H2]2'-deoxythymidine labeled duplex dodecamer [d-(CGCGAAT*T*CGC)]2 have been investigated by solid-state 2H NMR. Quadrupolar echo line shapes, spin-lattice relaxation, and quadrupolar echo decay times were obtained over hydration levels ranging from W = 0.0 to 25.2 (moles of H2O/mole of nucleotide). Variation of the line shape with changing hydration level was analyzed by using models employed in previous investigations of dodecamer base and sugar dynamics. Both fast local motions and a slower helix motion were present within the oligonucleotide. The fast motion was modeled as a four-site libration whose amplitude increased with hydration level. The root mean square amplitude of this librational model was 2-6 degrees larger than the amplitude observed in either the furanose ring or base labeled material for the entire range of hydration levels investigated. The observed line shape was inconsistent with a rapid three-site trans-gauche isomerization. A slow motion about the helix axis was observed at low water levels and increased in rate and amplitude with hydration. This motional model is in agreement with previous oligonucleotide studies.     

Automated screening for metabolites in complex mixtures using 2D COSY NMR spectroscopy

One of the greatest challenges in metabolomics is the rapid and unambiguous identification and quantification of metabolites in a biological sample. Although one-dimensional (1D) proton nuclear magnetic resonance (NMR) spectra can be acquired rapidly, they are complicated by severe peak overlap that can significantly hinder the automated identification and quantification of metabolites. Furthermore, it is currently not reasonable to assume that NMR spectra of pure metabolites are available a priori for every metabolite in a biological sample. In this paper we develop and report on tests of methods that assist in the automatic identification of metabolites using proton two-dimensional (2D) correlation spectroscopy (COSY) NMR. Given a database of 2D COSY spectra for the metabolites of interest, our methods provide a list sorted by a heuristic likelihood of the metabolites present in a sample that has been analyzed using 2D COSY NMR. Our models attempt to correct the displacement of the peaks that can occur from one sample to the next, due to pH, temperature and matrix effects, using a statistical and chemical model. The correction of one peak can result in an implied correction of others due to spin–spin coupling. Furthermore, these displacements are not independent: they depend on the relative position of functional groups in the molecule. We report experimental results using defined mixtures of amino acids as well as real complex biological samples that demonstrate that our methods can be very effective at automatically and rapidly identifying metabolites.     

Structural analysis of substance P using molecular dynamics and NMR spectroscopy.

The present work is a combined structural study, using Nuclear Magnetic Resonance (NMR) and Molecular Dynamics(MD), of the amidated and the free acid forms of substance P in water and methanol. The results obtained using both approaches were compared in order to characterize the structural features of both peptides in solution. From the NMR experiments it was derived that the free acid form adopts an extended conformation at the N-terminus and a helical conformation at the C-terminal segment of the peptide in both water and methanol; these structural features are in qualitative agreement with the results of the MD simulations. No significant differences in behavior were observed between the amidated and the free acid forms of the peptide in the simulations and in the experiments carried out in water, suggesting that the different activities of these analogs are due to their different mode of interaction with the receptor rather than to their structural preferences. Finally, we propose that the structure of substance P can be partially inferred from its sequence due to the presence of a Pro-X-Pro motif on the N-terminus and a Gly-Leu sequence on the C-terminus.