Structural study of copper(I)-bleomycin

Previous NMR studies on Cu(I)-bleomycin have suggested that this adduct has a geometry distinct from Fe(II)BLM. The coordination chemistry of this bleomycin derivative has been investigated through the extension of the NMR data reported previously, and the use of molecular dynamics calculations. The data collected from the NMR experiments support the coordination to the metal center of the primary and secondary amines in beta-aminoalanine and the pyrimidine ring. The detection in the NMR spectra of the signal derived from the amide hydrogen in beta-hydroxyhistidine indicates that this amide is protonated in Cu(I)-bleomycin, precluding participation of the pyrimidinyl carboxamide nitrogen in the coordination of Cu(I), as previously reported. Three-dimensional solution structures compatible with the NMR data have been assayed for Cu(I)-bleomycin for the first time by way of molecular dynamics calculations, and two models showing four and five coordination have been found to be those that better fit the experimental data. In both models the primary amine in beta-aminoalanine is coordinated such that it is located on the same side, with respect to the coordination cage, as the peptide linker fragment. This result seems important for the favored models to be compatible with either their possible oxidation to become one of the reported structures for Cu(II)BLM, or their transformation into Fe(II) adducts able to cause DNA damage.     

The TpG chelate of cis(diammineplatinum) forms two head-to-head rotamers in H2O solution

 d(TpG)^– reacts with cis-[Pt(NH₃)₂(H₂O)₂]²⁺ in two steps to yield the platinum chelate cis-[Pt(NH₃)₂{d(TpG)-N3(1),N7(2)}]. In the latter, hindered rotation of the bases leads to an equilibrium between two rotamers interconverting slowly on the NMR time scale. The structure of the two rotameric chelates was studied by means of ¹H NMR and molecular modeling techniques. The major and minor rotamers could be assigned unambiguously to the two head-to-head conformational domains which are characterized by syn/anti and anti/anti sugar-base orientations, respectively. Molecular models derived for both rotamers show that the orientations of the bases are mutually quasi-enantiomeric. The interconversion between the two rotamers (k ≈ 1 s^–1 at 293 K) is approximately 10⁴ times faster than the analogous rotamer interconversion observed in cis-[Pt(NH₃)₂{r(CpG)-N3(1),N7(2)}]⁺ [Girault J-P, Chottard G, Lallemand J-Y, Huguenin F, Chottard J-C (1984) J Am Chem Soc 106 : 7227–7232], suggesting that the steric clash of the exocyclic amino group of the platinum-bound cytosine with the ligands in cis position is more severe than that of the two thymine oxo groups. Received: 23 June 1997 / Accepted: 30 September 1997     

Molecular dynamics (MD) simulations for the prediction of chiral discrimination of N-acetylphenylalanine enantiomers by cyclomaltoheptaose (beta-cyclodextrin, beta-CD) based on the MM-PBSA (molecular mechanics-Poisson-Boltzmann surface area) approach

Molecular dynamics (MD) simulations were performed for the prediction of chiral discrimination of N-acetylphenylalanine enantiomers by cyclomaltoheptaose (beta-cyclodextrin, beta-CD). Binding free energies and various conformational properties were obtained using by the MM-PBSA (molecular mechanics Poisson-Boltzmann/surface area) approach. The calculated relative difference (DeltaDeltabinding) of binding free energy was in fine agreement with the experimentally determined value. The difference of rotameric distributions of guest N-acetylphenylalanine enantiomers complexed with the host, beta-CD, was observed after the conformational analyses, suggesting that the conformational changes of guest captured within host cavity would be a decisive factor for enantiodifferentiation at a molecular level.     

A computer-based protocol for semiautomated assignments and 3D structure determination of proteins

Summary A strategy is presented for the semiautomated assignment and 3D structure determination of proteins from heteronuclear multidimensional Nuclear Magnetic Resonance (NMR) data. This approach involves the computer-based assignment of the NMR signals, identification of distance restraints from nuclear Overhauser effects, and generation of 3D structures by using the NMR-derived restraints. The protocol is described in detail and illustrated on a resonance assignment and structure determination of the FK506 binding protein (FKBP, 107 amino acids) complexed to the immunosuppressant, ascomycin. The 3D structures produced from this automated protocol attained backbone and heavy atom rmsd of 1.17 and 1.69 Å, respectively. Although more highly resolved structures of the complex have been obtained by standard interpretation of NMR data (Meadows et al. (1993) Biochemistry, 32, 754–765), the structures generated with this automated protocol required minimal manual intervention during the spectral assignment and 3D structure calculations stages. Thus, the protocol may yield an approximate order of magnitude reduction in the time required for the generation of 3D structures of proteins from NMR data.     

