Calculation of protein ionization equilibria with conformational sampling: pK(a) of a model leucine zipper, GCN4 and barnase.

The use of conformational ensembles provided by nuclear magnetic resonance (NMR) experiments or generated by molecular dynamics (MD) simulations has been regarded as a useful approach to account for protein motions in the context of pK(a) calculations, yet the idea has been tested occasionally. This is the first report of systematic comparison of pK(a) estimates computed from long multiple MD simulations and NMR ensembles. As model systems, a synthetic leucine zipper, the naturally occurring coiled coil GCN4, and barnase were used. A variety of conformational averaging and titration curve-averaging techniques, or combination thereof, was adopted and/or modified to investigate the effect of extensive global conformational sampling on the accuracy of pK(a) calculations. Clustering of coordinates is proposed as an approach to reduce the vast diversity of MD ensembles to a few structures representative of the average electrostatic properties of the system in solution. Remarkable improvement of the accuracy of pK(a) predictions was achieved by the use of multiple MD simulations. By using multiple trajectories the absolute error in pK(a) predictions for the model leucine zipper was reduced to as low as approximately 0.25 pK(a) units. The validity, advantages, and limitations of explicit conformational sampling by MD, compared with the use of an average structure and a high internal protein dielectric value as means to improve the accuracy of pK(a) calculations, are discussed.     

Solution structure of a D,L-alternating oligonorleucine as a model of double-stranded antiparallel beta-helix

Conformational characteristics of alternating D,L linear peptides are of particular interest because of their capacity to form transmembrane channels with different transport properties, as some natural antibiotics do. Single- and double-stranded beta-helical structures are common for alternating D,L peptides. The stability of the beta-helix depends on several structural factors, such as the backbone peptide length, type and position of side chains, and nature of terminal groups. The NMR and molecular dynamics solution conformation of a synthetic alternating D,L-oligopeptide with 15 norleucines (XVMe) has been used as a model to get insight in to the conformational features of double-stranded beta-helix structures. The NH chemical shift values (delta(NH)) and long-range nuclear Overhauser effects (NOE) cross peaks, in particular interstrand connectivities, clearly point to an antiparallel double-stranded beta-helix for the XVMe major conformation in solution. An extensive set of distances (from NOE cross peaks) and H-bonds (from delta(NH)) has been included in the molecular dynamics calculations. The experimental NMR data and theoretical calculations clearly indicate that the most probable conformation of XVMe in solution is a double-strand antiparallel beta(5.6) increasing decreasing-helix structure.     

Simulated annealing with restrained molecular dynamics using CONGEN: energy refinement of the NMR solution structures of epidermal and type-alpha transforming growth factors.

The new functionality of the program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168; Bassolino-Klimas D et al., 1996, Protein Sci 5:593-603) has been applied for energy refinement of two previously determined solution NMR structures, murine epidermal growth factor (mEGF) and human type-alpha transforming growth factor (hTGF alpha). A summary of considerations used in converting experimental NMR data into distance constraints for CONGEN is presented. A general protocol for simulated annealing with restrained molecular dynamics is applied to generate NMR solution structures using CONGEN together with real experimental NMR data. A total of 730 NMR-derived constraints for mEGF and 424 NMR-derived constraints for hTGF alpha were used in these energy-refinement calculations. Different weighting schemes and starting conformations were studied to check and/or improve the sampling of the low-energy conformational space that is consistent with all constraints. The results demonstrate that loosened (i.e., "relaxed") sets of the EGF and hTGF alpha internuclear distance constraints allow molecules to overcome local minima in the search for a global minimum with respect to both distance restraints and conformational energy. The resulting energy-refined structures of mEGF and hTGF alpha are compared with structures determined previously and with structures of homologous proteins determined by NMR and X-ray crystallography.     

Automated NMR structure determination of stereo-array isotope labeled ubiquitin from minimal sets of spectra using the SAIL-FLYA system

