Computer modeling studies of the structure of a repressor

Due to advances in molecular biology the DNA sequences of structural genes coding for proteins are often known before a protein is characterized or even isolated. The function of a protein whose amino acid sequence has been deduced from a DNA sequence may not even be known. This has created greater interest in the development of methods to predict the tertiary structures of proteins. The a priori prediction of a protein's structure from its amino acid sequence is not yet possible. However, since proteins with similar amino acid sequences are observed to have similar three-dimensional structures, it is possible to use an analogy with a protein of known structure to draw some conclusions about the structure and properties of an uncharacterized protein. The process of predicting the tertiary structure of a protein relies very much upon computer modeling and analysis of the structure. The prediction of the structure of the bacteriophage 434 cro repressor is used as an example illustrating current procedures.     

Protein solution structure calculations in solution: Solvated molecular dynamics refinement of calbindin D9k

The three-dimensional solution structures of proteins determinedwith NMR-derived constraints are almost always calculated in vacuo. Thesolution structure of (Ca²⁺)_{_2}-calbindinD_9k has been redetermined by new restrained molecular dynamics(MD) calculations that include Ca²⁺ ions and explicit solventmolecules. Four parallel sets of MD refinements were run to provide accuratecomparisons of structures produced in vacuo, in vacuo withCa²⁺ ions, and with two different protocols in a solvent bathwith Ca²⁺ ions. The structural ensembles were analyzed interms of structural definition, molecular energies, packing density,solvent-accessible surface, hydrogen bonds, and the coordination of calciumions in the two binding loops. Refinement including Ca²⁺ ionsand explicit solvent results in significant improvements in the precisionand accuracy of the structure, particularly in the binding loops. Theseresults are consistent with results previously obtained in free MDsimulations of proteins in solution and show that the rMD refinedNMR-derived solution structures of proteins, especially metalloproteins, canbe significantly improved by these strategies.     

A high-resolution 1H NMR approach for structure determination of membrane peptides and proteins in non-deuterated detergent: Application to mastoparan X solubilized in n-octylglucoside

Summary Application of ¹H 2D NMR methods to solubilized membrane proteins and peptides has up to now required the use of selectively deuterated detergents. The unavailability of any of the common biochemical detergents in deuterated form has therefore limited to some extent the scope of this approach. Here a ¹H NMR method is described which allows structure determination of membrane peptides and small membrane proteins by ¹H 2D NMR in any type of non-deuterated detergent. The approach is based on regioselective excitation of protein resonances with DANTE-Z or spin-pinging pulse trains. It is shown that regioselective excitation of the amide-aromatic region of solubilized membrane proteins and peptides leads to an almost complete suppression of the two orders of magnitude higher contribution of the protonated detergent to the ¹H NMR spectrum. Consistently TOCSY, COSY and NOESY sequences incorporating such regioselective excitation in the F2 dimension yield protein ¹H 2D NMR spectra of quality comparable to those obtained in deuterated detergents. Regioselective TOCSY and NOESY spectra display all through-bond and through-space correlations within amide-aromatic protons and between these protons and aliphatic and -protons. Regioselective COSY spectra provide scalar coupling constants between amide and -protons. Application of the method to the membrane-active peptide mastoparan X, solubilized in n-octylglucoside, yields complete sequence-specific assignments and extensive secondary structure-related spatial proximities and coupling constants. It is shown that mastoparan adopts an -helical conformation when bound to nonionic detergent micelles. The present method is expected to increase the applicability of ¹H solution NMR methods to membrane proteins and peptides.Abbreviations 2D NMR two-dimensional NMR - COSY correlated spectroscopy - DANTE delays alternating nutations for tailored excitation - NOESY nuclear Overhauser enhancement spectroscopy - TOCSY total correlation spectroscopy     

How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation

Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from ¹⁵N relaxation and side-chain dynamics from ¹³C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.     

One simple step in the identification of the cofactors signals, one giant leap for the solution structure determination of multiheme proteins

Multiheme proteins play major roles in various biological systems. Structural information on these systems in solution is crucial to understand their functional mechanisms. However, the presence of numerous proton-containing groups in the heme cofactors and the magnetic properties of the heme iron, in particular in the oxidised state, complicates significantly the assignment of the NMR signals. Consequently, the multiheme proteins superfamily is extremely under-represented in structural databases, which constitutes a severe bottleneck in the elucidation of their structural-functional relationships. In this work, we present a strategy that simplifies the assignment of the NMR signals in multiheme proteins and, concomitantly, their solution structure determination, using the triheme cytochrome PpcA from the bacterium Geobacter sulfurreducens as a model. Cost-effective isotopic labeling was used to double label (¹³C/¹⁵N) the protein in its polypeptide chain, with the correct folding and heme post-translational modifications. The combined analysis of ¹H-¹³C HSQC NMR spectra obtained for labeled and unlabeled samples of PpcA allowed a straight discrimination between the heme cofactors and the polypeptide chain signals and their confident assignment. The results presented here will be the foundations to assist solution structure determination of multiheme proteins, which are still very scarce in the literature.     

