Conformation and hydration of aspartame

Conformational free energy calculations using an empirical potential (ECEPP/2) and the hydration shell model were carried out on the neutral, acidic, zwitterionic, and basic forms of aspartame in the hydrated state. The results indicate that as the molecule becomes more charged, the number of low energy conformations becomes smaller and the molecule becomes less flexible. The calculated free energies of hydration of charged aspartames show that hydration has a significant effect on the conformation in solution. Only two feasible conformations were found for the zwitterionic form, and these are consistent with the conformations deduced from NMR and X-ray diffraction experiments. The calculated free energy difference between these two conformations was 1.25 kcal/mol. The less favored of the two solvated conformations can be expected to be stabilized by hydrophobic interaction of the phenyl groups in the crystal.     

Preferred conformations of cyclic Ac-Cys-Pro-Xaa-Cys-NHMe peptides: a model for chain reversal and active site of disulfide oxidoreductase

The conformational study on cyclic Ac-Cys-Pro-Xaa-Cys-NHMe (Ac-CPXC-NHMe; X=Ala, Val, Leu, Aib, Gly, His, Phe, Tyr, Asn and Ser) peptides has been carried out using the Empirical Conformational Energy Program for Peptides, version 3 (ECEPP/3) force field and the hydration shell model in the unhydrated and hydrated states. This work has been undertaken to investigate structural implications of the CPXC sequence as the chain reversal for the initiation of protein folding and as the motif for active site of disulfide oxidoreductases. The backbone conformation DAAA is commonly the most feasible for cyclic CPXC peptides in the hydrated state, which has a type I beta-turn at the Pro-Xaa sequence. The proline residue and the hydrogen bond between backbones of two cystines as well as the formation of disulfide bond appear to play a role in stabilizing this preferred conformation of cyclic CPXC peptides. However, the distributions of backbone conformations and beta-turns may indicate that the cyclic CPXC peptide seems to exist as an ensemble of beta-turns and coiled conformations in aqueous solution. The intrinsic stability of the cyclic CPXC motif itself for the active conformation seems to play a role in determining electrochemical properties of disulfide oxidoreductases.     

Preferred conformations of a linear RGD tripeptide.

Conformational study of RGD tripeptides in the nonhydrated and hydrated states was carried out using an empirical potential function ECEPP/3 and the hydration shell model in order to investigate preferred conformations and factors responsible for their stability. RGD tripeptides in the nonhydrated and hydrated states can be interpreted as existing as an ensemble of feasible conformations rather than as a single dominant conformation from the analysis of distributions of backbone conformations, hydrogen bonds and beta-turns. The different distributions of conformations for the neutral and zwitterionic RGD tripeptides in both states may indicate that the conformation of the RGD tripeptide is liable to depend on solvent polarity and pH values. beta-Turn populations for the neutral tripeptide in both states are reasonably consistent with NMR measurements on linear RGD-containing peptides. The degradation of RGD tripeptide seems to be attributed mainly to the hydrogen bonds between the Asp side-chain and the backbone of Asp residue or C-terminal NHMe group, rather than to the flexible backbones of Gly and Asp residues.     

Conformational study of asialo-GM1 (GA1) ganglioside

In order to investigate the significance of preferred conformations of the saccharide for the steric orientation and recognition of glycosphingolipids at the membrane surface, the conformational free energy calculations were carried out on the asialo-GM1 [GA1; β-D -Gal(1 → 3) β-D -GalNAc(1 → 4) β-D -Gal(1 → 4) β-D -Glc-O-ceramide] using a new program CONCARB (CONformational study program for CARBohydrate) in the unhydrated and hydrated states. The overall backbone conformation of GA1 appears nearly to be extended with a little bent at the glycosidic II–III linkage, in which two pyranose rings of Gal(IV)-GalNAc-(III) moiety orient approximately perpendicular to those of Gal(II)-Glc(I) moiety. This is consistent with the structures deduced from high-sensitivity differential scanning calorimetry experiments and the nmr study on GA1. The calculated glycosidic torsion angles of the lowest free energy conformation of GA1 in the hydrated state are in accord with the structures of relevant oligosaccharides deduced from nmr experiments and hard sphere exoanomeric calculations. A comparison of the values of glycosidic torsion angles ϕ and π of GA1 and its constituent oligosaccharides indicates that the overall backbone conformation of each oligosaccharide is retained when the oligosaccharide chain becomes longer. This implies that the short-range interactions between the nearest-neighbored saccharides are of significant importance in stabilizing the overall backbone conformation of GA1 in both the unhydrated and hydrated states. The different orientation and hydrogen bonds of hydroxymethyl and hydroxyl groups from one oligosaccharide to another suggest that the medium- and long-range interactions are also of consequence. Hydration seems to affect significantly the conformation of these groups, but not to perturb remarkably the overall backbone conformation of GA1. © 1997 John Wiley & Sons, Inc. Biopoly 42: 19–35, 1997     

