Distance Geometry-simulated Annealing Research Articles

Chromosome3D: reconstructing three-dimensional chromosomal structures from Hi-C interaction frequency data using distance geometry simulated annealing.

BackgroundReconstructing three-dimensional structures of chromosomes is useful for visualizing their shapes in a cell and interpreting their function. In this work, we reconstruct chromosomal structures from Hi-C data by translating contact counts in Hi-C data into Euclidean distances between chromosomal regions and then satisfying these distances using a structure reconstruction method rigorously tested in the field of protein structure determination.ResultsWe first evaluate the robustness of the overall reconstruction algorithm on noisy simulated data at various levels of noise by comparing with some of the state-of-the-art reconstruction methods. Then, using simulated data, we validate that Spearman’s rank correlation coefficient between pairwise distances in the reconstructed chromosomal structures and the experimental chromosomal contact counts can be used to find optimum conversion rules for transforming interaction frequencies to wish distances. This strategy is then applied to real Hi-C data at chromosome level for optimal transformation of interaction frequencies to wish distances and for ranking and selecting structures. The chromosomal structures reconstructed from a real-world human Hi-C dataset by our method were validated by the known two-compartment feature of the human chromosome organization. We also show that our method is robust with respect to the change of the granularity of Hi-C data, and consistently produces similar structures at different chromosomal resolutions.ConclusionChromosome3D is a robust method of reconstructing chromosome three-dimensional models using distance restraints obtained from Hi-C interaction frequency data. It is available as a web application and as an open source tool at http://sysbio.rnet.missouri.edu/chromosome3d/.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3210-4) contains supplementary material, which is available to authorized users.

Conformations of a model cyclic hexapeptide, CYIQNC: 1H-NMR and molecular dynamics studies

Solution conformation of the cyclic hexapeptide sequence, [cyclo-S-Cys-Tyr-Ile-Gln-Asn-Cys-S] (CYIQNC) – a disulfide-linked fragment of a neurohypophyseal peptide hormone oxytocin (OT) – has been investigated by high-field one-dimensional (1D) and two-dimensional (2D) NMR spectroscopic methods and compared with the results obtained from computer simulation studies. 1H-NMR results based on temperature dependence of amide proton chemical shifts and nuclear Overhauser effect indicate that peptide in solution populates different conformations, characterized by two fused β-turns. The segment Ile3-Gln4-Asn5-Cys6 yields a preferred type-III β-turn at residues 4, 5 (HB, 3HN → 6CO), while the segment Cys6, Cys1-Tyr2-Ile3 exhibits inherently weaker, flexible β-turn either of type I/II’/III/half-turn at residues 1, 2 (HB, 6HN → 3CO). The computer simulation studies using a mixed protocol of distance geometry-simulated annealing followed by constrained minimization, restrained molecular dynamics, and energy minimization showed the possibility of existence of additional conformations with the hydrogen bonds, (a) 5HN → 3CO and (b) 2HN → 6CO. These results, therefore, indicate that the additional conformations obtained from both NMR and simulation studies can also be possible to the peptide. These additional conformations might have very small population in the solution and did not show their signatures in these conditions. These findings will be helpful in designing more analogs with modifications in the cyclic moiety of OT.

Abstract 5134: G-Quadruplex structure formation in the proximal promoter region of the human vascular endothelial growth factor gene

