Solution conformations of amphidinolide H

Solution conformations of amphidinolide H (1), a 26-membered macrolide exhibiting potent cytotoxic and antitumor activity, in CDCl3 and DMSO-d6 were investigated on the basis of NMR data, distance geometry calculation, and restrained energy minimization. Three-dimensional conformations in CDCl3 were suggested to be close to the X-ray structure of 1, while those in DMSO-d6 were indicated to be different from both those in CDCl3 and the X-ray structure.     

The task of determining the local structure of proteins by Overhauser effect nuclear spectroscopy addressed once more]

The testing of the earlier developed theoretical method for determining the backbone protein conformations (the local structure) on the basis of the two-dimensional nuclear Overhauser effect (NOE) spectroscopy has been fulfilled. The method approval has been carried out by the calculation (based upon spectral NOE parameters) of the local plastocyanin and bovine pancreatic trypsin inhibitor structures followed by the comparison of the received conformational parameters with the X-ray data. The comparison of the molecular conformations in solution and crystal has been implemented for different fragments of the polypeptide chain (beta-structures, alpha-helices, irregular segments) using the mathematical statistics methods. The verification of the "zero" hypothesis about the similarity of phi and psi variation rows which was carried out at the reliability level of 0.99 showed that in both cases there were no systematic deviations of dihedral angles of the compared conformations and that their dispersion differences were statistically indiscernible. It has been concluded that the approved method permits to determine the local structure of the conformationally rigid proteins (or their fragments) at the level close to that which provides the high resolution X-ray analysis.     

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.     

Comparative conformational analysis of cholesterol and ergosterol by molecular mechanics

A comparative conformational analysis of cholesterol and ergosterol has been carried out using molecular mechanics methods. These studies are aimed at giving a better understanding of the molecular nature of the interaction of these sterols with polyene macrolide antibiotics. Structures of cholesterol and ergosterol determined by X-ray methods have been used as initial geometries of these molecules for force field calculations. The calculation of steric energy has also been made for conformations which do not appear in the crystal. The latter conformers have different conformations of the side chain as well as different conformations of rings A and D. The rotational barriers around bonds C17–C20 and C20–C22 have also been calculated. The results obtained on differences and similarities in the conformations of cholesterol and ergosterol allow us to postulate a mechanism for differential interaction with the antibiotics. The relatively rigid side chain of ergosterol (stretched molecule) in comparison with the flexible side chain of cholesterol (bent molecule), allows better intermolecular contact of the first sterol molecule with a polyene macrolide and in consequence facilitates complex formation involving Van der Waal's forces.     

Assembly of polypeptide and protein backbone conformations from low energy ensembles of short fragments: development of strategies and construction of models for myoglobin, lysozyme, and thymosin beta 4.

Recently we developed methods for the construction of knowledge-based mean fields from a data base of known protein structures. As shown previously, this approach can be used to calculate ensembles of probable conformations for short fragments of polypeptide chains. Here we develop procedures for the assembly of short fragments to complete three-dimensional models of polypeptide chains. The amino acid sequence of a given protein is decomposed into all possible overlapping fragments of a given length, and an ensemble of probable conformations is calculated for each fragment. The fragments are assembled to a complete model by choosing appropriate conformations from the individual ensembles and by averaging over equivalent angles. Finally a consistent model is obtained by rebuilding the conformation from the average angles. From the average angles the local variability of the structure can be calculated, which is a useful criterion for the reliability of the model. The procedure is applied to the calculation of the local backbone conformations of myoglobin and lysozyme whose structures have been solved by X-ray analysis and thymosin beta 4, a polypeptide of 43 amino acid residues whose structure was recently investigated by NMR spectroscopy. We demonstrate that substantial fractions of the calculated local backbone conformations are similar to the experimentally determined structures.     

The X-ray investigation of 16 alpha,17 alpha-thiazolidine derivatives of delta 4- and delta 5-pregnanes and conformational analysis of their oxathiolane analog

