Conformational sampling of peptides in cellular environments

Biological systems provide a complex environment that can be understood in terms of its dielectric properties. High concentrations of macromolecules and cosolvents effectively reduce the dielectric constant of cellular environments, thereby affecting the conformational sampling of biomolecules. To examine this effect in more detail, the conformational preference of alanine dipeptide, poly-alanine, and melittin in different dielectric environments is studied with computer simulations based on recently developed generalized Born methodology. Results from these simulations suggest that extended conformations are favored over alpha-helical conformations at the dipeptide level at and below dielectric constants of 5-10. Furthermore, lower-dielectric environments begin to significantly stabilize helical structures in poly-alanine at epsilon = 20. In the more complex peptide melittin, different dielectric environments shift the equilibrium between two main conformations: a nearly fully extended helix that is most stable in low dielectrics and a compact, V-shaped conformation consisting of two helices that is preferred in higher dielectric environments. An additional conformation is only found to be significantly populated at intermediate dielectric constants. Good agreement with previous studies of different peptides in specific, less-polar solvent environments, suggest that helix stabilization and shifts in conformational preferences in such environments are primarily due to a reduced dielectric environment rather than specific molecular details. The findings presented here make predictions of how peptide sampling may be altered in dense cellular environments with reduced dielectric response.     

Conformational behaviour of the architectural units of peptides and proteins. Assessment of current understanding by ab initio quantum mechanical methods and refinement of the dipeptide energy surface.

The conformational energy surfaces of analogues of the dipeptide unit of polypeptides and proteins are calculated by ab initio methods using extended basis sets.The calculations are not particularly sensitive to the choice of (extended) basis set.The calculations are shown to support a particular empirical method parameterized with respect to crystal data. Non-hydrogen bonded conformations agree to within 3 kcal mol^?1, even for conformations in which quite considerable degrees of atomic overlap occur.Hydrogen bonded conformations, are, however, in less satisfactory agreement and it is the ab initio calculations which appear to be at fault.A simple correction is applied to the ab initio energy for hydrogen bonded conformations, and with the use of the empirical energy surface a full quantum mechanical conformational energy map is interpolated for the alanyl dipeptide.The effect of flexibility in the peptide backbone is taken into account, and supports recent empirical findings that distortions in valence angles must be considered in calculations of the conformational behaviour of peptides.     

Monte Carlo calculations on the conformations of models for the glycopeptide linkage of glycoproteins

In order to search for probable conformations of the peptide, the amino acid side chain, and the carbohydrate linkage in glycoproteins, conformational energy surfaces of glycopeptide model compounds were studied by Monte Carlo methods using the Metropolis algorithm. The potential energies were composed of empirical energy functions which include nonbonded interactions, electrostatics, hydrogen bonding, and torsional energies specified by parameters which have been used for peptides. Calculations were performed on 1-N-acetyl-2-acetamido-beta-D-glucopyranosyl amine and the glycosylated dipeptide N-acetyl-delta-N-(2-acetamido-beta-D-glucopyranosyl)-L-asparaginyl-N'-methyl amide as models for N-glycosylated peptides and on methyl-2-acetamido-alpha-D-galactopyranoside as well as the glycosylated dipeptides N-acetyl-gamma-O-(2-acetamido-alpha-D-galactopyranosyl)-L-threonyl-N'-methyl amide and its seryl analog as models for O-glycosylated glycoproteins. The probable conformations of these compounds were analyzed by single-angle probability tables and by two-dimensional conformation density maps projected from the Markov chains which contained up to six independently varied conformational dihedral angles. The presence of high barriers to rotation required the use of search strategies which resulted in a rather low acceptance rate for new conformations in the Metropolis algorithm in order to avoid trapping of the Markov chain in local energy minima. This problem contributed to the failure of these calculations to attain complete convergence to the thermodynamic limit for the glycosylated dipeptide models in which six dihedral angles were independently varied. Analysis of the results shows that the conformational space available to the highly crowded axial glycosides of the alpha-O-GalNAc type is much more restricted than that for the N-asparaginyl glycopeptides. The most probable conformation for the O-glycosylated peptides is is a beta-turn while N-glycosylated peptides may be either in a beta-turn or an extended conformation.     