Influence of hydration on the conformational stability and formation of bends in terminally blocked dipeptides

The conformations of 23 terminally blocked dipeptide sequences were examined by conformational energy calculations that included the effects of the aqueous solvent. Starting structures were derived from combinations of minimum-energy conformations of hydrated single residues. Their conformational energies were then minimized using the ECEPP potential (Empirical Conformational Energy Program for Peptides) with hydration included. Short-range interactions dominate in stabilizing the conformations of the hydrated dipeptides. Differences between conformational stabilities of hydrated and unhydrated dipeptides in many cases are due to the competition of solute–water and intramolecular hydrogen bonds. In other cases, perturbation of the hydration shell of the solute by close approach of solute atoms alters conformational preferences. Probabilities of formation of bends were calculated and compared to the corresponding quantities for unhydrated dipeptides and to those calculated from x-ray structures. For bends in dipeptides containing two nonpolar amino acids, computations omitting hydration yield better results. However, better agreement with experimental (x-ray) bend probabilities for dipeptides containing glycine or polar amino acids is obtained only in some sequences when hydration is included. The results are rationalized by the observation that, in proteins, bends containing nonpolar sequences occur on the inside, shielded from the solvent. Bends containing glycine or polar amino acids occur frequently on the surface of the protein, but they are not completely hydrated.     

Model for the conformational analysis of hydrated peptides. Effect of hydration on the conformational stability of the terminally blocked residues of the 20 naturally occurring amino acids

The effects of hydration are included in empirical conformational energy computations on oligopeptides by means of a modified hydration-shell model. Free energy terms are introduced to account for “specific hydration” due to water–solute hydrogen bonding and for “nonspecific hydration” describing the interaction of the solute with water molecules in a first-neighbor shell. The dielectric constant has been doubled (over the value used for calculations in the absence of water) to take into account the presence of solvent. Computations were carried out for the N-acetyl-N′-methylamides of the 20 naturally occurring amino acids. Conformational energy maps are compared with similar maps calculated in the absence of hydration. Minimum-energy conformations are located and compared with the corresponding minima for unhydrated peptides in terms of ordering with respect to potential energy, the dihedral angles at the minima, and the presence of intramolecular hydrogen bonds. The Boltzmann factors for various conformational regions are altered significantly on hydration in some cases. These changes can be explained in terms of differences in the hydration free energy terms for various conformations.     

Free energy landscapes of peptides by enhanced conformational sampling

The free energy landscapes of peptide conformations in water have been observed by the enhanced conformational sampling method, applying the selectively enhanced multicanonical molecular dynamics simulations. The conformations of the peptide dimers, -Gly-Gly-, -Gly-Ala-, -Gly-Ser-, -Ala-Gly-, -Asn-Gly-, -Pro-Gly-, -Pro-Ala-, and -Ala-Ala-, which were all blocked with N-terminal acetyl and C-terminal N-methyl groups, were individually sampled with the explicit TIP3P water molecules. From each simulation trajectory, we obtained the canonical ensemble at 300 K, from which the individual three-dimensional landscape was drawn by the potential of mean force using the three reaction coordinates: the backbone dihedral angle, psi, of the first amino acid, the backbone dihedral angle, phi, of the second amino acid, and the distance between the carbonyl oxygen of the N-terminal acetyl group and the C-terminal amide proton. The most stable state and several meta-stable states correspond to extended conformations and typical beta-turn conformations, and their free energy values were accounted for from the potentials of mean force at the states. In addition, the contributions from the intra-molecular energies of peptides and those from the hydration effects were analyzed. Consequently, the stable beta-turn conformations in the free energy landscape were consistent with the empirically preferred beta-turn types for each amino acid sequence. The thermodynamic values for the hydration effect were decomposed and they correlated well with the empirical values estimated from the solvent accessible surface area of each molecular conformation during the trajectories. The origin of the architecture of protein local fragments was analyzed from the viewpoint of the free energy and its decomposed factors.     