Abstract The human vascular endothelial growth factor gene, VEGF, is a key regulator of angiogenesis. It plays an important role in tumor survival, growth and metastasis. It is over-expressed in many types of human cancers including glioma and renal cell carcinoma, ovarian and pancreatic cell carcinomas. The VEGF proximal promoter region contains a poly G/poly C rich tract that is essential for basal and inducible VEGF expression. The guanine-rich (G-rich) strand on this tract is shown to form a specific secondary structure, the G-quadruplex. Here we have elucidated the solution structure of the major intramolecular G-quadruplex formed on the G-rich strand of this region in K+ solution by proton Nuclear Magnetic Resonance (NMR) spectroscopy followed by the structural calculation of the G-quadruplex using Distance Geometry Simulated Annealing (DGSA) module of XPLOR-NIH software. It is inferred that this strand adopts a well-defined parallel-stranded G-quadruplex structure, which contains three G-tetrads stabilized by Hoogsteen hydrogen bonding, as well as three double-chain reversal loops containing 1:4:1 bases, a feature important for the formation of the most thermodynamically stable conformation. The size of the 4-base loop plays a critical role in determining the most favored folded pattern of a G-quadruplex. These 4-base loop residues interact with the residues of the 3’ flanking region of VEGF sequence; also, the 1-base double-chain reversal loop on two edges of the tetrad are very stable, which overall contribute in stabilizing the 1:4:1 conformer. We have found that the two cytosines in this central 4-base loop are involved in hydrogen bonding with the 3’ terminal end residues of VEGF, which is a salient feature as this interaction is sequence-specific and differs from that in other parallel-stranded structures, e.g. c-myc, significantly. These results are also supported by the DMS footprinting results by Guo et al. (2008) which show that the guanines involved in the tetrad formation in this major G-quadruplex loop isomer 1:4:1 are well-protected against DMS methylation. Therefore, the knowledge about the secondary structure of this non B-DNA conformation of the VEGF promoter region is essential for using the G- quadruplex as a target for anticancer drugs, which could be a novel approach to anti-angiogenesis drug discovery in cancer therapy. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 5134. doi:10.1158/1538-7445.AM2011-5134

Isolation, gene expression and solution structure of a novel moricin analogue, antibacterial peptide from a lepidopteran insect, Spodoptera litura

An antibacterial peptide was isolated from a lepidopteran insect, Spodoptera litura. The molecular mass of this peptide was determined to be 4489.55 by matrix assisted laser desorption/ionization-time of flight mass (MALDI-TOF MS) spectrometry. The peptide consists of 42 amino acids and the sequence has 69–98% identity to those of moricin-related peptides, antibacterial peptides from lepidopetran insects. Thus, the peptide was designated S. litura (Sl) moricin. Sl moricin showed a broad antibacterial spectrum against Gram-positive and negative bacteria. Sl moricin gene was inducible by bacterial injection and expressed tissue—specifically in the fat body and hemocytes. Furthermore, the solution structure of Sl moricin was determined by two-dimensional (2D) 1H-nuclear magnetic resonance (NMR) spectroscopy and hybrid distance geometry-simulated annealing calculation. The tertiary structure revealed a long α-helix containing eight turns along nearly the full length of the peptide like that of moricin, confirming that Sl moricin is a new moricin-like antibacterial peptide. These results suggest that moricin is present not only in B. mori but also in other lepidopteran insects forming a gene family.

CSGT-DFT calculation of 13C and 15N NMR shielding of the backbone amide group, 13C α, and 13C β in ω-Conotoxin GVIA

Density functional theory (DFT) combined with the continuous set of gauge transformations (CSGT) method was employed to calculate chemical shielding of the carbonyl carbon ( 13C O) and nitrogen ( 15N) atoms of the backbone amide groups in ω-conotoxin GVIA, as well as alpha carbon ( 13C α) and beta carbon ( 13C β) shielding. Experimental 13C and 15N chemical shifts were determined using multidimensional NMR spectroscopy. The three-dimensional solution structure of ω-conotoxin GVIA was solved using a hybrid distance geometry-simulated annealing protocol. Peptide fragments with experimentally measured bond lengths, bond angles, and torsion angles were used as model structures in the CSGT-DFT calculations. Linear correlation plots were constructed to compare theoretical data with experiment.

Structural Basis for New Pattern of Conserved Amino Acid Residues Related to Chitin-binding in the Antifungal Peptide from the Coconut Rhinoceros Beetle Oryctes rhinoceros

Scarabaecin isolated from hemolymph of the coconut rhinoceros beetle Oryctes rhinoceros is a 36-residue polypeptide that has antifungal activity. The solution structure of scarabaecin has been determined from twodimensional 1H NMR spectroscopic data and hybrid distance geometry-simulated annealing protocol calculation. Based on 492 interproton and 10 hydrogen-bonding distance restraints and 36 dihedral angle restraints, we obtained 20 structures. The average backbone root-mean-square deviation for residues 4-35 is 0.728 +/- 0.217 A from the mean structure. The solution structure consists of a two-stranded antiparallel beta-sheet connected by a type-I beta-turn after a short helical turn. All secondary structures and a conserved disulfide bond are located in the C-terminal half of the peptide, residues 18-36. Overall folding is stabilized by a combination of a disulfide bond, seven hydrogen bonds, and numerous hydrophobic interactions. The structural motif of the C-terminal half shares a significant tertiary structural similarity with chitin-binding domains of plant and invertebrate chitin-binding proteins, even though scarabaecin has no overall sequence similarity to other peptide/polypeptides including chitin-binding proteins. The length of its primary structure, the number of disulfide bonds, and the pattern of conserved functional residues binding to chitin in scarabaecin differ from those of chitin-binding proteins in other invertebrates and plants, suggesting that scarabaecin does not share a common ancestor with them. These results are thought to provide further strong experimental evidence to the hypothesis that chitin-binding proteins of invertebrates and plants are correlated by a convergent evolution process.