Stereo-array isotope labeling (SAIL) has been combined with the fully automated NMR structure determination algorithm FLYA to determine the three-dimensional structure of the protein ubiquitin from different sets of input NMR spectra. SAIL provides a complete stereo- and regio-specific pattern of stable isotopes that results in sharper resonance lines and reduced signal overlap, without information loss. Here we show that as a result of the superior quality of the SAIL NMR spectra, reliable, fully automated analyses of the NMR spectra and structure calculations are possible using fewer input spectra than with conventional uniformly ¹³C/¹⁵N-labeled proteins. FLYA calculations with SAIL ubiquitin, using a single three-dimensional “through-bond” spectrum (and 2D HSQC spectra) in addition to the ¹³C-edited and ¹⁵N-edited NOESY spectra for conformational restraints, yielded structures with an accuracy of 0.83–1.15 Å for the backbone RMSD to the conventionally determined solution structure of SAIL ubiquitin. NMR structures can thus be determined almost exclusively from the NOESY spectra that yield the conformational restraints, without the need to record many spectra only for determining intermediate, auxiliary data of the chemical shift assignments. The FLYA calculations for this report resulted in 252 ubiquitin structure bundles, obtained with different input data but identical structure calculation and refinement methods. These structures cover the entire range from highly accurate structures to seriously, but not trivially, wrong structures, and thus constitute a valuable database for the substantiation of structure validation methods.     

Simulated annealing with restrained molecular dynamics using a flexible restraint potential: theory and evaluation with simulated NMR constraints.

A new functional representation of NMR-derived distance constraints, the flexible restraint potential, has been implemented in the program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168) for molecular structure generation. In addition, flat-bottomed restraint potentials for representing dihedral angle and vicinal scalar coupling constraints have been introduced into CONGEN. An effective simulated annealing (SA) protocol that combines both weight annealing and temperature annealing is described. Calculations have been performed using ideal simulated NMR constraints, in order to evaluate the use of restrained molecular dynamics (MD) with these target functions as implemented in CONGEN. In this benchmark study, internuclear distance, dihedral angle, and vicinal coupling constant constraints were calculated from the energy-minimized X-ray crystal structure of the 46-amino acid polypeptide crambin (ICRN). Three-dimensional structures of crambin that satisfy these simulated NMR constraints were generated using restrained MD and SA. Polypeptide structures with extended backbone and side-chain conformations were used as starting conformations. Dynamical annealing calculations using extended starting conformations and assignments of initial velocities taken randomly from a Maxwellian distribution were found to adequately sample the conformational space consistent with the constraints. These calculations also show that loosened internuclear constraints can allow molecules to overcome local minima in the search for a global minimum with respect to both the NMR-derived constraints and conformational energy. This protocol and the modified version of the CONGEN program described here are shown to be reliable and robust, and are applicable generally for protein structure determination by dynamical simulated annealing using NMR data.     

Membrane‐induced structure of novel human tachykinin hemokinin‐1 (hHK1)

