Beyond the nearest-neighbor Zimm-Bragg model for helix-coil transition in peptides

The nearest-neighbor (micro = 1) variant of the Zimm and Bragg (ZB) model has been extensively used to describe the helix-coil transition in biopolymers. In this work, we investigate the helix-coil transition for a 21-residue alanine peptide (AP) with the ZB model up to fourth nearest neighbor (micro = 1, 2, 3, and 4). We use a matrix approach that takes into account combinations of any number of helical stretches of any length and therefore gives the exact statistical weight of the chain within the assumptions of the ZB model. The parameters of the model are determined by fitting the temperature-dependent circular dichroism and Fourier transform infrared experimental spectra of the AP. All variants of the model fit the experimental data, thus giving similar results in terms of the macroscopic observables, such as temperature-dependent fractional helicity. However, the resulting microscopic parameters, such as distributions of the individual residue helical probabilities and free energy surfaces, vary significantly depending on the variant of the model. Overall, the mean residue enthalpy and entropy (in the absolute value) both increase with micro, but combined yield essentially the same "effective" value of the ZB propagation parameters for all micro. Greater helical probabilities for individual residues are predicted for larger micro, in particular, near the center of the sequence. The ZB nucleation parameters increase with increasing micro, which results in a lower free energy barrier to helix nucleation and lower apparent "cooperativity" of the transition. The significance of the long-range interactions for the predictions of ZB model for helix-coil transition, the calculated model parameters and the limitations of the model are discussed.     

Theoretical study of the partial molar volume change associated with the pressure-induced structural transition of ubiquitin

The partial molar volume (PMV) change associated with the pressure-induced structural transition of ubiquitin is analyzed by the three-dimensional reference interaction site model (3D-RISM) theory of molecular solvation. The theory predicts that the PMV decreases upon the structural transition, which is consistent with the experimental observation. The volume decomposition analysis demonstrates that the PMV reduction is primarily caused by the decrease in the volume of structural voids in the protein, which is partially canceled by the volume expansion due to the hydration effects. It is found from further analysis that the PMV reduction is ascribed substantially to the penetration of water molecules into a specific part of the protein. Based on the thermodynamic relation, this result implies that the water penetration causes the pressure-induced structural transition. It supports the water penetration model of pressure denaturation of proteins proposed earlier.     

Morphologies and conformation transition of lentinan in aqueous NaOH solution

Molecular morphologies and conformation transition of lentinan, a beta-(1-->3)-D-glucan from Lentinus edodes, were studied in aqueous NaOH solution by atomic force microscopy (AFM), viscometry, multiangle laser light scattering, and optical rotation measurements. The results revealed that lentinan exists as triple-helical chains and as single random-coil chains at NaOH concentration lower than 0.05M and higher than 0.08M, respectively. Moreover, the dramatic changes in weight-average molecular weight Mw, radius of gyration [s2](1/2), intrinsic viscosity [eta], as well as specific optical rotation at 589 nm [alpha]589 occurred in a narrow range of NaOH concentration between 0.05 and 0.08M NaOH, indicating that the helix-coil conformation transition of lentinan was carried out more easily than that of native schizophyllan and scleroglucan, and was irreversible. For the first time, we confirmed that the denatured lentinan molecule, which was dissolved in 0.15M NaOH to be disrupted into single coil chains, could be renatured as triple helical chain by dialyzing against abundant water in the regenerated cellulose tube at ambient temperature (15 degrees C). In view of the AFM image, lentinan in aqueous solution exhibited the linear, circular, and branched species of triple helix compared with native linear schizophyllan or scleroglucan.     

Temperature and pressure effects on conformational equilibria of alanine dipeptide in aqueous solution

We investigated the temperature and pressure effects on conformational equilibria of N-acetyl-L-alanine-N'-methylamide (AAlaMA) in aqueous solution by Raman spectroscopy. Scattering intensities in the skeletal stretching mode of AAlaMA in aqueous solution were decomposed into some component bands by the spectra analysis. Our results indicate that each component band for AAlaMA adopts not only the P(II) and alpha(R) conformations but also the C(7eq) conformation. From temperature and pressure dependencies of the band intensities, we determined the enthalpy differences and the volume differences between the conformers. The C(7eq) conformer is enthalpically most stable due to the intramolecular hydrogen bond. The partial molar volume of the C(7eq) conformer is the smallest through the solvent-exclusion effect rather than the solute-solvent electrostatic interaction effect.     