Preferred conformations of RGDX tetrapeptides to inhibit the binding of fibrinogen to platelets

The conformational study on Arg-Gly-Asp (RGD)-containing tetrapeptides in the unhydrated and hydrated states has been carried out using the force field ECEPP/3 and the hydration shell model. The tetrapeptides studied here are H-RGDX-OH (X = Trp, Tyr, Phe, Leu, Val, Cys, Gln, and Ser), which show the inhibitory activity for binding of fibrinogen to platelets in the order of RGDW approximately equal to RGDY approximately equal to RGDF approximately equal to RGDL > RGDV > or = RGDC > or = RGDQ > or = RGDS. The backbone conformations with two C(7) backbone-to-backbone hydrogen bonds between Asp and Arg residues and between Xaa and Gly residues are in common most probable for the RGD sequence of RGDX tetrapeptides in the hydrated state. The dominant beta-turns for RGDX are found to be the types V' and IV at Gly-Asp and Asp-Xaa sequences, respectively, which are quite similar to the types II' and I (or II), respectively. However, it cannot be ruled out that the extended conformations are also remarkably feasible for RGDX tetrapeptides in water by peering the distributions of backbone conformations. These calculated results are consistent with the experimental results on RGD-containing proteins and conformationally constrained RGD-containing peptides. The reason why the RGDX becomes more potent as the side chain of the X residue is more hydrophobic may be ascribed to that the more hydrophobic is the residue X, the more populated are beta-turn structures for the Gly-Asp sequence. The hydrophobic side chain of X residue exposed to water is likely to interact with the hydrophobic region of receptor easily.     

Empirical solvation models in the context of conformational energy searches: application to bovine pancreatic trypsin inhibitor.

Continuum solvation models that estimate free energies of solvation as a function of solvent accessible surface area are computationally simple enough to be useful for predicting protein conformation. The behavior of three such solvation models has been examined by applying them to the minimization of the conformational energy of bovine pancreatic trypsin inhibitor. The models differ only with regard to how the constants of proportionality between free energy and surface area were derived. Each model was derived by fitting to experimentally measured equilibrium solution properties. For two models, the solution property was free energy of hydration. For the third, the property was NMR coupling constants. The purpose of this study is to determine the effect of applying these solvation models to the nonequilibrium conformations of a protein arising in the course of global searches for conformational energy minima. Two approaches were used: (1) local energy minimization of an ensemble of conformations similar to the equilibrium conformation and (2) global search trajectories using Monte Carlo plus minimization starting from a single conformation similar to the equilibrium conformation. For the two models derived from free energy measurements, it was found that both the global searches and local minimizations yielded conformations more similar to the X-ray crystallographic structures than did searches or local minimizations carried out in the absence of a solvation component of the conformational energy. The model derived from NMR coupling constants behaved similarly to the other models in the context of a global search trajectory. For one of the models derived from measured free energies of hydration, it was found that minimization of an ensemble of near-equilibrium conformations yielded a new ensemble in which the conformation most similar to the X-ray determined structure PTI4 had the lowest total free energy. Despite the simplicity of the continuum solvation models, the final conformation generated in the trajectories for each of the models exhibited some of the characteristics that have been reported for conformations obtained from molecular dynamics simulations in the presence of a bath of explicit water molecules. They have smaller root mean square (rms) deviations from the experimentally determined conformation, fewer incorrect hydrogen bonds, and slightly larger radii of gyration than do conformations derived from search trajectories carried out in the absence of solvent.     

Effects of hydrated water on protein unfolding