PPT‐C encoded hemokinin‐1(hHK‐1) of Homo sapiens (TGKASQFFGLM) is a structurally distinct neuropeptide among the tachykinin family that participate in the NK‐1 receptor downstream signaling processes. Subsequently, signal transduction leads to execution of various effector functions which includes aging, immunological, and central nervous system (CNS) regulatory actions. However the conformational pattern of ligand receptor binding is unclear. The three‐dimensional structure of the hemokinin‐1 in aqueous and micellar environment has been studied by one and two‐dimensional proton nuclear magnetic resonance (2D 1H‐NMR spectroscopy) and distance geometry calculations. Data shows that hemokinin‐1 was unstructured in aqueous environment; anionic detergent SDS induces α‐helix formation. Proton NMR assignments have been carried out with the aid of correlation spectroscopy (gradient‐COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The inter proton distances and dihedral angle constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (CYANA) to generate a family of structures, which have been refined using restrained energy minimization and dynamics. Helical conformation is observed from residue K3‐M11. The conformational range of the peptide revealed by NMR studies has been analyzed in terms of characteristic secondary features. Observed conformational features have been compared to that of Substance P potent NK1 agonist. Thus the report provides a structural insight to study hHK‐1‐NK1 interaction that is essential for hHK1 based signaling events. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 702–710, 2015.     

Theoretical studies of the structure and molecular dynamics of a peptide crystal

Energy minimizations and molecular dynamics simulations have been performed on the cyclic peptide cyclo-(Ala-Pro-D-Phe)2 in both the isolated and crystal states. The results of these calculations have been analyzed, both to investigate our ability to reproduce experimental data (structure and vibrational and NMR spectra) and to investigate the effects of environment on the energy, structure, and dynamics of peptides. Comparison of the minimized and time-averaged crystal systems with the experimental peptide structure shows that the calculations have closely reproduced the experimental structure. Molecular dynamics of the isolated molecule has led to a new conformation, which is approximately equal to 8.5 kcal/mol more stable than the conformation that exists in the crystal, the latter conformation being stabilized by intermolecular (packing) forces. This illustrates the considerable effect that environment can have on the conformation of peptides. The crystal environment has also been shown to significantly reduce the dynamic conformational fluctuations seen for the isolated molecule. The behavior of the peptide during the isolated simulation also supports the experimental NMR observation of a symmetric structure that differs from the asymmetric, instantaneous structures which characterize the molecule during the dynamics. Calculations of vibrational frequencies of the peptide in the crystal and isolated states show the expected shifts in bond-stretching frequencies due to intermolecular interactions. Finally, we have calculated NMR coupling constants from the dynamics simulation of the isolated peptide and have compared these with the experimental values. This has led to a possible reinterpretation of the experimental data.     

Three-dimensional structure of neuropeptide k bound to dodecylphosphocholine micelles

Neuropeptide K (NPK), an N-terminally extended form of neurokinin A (NKA), represents the most potent and longest lasting vasodepressor and cardiomodulatory tachykinin reported thus far. NPK has been shown to have high selectivity for the NK2 receptor. Because the micelle-associated structure may be relevant to the NPK-receptor interaction, the three-dimensional structure of the NPK in aqueous and micellar environments has been studied by two-dimensional proton nuclear magnetic resonance (2D (1)H NMR spectroscopy) and distance geometry calculations. Proton NMR assignments have been carried out with the aid of correlation spectroscopy (DQF-COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The interproton distances and dihedral angle constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (DYANA) to generate a family of structures, which have been refined using restrained energy minimization and dynamics. The results show that in an aqueous environment NPK lacks a definite secondary structure, although some turn-like elements are present in the N terminus. The structure is well-defined in the presence of dodecylphosphocholine micelles. The global fold of NPK bound to DPC micelles consists of two well-defined helices from residues 9 to 18 and residues 27 to 33 connected by a noncanonical beta turn. The N terminus of the peptide is characterized by a 3(10) helix or a series of dynamic beta turns. The conformational range of the peptide revealed by NMR and circular dichroism (CD) studies has been analyzed in terms of characteristic secondary features. The observed conformational features have been further compared to a NKA and neuropeptide gamma (NPgamma) potent endogenous agonist for the NK2 receptor.     

Knowledge-based versus experimentally acquired distance and angle constraints for NMR structure refinement

Nuclear Overhauser effects (NOE) distance constraints and torsion angle constraints are major conformational constraints for nuclear magnetic resonance (NMR) structure refinement. In particular, the number of NOE constraints has been considered as an important determinant for the quality of NMR structures. Of course, the availability of torsion angle constraints is also critical for the formation of correct local conformations. In our recent work, we have shown how a set of knowledge-based short-range distance constraints can also be utilized for NMR structure refinement, as a complementary set of conformational constraints to the NOE and torsion angle constraints. In this paper, we show the results from a series of structure refinement experiments by using different types of conformational constraints--NOE, torsion angle, or knowledge-based constraints--or their combinations, and make a quantitative assessment on how the experimentally acquired constraints contribute to the quality of structural models and whether or not they can be combined with or substituted by the knowledge-based constraints. We have carried out the experiments on a small set of NMR structures. Our preliminary calculations have revealed that the torsion angle constraints contribute substantially to the quality of the structures, but require to be combined with the NOE constraints to be fully effective. The knowledge-based constraints can be functionally as crucial as the torsion angle constraints, although they are statistical constraints after all and are not meant to be able to replace the latter.     

The structure of d(CGCGAAT[]TCGCG) . d(CGCGAATTCGCG); the incorporation of a thymine photodimer into a B-DNA helix

In the light of the biological significance of thymine photodimers , studies of the energetics of the dodecanucleotide fragment d( CGCGAATTCGCG )2 have been carried out using the methods of molecular mechanics, with and without incorporation of a thymine dimer in the cis-syn configuration. The results of the calculations suggest that the thymine dimerized structures show no gross distortion in the double helix with the conformational changes relative to the normal B-DNA double helix restricted largely to the dimer region. The energetics of dTp[]dT reveal a number of conformers which are energetically almost equally favorable and are, as a group, qualitatively consistent with NMR studies on this molecule. The biological implications of the results of the conformational studies, reported here, have been examined vis-a-vis the currently available models for the recognition of DNA "damage" by repair enzymes.     

Conformational studies on somatostatin. II. The C-terminal hexapeptide fragment

The results of a conformational study on the C terminal hexapeptide of Somatostatin are presented. Semi-empirical energy calculations and high resolution NMR methods have been used to obtain information on the conformational properties of SRIF9-14 in [2H6]dimethylsulfoxide and 2H2O. It is concluded from the energy calculations that the peptide has an averaged conformation in which semi extended and folded structures are important. Only some of the folded conformations can explain the chemical shift differences between the amino acid residues Thr10 and Thr12 as a ring current shift by the Phe11 aromatic ring on Thr10. The nonequivalence is more pronounced in dimethyl-sulfoxide (0.23--0.15 ppm) where it decreases with increasing temperature towards the temperature independent value in 2H2O (0.03 ppm). This suggests that the folded conformations are somewhat predominant in dimethylsulfoxide solutions. In 2H2O the semi extended and folded structures are statistically equally important and the peptide is more flexible. A comparison with a study on the smaller fragments SRIF10-12 and SRIF10-13 which have similar conformational properties, demonstrates the usefulness of the fragment approach in conformational studies of peptides.     

Solution structure of tertiapin determined using nuclear magnetic resonance and distance geometry

The solution structure of tertiapin, a 21-residue bee venom peptide, has been characterized by circular dichroism (CD), two-dimensional nuclear magnetic resonance (NMR) spectroscopy, and distance geometry. A total of 21 lowest error structures were obtained from distance geometry calculations. Superimposition of these structures shows that the backbone of tertiapin is very well defined. One type-I reverse turn from residue 4 to 7 and an α-helix from residue 12 to 19 exist in the structure of tertiapin. The α-helical region is best defined from both conformational analysis and structural superimposition. The overall three-dimensional structure of tertiapin is highly compact resulting from side chain interactions. The structural information obtained from CD and NMR are compared for both tertiapin and apamin (ref. 3), another bee venom peptide. Tertiapin and apamin have some similar secondary structure, but display different tertiary structures. © 1993 Wiley-Liss, Inc.     

Correlation between changes in nuclear magnetic resonance order parameters and conformational entropy: molecular dynamics simulations of native and denatured staphylococcal nuclease

Recent work has suggested that changes in NMR order parameters may quantitatively reflect changes in the conformational entropy of a protein ensemble. The extent of the mathematical relationship between local entropy changes as seen by NMR order parameters and the full protein entropy change is a complex issue. As a step towards a fuller understanding of this problem, molecular dynamics calculations of both native and denatured staphylococcal nuclease were performed. The N-H bond vector motion, in both explicit and implicit solvent, was analyzed to estimate local and global entropy changes. The calculated N-H bond vector order parameters from simulation agreed on average with experimental values for both native and denatured structures. However, the inverted-U profile of order parameters versus residue number observed experimentally for denatured nuclease was only partially reproduced by simulation of compact denatured structures. Comparisons made across the full set of simulations revealed a correlation between the N-H order parameter-based conformational entropy change and the total quasiharmonic-based conformational entropy change between the native and denatured structures. The calculations showed that about 25% of the total entropy change was reflected by changes in simulated S2 values. This result suggests that NMR-derived order parameters may be used to provide a reasonable estimate of the total conformational entropy change on protein folding.     

Structural Characterization of Neurokinin-3 Receptor Selective Peptide Agonist Scyliorhinin II Bound to DPC Micelles

Abstract

Scyliorhinin II, a cyclic Tachykinin peptide, is a potent NK3 receptor agonist. The pharmacology of NK3 receptor is least characterized out of the three tachykinin receptor subtypes cloned and characterized for Tachykinins. To understand the structural basis of peptide-receptor interaction, the three-dimensional structure of the Scyliorhinin II in aqueous and micellar environments has been studied by two-dimensional proton nuclear magnetic resonance (2D ¹H-NMR spectroscopy) and distance geometry calculations. Proton NMR assignments have been carried out with the aid of correlation spectroscopy (gradient-COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The inter proton distances and dihedral angle constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (DYANA) to generate a family of structures, which have been refined using restrained energy minimization and dynamics. The results show that in an aqueous environment, Scyliorhinin II lacks a definite secondary structure. The structure is well-defined in presence of dodecyl phosphocholine micelles. The global fold of Scyliorhinin II bound to DPC micelles consists of a well-defined helix in the C-terminal region from residue 12–18 and a series of turns towards N-terminus. The structure is further stabilized by disulfide bond between Cys7 and Cys13. The conformational range of the peptide revealed by NMR and CD studies has been analyzed in terms of characteristic secondary features. Observed conformational features have been compared with those of Substance P, Neurokinin A and Neurokinin B, potent NK1, NK2, and NK3 agonists, respectively.     

Bioactive sesquiterpenes from Santolina rosmarinifolia subsp. Canescens. A conformational analysis of the germacrane ring.

The hexane extract of aerial parts of Santolina rosmarinifolia subsp. canescens afforded eight new sesquiterpenes in addition to known compounds. Their structures were determined by spectroscopic methods and chemical transformations. The conformational analysis of the germacrane constituents was carried out by spectroscopic methods, including NMR at varying temperature and by molecular mechanics calculations. The antifeedant, antibacterial and antitumoral activity of selected compounds has been tested.     

Limited sampling of conformational space by the distance geometry algorithm: implications for structures generated from NMR data

Calculations with a metric matrix distance geometry algorithm were performed that show that the standard implementation of the algorithm generally samples a very limited region of conformational space. This problem is most severe when only a small amount of distance information is used as input for the algorithm. Control calculations were performed on linear peptides, disulfide-linked peptides, and a double-stranded DNA decamer where only distances defining the covalent structures of the molecules (as well as the hydrogen bonds for the base pairs in the DNA) were included as input. Since the distance geometry algorithm is commonly used to generate structures of biopolymers from distance data obtained from NMR experiments, simulations were performed on the small globular protein basic pancreatic trypsin inhibitor (BPTI) that mimic calculations performed with actual NMR data. The results on BPTI and on the control peptides indicate that the standard implementation of the algorithm has two main problems: first, that it generates extended structures; second, that it has a tendency to consistently produce similar structures instead of sampling all structures consistent with the input distance information. These results also show that use of a simple root-mean-square deviation for evaluating the quality of the structures generated from NMR data may not be generally appropriate. The main sources of these problems are identified, and our results indicate that the problems are not a fundamental property of the distance geometry algorithm but arise from the implementations presently used to generate structures from NMR data. Several possible methods for alleviating these problems are discussed.     

Assessing the acid-base and conformational properties of histidine residues in human prion protein (125-228) by means of pK(a) calculations and molecular dynamics simulations

A thorough study of the acid-base behavior of the four histidines and the other titratable residues of the structured domain of human prion protein (125-228) is presented. By using multi-tautomer electrostatic calculations, average titration curves have been built for all titratable residues, using the whole bundles of NMR structures determined at pH 4.5 and 7.0. According to our results, (1) only histidine residues are likely to be involved in the first steps of the pH-driven conformational transition of prion protein; (2) the pK(a)'s of His140 and His177 are approximately 7.0, whereas those of His155 and His187 are < 5.5. 10-ns long molecular dynamics simulations have been performed on five different models, corresponding to the most significant combinations of histidine protonation states. A critical comparison between the available NMR structures and our computational results (1) confirms that His155 and His187 are the residues whose protonation is involved in the conformational rearrangement of huPrP in mildly acidic condition, and (2) shows how their protonation leads to the destructuration of the C-terminal part of HB and to the loss of the last turn of HA that represent the crucial microscopic steps of the rearrangement.     

Improving NMR protein structure quality by Rosetta refinement: a molecular replacement study

The structure of human protein HSPC034 has been determined by both solution nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography. Refinement of the NMR structure ensemble, using a Rosetta protocol in the absence of NMR restraints, resulted in significant improvements not only in structure quality, but also in molecular replacement (MR) performance with the raw X-ray diffraction data using MOLREP and Phaser. This method has recently been shown to be generally applicable with improved MR performance demonstrated for eight NMR structures refined using Rosetta (Qian et al., Nature 2007;450:259-264). Additionally, NMR structures of HSPC034 calculated by standard methods that include NMR restraints have improvements in the RMSD to the crystal structure and MR performance in the order DYANA, CYANA, XPLOR-NIH, and CNS with explicit water refinement (CNSw). Further Rosetta refinement of the CNSw structures, perhaps due to more thorough conformational sampling and/or a superior force field, was capable of finding alternative low energy protein conformations that were equally consistent with the NMR data according to the Recall, Precision, and F-measure (RPF) scores. On further examination, the additional MR-performance shortfall for NMR refined structures as compared with the X-ray structure were attributed, in part, to crystal-packing effects, real structural differences, and inferior hydrogen bonding in the NMR structures. A good correlation between a decrease in the number of buried unsatisfied hydrogen-bond donors and improved MR performance demonstrates the importance of hydrogen-bond terms in the force field for improving NMR structures. The superior hydrogen-bond network in Rosetta-refined structures demonstrates that correct identification of hydrogen bonds should be a critical goal of NMR structure refinement. Inclusion of nonbivalent hydrogen bonds identified from Rosetta structures as additional restraints in the structure calculation results in NMR structures with improved MR performance.     

Quantum chemical studies on the conformational structure of nucleic acids. IV. Calculation of backbone structure by CNDO method

Quantum chemical calculations using the CNDO/2 method, have been carried out to determine the energetically favoured ranges of the torsional angles (φ′, ω′, ω, φ, ψ) which fix the conformational structure of nucleic acid backbone. The two dimensional isoenergy maps have been constructed in the (ω′, ω) and (φ, ψ) hyperspaces. The variation of total energy with respect to φ′ has also been studied. The results show that the non-bonding interactions play a major role in the conformational stability of nucleic acids and polynucleotides. The theoretical predictions show good correspondence with the experimental data (X-ray and ¹³C NMR) as well as the other reported theoretical calculations (EHT, PCILO and classical potential functions). The most favoured structure has the conformational angles close to 240, 290, 290, 180 and 60° and these values lead to a helical structure with a pitch of 34 Å and about ten nucleotide units per turn of the helix. The proposed models of Watson &; Crick, DNA-B and DNA-C lie in high energy regions.     

Elucidation of the structure of constrained bicyclopeptides in solution by two-dimensional cross-relaxation spectroscopy: Amatoxin analogues

The evaluation of peptide structures in solution is made feasible by the combined use of two-dimensional NMR in the laboratory (NOESY) and rotating frames (ROESY), and by the use of molecular dynamics calculations. The present paper describes how both the NMR method and molecular dynamics calculations were applied to very rigid synthetic bicyclic peptides that are analogues of natural amatoxins. The NMR theory, which allows the estimate of interatomic distances between interacting nuclei, is briefly discussed. The experimental data were compared with those of known solid-state structures. Three amatoxin analogues have been examined. Of these, one is biologically active (S-deoxo γ[R] OH-Ile³-amaninamide) and its structure in the solid state has recently been worked out. The second and third analogues (S-deoxo-Ile³ -Ala⁵-amaninamide and S-deoxo-D -Ile³ -amaninamide, respectively) are inactive and their solid-state structures are unknown. The data presented confirm the authors' previous hypothesis that lack of biological activity of S-deoxo-Ile³-Ala⁵- amaninamide is due to the masking of the tryptophan ring by the methyl group of L -Ala and not to massive conformational changes of the analogue.