Nanosecond study of fluorescently labeled troponin C.

The time-resolved extrinsic fluorescence of rabbit skeletal troponin C was studied with the protein labeled at Cys-98 with N-(iodoacetyl)-N'-(5-sulfo-1-naphthyl)ethylenediamine. Both the intensity and anisotropy decays followed a biexponential decay law, regardless of the ionic condition, pH, viscosity or temperature. The lifetimes and their fractional amplitudes were insensitive to Mg2+, and the lifetimes were also insensitive to Ca2+. In response to Ca2+ binding to all four sites, the fractional amplitude (alpha 1) associated with the short lifetime (tau 1) decreased by a factor of two, thus increasing the ratio of the two amplitudes alpha 2/alpha 1 from 1.6 to 4.3. These amplitude changes suggest the existence of two conformational states of TnC-IAEDANS, with the conformation associated with the long-decay component (tau 2) being promoted by saturation of the two Ca(2+)-specific sites. At pH 5.2 the ratio alpha 2/alpha 1 for the apo-protein was 3.5 indicating different relative populations of the two decay components when compared with pH 7.2. In the presence of Ca2+ at the lower pH, alpha 2/alpha 1 decreased to 2.1, suggesting a shift of the conformations in favor of the short-decay component. Thus Ca2+ elicited different conformational changes in TnC at the two pH values. The recovered anisotropies suggest that there were fast molecular motions that were not resolved in the present experiments, and some of these motions were sensitive to Ca2+ binding to the specific sites. These results support the notion of communication between the N-domain and the C-terminal end of the central helix of troponin C.     

Conformational substates of calmodulin revealed by single-pair fluorescence resonance energy transfer: influence of solution conditions and oxidative modification

A calmodulin (CaM) mutant (T34,110C-CaM) doubly labeled with fluorescence probes AlexaFluor 488 and Texas Red in opposing domains (CaM-DA) has been used to examine conformational heterogeneity in CaM by single-pair fluorescence resonance energy transfer (spFRET). Burst-integrated FRET efficiencies of freely diffusing CaM-DA single molecules yielded distributions of distance between domains of CaM-DA. We recently reported distinct conformational substates of Ca(2+)-CaM-DA and apoCaM-DA, with peaks in the distance distributions centered at approximately 28 A, 34-38 A, and 55 A [Slaughter et al. (2004) J. Phys. Chem. B 108, 10388-10397]. In the present study, shifts in the amplitudes and center distances of the conformational substates were detected with variation in solution conditions. The amplitude of an extended conformation was observed to change as a function of Ca(2+) over a free Ca(2+) range that is consistent with binding to the high affinity, C-terminal Ca(2+) binding sites, suggesting the existence of communication between lobes of CaM. Lowering pH shifted the relative amplitudes of the conformations, with a marked increase in the presence of the compact conformations and an almost complete absence of the extended conformation. In addition, the single-molecule distance distribution of apoCaM-DA at reduced ionic strength was shifted to longer distance and showed evidence of an increase in conformational heterogeneity relative to apoCaM-DA at physiological ionic strength. Oxidation of methionine residues in CaM-DA produced a substantial increase in the amplitude of the extended conformation relative to the more compact conformation. The results are considered in light of a hypothesis that suggests that electrostatic interactions between charged amino acid side chains play an important role in determining the most stable CaM conformation under varying solution conditions.     

pH dependence of the urea and guanidine hydrochloride denaturation of ribonuclease A and ribonuclease T1

To investigate the pH dependence of the conformational stability of ribonucleases A and T1, urea and guanidine hydrochloride denaturation curves have been determined over the pH range 2-10. The maximum conformational stability of both proteins is about 9 kcal/mol and occurs near pH 4.5 for ribonuclease T1 and between pH 7 and 9 for ribonuclease A. The pH dependence suggests that electrostatic interactions among the charged groups make a relatively small contribution to the conformational stability of these proteins. The dependence of delta G on urea concentration increases from about 1200 cal mol-1 M-1 at high pH to about 2400 cal mol-1 M-1 at low pH for ribonuclease A. This suggests that the unfolded conformations of RNase A become more accessible to urea as the net charge on the molecule increases. For RNase T1, the dependence of delta G on urea concentration is minimal near pH 6 and increases at both higher and lower pH. An analysis of information of this type for several proteins in terms of a model developed by Tanford [Tanford, C. (1964) J. Am. Chem. Soc. 86, 2050-2059] suggests that the unfolded states of proteins in urea and GdnHCl solutions may differ significantly in the extent of their interaction with denaturants. Thus, the conformations assumed by unfolded proteins may depend to at least some extent on the amino acid sequence of the protein.     