The conformational stability of a protein in aqueous solution is described in terms of the thermodynamic properties such as unfolding Gibbs free energy, which is the difference in the free energy (Gibbs function) between the native and random conformations in solution. The properties are composed of two contributions, one from enthalpy due to intramolecular interactions among constituent atoms and chain entropy of the backbone and side chains, and the other from the hydrated water around a protein molecule. The hydration free energy and enthalpy at a given temperature for a protein of known three-dimensional structure can be calculated from the accessible surface areas of constituent atoms according to a method developed recently. Since the hydration free energy and enthalpy for random conformations are computed from those for an extended conformation, the thermodynamic properties of unfolding are evaluated quantitatively. The evaluated hydration properties for proteins of known transition temperature (Tm) and unfolding enthalpy (delta Hm) show an approximately linear dependence on the number of constituent heavy atoms. Since the unfolding free energy is zero at Tm, the enthalpy originating from interatomic interactions of a polypeptide chain and the chain entropy are evaluated from an experimental value of delta Hm and computed properties due to the hydrated water around the molecule at Tm. The chain enthalpy and entropy thus estimated are largely compensated by the hydration enthalpy and entropy, respectively, making the unfolding free energy and enthalpy relatively small. The computed temperature dependences of the unfolding free energy and enthalpy for RNase A, T4 lysozyme, and myoglobin showed a good agreement with the experimental ones.(ABSTRACT TRUNCATED AT 250 WORDS)     

New developments of the electrostatically driven monte carlo method: Test on the membrane-bound portion of melittin

The electrostatically driven Monte Carlo (EDMC) method has been greatly improved by adding a series of new features, including a procedure for cluster analysis of the accepted conformations. This information is used to guide the search for the global energy minimum. Alternative procedures for generating perturbed conformations to sample the conformational space were also included. These procedures enhance the efficiency of the method by generating a larger number of low-energy conformations. The improved EDMC method has been used to explore the conformational space of a 20-residue polypeptide chain whose sequence corresponds to the membrane-bound portion of melittin. The ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm was used to describe the conformational energy of the chain. After an exhaustive search involving 14 independent runs, the lowest energy conformation (LEC) (−91.0 kcal/mol) of the entire study was encountered in four of the runs, while conformations higher in energy by no more than 1.8 kcal/mol were found in the remaining runs with the exception of one of them (run 8). The LEC is identical to the conformation found recently by J. Lee, H.A. Scheraga, and S. Rackovsky [(1998) “Conformational Analysis of the 20-Residue Membrane-Bound Portion of Melittin by Conformational Space Annealing,” Biopolymers, Vol. 46, pp. 103–115] as the lowest energy conformation obtained in their study using the conformational space annealing method. These results suggest that this conformation corresponds to the global energy minimum of the ECEPP/3 potential function for this specific sequence; it also appears to be the conformation of lowest free energy. © 1998 John Wiley & Sons, Inc. Biopoly 46: 117–126, 1998     

Coupling Between Folding and Ionization Equilibria: Effects of pH on the Conformational Preferences of Polypeptides

A new approach to the conformational study of polypeptides is presented. It considers explicitly the coupling between the conformation of the molecule and the ionization equilibria at a given pH value. Calculations of the solvation free energy and free energy of ionization of a 17-residue polypeptide are carried out using a fast multigrid boundary element method (MBE). The MBE method uses an adaptive tessellation of the molecular surface by boundary elements with non-regular size to solve the Poisson equation rapidly, and with a high degree of accuracy. The MBE method is integrated into the ECEPP (Empirical Conformational Energy Program for Peptides) algorithm to compute the coupling between the ionization state and the conformation of the molecule.This approach has been applied to study the conformational preference of a short polypeptide for which the available NMR and CD experimental data indicate that conformations containing a right-handed α-helical segment are energetically more favorable at low values of pH. The results of calculations using the present method agree quite well with experiments, in contrast to previous applications with standard techniques (using pre-assigned charges at each pH) that were not able to reproduce the experimental findings. Also, it is shown how the coupling to the conformation leads to different degrees of ionization of a given type of residue, for example glutamic acid, at different positions in the amino acid sequence, at any given pH. The results of this study provide a sound basis to discuss the origin of the stability of polypeptide conformations, and its dependence on the environmental conditions.     

