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.     

Conformational and biological characterization of human parathyroid hormone hPTH(1-34) analogues containing beta-amino acid residues in positions 17-19

The N-terminal 1-34 fragment of parathyroid hormone (PTH) elicits the full spectrum of bone-related biological activities of the intact native sequences. It has been suggested that the structural elements essential for bioactivity are two helical segments located at the N-terminal and C-terminal sequences, connected by hinges or flexible points around positions 12 and 19. In order to assess the relevance of the local conformation around Gly(18) upon biological function, we synthesized and characterized the following human (h) PTH(1-34) analogues containing beta-amino acid residues: [analogues: see text]. Biological activity and binding affinity of analogue I are one order of magnitude lower than those of the parent compound. In analogue II, both binding affinity and biological activity are partially recovered. Analogues III and V have no binding affinity and very low biological activity. Both bioactivity and binding affinity are partially recovered in analogue IV. The conformational properties of the analogues in aqueous solution containing dodecylphosphocholine micelles were studied by CD, 2D-nuclear magnetic resonance and molecular dynamics calculations. The results confirmed the presence in all analogues of two helical segments located at the N-terminal and C-terminal sequences. The insertion of beta-amino acid residues around position 18 does not cause appreciable conformational differences in the five analogues. The differences in biological activity and binding affinity among the five analogues cannot be related to structural differences in the membrane mimetic environment reported in this study. Our results stress the importance of the side-chain functionalities in the sequence 17-19 for biological function.     

The dynamic properties of melittin in solution

Molecular dynamics simulations are described for the peptide melittin. The atomic trajectories are calculated both with normal potential energy functions and with additional distance restraints deduced from nuclear Overhauser effects observed in NMR experiments. The results are compared with NRM data on coupling constants and amide exchange rates and witt B-factors from X-ray crystallography. The observed correlations between experiment and molecular dynamics simulations suggest a relatively mobile C-terminus and relatively high flexibility around residue 11. It is noted that the high conformational variation around residue 11 is due in part to the presence of a proline at position 14 which results in a missing H-bond in the largely -helical structure. It is also noted that a proline is a common feature of many putative membrane spanning helices. A role for such prolines is suggested.     

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     

Structure-function relationship studies of bovine parathyroid hormone [bPTH(1-34)] analogues containing alpha-amino-iso-butyric acid (Aib) residues

The N-terminal 1-34 fragments of the parathyroid hormone (PTH) and parathyroid hormone-related protein (PTHrP) elicit the full spectrum of bone-related biological activities of the intact native sequences. It has been suggested that the structural elements essential for bioactivity are two helical segments located at the N-terminal and C-terminal sequences, connected by hinges or flexible points around positions 12 and 19. In order to assess the relevance of the local conformation around Gly(12) upon biological function, we synthesized and characterized the following PTH(1-34) analogues containing Aib residues: (I) A-V-S-E-I-Q-F-nL-H-N-Aib-G-K-H-L-S-S-nL-E-R-V-E-Nal-L-R-K-K-L-Q-D-V-H-N-Y-NH(2) ([Nle(8,18), Aib(11), Nal(23),Tyr(34)]bPTH(1-34)-NH(2)); (II) A-V-S-E-I-Q-F-nL-H-N-L-Aib-K-H-L-S-S-nL-E-R-V-E-Nal-L-R-K-K-L-Q-D-V-H-N-Y-NH(2) ([Nle(8,18), Aib(12),Nal(23),Tyr(34)]bPTH(1-34)-NH(2)); (III) A-V-S-E-I-Q-F-nL-H-N-L-G-Aib-H-L-S-S-nL-E-R-V-E-Nal-L-R-K-K-L-Q-D-V-H-N-Y-NH(2) ([Nle(8,18), Aib(13), Nal(23),Tyr(34)]bPTH(1-34)-NH(2)); (IV) A-V-S-E-I-Q-F-nL-H-N-Aib-Aib-K-H-L-S-S-nL-E-R-V-E-Nal-L-R-K-K-L-Q-D-V-H-N-YNH(2) ([Nle(8,18), Aib(11,12), Nal(23),Tyr(34)]bPTH(1-34)-NH(2)); (V) A-V-S-E-I-Q-F-nL-H-N-L-Aib-Aib-H-L-S-S-nL-E-R-V-E-Nal-L-R-K-K-L-Q-D-V-H-N-Y-NH(2) ([Nle(8,18), Aib(12,13),Nal(23),Tyr(34)]bPTH(1-34)-NH(2)). (nL= Nle; Nal= L-(2-naphthyl)-alanine; Aib= alpha-amino-isobutyric acid.) The introduction of Aib residues at position 11 in analogue I or at positions 11 and 12 in analogue IV resulted in a 5-20-fold lower efficacy and a substantial loss of binding affinity compared to the parent compound [Nle(8,18), Nal(23),Tyr(34)]bPTH(1-34)-NH(2). Both binding affinity and adenylyl cyclase stimulation activity are largely restored when the Aib residues are introduced at position 12 in analogue II, 13 in analogue III, and 12-13 in analogue V. The conformational properties of the analogues in aqueous solution containing dodecylphosphocholine micelles were studied by CD, two-dimensional (2D) NMR and computer simulations. The results indicated the presence of two helical segments in all analogues, located at the N-terminal and C-terminal sequences. Insertion of Aib residues at positions 12 and 13, or of Aib dyads at positions 11-12 and 12-13, enhances the stability of the N-terminal helix of all analogues. In all analogues the Aib residues are included in the helical segments. These results confirmed the importance of the helical structure in the N-terminal activation domain, as well as of the presence of the Leu(11) hydrophobic side chain in the native sequence, for PTH-like bioactivity.     

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.     