Solution structure and backbone dynamics of the catalytic domain of matrix metalloproteinase-2 complexed with a hydroxamic acid inhibitor

MMP-2 is a member of the matrix metalloproteinase family that has been implicated in tumor cell metastasis and angiogenesis. Here, we describe the solution structure of a catalytic domain of MMP-2 complexed with a hydroxamic acid inhibitor (SC-74020), determined by three-dimensional heteronuclear NMR spectroscopy. The catalytic domain, designated MMP-2C, has a short peptide linker replacing the internal fibronectin-domain insertion and is enzymatically active. Distance geometry-simulated annealing calculations yielded 14 converged structures with atomic root-mean-square deviations (r.m.s.d.) of 1.02 and 1.62 Å from the mean coordinate positions for the backbone and for all heavy atoms, respectively, when 11 residues at the N-terminus are excluded. The structure has the same global fold as observed for other MMP catalytic domains and is similar to previously solved crystal structures of MMP-2. Differences observed between the solution and the crystal structures, near the bottom of the S1′ specificity loop, appear to be induced by the large inhibitor present in the solution structure. The MMP-2C solution structure is compared with MMP-8 crystal structure bound to the same inhibitor to highlight the differences especially in the S1′ specificity loop. The finding provides a structural explanation for the selectivity between MMP-2 and MMP-8 that is achieved by large inhibitors.

Solution structure of the tumor necrosis factor receptor-1 death domain

Tumor necrosis factor receptor-1 death domain (TNFR-1 DD) is the intracellular functional domain responsible for the receptor signaling activities. The solution structure of the R347K mutant of TNFR-1 DD was solved by NMR spectroscopy. A total of 20 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 1167 distance constraints and 117 torsion angle constraints. The atomic rms distribution about the mean coordinate positions for the 20 structures for residues composing the secondary structure region is 0.40 Å for the backbone atoms and 1.09 Å for all atoms. The structure consists of six antiparallel α-helices arranged in a similar fashion to the other members of the death domain superfamily. The secondary structure and three-dimensional structure of R347K TNFR1-DD are very similar to the secondary structure and deduced topology of the R347A TNFR1-DD mutant. Mutagenesis studies identified critical residues located in α2 and part of α3 and α4 that are crucial for self-interaction and interaction with TRADD. Structural superposition with previously solved proteins in the death domain superfamily reveals that the major differences between the structures reside in α2, α3, and α4. Interestingly, these regions correspond to the binding sites of TNFR1-DD, providing a structural basis for the specificity of death domain interactions and its subsequent signaling event.

Solution Structure of B. subtilis Acyl Carrier Protein

Background: Acyl carrier protein (ACP) is a fundamental component of fatty acid biosynthesis in which the fatty acid chain is elongated by the fatty acid synthetase system while attached to the 4′-phosphopantetheine prosthetic group (4′-PP) of ACP. Activation of ACP is mediated by holo-acyl carrier protein synthase (ACPS) when ACPS transfers the 4′-PP moiety from coenzyme A (CoA) to Ser36 of apo-ACP. Both ACP and ACPS have been identified as essential for E. coli viability and potential targets for development of antibiotics. Results: The solution structure of B. subtilis ACP (9 kDa) has been determined using two-dimensional and three-dimensional heteronuclear NMR spectroscopy. A total of 22 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 1050 experimental NMR restraints. The atomic rmsd about the mean coordinate positions for the 22 structures is 0.45 ± 0.08 Å for the backbone atoms and 0.93 ± 0.07 Å for all atoms. The overall ACP structure consists of a four α-helical bundle in which 4′-PP is attached to the conserved Ser36 that is located in α helix II. Conclusions: Structural data were collected for both the apo and holo forms of ACP that suggest that the two forms of ACP are essentially identical. Comparison of the published structures for E. coli ACP and actinorhodin polyketide synthase acyl carrier protein ( act apo-ACP) from Streptomyces coelicolor A3(2) with B. subtilis ACP indicates similar secondary structure elements but an extremely large rmsd between the three ACP structures (>4.3 Å). The structural difference between B. subtilis ACP and both E. coli and act apo-ACP is not attributed to an inherent difference in the proteins, but is probably a result of a limitation in the methodology available for the analysis for E. coli and act apo-ACP. Comparison of the structure of free ACP with the bound form of ACP in the ACP–ACPS complex reveals a displacement of helix II in the vicinity of Ser36. The induced perturbation of ACP by ACPS positions Ser36 proximal to coenzyme A and aligns the dipole of helix II to initiate transfer of 4′-PP to ACP.