Structure determination of symmetric homo-oligomers by a complete search of symmetry configuration space, using NMR restraints and van der Waals packing

Structural studies of symmetric homo-oligomers provide mechanistic insights into their roles in essential biological processes, including cell signaling and cellular regulation. This paper presents a novel algorithm for homo-oligomeric structure determination, given the subunit structure, that is both complete, in that it evaluates all possible conformations, and data-driven, in that it evaluates conformations separately for consistency with experimental data and for quality of packing. Completeness ensures that the algorithm does not miss the native conformation, and being data-driven enables it to assess the structural precision possible from data alone. Our algorithm performs a branch-and-bound search in the symmetry configuration space, the space of symmetry axis parameters (positions and orientations) defining all possible C(n) homo-oligomeric complexes for a given subunit structure. It eliminates those symmetry axes inconsistent with intersubunit nuclear Overhauser effect (NOE) distance restraints and then identifies conformations representing any consistent, well-packed structure to within a user-defined similarity level. For the human phospholamban pentamer in dodecylphosphocholine micelles, using the structure of one subunit determined from a subset of the experimental NMR data, our algorithm identifies a diverse set of complex structures consistent with the nine intersubunit NOE restraints. The distribution of determined structures provides an objective characterization of structural uncertainty: backbone RMSD to the previously determined structure ranges from 1.07 to 8.85 A, and variance in backbone atomic coordinates is an average of 12.32 A(2). Incorporating vdW packing reduces structural diversity to a maximum backbone RMSD of 6.24 A and an average backbone variance of 6.80 A(2). By comparing data consistency and packing quality under different assumptions of oligomeric number, our algorithm identifies the pentamer as the most likely oligomeric state of phospholamban, demonstrating that it is possible to determine the oligomeric number directly from NMR data. Additional tests on a number of homo-oligomers, from dimer to heptamer, similarly demonstrate the power of our method to provide unbiased determination and evaluation of homo-oligomeric complex structures.     

Identification of key residues involved in fibril formation by the conserved N-terminal region of Plasmodium falciparum merozoite surface protein 2 (MSP2)

Merozoite surface protein 2 (MSP2) from the human malaria parasite Plasmodium falciparum is expressed as a GPI-anchored protein on the merozoite surface. MSP2 is assumed to have a role in erythrocyte invasion and is a leading vaccine candidate. Recombinant MSP2 forms amyloid-like fibrils upon storage, as do peptides corresponding to sequences in the conserved N-terminal region, which constitutes the structural core of fibrils formed by full-length MSP2. We have investigated the roles of individual residues in fibril formation and local ordered structure in two peptides, a recombinant 25-residue peptide corresponding to the entire N-terminal domain of mature MSP2 and an 8-residue peptide from the central region of this domain (residues 8–15). Both peptides formed fibrils that were similar to amyloid-like fibrils formed by full-length MSP2. Phe11 and Ile12 have important roles both in stabilising local structure in these peptides and promoting fibril formation; the F11A and I12A mutants of MSP2_8–15 were essentially unstructured in solution and fibril formation at pH 7.4 and 4.7 was markedly retarded. The T10A mutant showed intermediate behaviour, having a less well defined structure than wild-type and slower fibril formation at pH 7.4. The mutation of Phe11 and Ile12 in MSP2_1–25 significantly retarded but did not abolish fibril formation, indicating that these residues also play a key role in fibril formation by the entire N-terminal conserved region. These mutations had little effect on the aggregation of full-length MSP2, however, suggesting that regions outside the conserved N-terminus have unanticipated importance for fibril formation in the full-length protein.     

The role of dynamics in modulating ligand exchange in intracellular lipid binding proteins

Lipids are essential for many biological processes and crucial in the pathogenesis of several diseases. Intracellular lipid-binding proteins (iLBPs) provide mobile hydrophobic binding sites that allow hydrophobic or amphipathic lipid molecules to penetrate into and across aqueous layers. Thus iLBPs mediate the lipid transport within the cell and participate to a spectrum of tissue-specific pathways involved in lipid homeostasis. Structural studies have shown that iLBPs' binding sites are inaccessible from the bulk, implying that substrate binding should involve a conformational change able to produce a ligand entry portal. Many studies have been reported in the last two decades on iLBPs indicating that their dynamics play a pivotal role in regulating ligand binding and targeted release. The ensemble of reported data has not been reviewed until today. This review is thus intended to summarize and possibly generalize the results up to now described, providing a picture which could help to identify the missing notions necessary to improve our understanding of the role of dynamics in iLBPs' molecular recognition. Such notions would clarify the chemistry of lipid binding to iLBPs and set the basis for the development of new drugs.     