Three-dimensional structure of ectatomin from Ectatomma tuberculatum ant venom

Summary Two-dimensional ¹H NMR techniques were used to determine the spatial structure of ectatomin, a toxin from the venom of the ant Ectatomma tuberculatum. Nearly complete proton resonance assignments for two chains of ectatomin (37 and 34 amino acid residues, respectively) were obtained using 2D TOCSY, DQF-COSY and NOESY experiments. The cross-peak volumes in NOESY spectra were used to define the local structure of the protein and generate accurate proton-proton distance constraints employing the MARDIGRAS program. Disulfide bonds were located by analyzing the global fold of ectatomin, calculated with the distance geometry program DIANA. These data, combined with data on the rate of exchange of amide protons with deuterium, were used to obtain a final set of 20 structures by DIANA. These structures were refined by unrestrained energy minimization using the CHARMm program. The resulting rms deviations over 20 structures (excluding the mobile N- and C-termini of each chain) are 0.75 ? for backbone heavy atoms, and 1.25 ? for all heavy atoms. The conformations of the two chains are similar. Each chain consists of two α-helices and a hinge region of four residues; this forms a hairpin structure which is stabilized by disulfide bridges. The hinge regions of the two chains are connected together by a third disulfide bridge. Thus, ectatomin forms a four-α-helical bundle structure.     

The solution structure of a superpotent B-chain-shortened single-replacement insulin analogue

This paper reports on an insulin analogue with 12.5-fold receptor affinity, the highest increase observed for a single replacement, and on its solution structure, determined by NMR spectroscopy. The analogue is [D-AlaB26]des-(B27-B30)-tetrapeptide-insulin-B26-amide. C-terminal truncation of the B-chain by four (or five) residues is known not to affect the functional properties of insulin, provided the new carboxylate charge is neutralized. As opposed to the dramatic increase in receptor affinity caused by the substitution of D-Ala for the wild-type residue TyrB26 in the truncated molecule, this very substitution reduces it to only 18% of that of the wild-type hormone when the B-chain is present in full length. The insulin molecule in solution is visualized as an ensemble of conformers interrelated by a dynamic equilibrium. The question is whether the "active" conformation of the hormone, sought after in innumerable structure/function studies, is or is not included in the accessible conformational space, so that it could be adopted also in the absence of the receptor. If there were any chance for the active conformation, or at least a predisposed state to be populated to a detectable extent, this chance should be best in the case of a superpotent analogue. This was the motivation for the determination of the three-dimensional structure of [D-AlaB26]des-(B27-B30)-tetrapeptide-insulin-B26-amide. However, neither the NMR data nor CD spectroscopic comparison of a number of related analogues provided a clue concerning structural features predisposing insulin to high receptor affinity. After the present study it seems more likely than before that insulin will adopt its active conformation only when exposed to the force field of the receptor surface.     

Proton nuclear magnetic resonance assignments and secondary structure determination of the ColE1 rop (rom) protein

The complete resonance assignment of the ColE1 rop (rom) protein at pH 2.3 was obtained by two-dimensional (2D) proton nuclear magnetic resonance spectroscopy (1H NMR) at 500 and 600 MHz using through-bond and through-space connectivities. Sequential assignments and elements of regular secondary structure were deduced by analysis of nuclear Overhauser enhancement spectroscopy (NOESY) experiments and 3JHN alpha coupling constants. One 7.2-kDa monomer of the homodimer consists of two antiparallel helices connected by a hairpin loop at residue 31. The C-terminal peptide consisting of amino acids 59-63 shows no stable conformation. The dimer forms a four-helix bundle with opposite polarization of neighboring elements in agreement with the X-ray structure.     

1H and 15N NMR resonance assignments and solution secondary structure of oxidized Desulfovibrio desulfuricans flavodoxin

Summary Sequence-specific ¹H and ¹⁵N resonance assignments have been made for 137 of the 146 nonprolyl residues in oxidized Desulfovibrio desulfuricans [Essex 6] flavodoxin. Assignments were obtained by a concerted analysis of the heteronuclear three-dimensional ¹H-¹⁵N NOESY-HMQC and TOCSY-HMQC data sets, recorded on uniformly ¹⁵N-enriched protein at 300 K. Numerous side-chain resonances have been partially or fully assigned. Residues with overlapping ¹H^N chemical shifts were resolved by a three-dimensional ¹H-¹⁵N HMQC-NOESY-HMQC spectrum. Medium-and long-range NOEs, ³J_NH coupling constants, and ¹H^N exchange data indicate a secondary structure consisting of five parallel -strands and four -helices with a topology similar to that of Desulfovibrio vulgaris [Hidenborough] flavodoxin. Prolines at positions 106 and 134, which are not conserved in D. vulgaris flavodoxin, contort the two C-terminal -helices.Abbreviations CSI chemical shift index - DQF-COSY double-quantum-filtered correlation spectroscopy - DIPSI decoupling in the presence of scalar interactions - FMN flavin mononucleotide - GARP globally optimized alternating phase rectangular pulse - HMQC heteronuclear multiple-quantum coherence - HSQC heteronuclear single-quantum coherence - NOE nuclear Overhauser effect - NOESY nuclear Overhauser enhancement spectroscopy - TOCSY total correlation spectroscopy - TPPI time-proportional phase increments - TSP 3-(trimethylsilyl)propionic-2,2,3,3-d ₄ acid, sodium salt     

