Local structure of cytochrome c from horse heart in solution. Conformational analysis using data of two-dimensional nuclear Overhauser effect spectroscopy]

Using the earlier suggested method the calculation of the backbone conformations of horse heart cytochrome c in oxidized (ferricytochrome c) and reduced (ferrocytochrome c) states has been performed by the two-dimensional nuclear Overhauser effect spectroscopy data. For both protein forms the secondary structure elements have been revealed and the conformations of the irregular polypeptide chain segments have been analysed. The similarity of the secondary structures of ferri- and ferrocytochrome c in solution was established from the comparison of their conformations. Small differences between the conformations of two molecule forms are shown to be localized within the polypeptide chain fragments situated in the spatial structure near the heme crevice. The comparison of the dihedral phi and psi angles in the calculated conformations of horse cytochrome C with the corresponding characteristics of X-ray structures of tuna ferri- and ferrocytochrome c made for the oxidized and reduced protein forms using the quantitative criteria testifies the similarity of their conformations in solution and crystal. In is shown that the conformational changes of the separate amino acid residues which take place as the result of the "solution-to-crystal" transition occur on the surface fragments of protein globule and do not lead to essential alterations of the secondary molecule structure.     

Flexible ligand docking using conformational ensembles.

Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.     

Conformational features responsible for the binding of cyclic analogues of enkephalin to opioid receptors. III. Probable binding conformations of mu-agonists with phenylalanine in position 3.

Theoretical conformational analysis was carried out for several tetrapeptide analogues of beta-casomorphin and dermorphin containing a Phe residue in position 3. Sets of low-energy backbone structures of the mu-selective peptides [N-Me-Phe3, D-Pro4]-morphiceptin and Tyr-D-Orn-Phe-Asp-NH2 were obtained. These sets of structures were compared for geometrical similarity between themselves and with the low-energy conformations found for the delta-selective peptide Tyr-D-Cys-Phe-D-Pen-OH and nonactive peptide Tyr-Orn-Phe-Asp-NH2. Two pairs of geometrically similar conformations of mu-selective peptides, sharing no similarity with the conformations of peptides showing low affinity to the mu-receptor, were selected as two alternative models of probable mu-receptor-bound backbone conformations. Both models share geometrical similarity with the low-energy structures of the linear mu-selective peptide Tyr-D-Ala-Phe-Phe-NH2. Putative binding conformations of Tyr1 and Phe3 side chains are also discussed.     

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.     

Conformational similarity among amino acid residues: 1. Analysis of protein crystal structure data

The probability distribution in the (?,ψ)-plane obtained for each amino acid residue from cyrstal structure data of globular proteins is compared. This has shown amino acid residues. Pro and Gly to be conformationally unique. Conformational similarity in the (?,ψ)-plane of amino acid reced does not necessarily mean that they will have the same chemical or biochemical properties or similar secondary structures. A set of amino acid residues are given which can adopt the conformations of other amino acid residues without much difficulty either in the whole (?,ψ)-plane or in regions, where the observed conformations are maximum.     

Crystal structure of a human rhinovirus that displays part of the HIV-1 V3 loop and induces neutralizing antibodies against HIV-1

We report the 2.7 A resolution structure of a chimeric rhinovirus, MN-III-2, that displays part of the HIV-1 gp120 V3 loop and elicits HIV-neutralizing antibodies. The V3 loop insert is dominated by two type I beta turns. The structures of two adjacent tripeptides resemble those of analogous segments in three Fab/V3 loop peptide complexes. Although two of the three corresponding antibodies bind and neutralize MN-III-2 well, only one of the three can bind without significant rearrangement. These results suggest that the V3 loop insert: (1) can share some local conformational similarity to V3 loop sequences presented on different structural frameworks; (2) must be able to adopt multiple conformations, even in a relatively constrained environment; and (3) may mimic the conformational variability of the epitope on HIV-1, increasing the likelihood of eliciting appropriate neutralizing immune responses.     

NMR studies of defensin antimicrobial peptides. 1. Resonance assignment and secondary structure determination of rabbit NP-2 and human HNP-1.