Detergents for the stabilization and crystallization of membrane proteins

The use of detergents for the structural study of membrane proteins is discussed with an emphasis on practical issues relating to membrane solubilization, protein aggregation, detergent purity and detergent quantitation. Detergents are useful reagents as mimics of lipid bilayers because of their self-assembling properties, but as a result, they have complex properties in solution. It can be difficult to maintain a solubilized membrane protein in a native conformational state, and the non-specific aggregation of detergent-solubilized proteins is a common problem. Empirical "stability screens" can be helpful in choosing which detergents, and which detergent concentrations, may be optimal for a given system.     

Solution structure of epiregulin and the effect of its C-terminal domain for receptor binding affinity

Epiregulin (EPR), a novel member of epidermal growth factor (EGF) family, is a ligand for ErbB-1 and ErbB-4 receptors. The binding affinity of EPR for the receptors is lower than those of other EGF-family ligands. The solution structure of EPR was determined using two-dimensional nuclear magnetic resonance spectroscopy. The secondary structure in the C-terminal domain of EPR is different from other EGF-family ligands because of the lack of hydrogen bonds. The structural difference in the C-terminal domain may provide an explanation for the reduced binding affinity of EPR to the ErbB receptors.     

Fixing the conformations of diamineplatinum(II)-GpG chelates: NMR and CD signatures of individual rotamers

The bulky, asymmetric analog of the antitumor drug cisplatin, [PtCl(2)(tmen)] (tmen = N,N,N'-trimethylethylenediamine), was used to produce crosslinks with the dinucleotide d(GpG), modeling the most frequent lesions that cisplatin and its analogs cause to DNA. The ligand tmen was chosen because it is expected to constrain the guanine cis to the NMe(2) group in the adduct [Pt(tmen){d(GpG)}](+) to an orientation perpendicular to the coordination plane and to stabilize the other guanine in an oblique orientation, thus maintaining a head-to-head geometry typical of cisplatin-d(GpG) crosslinks within single- and double-stranded DNA. Of the four possible combinations of tmen chirality (R or S symmetry of the coordinated NHMe group) and crosslink direction (5'-G bound cis to the secondary or the tertiary amino group of tmen), two isomers were preponderantly formed, [Pt(R-tmen){d(GpG)}](+) with 5'-G bound cis to NMe(2) and [Pt(S-tmen){d(GpG)}](+) with 5'-G bound cis to NHMe. The former was shown to have a right-handed R2 orientation of guanines similar to that found in duplex DNA, whereas the latter had a left-handed L1 orientation that modeled cisplatin-d(GpG) adducts within single-stranded DNA. The R2 rotamer was found to be in an equilibrium (as observed using EXSY spectroscopy) with a minor fraction (< or =4%) of a Delta-HT rotamer related to R2 by rotation of the 3'-G about the Pt-N7 bond. The major rotamers R2 and L1 were isolated using reverse-phase HPLC, and their NMR and CD signatures were compared to those of the corresponding rotamers of the less hindered adduct [Pt(dmen)(GpG)](+) (dmen = N,N-dimethylethylenediamine). From this and other comparisons with previously reported platinum dinucleotide complexes, and from molecular modeling, it could be concluded that both steric repulsion between guanine and substituents of the cis amino group and N-H...O6 hydrogen bonding are significant effects favoring the oblique orientation of one guanine base typical of the HH rotamers of [Pt(diamine){d(GpG)}](+) and [Pt(diamine)(GpG)](+) complexes.     

Structure and dynamics of the von Willebrand Factor C6 domain

Von Willebrand disease (VWD) is a bleeding disorder with different levels of severity. VWD-associated mutations are located in the von Willebrand factor (VWF) gene, coding for the large multidomain plasma protein VWF with essential roles in hemostasis and thrombosis. On the one hand, a variety of mutations in the C-domains of VWF are associated with increased bleeding upon vascular injury. On the other hand, VWF gain-of-function (GOF) mutations in the C4 domain have recently been identified, which induce an increased risk of myocardial infarction. Mechanistic insights into how these mutations affect the molecular behavior of VWF are scarce and holistic approaches are challenging due to the multidomain and multimeric character of this large protein. Here, we determine the structure and dynamics of the C6 domain and the single nucleotide polymorphism (SNP) variant G2705R in C6 by combining nuclear magnetic resonance spectroscopy, molecular dynamics simulations and aggregometry. Our findings indicate that this mutation mostly destabilizes VWF by leading to a more pronounced hinging between both subdomains of C6. Hemostatic parameters of variant G2705R are close to normal under static conditions, but the missense mutation results in a gain-of-function under flow conditions, due to decreased VWF stem stability. Together with the fact that two C4 variants also exhibit GOF characteristics, our data underline the importance of the VWF stem region in VWF’s hemostatic activity and the risk of mutation-associated prothrombotic properties in VWF C-domain variants due to altered stem dynamics.     

