Interpretation of some Conformational Data in Myoglobin and Lysozyme

THE conformations of the polypeptide chains of myoglobin¹ and lysozyme² have been successfully simulated with the aid of computed Van der Waals contact and energy maps of the theoretical independent peptide unit (IPU)^3–5. The non-glycyl experimental points plotted on an alanyl IPU are rather scattered on the allowed conformational regions of the map⁶, especially in the case of lysozyme. By contrast, well defined clusters of points can be observed when only the amino-acid residues in segments of the helical secondary structure (mainly α and β chains) are plotted. In addition, clusters of points, albeit less well defined, can be observed by plotting the points relative to the experimental conformations of the first non-helical amino-acid residue next to a more or less folded segment of that α-helical type so frequently present in globular proteins (Fig. 1).     

Alteration of the DNA double helix conformation upon incorporation of mispairs as revealed by energy computations and pathways of point mutations.

To explain biochemical and genetic data on spontaneous nucleotide replacements in nucleic acid biosynthesis all the 8 mispairs in normal tautomeric forms have been considered. Possible B-conformations of DNA fragments containing each of such mispairs incorporated between Watson-Crick pairs have been found using computations of the energy of non-bonded interactions via classical potential functions. These conformations have no reduced interatomic contacts. The values of each dihedral angle of the sugar-phosphate backbone fall within the limits of those of double-helical fragments of B-DNA in crystals. These values differ from those of the corresponding angles for the low-energy polynucleotide conformations consisting of canonical pairs by no more than 30 degrees (except for the fragment with the U:U pair for which the C4'-C3'-O-P angle differs by about 50 degrees). The difference in experimentally observed frequencies of various nucleotide replacements in DNA biosynthesis correlates with the difference in the energy of non-bonded interactions and with the extent of the sugar-phosphate backbone distortion for the fragments containing the mispairs which serve as intermediates for the replacements.     

Conformational behavior of chondroitin and chondroitin sulfate in relation to their physical properties as inferred by molecular modeling

Chondroitin and chondroitin sulfates belong to the family of glycosaminoglycans. They are most widely distributed in animal tissues, where they are involved in structural functions and in cell-cell communication. Their basic structures consist of a disaccharidic repeating unit of beta-D-glucuronic acid (GlcA) and 2-acetamido-2-deoxy-beta-D-galactose (GalNAc), this latter being sulfated at different positions. Molecular mechanics has been applied to calculate the adiabatic energy maps for each of the constituting disaccharides of chondroitin, chondroitin 4-sulfate, and chondroitin 6-sulfate using the MM3 force field. Based on these maps, higher levels of structural organization have been simulated. On one hand, the disordered state is studied through a Metropolis-based algorithm; the resulting chains present a behavior of semirigid polymers, with an order of stiffness: chondroitin 4-sulfate > chondroitin > chondroitin 6-sulfate. On the other hand, the exploration of the stable ordered forms leads to numerous helical conformations of comparable energies. Several of these conformations correspond to the experimentally observed ones. The ability of coordination with cations has also been explored, resulting in a preferential stereospecificity for calcium ions when compared to sodium ions.     

Native and nonnative conformational preferences in the urea-unfolded state of barstar

The refolding of barstar from its urea-unfolded state has been studied extensively using various spectroscopic probes and real-time NMR, which provide global and residue-specific information, respectively, about the folding process. Here, a preliminary structural characterization by NMR of barstar in 8 M urea has been carried out at pH 6.5 and 25 degrees C. Complete backbone resonance assignments of the urea-unfolded protein were obtained using the recently developed three-dimensional NMR techniques of HNN and HN(C)N. The conformational propensities of the polypeptide backbone in the presence of 8 M urea have been estimated by examining deviations of secondary chemical shifts from random coil values. For some residues that belong to helices in native barstar, 13C(alpha) and 13CO secondary shifts show positive deviations in the urea-unfolded state, indicating that these residues have propensities toward helical conformations. These residues are, however, juxtaposed by residues that display negative deviations indicative of propensities toward extended conformations. Thus, segments that are helical in native barstar are unlikely to preferentially populate the helical conformation in the unfolded state. Similarly, residues belonging to beta-strands 1 and 2 of native barstar do not appear to show any conformational preferences in the unfolded state. On the other hand, residues belonging to the beta-strand 3 segment show weak nonnative helical conformational preferences in the unfolded state, indicating that this segment may possess a weak preference for populating a helical conformation in the unfolded state.     