Molecular dynamics studies of the conformation of sorbitol

Molecular dynamics simulations of a 3 molal aqueous solution of d-sorbitol (also called d-glucitol) have been performed at 300 K, as well as at two elevated temperatures to promote conformational transitions. In principle, sorbitol is more flexible than glucose since it does not contain a constraining ring. However, a conformational analysis revealed that the sorbitol chain remains extended in solution, in contrast to the bent conformation found experimentally in the crystalline form. While there are 243 staggered conformations of the backbone possible for this open-chain polyol, only a very limited number were found to be stable in the simulations. Although many conformers were briefly sampled, only eight were significantly populated in the simulation. The carbon backbones of all but two of these eight conformers were completely extended, unlike the bent crystal conformation. These extended conformers were stabilized by a quite persistent intramolecular hydrogen bond between the hydroxyl groups of carbon C-2 and C-4. The conformational populations were found to be in good agreement with the limited available NMR data except for the C-2–C-3 torsion (spanned by the O-2–O-4 hydrogen bond), where the NMR data support a more bent structure.     

Contribution of histidine residues to the conformational stability of ribonuclease T1 and mutant Glu-58----Ala

The pK values of the histidine residues in ribonuclease T1 (RNase T1) are unusually high: 7.8 (His-92), 7.9 (His-40), and 7.3 (His-27) [Inagaki et al. (1981) J. Biochem. 89, 1185-1195]. In the RNase T1 mutant Glu-58----Ala, the first two pK values are reduced to 7.4 (His-92) and 7.1 (His-40). These lower pKs were expected since His-92 (5.5 A) and His-40 (3.7 A) are in close proximity to Glu-58 at the active site. The conformational stability of RNase T1 increases by over 4 kcal/mol between pH 9 and 5, and this can be entirely accounted for by the greater affinity for protons by the His residues in the folded protein (average pK = 7.6) than in the unfolded protein (pk approximately 6.6). Thus, almost half of the net conformational stability of RNase T1 results from a difference between the pK values of the histidine residues in the folded and unfolded conformations. In the Glu-58----Ala mutant, the increase in stability between pH 9 and 5 is halved (approximately 2 kcal/mol), as expected on the basis of the lower pK values for the His residues in the folded protein (average pK = 7.1). As a consequence, RNase T1 is more stable than the mutant below pH 7.5, and less stable above pH 7.5. These results emphasize the importance of measuring the conformational stability as a function of pH when comparing proteins differing in structure.     

Reversible stages of the low-pH-triggered conformational change in influenza virus hemagglutinin

The refolding of the prototypic fusogenic protein hemagglutinin (HA) at the pH of fusion is considered to be a concerted and irreversible discharge of a loaded spring, with no distinct intermediates between the initial and final conformations. Here, we show that HA refolding involves reversible conformations with a lifetime of minutes. After reneutralization, low pH-activated HA returns from the conformations wherein both the fusion peptide and the kinked loop of the HA2 subunit are exposed, but the HA1 subunits have not yet dissociated, to a structure indistinguishable from the initial one in functional, biochemical and immunological characteristics. The rate of the transition from reversible conformations to irreversible refolding depends on the pH and on the presence of target membrane. Importantly, recovery of the initial conformation is blocked by the interactions between adjacent HA trimers. The existence of the identified reversible stage of refolding can be crucial for allowing multiple copies of HA to synchronize their release of conformational energy, as required for fusion.     

Conformations of disulfide bridges in proteins