Two-dimensional nuclear magnetic resonance spectroscopy has been used to make resonance assignments of the proton spectra of two defensin antimicrobial peptides, human neutrophil peptide HNP-1 and rabbit neutrophil peptide NP-2. The secondary structures of these peptides were determined from analysis of the proton-proton NOEs and from the positions of slowly exchanging amide protons. Both peptides contain a long stretch of a double-stranded antiparallel beta-sheet in a hairpin conformation that contains a beta-bulge, a short region of triple-stranded beta-sheet, and several tight turns. The NMR results clearly show that HNP-1 forms a dimer or higher order aggregate in solution and that Pro8 exists as a cis peptide bond. The NMR data on these peptides are compared with NMR data for a homologous peptide NP-5 [Bach, A. C., Selsted, M. E., & Pardi, A. (1987) Biochemistry 26, 4389-4397]. Analysis of the conformation-dependent proton chemical shifts shows that it is not possible to confidently judge the structural similarity of the three defensins from chemical shift data alone. However, comparison of the 3JHN alpha coupling constants in NP-2 and NP-5 indicates that the backbone conformations for these peptides are very similar. A more detailed comparison of the solution conformations of the defensins peptides is made in the following paper in this issue where the NMR data are used as input for distance geometry and molecular dynamics calculations to determine the three-dimensional structures of HNP-1 and NP-2.     

Computational ligand-based rational design: Role of conformational sampling and force fields in model development

A significant number of drug discovery efforts are based on natural products or high throughput screens from which compounds showing potential therapeutic effects are identified without knowledge of the target molecule or its 3D structure. In such cases computational ligand-based drug design (LBDD) can accelerate the drug discovery processes. LBDD is a general approach to elucidate the relationship of a compound's structure and physicochemical attributes to its biological activity. The resulting structure-activity relationship (SAR) may then act as the basis for the prediction of compounds with improved biological attributes. LBDD methods range from pharmacophore models identifying essential features of ligands responsible for their activity, quantitative structure-activity relationships (QSAR) yielding quantitative estimates of activities based on physiochemical properties, and to similarity searching, which explores compounds with similar properties as well as various combinations of the above. A number of recent LBDD approaches involve the use of multiple conformations of the ligands being studied. One of the basic components to generate multiple conformations in LBDD is molecular mechanics (MM), which apply an empirical energy function to relate conformation to energies and forces. The collection of conformations for ligands is then combined with functional data using methods ranging from regression analysis to neural networks, from which the SAR is determined. Accordingly, for effective application of LBDD for SAR determinations it is important that the compounds be accurately modelled such that the appropriate range of conformations accessible to the ligands is identified. Such accurate modelling is largely based on use of the appropriate empirical force field for the molecules being investigated and the approaches used to generate the conformations. The present chapter includes a brief overview of currently used SAR methods in LBDD followed by a more detailed presentation of issues and limitations associated with empirical energy functions and conformational sampling methods.     

Structural characterization of manganese(II)-nucleotide complexes bound to yeast 3-phosphoglycerate kinase: 13C relaxation measurements using [U-13C]ATP and [U-13C]ADP

A complete characterization of the conformations of Mn.ADP and Mn.ATP bound to the active site of yeast 3-P-glycerate kinase is presented. These conformations have been deduced on the basis of paramagnetic effects on 13C spin-lattice relaxation rates in [U-13C]nucleotides due to Mn(II), used as a substituent activating cation. The 13C relaxation measurements were performed on exclusively enzyme-bound complexes E.Mn.[U-13C]ATP and E.Mn.[U-13C]ADP at three distinct 13C NMR frequencies: 75.4, 125.7, and 181 MHz. The frequency dependence of the relaxation data has been analyzed in an effort to evaluate distances from the cation for all 10 13C nuclei in the adenosine moieties of E.Mn.ATP and E.Mn.ADP. These distance data, taken along with previously published cation-31P distances, have been used as constraints in the molecular modeling program Quanta, in which molecular dynamics simulations and energy minimization have been performed to determine the conformations that are compatible with the distance data. It was possible to model the distances on the basis of a single enzyme-bound conformation for each of the nucleotides. The details of the enzyme-bound Mn.ATP and Mn.ADP conformations are distinguishably different from each other, indicating that structural alterations occur in the enzyme-bound reaction complex as the enzyme turns over. For example, when the adenosine moieties in the bound structures of Mn.ATP and Mn.ADP are superposed, the cation is found to be displaced by approximately 2.4 A between the two conformations, suggesting that these structural changes may involve movements with significant amplitudes. Furthermore, the NMR-determined structures of enzyme-bound Mn.ATP and Mn.ADP are significantly different from those in published X-ray crystal structures of the enzyme-nucleotide complexes.     