An X-ray study of 3,20-dioxo-4-pregnene-[16 alpha,17 alpha-d]--2',2'-dimethylthiazolidine (I) and 3 beta-hydroxy-20-oxo-5--pregnene-[16 alpha,17 alpha-d]-2',2'- dimethylthiazolidine (II) has been carried out. Two independent molecules in crystal II have significantly different conformations of the D and E rings, although according to the atom-atom potential calculations the energy of interaction of these molecules with their neighbors in crystal is the same. The calculation of conformational energy of 3,20-dioxo-4-pregnene-[17 alpha,16 alpha-d]-2',2'--dimethyloxathiolane (III) by the molecular mechanics method (MMM) indicates a possibility of existence of two similar conformers also for this molecule. The MMM calculation shows also that the conformation of molecule III (as well as progesterone) with the 17 beta-acetyl group torsion angle C(16)C(17)C(20)0(20) close to -120 degrees is possible.     

From determination of amino acid residue conformation to reconstruction of the spatial structure of proteins (from the data of nuclear Overhauser enhancement spectroscopy

A new approach to the calculation of the spatial protein structure based on the joint utilization of the theoretical conformational analysis method and nuclear Overhauser enhancement (NOE) spectroscopy data is proposed and verified. The quality in determining various molecule structural parameters is estimated in terms of the expected NOE spectral parameters derived from the X-ray analysis data of the avian pancreatic polypeptide. The proposed approach is shown to correctly determine such structural parameters of protein molecules as local amino acid residue conformations, reciprocal spatial orientation of the C alpha atoms neighbouring along amino acid sequence and reapproached segments of the polypeptide chain. Spatially remote molecule fragments are mainly responsible for the error in determining structural parameters.     

X-ray fibre diffraction of cartilage proteoglycan aggregates containing hyaluronic acid (Short Communication)

Ordered conformations of proteoglycan-hyaluronic acid aggregates in the intercellular matrix in cartilage were observed by X-ray diffraction. The sodium salt form of three samples, (a) aggregated proteoglycan, (b) disaggregated proteoglycan and (c) reconstituted disaggregated proteoglycan, give essentially similar X-ray fibre-type diffraction photographs. The patterns correlate with the chondroitin 4-sulphate component and can be interpreted as twofold helical conformations, similar to that observed previously for the free acid form of chondroitin 4-sulphate (Isaac & Atkins, 1973). The information takes us one step nearer the situation found in cartilage in vivo.     

Bayesian estimation of NMR restraint potential and weight: a validation on a representative set of protein structures

The classical procedure for nuclear magnetic resonance structure calculation allocates empirical distance ranges and uses historical values for weighting factors. However, Bayesian analysis suggests that there are more optimal choices for potential shape (bounds-free log-harmonic shape) and restraints weights. We compare the classical protocol with the Bayesian approach for more than 300 protein structures. We analyze the conformation similarity to the corresponding X-ray crystal structure, the distribution of the conformations around their average, and independent validation criteria. On average, the log-harmonic potential reduces the difference to the X-ray crystal structure. If the log-harmonic potential is used, the constant weighting tightens the distribution around the average conformation, with respect to the distributions obtained with Bayesian weighting. Conversely, the structure quality is improved by the Bayesian weighting over the classical procedure, whereas constant weighting worsens some criteria. The quality improvement obtained with the log-harmonic potential coupled to Bayesian weighting validates this approach on a representative set of protein structures.     

Classifying protein kinase structures guides use of ligand-selectivity profiles to predict inactive conformations: structure of lck/imatinib complex

We report a clustering of public human protein kinase structures based on the conformations of two structural elements, the activation segment and the C-helix, revealing three discrete clusters. One cluster includes kinases in catalytically active conformations. Each of the other clusters contains a distinct inactive conformation. Typically, kinases adopt at most one of the inactive conformations in available X-ray structures, implying that one of the conformations is preferred for many kinases. The classification is consistent with selectivity profiles of several well-characterized kinase inhibitors. We show further that inhibitor selectivity profiles guide kinase classification. For example, selective inhibition of lck among src-family kinases by imatinib (Gleevec) suggests that the relative stabilities of inactive conformations of lck are different from other src-family kinases. We report the X-ray structure of the lck/imatinib complex, confirming that the conformation adopted by lck is distinct from other structurally-characterized src-family kinases and instead resembles kinases abl1 and kit in complex with imatinib. Our classification creates new paths for designing small-molecule inhibitors.     

Molecular orbital studies on the conformation of phospholipids. II. Preferred conformations of hydrocarbon chains and molecular organization in biomembranes.

The preferred conformations of the nonpolar β and γ (hydrocarbon) chains in phospholipids have been derived using EHT and CNDO calculations. These calculations indicate that the possible conformations of phospholipids are highly restricted. The calculations find support from X-ray diffraction studies and NMR measurements on model compounds. When considering conformations relevant to structures in cell membranes, a further selection is possible because of the fact that in aqueous solutions hydrophobic interactions stabilize an arrangement where the hydrocarbon chains (β and γ) are stacked almost parallel to one another, leading to a bilayer structure. The various models for β and γ-chains which satisfy this condition have been compared and it has been shown that of these only four are favoured by energy considerations. These arrangements differ from one another in the orientation of the β-chain and γ-chains in the interior of the bilayer structure. A low energy pathway connects these conformations and thus the molecule can easily flip from one stable bilayer arrangement to another. The possible conformations of the polar group (α) are likewise restricted. The proposed model provides explanations to a number of dynamic and static properties of phospholipids, in particular to the observed NMR coupling constants, ¹H and ¹³C relaxation times, studies based on ESR spin labels and the observed X-ray diffraction results on model compounds.     

Temperature and solvent dependence of the dynamical landscape of tau protein conformations

We report the variation with temperature of the ensemble distribution of conformations spanned by the tau protein in its dynamical states measured by small-angle X-ray scattering (SAXS) using synchrotron radiation. The SAXS data show a clear temperature variation of the distribution of occupied protein conformations from 293 to 318 K. More conformations with a smaller radius of gyration are occupied at higher temperature. The protein–solvent interactions are shown by computer simulation to be essential for controlling the dynamics of protein conformations, providing evidence for the key role of water solvent in the protein dynamics, as proposed by Giorgio Careri.     

Calculation of protein conformation by the build-up procedure. Application to bovine pancreatic trypsin inhibitor using limited simulated nuclear magnetic resonance data

Low-energy conformations of a set of tetrapeptides derived from the small protein bovine pancreatic trypsin inhibitor (BPTI) were generated by a build-up procedure from the low-energy conformations of single amino acid residues. At each stage, various-size fragments were built up from all combinations of smaller ones, the total energies were then minimized, and the low-energy conformations were retained for the next stage. The energies of the tetrapeptides were re-ordered by including the effects of hydration. No information other than the amino acid sequence was used to obtain the low-energy conformations of the hydrated tetrapeptides. The latter were then supplemented with a limited set of simulated NMR distance information, derived from the X-ray structure of BPTI, to provide a basis for building the rest of the whole protein molecule by the same procedure. A total of 189 upper bounds, plus 12 pairs of upper and lower bounds pertaining to the location of the three disulfide bonds in this molecule, were used. Four sets of conformations of the entire molecule were generated by utilizing different combinations of smaller fragments. It was possible to obtain low-energy conformations with small rms deviations, 1.1 to 1.4 A for the alpha-carbons, from the structure derived by X-ray diffraction. The average deviations of the backbone dihedral angles were also low, viz. 23 degrees to 26 degrees.     

Calculation of conformational ensembles from potentials of mean force. An approach to the knowledge-based prediction of local structures in globular proteins

We present a prototype of a new approach to the folding problem of polypeptide chains. This approach is based on the analysis of known protein structures. It derives the energy potentials for the atomic interactions of all amino acid residue pairs as a function of the distance between the involved atoms. These potentials are then used to calculate the energies of all conformations that exist in the data base with respect to a given sequence. Then, by using only the most stable conformations, clusters of the most probable conformations for the given sequence are obtained. To discuss the results properly we introduce a new classification of segments based on their conformational stability. Special care is taken to allow for sparse data sets. The use of the method is demonstrated in the discussion of the identical oligopeptide sequences found in different conformations in unrelated proteins. VNTFV, for example, adopts a beta-strand in ribonuclease but it is found in an alpha-helical conformation in erythrocruorin. In the case of VNTFV the ensemble obtained consists of a single cluster of beta-strand conformations, indicating that this may be the preferred conformation for the pentapeptide. When the flanking residues are included in the calculation the hepapeptide P-VNTFV-H (ribonuclease) again yields an ensemble of beta-strands. However, in the ensemble of D-VNTFV-A (erythrocruorin) the major cluster is of alpha-helical type. In the present study we concentrate on the local aspects of protein conformations. However, the theory presented is quite general and not restricted to oligopeptides. We indicate extensions of the approach to the calculation of global conformations of proteins as well as conceivable applications to a number of molecular systems.     

Homology modeling of histidine-containing phosphocarrier protein and eosinophil-derived neurotoxin: Construction of models and comparison with experiment

Homology modeling methods have been used to construct models of two proteins—the histidine-containing phosphocarrier protein (HPr) from Mycoplasma capricolum and human eosinophil-derived neurotoxin (EDN). Comparison of the models with the subsequently determined X-ray crystal structures indicates that the core regions of both proteins are reasonably well reproduced, although the template structures are closer to the X-ray structures in these regions—possible enhancements are discussed. The conformations of most of the side chains in the core of HPr are well reproduced in the modeled structure. As expected, the conformations of surface side chains in this protein differ significantly from the X-ray structure. The loop regions of EDN were incorrectly modeled—reasons for this and possible enhancements are discussed. © 1995 Wiley-Liss, Inc.     

Variable-Target-Function and Build-up Procedures for the Calculation of Protein Conformation. Application to Bovine Pancreatic Trypsin Inhibitor Using Limited Simulated Nuclear Magnetic Resonance Data

Abstract

An implementation of the variable-target-function procedure, first introduced by Braun and Gō [W. Braun and N. Gō, J. Mol. Biol. 186, 611–626 (1985)], has been used to generate conformations of the small protein bovine pancreatic trypsin inhibitor (BPTI), given a limited set of simulated data that could be obtained by nuclear magnetic resonance (NMR) techniques. A hybrid strategy was also used to calculate conformations of BPTI, given the same information. In the hybrid strategy, low-energy structures of medium-size fragments (decapeptides) of BPTI were generated using the variable-target-function method, followed by restrained energy optimization. The low-energy conformations were used as a basis to build up the complete fifty-eight-residue BPTI molecule. By using the variable-target-function approach, in which energy considerations were not introduced until full conformations of the entire BPTI molecule had been generated, it was not possible to obtain calculated structures with rms deviations from the X-ray conformation of less than 1.6 Å for the α-carbons. On the other hand, with the hybrid strategy, which involved the consideration of realistic energy terms in the early stages of the calculations, it was possible to calculate low-energy conformations of BPTI with rms deviations from the X-ray structure of 1.06 to 1.50 Å for the α-carbons. When the rms deviations were computed along the amino acid sequence, it was found that there was a good correlation between deviations among the calculated structures and deviations from the X-ray structure.     

Conformational analysis of protein side chains using two-dimensional Overhauser nuclear effect spectroscopy]