A method for determining overall protein fold from NMR distance restraints

Summary We describe a simple method for determining the overall fold of a polypeptide chain from NOE-derived distance restraints. The method uses a reduced representation consisting of two particles per residue, and a force field containing pseudo-bond and pseudo-angle terms, an electrostatic term, but no van der Waals or hard shell repulsive terms. The method is fast and robust, requiring relatively few distance restraints to approximate the correct fold, and the correct mirror image is readily determined. The method is easily implemented using commercially available molecular modeling software.     

SESAME: A least-squares approach to the evaluation of protein structures computed from NMR data

Summary A method is proposed for defining a probability distribution on an ensemble of protein conformations from a 2D NOE spectrum, while at the same time back-calculating the experimental spectrum from the ensemble. This enables one to assess the relative quality and significance of the conformations, and to test the consistency of the ensemble as a whole with the experimental spectrum. The method eliminates the need to integrate the cross-peak intensities and is surprisingly insensitive to random noise in the spectrum. In this communication, these advantages are demonstrated by applying the method to simulated data, for which the correct result is already known.Abbreviations NOE nuclear Overhauser effect - BPTI basic pancreatic trypsin inhibitor - rmsd root-mean-square deviation     

Effect of glucose-6-P on the catalytic and structural properties of glycogen phosphorylase a

A monoclonal antibody to porcine beta-lipotropin has been produced which binds to the N-terminal (gamma-lipotropin) portion of the molecule. The antibody can be used to detect beta-lipotropin as well as other beta-endorphin precursors (predominantly a Mr 38 000 polypeptide) using radiobinding assay or the immunoblotting technique. Purification of the peptides can be readily achieved by affinity chromatography using the monoclonal antibody covalently bound to Sepharose 4B. As the antibody recognises the N-terminal part of beta-lipotropin, it can be used to detect and purify beta-lipotropin and other beta-endorphin precursors in the presence of beta-endorphin.     

The new program OPAL for molecular dynamics simulations and energy refinements of biological macromolecules

Summary A new program for molecular dynamics (MD) simulation and energy refinement of biological macromolecules, OPAL, is introduced. Combined with the supporting program TRAJEC for the analysis of MD trajectories, OPAL affords high efficiency and flexibility for work with diferent force fields, and offers a user-friendly interface and extensive trajectory analysis capabilities. Salient features are computational speeds of up to 1.5 GFlops on vector supercomputers such as the NEC SX-3, ellipsoidal boundaries to reduce the system size for studies in explicit solvents, and natural treatment of the hydrostatic pressure. Practical applications of OPAL are illustrated with MD simulations of pure water, energy minimization of the NMR structure of the mixed disulfide of a mutant E. coli glutaredoxin with glutathione in different solvent models, and MD simulations of a small protein, pheromone Er-2, using either instantaneous or time-averaged NMR restraints, or no restraints.Abbreviations D diffusion constant in cm²/s - Er-2 pheromone 2 from Euplotes raikovi - GFlop one billion floating point operations per second - Grx(C14S)-SG mixed disulfide between a mutant E. coli glutaredoxin, with Cys¹⁴ replaced by Ser, and glutathione - MD molecular dynamics - NOE nuclear Overhauser enhancement - rmsd root-mean-square deviation - density in g/cm³     

Coarse-grained modelling of urea-adamantyl functionalised poly(propylene imine) dendrimers

To investigate the behaviour of poly(propylene imine) dendrimers – and urea–adamantyl functionalised ones – in solution using molecular dynamics simulations, we developed a coarse-grained model to tackle the relatively large system sizes and time scales needed. Harmonic bond and angle potentials were derived from atomistic simulations using an iterative Boltzmann inversion scheme, modified to incorporate Gaussian fits of the bond and angle distributions. With the coarse-grained model and accompanying force field simulations of generations 1–7 of both dendrimer types in water were performed. They compare favourably with atomistic simulations and experimental results on the basis of size, shape, monomer density, spacer back-folding and atomic form factor measurements. These results show that the structural dynamics of these dendrimers originate from flexible chains constrained by configurational and spatial requirements. Large dendrimers are more rigid and spherical, while small ones are flexible, alternatively rod-like and globular.     