Unconstrained peptoid tetramer exhibits a predominant conformation in aqueous solution

Conformational control in peptoids, N-substituted glycines, is crucial for the design and synthesis of biologically-active compounds and atomically-defined nanomaterials. While there are a growing number of structural studies in solution, most have been performed with conformationally-constrained short sequences (e.g., sterically-hindered sidechains or macrocyclization). Thus, the inherent degree of heterogeneity of unconstrained peptoids in solution remains largely unstudied. Here, we explored the folding landscape of a series of simple peptoid tetramers in aqueous solution by NMR spectroscopy. By incorporating specific ¹³C-probes into the backbone using bromoacetic acid-2-¹³C as a submonomer, we developed a new technique for sequential backbone assignment of peptoids based on the 1,n-Adequate pulse sequence. Unexpectedly, two of the tetramers, containing an N-(2-aminoethyl)glycine residue (Nae), had preferred conformations. NMR and molecular dynamics studies on one of the tetramers showed that the preferred conformer (52%) had a trans-cis-trans configuration about the three amide bonds. Moreover, >80% of the ensemble contained a cis amide bond at the central amide. The backbone dihedral angles observed fall directly within the expected minima in the peptoid Ramachandran plot. Analysis of this compound against similar peptoid analogs suggests that the commonly used Nae monomer plays a key role in the stabilization of peptoid structure via a side-chain-to-main-chain interaction. This discovery may offer a simple, synthetically high-yielding approach to control peptoid structure, and suggests that peptoids have strong intrinsic conformational preferences in solution. These findings should facilitate the predictive design of folded peptoid structures, and accelerate application in areas ranging from drug discovery to biomimetic nanoscience.     

Knowledge-based Homology Modeling and Experimental Determination of Amino Acid Side Chain Accessibility by the Laser Photo CIDNP (Chemically Induced Dynamic Nuclear Polarization) Approach in Solution: Lessons from the Small Sialidase of Clostridium perfringens

The success of knowledge-based homology modelling is critically dependent on the predictive potency of the program structure-based calculations, which attempt to translate homologous sequences into three-dimensional structures, and on the actual relevance of the crystal structure for the protein topology. As quality control, experimental data for selected parameters of the proteins conformation are required. Using the crystal structure of the sialidase of Salmonella typhimurium as framework for model building of the homologous enzyme from Clostridium perfringens, a set of energy-minimised conformers is derived. These proteins present e.g. Tyr, Trp and His residues with an assessable area on the surface, since the side chains of these amino acid residues are responsive to chemically induced dynamic nuclear polarization (CIDNP), monitored by NMR. Hence, as first lesson, a comparative analysis for model-derived and experimentally determined values can be performed. The second lesson of this study concerns the notable impact of single amino acid substitutions (Tyr/Phe, Cys/Ser) on the surface accessibility of the CIDNP-reactive amino acid side chains in mutant forms of the sialidase. Corroborating the predictions from the theoretical calculations, the spectra of the engineered mutants reveal marked and non-uniform alterations. Thus, the effect of apparently rather conservative amino acid substitutions on a distinct conformational aspect of this protein, even at distant sites, should not be underestimated.     

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.     

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.     

Structure of Bombyx mori chemosensory protein 1 in solution

Chemosensory Proteins (CSPs) represent a family of conserved proteins found in insects that may be involved in chemosensory functions. BmorCSP1 is expressed mainly in antennae and legs of the silkworm moth Bombyx mori and was cloned from antennal cDNA. Here we report the determination of the structure of Bombyx mori CSP1 (BmorCSP1) by NMR. The overall fold of BmorCSP1 is globular and comprises six alpha-helices. These helices span residues 10-14, 17-27, 35-49, 57-72, 75-85, and 92-100. The internal hydrophobic sides of the helices are formed mostly by leucine and isoleucine residues and, therefore, well suited to constitute a binding site for hydrophobic ligands.     

Conformational dynamics of wild-type Lys-48-linked diubiquitin in solution

Proteasomal degradation is mediated through modification of target proteins by Lys-48-linked polyubiquitin (polyUb) chain, which interacts with several binding partners in this pathway through hydrophobic surfaces on individual Ub units. However, the previously reported crystal structures of Lys-48-linked diUb exhibit a closed conformation with sequestered hydrophobic surfaces. NMR studies on mutated Lys-48-linked diUb indicated a pH-dependent conformational equilibrium between closed and open states with the predominance of the former under neutral conditions (90% at pH 6.8). To address the question of how Ub-binding proteins can efficiently access the sequestered hydrophobic surfaces of Ub chains, we revisited the conformational dynamics of Lys-48-linked diUb in solution using wild-type diUb and cyclic forms of diUb in which the Ub units are connected through two Lys-48-mediated isopeptide bonds. Our newly determined crystal structure of wild-type diUb showed an open conformation, whereas NMR analyses of cyclic Lys-48-linked diUb in solution revealed that its structure resembled the closed conformation observed in previous crystal structures. Comparison of a chemical shift of wild-type diUb with that of monomeric Ub and cyclic diUb, which mimic the open and closed states, respectively, with regard to the exposure of hydrophobic surfaces to the solvent indicates that wild-type Lys-48-linked diUb in solution predominantly exhibits the open conformation (75% at pH 7.0), which becomes more populated upon lowering pH. The intrinsic properties of Lys-48-linked Ub chains to adopt the open conformation may be advantageous for interacting with Ub-binding proteins.     

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.