The method has been proposed to determine the conformations of protein side-chains (dihedral angles chi 1) using two-dimensional nuclear Overhauser effect spectroscopy data. This method is grounded of the algorithm prepared on the basis of joint consideration of proton-proton distance dependences in dipeptide units of L-amino acid residues on the dihedral angles phi, psi and chi 1 with the accounting of the local sterical conditions of the polypeptide chain. The obtained results gave the possibility to bring the different regions of space (phi, psi) of amino acid residues into the line with the specific sets of nuclear Overhauser effect spectral parameters which unambiguously characterize in most cases the conformational states of their side-chains. The method efficiency was displayed on the test calculation with the utilization as the experimental data of the "model" nuclear Overhauser effect contacts derived from the X-ray atomic coordinates of the bovine pancreatic trypsin inhibitor molecule.     

The energetics of off-rotamer protein side-chain conformations

Non-rotameric ("off-rotamer") conformations are commonly observed for the side-chains of protein crystal structures. This study examines whether such conformations are real or artifactual by comparing the energetics of on and off-rotamer side-chain conformations calculated with the CHARMM energy function. Energy-based predictions of side-chain orientation are carried out by rigid-geometry mapping in the presence of the fixed protein environment for 1709 non-polar side-chains in 24 proteins for which high-resolution (2.0 A or better) structures are available. For on-rotamer conformations, 97.6 % are correctly predicted; i.e. they correspond to the absolute minima of their local side-chain energy maps (generally to within 10 degrees or less). By contrast, for the observed off-rotamer side-chain conformations, 63.8 % are predicted correctly. This difference is statistically significant (P<0.001) and suggests that while most of the observed off-rotamer conformations are real, many of the erroneously predicted ones are likely to be artifacts of the X-ray refinements. Probabilities for off-rotamer conformations of the non-polar side-chains are calculated to be 5.0-6.1 % by adaptive umbrella-sampled molecular dynamics trajectories of individual amino acid residues in vacuum and in the presence of an average protein or aqueous dielectric environment. These results correspond closely to the 5.7 % off-rotamer fraction predicted by the rigid-geometry mapping studies. Since these values are about one-half of the 10.2 % off-rotamer fraction observed in the X-ray structures, they support the conclusion that many of the latter are artifacts. In both the rigid-geometry mapping and the molecular dynamics studies, the discrepancies between the predicted and observed fractions of off-rotamer conformations are largest for leucine residues (approximately 6 % versus 16.6 %). The simulations for the isolated amino acid residues indicate that the real off-rotamer frequency of 5-6 % is consistent with the internal side-chain and local side-chain-backbone energetics and does not originate from shifts due to the protein. The present results suggest that energy-based rotation maps can be used to find side-chain positional artifacts that appear in crystal structures based on refinements in the 2 A resolution range.     