The utility of residual dipolar couplings in detecting motion in carbohydrates: application to sucrose

The solution structure and dynamics of sucrose are examined using a combination of NMR residual dipolar coupling and molecular mechanics force fields. It is found that the alignment tensors of the individual rings are different, and that fitting 35 measured residual dipolar couplings to structures with specific phi, psi values indicates the presence of three major conformations: phi, psi=(120 degrees ,270 degrees), (45 degrees, 300 degrees) and (90 degrees ,180 degrees). Furthermore, fitting two structures simultaneously to the 35 residual dipolar couplings results in a substantial improvement in the fits. The existence of multiple conformations having similar stabilities is a strong indication of motion, due to the interconversion among these states. Results from four molecular mechanics force fields are in general agreement with the experimental results. However, there are major disagreements between force fields. Because fits of residual dipolar couplings to structures are dependent on the force field used to calculate the structures, multiple force fields were used to interpret NMR data. It is demonstrated that the pucker of the fructofuranosyl ring affects the calculated potential energy surface, and the fit to the residual dipolar couplings data. Previously published 13C nuclear relaxation results suggesting that sucrose is rigid are not inconsistent with the present results when motional timescales are considered.     

Raman optical activity characterization of native and molten globule states of equine lysozyme: comparison with hen lysozyme and bovine alpha-lactalbumin

Vibrational Raman optical activity (ROA) spectra of the calcium-binding lysozyme from equine milk in native and nonnative states are measured and compared with those of the homologous proteins hen egg white lysozyme and bovine alpha-lactalbumin. The ROA spectrum of holo equine lysozyme at pH 4.6 and 22 degrees C closely resembles that of hen lysozyme in regions sensitive to backbone and side chain conformations, indicating similarity of the overall secondary and tertiary structures. However, the intensity of a strong positive ROA band at approximately 1340 cm(-1), which is assigned to a hydrated form of alpha helix, is more similar to that in the ROA spectrum of bovine alpha-lactalbumin than hen lysozyme and may be associated with the greater flexibility and calcium-binding ability of equine lysozyme and bovine alpha-lactalbumin compared with hen lysozyme. In place of a strong sharp positive ROA band at approximately 1300 cm(-1) in hen lysozyme that is assigned to an alpha helix in a more hydrophobic environment, equine lysozyme shows a broader band centered at approximately 1305 cm(-1), which may reflect greater heterogeneity in some alpha-helical sequences. The ROA spectrum of apo equine lysozyme at pH 4.6 and 22 degrees C is almost identical to that of the holo protein, which indicates that loss of calcium has little influence on the backbone and side chain conformations, including the calcium-binding loop. From the similarity of their ROA spectra, the A state at pH 1.9 and both 2 and 22 degrees C and the apo form at pH 4.5 and 48 degrees C, which are partially folded denatured (molten globule or state A) forms of equine lysozyme, have similar structures that the ROA suggests contain much hydrated alpha helix. The A state of equine lysozyme is shown by these results to be more highly ordered than that of bovine alpha-lactalbumin, the ROA spectrum of which has more features characteristic of disordered states. A positive tryptophan ROA band at approximately 1551 cm(-1) in the native holo protein disappears in the A state, which is probably due to the presence of nonnative conformations of the tryptophans associated with a previously identified cluster of hydrophobic residues.     

Characterizing and predicting the functional and conformational diversity of seven-transmembrane proteins