Positional preference of proline in alpha-helices.

Conformational free energy calculations have been carried out for proline-containing alanine-based pentadecapeptides with the sequence Ac-(Ala)n-Pro-(Ala)m-NHMe, where n + m = 14, to figure out the positional preference of proline in alpha-helices. The relative free energy of each peptide was calculated by subtracting the free energy of the extended conformation from that of the alpha-helical one, which is used here as a measure of preference. The highest propensity is found for the peptide with proline at the N-terminus (i.e., Ncap + 1 position), and the next propensities are found at Ncap, N' (Ncap - 1), and C' (Ccap + 1) positions. These computed results are reasonably consistent with the positional propensities estimated from X-ray structures of proteins. The breaking in hydrogen bonds around proline is found to play a role in destabilizing alpha-helical conformations, which, however, provides the favored hydration of the corresponding N-H and C=O groups. The highest preference of proline at the beginning of alpha-helix appears to be due to the favored electrostatic and nonbonded energies between two residues preceding proline and the intrinsic stability of alpha-helical conformation of the proline residue itself as well as no disturbance in hydrogen bonds of alpha-helix by proline. The average free energy change for the substitution of Ala by Pro in a alpha-helix is computed to be 4.6 kcal/mol, which is in good agreement with the experimental value of approximately 4 kcal/mol estimated for an oligopeptide dimer and proteins of barnase and T4 lysozyme.     

Prediction of conformation of rat galanin in the presence and absence of water with the use of monte Carlo methods and the ECEPP/3 force field

The conformation of the 29-residue rat galanin neuropeptide was studied using the Monte Carlo with energy minimization (MCM) and electrostatically driven Monte Carlo (EDMC) methods. According to a previously elaborated procedure, the polypeptide chain was first treated in a united-residue approximation, in order to enable extensive exploration of the conformational space to be carried out (with the use of MCM), Then the low-energy united-residue conformations were converted to the all-atom representations, and EDMC simulations were carried out for the all-atom polypeptide chains, using the ECEPP/3 force field with hydration included. In order to estimate the effect of environment on galanin conformation, the low-energy conformations obtained as a result of these simulations were taken as starting structures for further EDMC runs that did not include hydration. The lowest-energy conformation obtained in aqueous solution calculations had a nonhelical N-terminal part packed against the nonpolar face of a residual helix that extended from Pro¹³ toward the C-terminus. One next lowest-energy structure was a nearly-all-helical conformation, but with a markedly higher energy. In contrast, all of the low-energy conformations in the absence of water were all-helical differing only by the extent to which the helix was kinked around Pro¹³. These results are in qualitative agreement with the available NMR and CD data of galanin in aqueous and nonaqueous solvents.     

Conformational Transitions in Poly d(CGCGCGTTAATT)

Abstract

Conformational studies on poly d(CGCGCGTTAATT) in solution by circular dichroism spectroscopy are reported. The polynucleotide exhibits B conformation in sodium chloride solution and on addition of NiCl₂ a B-Z transition is observed. NiCl₂ titrations carried out in the presence of 5M NaCl show a midpoint of transition at 2.25 mM NiCl₂ and a complete (maximum conversion to Z form) transition at 16 mM NiCl₂. In 60% alcohol the polynucleotide remains in B conformation. The polynucleotide isomerizes into ψ and A conformations in the presence of spermidine and spermine respectively. The thermodynamic parameters calculated from the melting profiles using a two state model show that the polynucleotide is almost equally stable in its B and Z conformations.     

Theoretical conformational analysis of methylamide N-acetyl-L-arginine]

The spatial structure of methylamide N-acetyl-L-argine was studied taking into account the non-valent and electrostati interactions, the torsion energy, and the distorsion of valency angles. Calculation of the favourable conformations of the molecule was carried out with the use of all the combinations of angles phi, psi, chi1 divided by chi4 as an intital approximation. These correspond to the low energy forms of the main chain and to the minima of the torsion potentials of the side chain. Conformational possibilities of arginine and lysine were compared. The calculated stable conformation of N-acetyl-L-arginine-methylamide are compared with the geometry of arginine residues in the proteins with known structure.     