Theoretical study for volume changes associated with the helix-coil transition of peptides.

We calculate the partial molar volumes and their changes associated with the coil(extended)-to-helix transition of two types of peptide, glycine-oligomer and glutamic acid-oligomer, in aqueous solutions by using the Kirkwood-Buff solution theory coupled with the three-dimensional reference interaction site model (3D-RISM) theory. The volume changes associated with the transition are small and positive. The volume is analyzed by decomposing it into five contributions following the procedure proposed by Chalikian and Breslauer: the ideal volume, the van der Waals volume, the void volume, the thermal volume, and the interaction volume. The ideal volumes and the van der Waals volumes do not change appreciably upon the transition. In the both cases of glycine-peptide and glutamic acid-peptide, the changes in the void volumes are positive, while those in the thermal volumes are negative, and tend to balance those in the void volumes. The change in the interaction volume of glycine-peptide does not significantly contribute, while that of glutamic acid-peptide makes a negative contribution.     

Molecular dynamics simulation study on the high-pressure behaviour of an AK16 peptide

While most proteins unfold under high-pressure conditions, some high-pressure experiments suggest that an AK16 peptide forms more helical structures. In order to understand this abnormality, molecular dynamics simulations with the simulated tempering method for the isobaric–isothermal ensemble were performed in a wide pressure range from 0.1 MPa to 1.4 GPa. It was found that the fraction of the folded state decreases once and increases after that with increasing pressure. The partial molar volume change from the folded state to unfolded state increases monotonically from a negative value to a positive value with pressure. The behaviour under high-pressure conditions is consistent with the experimental results. The radius of gyration of highly helical structures decreases with increasing pressure. Moreover, interatomic distances of AK16 become shorter at high pressure than at low pressure. These behaviours indicate that the helical structures are squeezed by high pressure.     

Effect of pressure on helix‐coil transition of an alanine‐based peptide: An FTIR study

Protein–protein interactions are mediated by complementary amino acids defining complementary surfaces. Typically not all members of a family of related proteins interact equally well with all members of a partner family; thus analysis of the sequence record can reveal the complementary amino acid partners that confer interaction specificity. This article develops methods for learning and using probabilistic graphical models of such residue “cross‐coupling” constraints between interacting protein families, based on multiple sequence alignments and information about which pairs of proteins are known to interact. Our models generalize traditional consensus sequence binding motifs, and provide a probabilistic semantics enabling sound evaluation of the plausibility of new possible interactions. Furthermore, predictions made by the models can be explained in terms of the underlying residue interactions. Our approach supports different levels of prior knowledge regarding interactions, including both one‐to‐one (e.g., pairs of proteins from the same organism) and many‐to‐many (e.g., experimentally identified interactions), and we present a technique to account for possible bias in the represented interactions. We apply our approach in studies of PDZ domains and their ligands, fundamental building blocks in a number of protein assemblies. Our algorithms are able to identify biologically interesting cross‐coupling constraints, to successfully identify known interactions, and to make explainable predictions about novel interactions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Thermally induced coil-to-helix transition of sodium gellan gum with different molar masses in aqueous salt solutions

Using 5 samples of well-purified Na-gellans (Na-gellans G1-G5, weight-average molar mass M(w) = 120 x 10(3)-32 x 10(3) at 40 degrees C), the effects of molar mass on the coil-to-double-helix transition in aqueous solutions with 25 mM NaCl were studied by light scattering and circular dichroism (CD) measurements, viscometry, and differential scanning calorimetry (DSC). From the temperature dependence of M(w), molar ellipticity at 201 nm [theta]201, intrinsic viscosity [eta], and DSC exothermic curves, it was found that the coil-to-double-helix transitions for G1-G5 samples took place at almost the same temperature. The [eta] and M(w) obtained in the temperature range from 40 to 25 degrees C can be explained by a simple coil/double-helix equilibrium model using the double-helix contents determined from CD data. The van't Hoff's transition enthalpy deltaH(vH) of Na-gellans depended on M(w). It is concluded that the coil-to-double-helix transitions of Na-gellans are all-or-none type transitions, and are accelerated with increasing M(w).     