High-resolution solution structure of the catalytic fragment of human collagenase-3 (MMP-13) complexed with a hydroxamic acid inhibitor

The high-resolution solution structure of the catalytic fragment of human collagenase-3 (MMP-13) complexed with a sulfonamide derivative of a hydroxamic acid compound (WAY-151693) has been determined by multidimensional heteronuclear NMR. A total of 30 structures were calculated for residues 7–164 by means of hybrid distance geometry-simulated annealing using a total of 3280 experimental NMR restraints. The atomic rms distribution about the mean coordinate positions for the 30 structures is 0.43(±0.05) Å for the backbone atoms, 0.80(±0.09) Å for all atoms, and 0.47(±0.04) Å for all atoms excluding disordered side-chains. The overall structure of MMP-13 is composed of a β-sheet consisting of five β-strands in a mixed parallel and anti-parallel arrangement and three α-helices where its overall fold is consistent with previously solved MMP structures. A comparison of the NMR structure of MMP-13 with the published 1.6 Å resolution X-ray structure indicates that the major differences between the structures is associated with loop dynamics and crystal-packing interactions. The side-chains of some active-site residues for the NMR and X-ray structures of MMP-13 adopt distinct conformations. This is attributed to the presence of unique inhibitors in the two structures that encounter distinct interactions with MMP-13. The major structural difference observed between the MMP-13 and MMP-1 NMR structures is the relative size and shape of the S1′ pocket where this pocket is significantly longer for MMP-13, nearly reaching the surface of the protein. Additionally, MMP-1 and MMP-13 exhibit different dynamic properties for the active-site loop and the structural Zn-binding region. The inhibitor WAY-151693 is well defined in the MMP-13 active-site based on a total of 52 distance restraints. The binding motif of WAY-151693 in the MMP-13 complex is consistent with our previously reported MMP-1:CGS-27023A NMR structure and is similar to the MMP-13: RS-130830 X-ray structure.

Solution structure of ZipA, a crucial component of Escherichia coli cell division.

ZipA, an essential component of cell division in Escherichia coli, interacts with the FtsZ protein at the midcell in one of the initial steps of septum formation. The high-resolution solution structure of the 144-residue C-terminal domain of E. coli ZipA (ZipA(185)(-)(328)) has been determined by multidimensional heteronuclear NMR. A total of 30 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 2758 experimental NMR restraints. The atomic root means square distribution about the mean coordinate positions for residues 6-142 for the 30 structures is 0.37 +/- 0.04 A for the backbone atoms, 0. 78 +/- 0.05 A for all atoms, and 0.45 +/- 0.04 A for all atoms excluding disordered side chains. The NMR solution structure of ZipA(185)(-)(328) is composed of three alpha-helices and a beta-sheet consisting of six antiparallel beta-strands where the alpha-helices and the beta-sheet form surfaces directly opposite each other. A C-terminal peptide from FtsZ has been shown to bind ZipA(185)(-)(328) in a hydrophobic channel formed by the beta-sheet providing insight into the ZipA-FtsZ interaction. An unexpected similarity between the ZipA(185)(-)(328) fold and the split beta-alpha-beta fold observed in many RNA binding proteins may further our understanding of the critical ZipA-FtsZ interaction.

NMR structure of free RGS4 reveals an induced conformational change upon binding Galpha.