The activation of seven-transmembrane receptors (7TMRs) allows cells to sense their environment and convert extracellular signals (like hormone binding) into intracellular signals (through G protein-coupled and/or β arrestin-coupled pathways). A single 7TMR is capable of transducing a wide spectrum of physiological responses inside a cell by coupling to these pathways. This intracellular pleiotropic action is enabled by multiple conformations exhibited by these receptors. Developments in membrane protein structure determination technologies have led to a rapid increase in crystal structures for many 7TMRs. Majority of these receptors have been crystallized in their inactive conformation and, for some, one of the many active conformations has also been crystallized. Given the topological constraints of a lipid bilayer that results in a single fold of seven almost parallel TM helices connected by mostly unstructured loops, these structures exhibit a diversity of conformations not only across the receptors but also across the different functional forms for receptors with structures for one of the functionally active conformations. Here we present a method to characterize this conformational diversity in terms of transmembrane helix topology (TMHTOP) parameters and how to use these helix orientation parameters to predict functionally-distinct multiple conformations for these receptors. The TMHTOP parameters enable a quantification of the structural changes that underlie 7TMR activation and also sheds a unique mechanistic light on the pleiotropic nature of these receptors. It provides a common language to describe the 7TMR activation mechanisms as well as differences across many receptors in terms of visually intuitive structural parameters. Protein structure prediction methods can use these parameters to describe 7TMR conformational ensembles, which coupled to experimental data can be used to develop testable hypotheses for the structural basis of 7TMR functions.     

Nuclear magnetic resonance solution and X-ray structures of squash trypsin inhibitor exhibit the same conformation of the proteinase binding loop

A comparison of the solution nuclear magnetic resonance (n.m.r.) structures of squash trypsin inhibitor from seeds of the squash Cucurbita maxima with the X-ray structure of a trypsin complex of the inhibitor shows that the n.m.r. and X-ray structures are similar in terms of the global folding and secondary structure. The average atomic root-mean-square difference between the 36 n.m.r. structures on the one hand and the X-ray structure is 0.96 A for the backbone atoms and 1.95 A for all heavy atoms. The n.m.r. and X-ray structures exhibit extremely similar conformations of the primary proteinase binding loop. Despite the overall similarity, there are small differences between the mean computed structure and the X-ray structure. The n.m.r. structures have slightly different positions of the segments from residues 16 to 18, and 24 and 25. The n.m.r. results show that the inclusion of stereospecific assignments and precise distance constraints results in a significant improvement in the definition of the n.m.r. structure, making possible a detailed analysis of the local conformations in the protein.     

Conformational study of cyclolinopeptide A. A distance geometry and molecular dynamics approach

The conformation of cyclolinopeptide A, c(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val), a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of distance geometry calculations and molecular dynamics simulations. The starting points for all the calculations were an X-ray structure and other structures obtained from distance geometry calculations based on NMR data. Restrained and unrestrained molecular dynamics simulations are reported in vacuo and in CCl4. Structural and dynamic properties are investigated and compared with those experimentally determined. The conformation obtained from the MD simulations which best reproduces the NMR parameters is at the same time one of the most stable ones and is also fairly similar to the crystal structure. An explanation for the occurrence of multiple conformations in solution at room temperature is given.     

Gaussian-weighted RMSD superposition of proteins: a structural comparison for flexible proteins and predicted protein structures

Many proteins contain flexible structures such as loops and hinged domains. A simple root mean square deviation (RMSD) alignment of two different conformations of the same protein can be skewed by the difference between the mobile regions. To overcome this problem, we have developed a novel method to overlay two protein conformations by their atomic coordinates using a Gaussian-weighted RMSD (wRMSD) fit. The algorithm is based on the Kabsch least-squares method and determines an optimal transformation between two molecules by calculating the minimal weighted deviation between the two coordinate sets. Unlike other techniques that choose subsets of residues to overlay, all atoms are included in the wRMSD overlay. Atoms that barely move between the two conformations will have a greater weighting than those that have a large displacement. Our superposition tool has produced successful alignments when applied to proteins for which two conformations are known. The transformation calculation is heavily weighted by the coordinates of the static region of the two conformations, highlighting the range of flexibility in the overlaid structures. Lastly, we show how wRMSD fits can be used to evaluate predicted protein structures. Comparing a predicted fold to its experimentally determined target structure is another case of comparing two protein conformations of the same sequence, and the degree of alignment directly reflects the quality of the prediction.     

The influenza fusion peptide adopts a flexible flat V conformation in membranes