Role of hydrophobicity and solvent-mediated charge-charge interactions in stabilizing alpha-helices.

A theoretical study to identify the conformational preferences of lysine-based oligopeptides has been carried out. The solvation free energy and free energy of ionization of the oligopeptides have been calculated by using a fast multigrid boundary element method that considers the coupling between the conformation of the molecule and the ionization equilibria explicitly, at a given pH value. It has been found experimentally that isolated alanine and lysine residues have somewhat small intrinsic helix-forming tendencies; however, results from these simulations indicate that conformations containing right-handed alpha-helical turns are energetically favorable at low values of pH for lysine-based oligopeptides. Also, unusual patterns of interactions among lysine side chains with large hydrophobic contacts and close proximity (5-6 A) between charged NH3+ groups are observed. Similar arrangements of charged groups have been seen for lysine and arginine residues in experimentally determined structures of proteins available from the Protein Data Bank. The lowest-free-energy conformation of the sequence Ac-(LYS)6-NMe from these simulations showed large pKalpha shifts for some of the NH3+ groups of the lysine residues. Such large effects are not observed in the lowest-energy conformations of oligopeptide sequences with two, three, or four lysine residues. Calculations on the sequence Ac-LYS-(ALA)4-LYS-NMe also reveal low-energy alpha-helical conformations with interactions of one of the LYS side chains with the helix backbone in an arrangement quite similar to the one described recently by (Proc. Natl. Acad. Sci. U.S.A. 93:4025-4029). The results of this study provide a sound basis with which to discuss the nature of the interactions, such as hydrophobicity, charge-charge interaction, and solvent polarization effects, that stabilize right-handed alpha-helical conformations.     

Monte carlo simulations of protein-like heteropolymers.

Properties of a simple model of polypeptide chains were studied by the means of the Monte Carlo method. The chains were built on the (310) hybrid lattice. The residues interacted with long-range potential. There were two kinds of residues: hydrophobic and hydrophilic forming a typical helical pattern -HHPPHPP-. Short range potential was used to prefer helical conformations of the chain. It was found that at low temperatures the model chain formes dense and partially ordered structures (non-unique). The presence of the local potential led to an increase of helicity. The effect of the interplay between the two potentials was studied. After the collapse of the chain further annealing caused rearrangement of helical structures. Dynamic properties of the chain at low temperature depended strongly on the local chain ordering.     

Evaluation of the conformational free energies of loops in proteins

In this paper we discuss the problem of including solvation free energies in evaluating the relative stabilities of loops in proteins. A conformational search based on a gas-phase potential function is used to generate a large number of trial conformations. As has been found previously, the energy minimization step in this process tends to pack charged and polar side chains against the protein surface, resulting in conformations which are unstable in the aqueous phase. Various solvation models can easily identify such structures. In order to provide a more severe test of solvation models, gas phase conformations were generated in which side chains were kept extended so as to maximize their interaction with the solvent. The free energies of these conformations were compared to that calculated for the crystal structure in three loops of the protein E. coli RNase H, with lengths of 7, 8, and 9 residues. Free energies were evaluated with a finite difference Poisson-Boltzmann (FDPB) calculation for electrostatics and a surface area-based term for nonpolar contributions. These were added to a gas-phase potential function. A free energy function based on atomic solvation parameters was also tested. Both functions were quite successful in selecting, based on a free energy criterion, conformations quite close to the crystal structure for two of the three loops. For one loop, which is involved in crystal contacts, conformations that are quite different from the crystal structure were also selected. A method to avoid precision problems associated with using the FDPB method to evaluate conformational free energies in proteins is described. © 1994 John Wiley & Sons, Inc.     

A stereochemical model for phospholipids: I. Conformational energy refinement and molecular packing of l-α-Dipalmitoyl lecithin (DPL)