Heterotrimeric guanine nucleotide-binding proteins (G-proteins) are transducers in many cellular transmembrane signaling systems where regulators of G-protein signaling (RGS) act as attenuators of the G-protein signal cascade by binding to the Galpha subunit of G-proteins (G(i)(alpha)(1)) and increasing the rate of GTP hydrolysis. The high-resolution solution structure of free RGS4 has been determined using two-dimensional and three-dimensional heteronuclear NMR spectroscopy. A total of 30 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 2871 experimental NMR restraints. The atomic rms distribution about the mean coordinate positions for residues 5-134 for the 30 structures is 0.47 +/- 0.05 A for the backbone atoms, 0. 86 +/- 0.05 A for all atoms, and 0.56 +/- 0.04 A for all atoms excluding disordered side chains. The NMR solution structure of free RGS4 suggests a significant conformational change upon binding G(i)(alpha)(1) as evident by the backbone atomic rms difference of 1. 94 A between the free and bound forms of RGS4. The underlying cause of this structural change is a perturbation in the secondary structure elements in the vicinity of the G(i)(alpha)(1) binding site. A kink in the helix between residues K116-Y119 is more pronounced in the RGS4-G(i)(alpha)(1) X-ray structure relative to the free RGS4 NMR structure, resulting in a reorganization of the packing of the N-terminal and C-terminal helices. The presence of the helical disruption in the RGS4-G(i)(alpha)(1) X-ray structure allows for the formation of a hydrogen-bonding network within the binding pocket for G(i)(alpha)(1) on RGS4, where RGS4 residues D117, S118, and R121 interact with residue T182 from G(i)(alpha)(1). The binding pocket for G(i)(alpha)(1) on RGS4 is larger and more accessible in the free RGS4 NMR structure and does not present the preformed binding site observed in the RGS4-G(i)(alpha)(1) X-ray structure. This observation implies that the successful complex formation between RGS4 and G(i)(alpha)(1) is dependent on both the formation of the bound RGS4 conformation and the proper orientation of T182 from G(i)(alpha)(1). The observed changes for the free RGS4 NMR structure suggest a mechanism for its selectivity for the Galpha-GTP-Mg(2+) complex and a means to facilitate the GTPase cycle.

NMR Solution Structure of Domain 1 of Human Annexin I Shows an Autonomous Folding Unit

Annexins are excellent models for studying the folding mechanisms of multidomain proteins because they have four-eight homologous helical domains with low identity in sequence but high similarity in folding. The structure of an isolated domain 1 of human annexin I has been determined by NMR spectroscopy. The sequential assignments of the 1H, 13C, and 15N resonances of the isolated domain 1 were established by multinuclear, multidimensional NMR spectroscopy. The solution structure of the isolated domain 1 was derived from 1,099 experimental NMR restraints using a hybrid distance geometry-simulated annealing protocol. The root mean square deviation of the ensemble of 20 refined conformers that represent the structure from the mean coordinate set derived from them was 0. 57 +/- 0.14 A and 1.11 +/- 0.19 A for the backbone atoms and all heavy atoms, respectively. The NMR structure of the isolated domain 1 could be superimposed with a root mean square deviation of 1.36 A for all backbone atoms with the corresponding part of the crystal structure of a truncated human annexin I containing all four domains, indicating that the structure of the isolated domain 1 is highly similar to that when it folded together with the other three domains. The result suggests that in contrast to isolated domain 2, which is largely unfolded in solution, isolated domain 1 constitutes an autonomous folding unit and interdomain interactions may play critical roles in the folding of annexin I.

The three-dimensional solution structure of the SRC homology domain-2 of the growth factor receptor-bound protein-2.

A set of high-resolution three-dimensional solution structures of the Src homology region-2 (SH2) domain of the growth factor receptor-bound protein-2 was determined using heteronuclear NMR spectroscopy. The NMR data used in this study were collected on a stable monomeric protein solution that was free of protein aggregates and proteolysis. The solution structure was determined based upon a total of 1439 constraints, which included 1326 nuclear Overhauser effect distance constraints, 70 hydrogen bond constraints, and 43 dihedral angle constraints. Distance geometry-simulated annealing calculations followed by energy minimization yielded a family of 18 structures that converged to a root-mean-square deviation of 1.09 A for all backbone atoms and 0.40 A for the backbone atoms of the central beta-sheet. The core structure of the SH2 domain contains an antiparallel beta-sheet flanked by two parallel alpha-helices displaying an overall architecture that is similar to other known SH2 domain structures. This family of NMR structures is compared to the X-ray structure and to another family of NMR solution structures determined under different solution conditions.