Knowledge about the influenza fusion peptide (FP) membrane insertion mode is crucial for understanding its fusogenic mechanism. NMR and electron paramagnetic resonance experiments showed that in micelles, the FP inserted as a fixed-angle inverted V. In membranes, however, it was shown to insert as a straight α-helix (by molecular-dynamics simulations) and to adopt multiple kinked conformations (by solid-state NMR). In this work we performed explicit-solvent molecular-dynamics simulations of the influenza FP, and its F9A and W14A mutants, in POPC membranes. The Hα1 chemical shifts predicted from the molecular-dynamics structures are in excellent agreement with the experimental values obtained for the three peptides. The peptide orientation and conformations observed from the simulations lead to a flexible flat-V model in which the peptide lies almost flat on the membrane surface and alternates between kinked and straight-helix conformations.     

Multiple helical conformations of the helix‐turn‐helix region revealed by NOE‐restrained MD simulations of tryptophan aporepressor,TrpR

The nature of flexibility in the helix‐turn‐helix region of E. coli trp aporepressor has been unexplained for many years. The original ensemble of nuclear magnetic resonance (NMR structures showed apparent disorder, but chemical shift and relaxation measurements indicated a helical region. Nuclear Overhauser effect (NOE) data for a temperature‐sensitive mutant showed more helical character in its helix‐turn‐helix region, but nevertheless also led to an apparently disordered ensemble. However, conventional NMR structure determination methods require all structures in the ensemble to be consistent with every NOE simultaneously. This work uses an alternative approach in which some structures of the ensemble are allowed to violate some NOEs to permit modeling of multiple conformational states that are in dynamic equilibrium. Newly measured NOE data for wild‐type aporepressor are used as time‐averaged distance restraints in molecular dynamics simulations to generate an ensemble of helical conformations that is more consistent with the observed NMR data than the apparent disorder in the previously reported NMR structures. The results indicate the presence of alternating helical conformations that provide a better explanation for the flexibility of the helix‐turn‐helix region of trp aporepressor. Structures representing these conformations have been deposited with PDB ID: 5TM0. Proteins 2017; 85:731–740. © 2016 Wiley Periodicals, Inc.     

A catalogue of phi, psi and omega torsion angles for tetrapeptides in proteins

This work examined the degree of consistency, over multiple occurrences in several proteins, of conformations of tetrapeptides. A data base of 114 proteins was taken from the Brookhaven Data Bank. The number of occurrences of every tetrapeptide in these proteins was counted. For each tetrapeptide which occurred more than 6 times all phi, psi and omega angles were calculated and averaged over the occurrences. Due to the limitations of the data base there were only 24 tetrapeptides in this category. This does not provide sufficient information for the extensive prediction of structures which are not homologous to those already in the data base but the effectiveness of structure prediction can be expected to increase with the size of the data base.     

Detection of native-like models for amino acid sequences of unknown three-dimensional structure in a data base of known protein conformations.

We present an approach which can be used to identify native-like folds in a data base of protein conformations in the absence of any sequence homology to proteins in the data base. The method is based on a knowledge-based force field derived from a set of known protein conformations. A given sequence is mounted on all conformations in the data base and the associated energies are calculated. Using several conformations and sequences from the globin family we show that the native conformation is identified correctly. In fact the resolution of the force field is high enough to discriminate between a native fold and several closely related conformations. We then apply the procedure to several globins of known sequence but unknown three dimensional structure. The homology of these sequences to globins of known structures in the data base ranges from 49 to 17%. With one exception we find that for all globin sequences one of the known globin folds is identified as the most favorable conformation. These results are obtained using a force field derived from a data base devoid of globins of known structure. We briefly discuss useful applications in protein structural research and future development of our approach.     

Serpins in the Caenorhabditis elegans genome.

Data mining in genome sequences can identify distant homologues of known protein families, and is most powerful if solved structures are available to reveal the three-dimensional implications of very dissimilar sequences. Here we describe putative serpin sequences identified with very high statistical significance in the Caenorhabditis elegans genome. When mapped onto vertebrate serpins such as alpha1-antitrypsin, they suggest novel structural features. Some appear complete, some show extensive deletions, and others appear to contain only the C-terminal part of the known serpin fold, probably in partnership with N-terminal regions that have conformations unlike those of known serpins. The observation of such striking sequence similarity, in proteins that must have significantly different overall structures, substantially extends the structural characteristics of the serpin family of proteins.