Practical applications of time-averaged restrained molecular dynamics to ligand-receptor systems: FK506 bound to the Q50R,A95H,K98I triple mutant of FKBP-13

Summary The ability of time-averaged restrained molecular dynamics (TARMD) to escape local low-energy conformations and explore conformational space is compared with conventional simulated-annealing methods. Practical suggestions are offered for performing TARMD calculations with ligand-receptor systems, and are illustrated for the complex of the immunosuppressant FK506 bound to Q50R,A95H,K98I triple mutant FKBP-13. The structure of ¹³C-labeled FK506 bound to triple-mutant FKBP-13 was determined using a set of 87 NOE distance restraints derived from HSQC-NOESY experiments. TARMD was found to be superior to conventional simulated-annealing methods, and produced structures that were conformationally similar to FK506 bound to wild-type FKBP-12. The individual and combined effects of varying the NOE restraint force constant, using an explicit model for the protein binding pocket, and starting the calculations from different ligand conformations were explored in detail.Abbreviations DG distance geometry - dmFKBP-12 double-mutant (R42K,H87V) FKBP-12 - FKBP-12 FK506-binding protein (12 kDa) - FKBP-13 FK506-binding protein (13 kDa) - HSQC heteronuclear single-quantum coherence - K_NOE force constant (penalty) for NOE-derived distance restraints - MD molecular dynamics - NOE nuclear Overhauser effect - SA simulated annealing - TARMD molecular dynamics with time-averaged restraints - tmFKBP-13 triple-mutant (Q50R,A95H,K98I) FKBP-13 - wtFKBP-12 wild-type FKBP-12     

Nuclear magnetic resonance analysis of carotenoids from the burgundy plumage of the Pompadour Cotinga (Xipholena punicea)

Previous analysis of carotenoids extracted from the burgundy plumage of the Pompadour Cotinga (Xipholena punicea) revealed six novel keto-carotenoid pigments with methoxyl groups in the C3-position of one or both β-rings. High performance liquid chromatography (HPLC), mass spectrometry, chemical analysis and, in some instances ¹H NMR spectroscopy were employed to determine the structures of the molecules. Further analysis by NMR was precluded due to lack of material. The recent acquisition of multiple feathers from X. punicea specimens has made it possible to complete this work using correlated homonuclear spectroscopy (COSY), nuclear overhauser effect spectroscopy (NOESY) and ¹H NMR. These new data conclusively confirm the structures of the six methoxy-carotenoids suggested by the earlier work. In addition, the resonance positions of the protons from the novel 3-methoxy-4-keto-β-ring and 2,3-didehydro-3-methoxy-4-keto-β-ring moieties are reported here for the first time.     

Identification of low micromolar dual inhibitors for aldose reductase (ALR2) and poly (ADP-ribose) polymerase (PARP-1) using structure based design approach

Clinical studies have revealed that diabetic retinopathy is a multifactorial disorder. Moreover, studies also suggest that ALR2 and PARP-1 co-occur in retinal cells, making them appropriate targets for the treatment of diabetic retinopathy. To find the dual inhibitors of ALR2 and PARP-1, the structure based design was carried out in parallel for both the target proteins. A series of novel thiazolidine-2,4-dione (TZD) derivatives were therefore rationally designed, synthesized and their in vitro inhibitory activities against ALR2 and PARP-1 were evaluated. The experimental results showed that compounds 5b and 5f, with 2-chloro and 4-fluoro substitutions, showed biochemical activities in micromolar and submicromolar range (IC₅₀ 1.34–5.03 μM) against both the targeted enzymes. The structure-activity relationship elucidated for these novel inhibitors against both the enzymes provide new insight into the binding mode of the inhibitors to the active sites of enzymes. The positive results of the biochemical assay suggest that these compounds may be further optimized and utilized for the treatment of diabetic retinopathy.     

H assignment and secondary structure determination of human melanoma growth stimulating activity (MGSA) by NMR spectroscopy

The solution structure of melanoma growth stimulating activity (MGSA) has been investigated using proton NMR spectroscopy. Sequential resonance assignments have been carried out, and elements of secondary structure have been identified on the basis of NOE, coupling constant, chemical shift, and amide proton exchange data. Long-range NOEs have established that MGSA is a dimer in solution. The secondary structure and dimer interface of MGSA appear to be similar to those found previously for the homologous chemokine interleukin-8 [Clore et al. (1990) Biochemistry 29, 1689-1696]. The MGSA monomer contains a three stranded anti-parallel β-sheet arranged in a ‘Greek-key’ conformation, and a C-terininal -helix (residues 58 69).