The conformational characteristics of disulfide bridges in proteins have been analyzed using a dataset of 22 protein structures, available at a resolution of less than or equal to 2.0 A, containing a total of 72 disulfide crosslinks. The parameters used in the analysis include (phi, psi) values at Cys residues, bridge dihedral angles chi ss, chi i1, chi j1, chi i2, and chi j2, the distances C alpha i-C alpha j and C beta i-C beta j between the C alpha and C beta atoms of Cys(i) and Cys(j). Eight families of bridge conformations with three or more occurrences have been identified on the basis of these stereochemical parameters. The most populated family corresponds to the "left handed spiral" identified earlier by Richardson [1981) Adv. Protein Chem. 34, 167-330). Disulfide bridging across antiparallel extended strands is observed in alpha-lytic protease, crambin, and beta-trypsin and this structure is shown to be very similar to those obtained in small cystine peptides. Solvent accessible surface area calculations show that the overwhelming majority of disulfide bridges are inaccessible to solvent.     

Conformation of glycopeptides with L- and D-stereoisomers of Ala and Glu amino acids

On the base of the semi-empirical method of conformational analysis low-energetic conformations have been found for the glycopeptides with the O-beta-N-acetyl-D-muramic acid or the O-beta-N-acetyl--D-glucosaminyl-(1-4)-N-acetyl-D-muramic acid in the glycan part and L-Ala--D-GluNH2, L-Ala--L-GluNH2, D-Ala--D-GluNH2 in the peptide part. It was shown that for either molecule preferentiability of any conformation is determined mainly by interactions between dipeptide and muramic acid residue.     

Characterization of pH-dependent conformational heterogeneity in Rhodospirillum rubrum cytochrome c2 using 15N and 1H NMR

The 15N-enriched ferricytochrome c2 from Rhodospirillum rubrum has been studied by 15N and 1H NMR spectroscopy as a function of pH. The 15N resonances of the heme and ligand tau nitrogen are broadened beyond detection because of paramagnetic relaxation. The 15N resonance of the ligand histidine phi nitrogen was unambiguously identified at 184 ppm (pH 5.6). The 15N resonances of the single nonligand histidine are observed only at low pH, as in the ferrocytochrome because of the severe broadening caused by tautomerization. The dependence of the 15N and 1H spectra of the ferricytochrome on pH indicated that the ligand histidine tau NH does not dissociate in the neutral pH range and is involved in a hydrogen bond, similar to that in the reduced state. Because neither deprotonated nor non-hydrogen-bonded forms of the ligand histidine are observed in the spectra of either oxidation state, the participation of such forms in producing heterogeneous populations having different electronic g tensors is ruled out. Transitions having pKa's of 6.2, 8.6, and 9.2 are observed in the ferricytochrome. The localized conformational change around the omega loops is observed in the neutral pH range, as in the ferrocytochrome. Structural heterogeneity leads to multiple resonances of the heme ring methyl at position 8. The exchange rate between the conformations is temperature dependent. The transition with a pKa of 6.2 is assigned to the His-42 imidazole group. The displacement of the ligand methionine, which occurs with a pKa of 9.2, causes gross conformational change near the heme center.(ABSTRACT TRUNCATED AT 250 WORDS)     

pH-Dependent Conformational Changes in Bacterial Hsp90 Reveal a Grp94-Like Conformation at pH 6 That Is Highly Active in Suppression of Citrate Synthase Aggregation

The molecular chaperone Hsp90 depends upon large conformational rearrangements for its function. One driving force for these rearrangements is the intrinsic ATPase activity of Hsp90, as seen with other chaperones. However, unlike other chaperones, structural and kinetic studies have shown that the ATPase cycle of Hsp90 is not conformationally deterministic. That is, rather than dictating the conformational state, ATP binding and hydrolysis shift the equilibrium between a preexisting set of conformational states in an organism-dependent manner. While many conformations of Hsp90 have been described, little is known about how they relate to chaperone function. In this study, we show that the conformational equilibrium of the bacterial Hsp90, HtpG, can be shifted with pH. Using small-angle X-ray scattering, we identify a two-state pH-dependent conformational equilibrium for apo HtpG. Our structural modeling reveals that this equilibrium is observed between the previously observed extended state and a second state that is strikingly similar to the recently solved Grp94 crystal structure. In the presence of nonhydrolyzable 5′-adenylyl-β,γ-imidodiphosphate, a third state, which is identical with the solved AMPPNP-bound structure from yeast Hsp90, is populated. Electron microscopy confirmed the observed conformational equilibria. We also identify key histidine residues that control this pH-dependent equilibrium; using mutagenesis, we successfully modulate the conformational equilibrium at neutral pH. Using these mutations, we show that the Grp94-like state provides stronger aggregation protection compared to the extended apo conformation in the context of a citrate synthase aggregation assay. These studies provide a more detailed view of HtpG's conformational dynamics and provide the first linkage between a specific conformation and chaperone function.     