The preferred conformations of L-α-dipalmitoyl-lecithin (DPL) have been refined using a steepest descent procedure throughout non-bonded potential energy calculations. The results indicate that energy differences between the conformers is very low and the energy parameters are quite constant around the minimum, suggesting a large degree of flexibility of this molecule.The molecular packing energy calculations have been performed by separating the two hydrocarbon chains from the polar head groups of the molecule. From the energy parameters for the packing of the aliphatic chains it results that for distances between adjacent chains up to 4·4 Å the intermolecular forces allow the maximum degree of freedom. This suggests that hydrocarbon chains do not play the main role in the packing process of DPL molecule. Therefore the energy parameters of the polar segment have been calculated, assuming that the C₂ asymmetric carbon atom represents the points of the hexagonal lattice and the rotation centre for each molecule. For the internal symmetry of this segment of the molecule two non-equivalent conformers have been selected over all sets of allowed conformations (the GGG and GGG¹ for the α₂, α₃ and α₅ torsion angles). The energy packing calculation has been carried out for two independent sets of data, with and without the electrostatic contributions. In the first case a unique topological situation is allowed with the P-N dipole lying parallel to the lattice plane. In the second case different situations including that with the P-N dipole lying orthogonal to the plane are allowed. These data are discussed in relation to different physical conditions.     

Conformational Studies of β-glucans

Steric and energy contour diagrams have been plotted for disaccharide-like and for helical structures of linear β-D -glucans having (1 → 2), (1 → 3) and (1 → 4) linkages. The allowed conformations constitute only about. 4% of the total conformations, indicating that the freedom of rotation of glucose residues is highly restricted in all the three polysaccharides. The additional restrictions of the monomer unit, as one passes from disaccharide to polysaccaride structures, are severe in the case of (1 → 2) and (1 → 3) linked polysaccharides but not in (1 → 4) linked polysaccharide. The difference in the nature of linkages also has shown to affect the energetically preferred conformations: (1 → 2) linkages lead only to left handed helical conformations; (1 → 3) linkages lead to both right and left handed wide and extended helical conformations, (1 → 4) linkages lead to both right and left handed extended helical conformations. The possible hydrogen bonds between adjacent residues are also dependent on the nature of linkage.     

Model building of DNA/RNA triple helix containing L-deoxyadenosine.

A three-dimensional model of DNA/RNA triple helix that contains a poly(L-deoxyadenosine) (L-dA) chain is proposed based on computer-assisted model building and energy calculations. The model building was performed by a new method that systematically searches possible conformations of nucleotide units in the helical chains. Two possible orientations of sugar-phosphate chains, in which two homopyrimidine strands are parallel or antiparallel with each other, were considered in the systematic search. Several possible base-pairing models, in which there are one Watson-Crick base pair and one other base pair, were also considered. Many possible models selected by the systematic search were further refined through molecular mechanics calculation incorporating a helical boundary condition. The preferred model, which was selected on the basis of potential energy, was the one with Watson-Crick and Hoogsteen base pairs and with its two polypyrimidine chains in the antiparallel orientation. The model can explain the experimental observation that poly(L-dA) forms a stable triple helix with poly(uridylic acid) (U) but not with poly(deoxythymidylic acid) (dT).     

Conformational analysis of a segment in bacterioopsin by two-dimensional (1)H-NMR spectroscopy]

The spatial structure of a synthetic peptide, an analogue of the membrane spanning segment B (residues 34-65) of bacterioopsin from Halobacterium halobium, has been refined. Backbone torsion angles were derived from intensities of short-range interproton NOEs. These, together with a complete set of the NOEs integral intensities formed the basis for the three-dimensional structure refinement by the energy minimization with consideration of NOE penalty functions. Analysis indicates the right-handed alpha-helical conformation of segment B extending from Asp-38 to Tyr-64 with a kink of the helical axis (27 degrees) at Pro-50. The most stable region with an average root-mean-square deviation of 0.43 A between the backbone atoms includes residues 42-60 in six energy refined structures. The N-terminal part of segment B (residues 34-37) has no ordered conformation. The inferred structure is in close agreement with the electron cryomicroscopy structure of bacteriorhodopsin, differing from it in conformations of most of the side chains.     

Nonrandom structure in the urea-unfolded Escherichia coli outer membrane protein X (OmpX)