Solution structure of the spectrin repeat: a left-handed antiparallel triple-helical coiled-coil

Cytoskeletal proteins belonging to the spectrin family have an elongated structure composed of repetitive units. The three-dimensional solution structure of the 16th repeat from chicken brain α-spectrin (R16) has been determined by NMR spectroscopy and distance geometry-simulated annealing calculations. We used a total of 1035 distance restraints, which included 719 NOE-based values obtained by applying the ambiguous restraints for iterative assignment (ARIA) method. In addition, we performed a direct refinement against 1H-chemical shifts. The final ensemble of 20 structures shows an average RMSD of 1.52 Å from the mean for the backbone atoms, excluding loops and N and C termini. R16 is made up of three antiparallel α-helices separated by two loops, and folds into a left-handed coiled-coil.The basic unit of spectrin is an antiparallel heterodimer composed of two homologous chains, β and α. These assemble a tetramer via a mechanism that relies on the completion of a single repeat by association of the partial repeats located at the C terminus of the β-chain (two helices) and at the N terminus of the α-chain (one helix). This tetramer is the assemblage able to cross-link actin filaments. Model building by homology of the “tetramerization” repeat from human erythrocyte spectrin illuminates the possible role of point mutations which cause hemolytic anemias.

NMR three‐dimensional solution structure of the serine protease inhibitor cyclotheonamide A

The nmr solution conformation of cyclotheonamide A (CtA) was determined in aqueous media. The data produced 15 distance and 10 torsional constraints which were used to generate conformations using restrained simulated annealing (SA) and distance geometry/simulated annealing (DG/SA) calculations. Two different calculation protocols were performed to ensure proper sampling of conformational space and even though the torsional restraints were input differently, both calculation methods yielded the same conformation of CtA. In the structure calculations, all solutions of the Karplus equation were sampled simultaneously using the restrained SA protocol and large ranges were used for the dihedral restraints in the DG/SA protocol because all solutions to the Karplus equation could not be sampled simultaneously. The solution conformation was also compared to the solid state x-ray conformations of CtA bound to thrombin and trypsin. The conformation of the residues important for active site binding (d-Phe, h-Arg, and Pro) are nearly identical in aqueous solution and solid state with largest differences at the a-Ala and v-Tyr residues. CtA appears to be preordered in structure and does not undergo a significant conformational change upon binding to the enzyme active site. © 1997 John Wiley & Sons, Inc.

NMR three-dimensional solution structure of the serine protease inhibitor cyclotheonamide A

The nmr solution conformation of cyclotheonamide A (CtA) was determined in aqueous media. The data produced 15 distance and 10 torsional constraints which were used to generate conformations using restrained simulated annealing (SA) and distance geometry/simulated annealing (DG/SA) calculations. Two different calculation protocols were performed to ensure proper sampling of conformational space and even though the torsional restraints were input differently, both calculation methods yielded the same conformation of CtA. In the structure calculations, all solutions of the Karplus equation were sampled simultaneously using the restrained SA protocol and large ranges were used for the dihedral restraints in the DG/SA protocol because all solutions to the Karplus equation could not be sampled simultaneously. The solution conformation was also compared to the solid state x-ray conformations of CtA bound to thrombin and trypsin. The conformation of the residues important for active site binding (d-Phe, h-Arg, and Pro) are nearly identical in aqueous solution and solid state with largest differences at the a-Ala and v-Tyr residues. CtA appears to be preordered in structure and does not undergo a significant conformational change upon binding to the enzyme active site. © 1997 John Wiley & Sons, Inc.

Structure and Stability of an Immunoglobulin Superfamily Domain from Twitchin, a Muscle Protein of the Nematode Caenorhabditis elegans

The NMR solution structure of an immunoglobulin superfamily module of twitchin (Ig 18′) has been determined and the kinetic and equilibrium folding behaviour characterised. Thirty molecular coordinates were calculated using a hybrid distance geometry-simulated annealing protocol based on 1207 distance and 48 dihedral restraints. The atomic rms distributions about the mean coordinate for the ensemble of structures is 0.55(±0.09) Å for backbone atoms and 1.10(±0.08) Å for all heavy atoms. The protein has a topology very similar to that of telokin and the titin Ig domains and thus it falls into the I set of the immunoglobulin superfamily. The close agreement between the predicted and observed structures of Ig 18′ demonstrates clearly that the I set profile can be applied in the structure prediction of immunoglobulin-like domains of diverse modular proteins. Folding studies reveal that the protein has relatively low thermodynamic stability, Δ G H 2O U–F=4.0 kcal mol −1at physiological pH. Unfolding studies suggest that the protein has considerable kinetic stability, the half life of the unfolding is greater than 40 minutes in the absence of denaturant.