NMR structures of the histidine-rich peptide LAH4 in micellar environments: membrane insertion, pH-dependent mode of antimicrobial action, and DNA transfection

The LAH4 family of histidine-rich peptides exhibits potent antimicrobial and DNA transfection activities, both of which require interactions with cellular membranes. The bilayer association of the peptides has been shown to be strongly pH-dependent, with in-planar alignments under acidic conditions and transmembrane orientations when the histidines are discharged. Therefore, we investigated the pH- and temperature-dependent conformations of LAH4 in DPC micellar solutions and in a TFE/PBS solvent mixture. In the presence of detergent and at pH 4.1, LAH4 adopts helical conformations between residues 9 and 24 concomitantly with a high hydrophobic moment. At pH 6.1, a helix-loop-helix structure forms with a hinge encompassing residues His¹⁰-Ala¹³. The data suggest that the high density of histidine residues and the resulting electrostatic repulsion lead to both a decrease in the pK values of the histidines and a less stable α-helical conformation of this region. The hinged structure at pH 6.1 facilitates membrane anchoring and insertion. At pH 7.8, the histidines are uncharged and an extended helical conformation including residues 4-21 is again obtained. LAH4 thus exhibits a high degree of conformational plasticity. The structures provide a stroboscopic view of the conformational changes that occur during membrane insertion, and are discussed in the context of antimicrobial activity and DNA transfection.     

Single molecule analyses of the conformational substates of calmodulin bound to the phosphorylase kinase complex

The four integral delta subunits of the phosphorylase kinase (PhK) complex are identical to calmodulin (CaM) and confer Ca(2+) sensitivity to the enzyme, but bind independently of Ca(2+). In addition to binding Ca(2+), an obligatory activator of PhK's phosphoryltransferase activity, the delta subunits transmit allosteric signals to PhK's remaining alpha, beta, and gamma subunits in activating the enzyme. Under mild conditions about 10% of the delta subunits can be exchanged for exogenous CaM. In this study, a CaM double-mutant derivatized with a fluorescent donor-acceptor pair (CaM-DA) was exchanged for delta to assess the conformational substates of PhKdelta by single molecule fluorescence resonance energy transfer (FRET) +/-Ca(2+). The exchanged subunits were determined to occupy distinct conformations, depending on the absence or presence of Ca(2+), as observed by alterations of the compact, mid-length, and extended populations of their FRET distance distributions. Specifically, the combined predominant mid-length and less common compact conformations of PhKdelta became less abundant in the presence of Ca(2+), with the delta subunits assuming more extended conformations. This behavior is in contrast to the compact forms commonly observed for many of CaM's Ca(2+)-dependent interactions with other proteins. In addition, the conformational distributions of the exchanged PhKdelta subunits were distinct from those of CaM-DA free in solution, +/-Ca(2+), as well as from exogenous CaM bound to the PhK complex as delta'. The distinction between delta and delta' is that the latter binds only in the presence of Ca(2+), but stoichiometrically and at a different location in the complex than delta.     

Conformation of di-n-propylglycine residues (Dpg) in peptides: crystal structures of a type I' beta-turn forming tetrapeptide and an alpha-helical tetradecapeptide.