Leap-dynamics: efficient sampling of conformational space of proteins and peptides in solution

A molecular simulation scheme, called Leap-dynamics, that provides efficient sampling of protein conformational space in solution is presented. The scheme is a combined approach using a fast sampling method, imposing conformational 'leaps' to force the system over energy barriers, and molecular dynamics (MD) for refinement. The presence of solvent is approximated by a potential of mean force depending on the solvent accessible surface area. The method has been successfully applied to N-acetyl-L-alanine-N-methylamide (alanine dipeptide), sampling experimentally observed conformations inaccessible to MD alone under the chosen conditions. The method predicts correctly the increased partial flexibility of the mutant Y35G compared to native bovine pancreatic trypsin inhibitor. In particular, the improvement over MD consists of the detection of conformational flexibility that corresponds closely to slow motions identified by nuclear magnetic resonance techniques.     

Molecular modelling of miraculin: Structural analyses and functional hypotheses

Miraculin is a plant protein that displays the peculiar property of modifying taste by swiching sour into a sweet taste. Its monomer is flavourless at all pH as well as at high concentration; the dimer form elicits its taste-modifying activity at acidic pH; a tetrameric form is also reported as active. Two histidine residues, located in exposed regions, are the main responsible of miraculin activity, as demonstrated by mutagenesis studies. Since structural data of miraculin are not available, we have predicted its three-dimensional structure and simulated both its dimer and tetramer forms by comparative modelling and molecular docking techniques. Finally, molecular dynamics simulations at different pH conditions have indicated that at acidic pH the dimer assumes a widely open conformation, in agreement with the hypotheses coming from other studies.     

1H, 13C and 15N assignments and chemical shift-derived secondary structure of intestinal fatty acid-binding protein

Summary Sequence-specific ¹H, ¹³C and ¹⁵N resonance assignments have been established for rat intestinal fatty acid-binding protein complexed with palmitate (15.4 kDa) at pH 7.2 and 37°C. The resonance assignment strategy involved the concerted use of seven 3D triple-resonance expriments (CC-TOCSY, HCCH-TOCSY, HNCO, HNCA, ¹⁵N-TOCSY-HMQC, HCACO and HCA(CO)N). A central feature of this strategy was the concurrent assignment of both backbone and side-chain aliphatic atoms, which was critical for overcoming ambiguities in the assignment process. The CC-TOCSY experiment provided the unambiguous links between the side-chain spin systems observed in HCCH-TOCSY and the backbone correlations observed in the other experiments. Assignments were established for 124 of the 131 residues, although 6 of the 124 had missing amide ¹H resonances, presumably due to rapid exchange with solvent under these experimental conditions. The assignment database was used to determine the solution secondary structure of the complex, based on chemical shift indices for the ¹H, ¹³C, ¹³C and ¹³CO atoms. Overall, the secondary structure agreed well with that determined by X-ray crystallography [Sacchettini et al. (1989) J. Mol. Biol., 208, 327–339], although minor differences were observed at the edges of secondary structure elements.     

寡聚脱氧核苷酸d（CCGTACGG）质子共振谱线归属和溶液物象表征

寡聚脱氧核苷酸ｄ（ＣＣＧＴＡＣＧＧ）质子共振谱线归属和溶液物象表征王萍,石根斌,宋国强,陈凯先,嵇汝运（中国科学院上海药物研究所,２０００３１）关键词寡聚脱氧核苷酸;２ＤＮＭＲ;溶液构象石蒜内铵是一种新型ＤＮＡ嵌合剂,它可以显著改变ＤＮＡ螺旋的构象。...     

NMR-based homology model for the solution structure of the C-terminal globular domain of EMILIN1

EMILIN1 is a glycoprotein of elastic tissues that has been recently linked to the pathogenesis of hypertension. The protein is formed by different independently folded structural domains whose role has been partially elucidated. In this paper the solution structure, inferred from NMR-based homology modelling of the C-terminal trimeric globular C1q domain (gC1q) of EMILIN1, is reported. The high molecular weight and the homotrimeric structure of the protein required the combined use of highly deuterated ¹⁵N, ¹³C-labelled samples and TROSY experiments. Starting from a homology model, the protein structure was refined using heteronuclear residual dipolar couplings, chemical shift patterns, NOEs and H-exchange data. Analysis of the gC1q domain structure of EMILIN1 shows that each protomer of the trimer adopts a nine-stranded β sandwich folding topology which is related to the conformation observed for other proteins of the family. Distinguishing features, however, include a missing edge-strand and an unstructured 19-residue loop. Although the current data do not allow this loop to be precisely defined, the available evidence is consistent with a flexible segment that protrudes from each subunit of the globular trimeric assembly and plays a key role in inter-molecular interactions between the EMILIN1 gC1q homotrimer and its integrin receptor α4β1. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.