Conformational study of [Met5]enkephalin-Arg-Phe in the presence of phosphatidylserine vesicles.

The interaction of [Met5]enkephalin-Arg.Phe with phosphatidylserine (PtdSer) was studied by circular dichroism (CD), two-dimensional nuclear magnetic resonance spectroscopy, hybrid distance geometry simulated annealing (DG-SA) and molecular dynamics (MD) calculations. The very low solubility of [Met5]enkephalin-Arg-Phe and the instability of the solution containing PtdSer vesicles at low pH values did not allow us to observe the amide proton resonances in the usual two-dimensional NMR work. NOESY cross-peaks of protons of side chains from two-dimensional NMR were converted into distances which were used as restraints for modelling with DG-SA and MD. Our results indicate that, in aqueous solutions at pH 7.68 [Met5]enkephalin-Arg-Phe exists in the absence of PtdSer as a random distribution of conformers, whereas in the presence of PtdSer it adopts conformations containing a common orientation of the bonds of C alpha 2, C alpha 3, C alpha 4, and C alpha 5, although different orientations of the peptide planes are consistent with the results. Two of the reported conformers from MD simulations are characterized by the presence of a 2<--4 gamma and inverse gamma turns centered on Gly3. A gradual decline of order was observed when moving from the central moiety of the peptide to both the N-terminus and C-terminus. Finally, the DG-SA and MD calculations resulted in a structure such that the orientation of the Phe4 and Met5 side chains favours hydrophobic interactions with the apolar portion of the PtdSer vesicle to form a hydrophobic cluster. These data support the hypothesis of a role of lipids to modify the conformation of [Met5]enkephalin-Arg-Phe to permit the interactions with the receptor site.

DNA Bending and Unwinding Associated with Actinomycin D Antibiotics Bound to Partially Overlapping Sites on DNA

Actinomycin D (ActD) is a potent anti-tumor antibiotic, that preferentially targets (G-C)·(G-C) steps on duplex DNA. We have reported on the solution structure of the ActD – d(A-A-A-G-C-T-T-T) complex (one drug per duplex) based on a combined application of NMR and molecular dynamics calculations. This study established that ActD binds to DNA through intercalation of the phenoxazone chromophore between (G-C)·(G-C) steps with the benzenoid and quinonoid-linked cyclic pentapeptide lactone rings spanning two base-pairs in opposite directions in the minor groove of the helix. This research is now extended to the binding of two ActD molecules to adjacent complexation sites within a (G-C-G-C)·(G-C-G-C) segment in the self-complementary d(A1-A2-G3-C4-G5-C6-T7-T8) duplex. The occupation of the central (C4-G5)·(C4-G5) segment between the two intercalation sites by the inwardly pointing cyclic pentapeptide lactone rings from adjacent bound ActD molecules should result in a potential steric clash in the center of the helix. The NMR data and its analysis on the ActD – d(A-A-G-C-G-C-T-T) complex (two drugs per duplex) establish that two ActD molecules intercalate into symmetry-related (G3-C4)·(G5-C6) steps with their attached benzenoid and quinonoid cyclic pentapeptide lactone rings positioned in the minor groove and directed towards the center and the ends of the helix, respectively. The solution structure of the complex was solved by using NMR restraints to guide distance geometry-simulated annealing and restrained molecular dynamics calculations including intensity-based refinements. The DNA helix exhibits a pronounced kink and is fully unwound at the central (C4-G5)·(C4-G5) step which results in an opening and widening of the minor groove to generate additional space for accommodation of the inwardly pointing benzenoid cyclic pentapeptide lactone rings in the complex. The outwardly and inwardly pointing cyclic pentapeptide lactone rings of symmetry- related ActD molecules retain similar conformations with the largest difference observed for the L-MeVal residues in the complex. The present study defines how structural changes primarily in the DNA associated with the directional bending of the helix towards the major groove and away from the bound drug opens up and widens the minor groove to accommodate two intercalated ActD molecules bound at partially overlapping sites on the DNA.