Conformational analysis of the polar head group in phosphatidylcholine bilayers: a structural change induced by cations

The conformation of the polar head group of phosphatidylcholine in a bilayer in the liquid-crystalline state was deduced by analyzing the deuterium quadrupole splittings of the choline group and the phosphorus chemical shift anisotropy of the phosphate group in combination with the restriction of the choline conformation determined in laser Raman studies. The latter efficiently reduced the number of candidates for the actual conformation. A family of conformations was obtained for both the dynamic-structure and rigid-structure models, respectively. The polar head group is oriented roughly parallel to the membrane surface in both models. Furthermore, they are close to conformation A of the crystal structure of 1,2-dimyristoyl-sn-glycero-3-phosphocholine. The dynamic-structure model was concluded to be more reasonable in view of the fact that the polar head-group structures in most crystals comprise two conformations, which are nearly mirror images of each other. Conformational analysis was also carried out for the polar head group in the presence of multivalent cations. A possible conformational change of the polar head group induced by cations is discussed in the light of the present results.     

Molecular simulation uncovers the conformational space of the λ Cro dimer in solution

The significant variation among solved structures of the λ Cro dimer suggests its flexibility. However, contacts in the crystal lattice could have stabilized a conformation which is unrepresentative of its dominant solution form. Here we report on the conformational space of the Cro dimer in solution using replica exchange molecular dynamics in explicit solvent. The simulated ensemble shows remarkable correlation with available x-ray structures. Network analysis and a free energy surface reveal the predominance of closed and semi-open dimers, with a modest barrier separating these two states. The fully open conformation lies higher in free energy, indicating that it requires stabilization by DNA or crystal contacts. Most NMR models are found to be unstable conformations in solution. Intersubunit salt bridging between Arg⁴ and Glu⁵³ during simulation stabilizes closed conformations. Because a semi-open state is among the low-energy conformations sampled in simulation, we propose that Cro-DNA binding may not entail a large conformational change relative to the dominant dimer forms in solution.     

Dependence of protein stability on the structure of the denatured state: free energy calculations of I56V mutation in human lysozyme.

Free energy calculations were carried out to understand the effect of the I56V mutation of human lysozyme on its thermal stability. In the simulation of the denatured state, a short peptide including the mutation site in the middle is employed. To study the dependence of the stability on the denatured-state structure, five different initial conformations, native-like, extended, and three random-coil-like conformations, were examined. We found that the calculated free energy difference, DeltaDeltaGcal, depends significantly on the structure of the denatured state. When native-like structure is employed, DeltaDeltaGcal is in good agreement with the experimental free energy difference, DeltaDeltaGexp, whereas in the other four models, DeltaDeltaGcal differs sharply from DeltaDeltaGexp. It is therefore strongly suggested that the structure around the mutation site takes a native-like conformation rather than an extended or random-coil conformation. From the free energy component analysis, it has been shown that free energy components originating from Lennard-Jones and covalent interactions dominantly determine DeltaDeltaGcal. The contribution of protein-protein interactions to the nonbonded component of DeltaDeltaGcal is about the same as that from protein-water interactions. The residues that are located in a hydrophobic core (F3, L8, Y38, N39, T40, and I89) contribute significantly to the nonbonded free energy component of DeltaDeltaGcal. We also propose a general computational strategy for the study of protein stability that is equally conscious of the denatured and native states.     

Conformational flexibility of dinucleoside dimers during unwinding from the B-form to an intercalation structure.

Intra-dinucleoside dimer conformational energy calculations were performed on a series of conformations corresponding to the unwinding of a dinucleoside dimer from the B-form type structure. Conformational states were characterized in terms of the distance, d, between parallel bases. All 16 possible sequences for (G) and (C) and/or (A) and (T) bases were considered. Both free space, and free space plus aqueous solvation energies were computed. It was found that the B-form like conformer is the most stable structures in free space and in an aqueous medium. However, the energy as function of d is base pair sequence dependent. This sequence dependence suggests a possible inherent specificity for intercalation of the dinucleoside dimer with a drug molecule.