On the basis of sequence-specific resonance assignments for the complete polypeptide backbone and most of the amino acid side chains by heteronuclear nuclear magnetic resonance (NMR) spectroscopy, the urea-unfolded form of the outer membrane protein X (OmpX) from Escherichia coli has been structurally characterized. (1)H-(1)H nuclear Overhauser effects (NOEs), dispersion of the chemical shifts, amide proton chemical shift temperature coefficients, amide proton exchange rates, and (15)N[(1)H]-NOEs show that OmpX in 8 M urea at pH 6.5 is globally unfolded, but adopts local nonrandom conformations in the polypeptide segments of residues 73-82 and 137-145. For these two regions, numerous medium-range and longer-range NOEs were observed, which were used as the input for structure calculations of these polypeptide segments with the program DYANA. The segment 73-82 forms a quite regular helical structure, with only loosely constrained amino acid side chains. In the segment 137-145, the tryptophan residue 140 forms the core of a small hydrophobic cluster. Both nonrandom structures are present with an abundance of about 25% of the protein molecules. The sequence-specific NMR assignment and the physicochemical characterization of urea-denatured OmpX presented in this paper are currently used as a platform for investigations of the folding mechanism of this integral membrane protein.     

Conformational analysis of xyloglucans

Xyloglucan isolated from the elongating regions of pea stems was examined using X-ray diffraction and energy calculations. The X-ray fibre pattern suggested that the backbone (1----4)-beta-D-glucan takes an extended two-fold helix similar to common cellulose. In order to study side chains (xylosyl or fucosyl-galactosyl-xylosyl residues) of the polysaccharide, energetically preferable conformations were searched by calculation of interactions between non-bonded atom pairs. A stepwise calculation for the conformation of fucosyl-galactosyl-xylosyl residue gave 10 allowed area (phi-psi) maps which are useful to deduce xyloglucan conformations of both monocotyledons and dicotyledons in the walls of growing plant cells.     

11-cis-retinal protonated Schiff base: influence of the protein environment on the geometry of the rhodopsin chromophore

Density functional theory (DFT) calculations based on the self-consistent-charge tight-binding approximation have been performed to study the influence of the protein pocket on the 3-dimensional structure of the 11-cis-retinal Schiff base (SB) chromophore. Starting with an effectively planar chromophore embedded in a protein pocket consisting of the 27 next-nearest amino acids, the relaxed chromophore geometry resulting from energy optimization and molecular dynamics (MD) simulations has yielded novel insights with respect to the following questions: (i) The conformation of the beta-ionone ring. The protein pocket tolerates both conformations, 6-s-cis and 6-s-trans, with a total energy difference of 0.7 kcal/mol in favor of the former. Of the two possible 6-s-cis conformations, the one with a negative twist angle (optimized value: -35 degrees ) is strongly favored, by 3.6 kcal/mol, relative to the one in which the dihedral is positive. (ii) Out-of-plane twist of the chromophore. The environment induces a nonplanar helical deformation of the chromophore, with the distortions concentrated in the central region of the chromophore, from C10 to C13. The dihedral angle between the planes formed by the bonds from C7 to C10 and from C13 to C15 is 42 degrees. (iii) The absolute configuration of the chromophore. The dihedral angle about the C12-C13 bond is +170 degrees from planar s-cis, which imparts a positive helicity on the chromophore, in agreement with earlier considerations based on theoretical and spectroscopic evidence.     

The effect of sequence patterns and local conformational preferences on the structure of collapsed polypeptide

We studied the structure of a polypeptide model by means of the Monte Carlo method. The model chain consisting of two kinds of residues (hydrophobic and hydrophilic) was confined on the (310) hybrid lattice. The residues interacted with the long-range contact potential. The short-range potential was also used by introducing the preferences of conformations corresponding to the helical structure. Simulations of the coil-to-globule collapse were done by an annealing process starting from high-temperature structures and then gradual cooling of the system. The parameters describing the behavior of the system were monitored. It has been found that in a case of a helical pattern -HHPPHPP- the collapsed chains consisted of helical fragments. The resulting structures were formed in such way that the polar residues were located in the outer shell of the molecule since the hydrophobic residues filled the inner part of molecule. The results showed that the proper balance between the magnitude of the potentials used in a model is important and its influence on final structures of molecules was discussed.     

An algorithm for determining the conformation of polypeptide segments in proteins by systematic search