The crystal structures of two oligopeptides containing di-n-propylglycine (Dpg) residues, Boc-Gly-Dpg-Gly-Leu-OMe (1) and Boc-Val-Ala-Leu-Dpg-Val-Ala-Leu-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe (2) are presented. Peptide 1 adopts a type I'beta-turn conformation with Dpg(2)-Gly(3) at the corner positions. The 14-residue peptide 2 crystallizes with two molecules in the asymmetric unit, both of which adopt alpha-helical conformations stabilized by 11 successive 5 --> 1 hydrogen bonds. In addition, a single 4 --> 1 hydrogen bond is also observed at the N-terminus. All five Dpg residues adopt backbone torsion angles (phi, psi) in the helical region of conformational space. Evaluation of the available structural data on Dpg peptides confirm the correlation between backbone bond angle N-C(alpha)-C' (tau) and the observed backbone phi,psi values. For tau > 106 degrees, helices are observed, while fully extended structures are characterized by tau < 106 degrees. The mean tau values for extended and folded conformations for the Dpg residue are 103.6 degrees +/- 1.7 degrees and 109.9 degrees +/- 2.6 degrees, respectively.     

江浙蝮蛇毒酸性与碱性磷脂酶A_2溶液构象变化的差示扫描量热法研究

运用差示扫描量热法,在不同ｐＨ值的缓冲溶液内和各种浓度的碱土族氯化物溶液内,研究了来自江浙蝮蛇（ＡｇｋｉｓｔｒｏｄｏｎｈａｌｙｓＰａｌｌａｓ）毒的酸性与碱性磷脂酶外Ａ２（ＰＬＡ２）的热变性过程。得到表征这两种酶溶液构象变化的热力学参数。依据这些参数研究了两者的溶液构象及其变化。在ｐＨ４．５以下,分子净荷正电的这两种酶在溶液中不形成可热致伸展的有序构象;ｐＨ高于４．５时,Ａｓｐ和Ｇｌｕ的侧链羧基以负离子形式存在有利于有序构象的稳定。Ｈｉｓ是决定ＰＬＡ２活力和热稳定性的重要残基。磷酸根离子和这两种酶有结合作用而降低有序构象的热稳定性。碱土族阳离子除和这两种酶结合外,还以依赖于离子强度的方式复杂地影响酶的溶液构象,但其作用不完全是静电性的,是或多或少地随离子的不同而不同的。计算给出酸性ＰＬＡ２的△Ｈｃｄ．     

Novel small renin inhibitors containing 4,5- or 3,5-dihydroxy-2-substituted-6-phenylhexanamide replacements at the P2P3 sites

Renin inhibitors containing a 4,5- or a 3,5-dihydroxy-2-substituted-6-phenylhexanamide fragment at the P₂---P₃ sites have been prepared and evaluated. The four possible diastereomeric diols of the two series of inhibitors were synthesized to determine the optimal configuration of the carbinol centers for these replacements. The most potent inhibitors of each series, 1a and 2c have a molecular weight of only 503 and IC₅₀ values of 23 and 20 nM in a human plasma renin assay at pH 6.0. Their very low aqueous solubility limited their further evaluation. The efficacy of these P₂---P₃ replacements is a result of their ability to maintain the important hydrogen-bonds with the enzyme. Due to conformational differences with the dipeptide, adjustment at the P₂ side chain was required. These 4,5- and 3,5-dihydroxyhexanamide segments could be seen as novel N-terminal dipeptide replacements.     

Urea denaturation of barnase: pH dependence and characterization of the unfolded state.

To investigate the pH dependence of the conformational stability of barnase, urea denaturation curves were determined over the pH range 2-10. The maximum conformational stability of barnase is 9 kcal mol-1 and occurs between pH 5 and 6. The dependence of delta G on urea concentration increases from 1850 cal mol-1 M-1 at high pH to about 3000 cal mol-1 M-1 near pH 3. This suggests that the unfolded conformations of barnase become more accessible to urea as the net charge on the molecule increases. Previous studies suggested that in 8 M urea barnase unfolds more completely than ribonuclease T1, even with the disulfide bonds broken [Pace, C.N., Laurents, D. V., & Thomson, J.A. (1990) Biochemistry 29, 2564-2572]. In support of this, solvent perturbation difference spectroscopy showed that in 8 M urea the Trp and Tyr residues in barnase are more accessible to perturbation by dimethyl sulfoxide than in ribonuclease T1 with the disulfide bonds broken.