Variable-target-function and build-up procedures for the calculation of protein conformation. Application to bovine pancreatic trypsin inhibitor using limited simulated nuclear magnetic resonance data

An implementation of the variable-target-function procedure, first introduced by Braun and Go [W. Braun and N. Go, J. Mol. Biol. 186, 611-626 (1985)], has been used to generate conformations of the small protein bovine pancreatic trypsin inhibitor (BPTI), given a limited set of simulated data that could be obtained by nuclear magnetic resonance (NMR) techniques. A hybrid strategy was also used to calculate conformations of BPTI, given the same information. In the hybrid strategy, low-energy structures of medium-size fragments (decapeptides) of BPTI were generated using the variable-target-function method, followed by restrained energy optimization. The low-energy conformations were used as a basis to build up the complete fifty-eight-residue BPTI molecule. By using the variable-target-function approach, in which energy considerations were not introduced until full conformations of the entire BPTI molecule had been generated, it was not possible to obtain calculated structures with rms deviations from the X-ray conformation of less than 1.6 A for the alpha-carbons. On the other hand, with the hybrid strategy, which involved the consideration of realistic energy terms in the early stages of the calculations, it was possible to calculate low-energy conformations of BPTI with rms deviations from the X-ray structure of 1.06 to 1.50 A for the alpha-carbons. When the rms deviations were computed along the amino acid sequence, it was found that there was a good correlation between deviations among the calculated structures and deviations from the X-ray structure.