The feasibility of determining the conformation of segments of a polypeptide chain up to six residues in length in globular proteins by means of a systematic search through the possible conformations has been investigated. Trial conformations are generated by using representative sets of phi, psi, and chi angles that have been derived from an examination of the distributions of these angles in refined protein structures. A set of filters based on simple rules that protein structures obey is used to reduce the number of conformations to a manageable total. The most important filters are the maintenance of chain integrity and the avoidance of too-short van der Waals contacts with the rest of the protein and with other portions of the segment under construction. The procedure is intended to be used with approximate models so that allowance is made throughout for errors in the rest of the structure. All possible main chains are first constructed and then all possible side-chain conformations are built onto each of these. The electrostatic energy, including a solvent screening term, and the exposed hydrophobic area are evaluated for each accepted conformation. The method has been tested on two segments of chain in the trypsin like enzyme from Streptomyces griseus. It is found that there is a wide spread of energies among the accepted conformations, and the lowest energy ones have satisfactorily small root mean square deviations from the X-ray structure.     

Conformational behavior of hyaluronan in relation to its physical properties as probed by molecular modeling

Hyaluronan (HA) is a linear charged polysaccharide whose structure is made up of repeating disaccharide units. Apparently conflicting reports have been published about the nature of the helical structure of HA in the solid state. Recent developments in the field of molecular modeling of polysaccharides offer new opportunities to reexamine the structural basis underlying the formation and stabilization of ordered structures and their interactions with counterions. The conformational spaces available and the low energy conformations for the disaccharide, trisaccharide, and tetrasaccharide segments of HA were investigated via molecular mechanics calculations using the MM3 force field. First, the results were used to access the configurational statistics of the corresponding polysaccharide. A disordered chain having a persistence length of 75 A at 25 degrees C is predicted. Then, the exploration of the stable ordered forms of HA led to numerous helical conformations, both left- and right-handed, having comparable energies. Several of these conformations correspond to the experimentally observed ones and illustrate the versatility of the polysaccharide. The double stranded helical forms have also been explored and theoretical structures have been compared to experimentally derived ones.     

Contributions of left-handed helical residues to the structure and stability of bacteriophage T4 lysozyme

Non-glycine residues in proteins are rarely observed to have "left-handed helical" conformations. For glycine, however, this conformation is common. To determine the contributions of left-handed helical residues to the stability of a protein, two such residues in phage T4 lysozyme, Asn55 and Lys124, were replaced with glycine. The mutant proteins fold normally and are fully active, showing that left-handed non-glycine residues, although rare, do not have an indispensable role in the folding of the protein or in its activity. The thermodynamic stability of the Lys124 to Gly variant is essentially identical with that of wild-type lysozyme. The Asn55 to Gly mutant protein is marginally less stable (0.5 kcal/mol). These results indicate that the conformational energy of a glycine and a non-glycine residue in the left-handed helical conformation are very similar. This is consistent with some theoretical energy distributions, but is inconsistent with others, which suggest that replacements of the sort described here might increase the stability of the protein by up to 5 kcal/mol. Crystallographic analysis of the mutant proteins shows that the backbone conformation of the Lys124 to Gly variant is essentially identical with that of the wild-type structure. In the case of the Asn55 to Gly replacement, however, the (phi, psi) values of residue 55 change by about 20 degrees. This suggests that the energy minimum for left-handed glycine residues is not the same as that for non-glycine residues. This is strongly indicated also by a survey of accurately determined protein crystal structures, which suggests that the energy minimum for left-handed glycine residues is near (phi = 90 degrees, psi = 0 degrees), whereas that for non-glycine residues is close to (phi = 60 degrees, psi = 30 degrees). This apparent energy minimum for glycine is not clearly predicted by any of the theoretical (phi, psi) energy contour maps.     

Biophysical studies on molecular mechanism of abortifacient action of prostaglandins. VI. Conformation energy calculation on PGE2, PGF2 alpha and 15-(s)-methyl PGF2 alpha

This paper reports the conformation energy (CE) calculations on PGE2, PGE2 alpha and 15-(s)-methyl PGE2 alpha on the basis of empirical potential energy functions for the simultaneous rotations around C7-C8 (psi), C12-C13 (theta) and C14-C15 (phi) bonds. The variation of the minimum conformation energy E for each isoenergy map in the psi theta plane with respect to phi gives two minima around 90 degrees and 240 degrees in PGE2, 60 degrees and 245 degrees in PGF2 alpha, and 60 degrees and 150 degrees in 15-(s)-methyl PGF2 alpha. The latter two forms also have a small dip at 270 degrees. The pattern of allowed low energy conformations of PGF2 alpha and 15-(s)-methyl PGF2 alpha is quite similar and is characterized by the existence of six low energy regions.