The effect of temperature and lipid on the conformational transition of gramicidin A in lipid vesicles

The present study investigated the effect of temperature and lipid/peptide molar ratio on the conformational changes of the membrane peptide gramicidin A from a double-stranded helix to a single-stranded helical dimmer in 1,2-dimyristoyl-glycerol-3-phosphochloine (DMPC) vesicles. Tryptophan fluorescence spectroscopy results suggested that the conformational transition fitted a three-state (two-step) "folding" model. Rate constants, k(1) and k(2), were determined for each of the two steps. Since k(1) and k(2) increased with an increase in temperature, we hypothesized that the process corresponded to the breakage and formation of the backbone hydrogen bonds. The k(1) was from 10 to 45 folds faster than k(2), except for lipid/peptide molar ratios above 89.21, where k(2) increased rapidly. At molar ratios below 89.21, k(2) was insensitive to changes in lipid concentration. To account for this phenomenon, we proposed that while the driving interaction at high molar ratios is between the indole rings of the tryptophan residues and the lipid head groups, at low molar ratios there may be an intermolecular interaction between the tryptophan residues that causes gramicidin A to form an organized aggregated network. This aggregated network, caused by the tryptophan-tryptophan interaction, may be the main effect responsible for the slow down of the conformation change.     

Breaking non-native hydrophobic clusters is the rate-limiting step in the folding of an alanine-based peptide

The formation mechanism of an alanine-based peptide has been studied by all-atom molecular dynamics simulations with a recently developed all-atom point-charge force field and the Generalize Born continuum solvent model at an effective salt concentration of 0.2M. Thirty-two simulations were conducted. Each simulation was performed for 100 ns. A surprisingly complex folding process was observed. The development of the helical content can be divided into three phases with time constants of 0.06-0.08, 1.4-2.3, and 12-13 ns, respectively. Helices initiate extreme rapidly in the first phase similar to that estimated from explicit solvent simulations. Hydrophobic collapse also takes place in this phase. A folding intermediate state develops in the second phase and is unfolded to allow the peptide to reach the transition state in the third phase. The folding intermediate states are characterized by the two-turn short helices and the transition states are helix-turn-helix motifs-both of which are stabilized by hydrophobic clusters. The equilibrium helical content, calculated by both the main-chain Phi-Psi torsion angles and the main-chain hydrogen bonds, is 64-66%, which is in remarkable agreement with experiments. After corrected for the solvent viscosity effect, an extrapolated folding time of 16-20 ns is obtained that is in qualitative agreement with experiments. Contrary to the prevailing opinion, neither initiation nor growth of the helix is the rate-limiting step. Instead, the rate-limiting step for this peptide is breaking the non-native hydrophobic clusters in order to reach the transition state. The implication to the folding mechanisms of proteins is also discussed.     

Engineered variants of InlB with an additional leucine‐rich repeat discriminate between physiologically relevant and packing contacts in crystal structures of the InlB:MET complex

The physiological relevance of contacts in crystal lattices often remains elusive. This was also the case for the complex between the invasion protein internalin B (InlB) from Listeria monocytogenes and its host cell receptor, the human receptor tyrosine kinase (RTK) MET. InlB is a MET agonist and induces bacterial host cell invasion. Activation of RTKs generally involves ligand‐induced dimerization of the receptor ectodomain. The two currently available crystal structures of the InlB:MET complex show the same arrangement of InlB and MET in a 1:1 complex, but different dimeric 2:2 assemblies. Only one of these 2:2 assemblies is predicted to be stable by a computational procedure. This assembly is mainly stabilized by a contact between the Cap domain of InlB from one and the Sema domain of MET from another 1:1 complex. Here, we probe the physiological relevance of this interaction. We generated variants of the leucine‐rich repeat (LRR) protein InlB by inserting an additional repeat between the first and the second LRR. This should allow formation of the 1:1 complex but disrupt the potential 2:2 complex involving the Cap‐Sema contact due to steric distortions. A crystal structure of one of the engineered proteins showed that it folded properly. Binding affinity to MET was comparable to that of wild‐type InlB. The InlB variant induced MET phosphorylation and cell scatter like wild‐type InlB. These results suggest that the Cap‐Sema interaction is not physiologically relevant and support the previously proposed assembly, in which a 2:2 InlB:MET complex is built around a ligand dimer.