Temperature-sensitive mutations of bacteriophage T4 lysozyme occur at sites with low mobility and low solvent accessibility in the folded protein

Twenty-five different temperature-sensitive point mutations at 20 sites in the lysozyme gene of bacteriophage T4 have been identified. All of the mutations alter amino acid side chains that have lower than average crystallographic thermal factors and reduced solvent accessibility in the folded protein. This suggests that the amino acids with well-defined conformations can form specific intramolecular interactions that make relatively large contributions to the thermal stability of the protein. Residues with high mobility or high solvent accessibility are much less susceptible to destabilizing substitutions, suggesting that, in general, such amino acids contribute less to protein stability. The pattern of the sites of ts substitutions observed in the folded conformation of T4 lysozyme suggests that severe destabilizing mutations that primarily affect the free energy of the unfolded state are rare. These results indicate that proteins can be stabilized by adding new interactions to regions that are rigid or buried in the folded conformation.     

Limited Proteolysis of Disulfide-reduced Ovalbumin by Subtilisin

Our previous study [Takahashiet al., J. Biochem., 109, 846–851 (1991)] has shown that the disulfide-reduced form of ovalbumin was proteolyzed by subtilisin into three major fragments. It was investigated whether or not these three fragments would be folded into one molecule. Gel permeation and ion-exchange chromatography indicated that the three fragments were eluted in a single peak. The proteolyzed protein had a CD spectrum that was almost indistinguishable from the disulfide-reduced, non-proteolyzed, form of ovalbumin. Differential scanning calorimetry, however, revealed, that the proteolyzed ovalbumin was denatured at a lower temperature than that of the disulfide-reduced, non-proteolyzed. protein. Thus, it is concluded that the three fragments were folded into a native-like conformation with decreased stability. Chemical analyses of the fragments purified by reverse-phase HPLC revealed that there was a cleavage site in the disulfide-reduced form of ovalbumin, at least at the amino-terminal side of Cys⁷³, in addition to the well-known cleavage sites in plakalbumin.     

Bioencapsulation of apomyoglobin in nanoporous organosilica sol–gel glasses: Influence of the siloxane network on the conformation and stability of a model protein

Nanoporous sol–gel glasses were used as host materials for the encapsulation of apomyoglobin, a model protein employed to probe in a rational manner the important factors that influence the protein conformation and stability in silica‐based materials. The transparent glasses were prepared from tetramethoxysilane (TMOS) and modified with a series of mono‐, di‐ and tri‐substituted alkoxysilanes, R_nSi(OCH₃)_4?n (R = methyl‐, n = 1; 2; 3) of different molar content (5, 10, 15%) to obtain the decrease of the siloxane linkage (? Si? O? Si? ). The conformation and thermal stability of apomyoglobin characterized by circular dichroism spectroscopy (CD) was related to the structure of the silica host matrix characterized by ²⁹Si MAS NMR and N₂ adsorption. We observed that the protein transits from an unfolded state in unmodified glass (TMOS) to a native‐like helical state in the organically modified glasses, but also that the secondary structure of the protein was enhanced by the decrease of the siloxane network with the methyl modification (n = 0 < n = 1 < n = 2 < n = 3; 0 < 5 < 10 < 15 mol %). In 15% trimethyl‐modified glass, the protein even reached a maximum molar helicity (?24,000 deg. cm² mol^?1) comparable to the stable folded heme‐bound holoprotein in solution. The protein conformation and stability induced by the change of its microlocal environment (surface hydration, crowding effects, microstructure of the host matrix) were discussed owing to this trend dependency. These results can have an important impact for the design of new efficient biomaterials (sensors or implanted devices) in which properly folded protein is necessary. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 895–906, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Possible mechanism for the dynamic stabilization of protein structure

There is an apparent contradiction between the long lifetime of the metastable structure of native proteins and the high rate of structural fluctuations, which result from the small activation energy required to change the native conformation.In this paper we point out that the observed stability of proteins is not a consequence of large potential barriers, but a result of the continuous reconstitution of the degraded structure by chain propagation. Polypeptide chains of proteins having naturally selected amino-acid sequences have regenerative ability which ensures the long lifetime of the native structure by making most of the fluctuations reversible.A simple calculation shows that in a certain fluctuation of an average protein molecule the probability of denaturation is less than 10⁻²⁵, therefore even the most rapid, picosecond time scale fluctuations cause spontaneous denaturation only in million year time scale. Hence, the generally observed spontaneous denaturation in vitro is rather a consequence of covalent structure modification or intermolecular interactions than a result of an intramolecular interconversion from the native conformation to another conformation.     

Polar and nonpolar atomic environments in the protein core: Implications for folding and binding

Hydrophobic interactions are believed to play an important role in protein folding and stability. Semi-empirical attempts to estimate these interactions are usually based on a model of solvation, whose contribution to the stability of proteins is assumed to be proportional to the surface area buried upon folding. Here we propose an extension of this idea by defining an environment free energy that characterizes the environment of each atom of the protein, including solvent, polar or nonpolar atoms of the same protein or of another molecule that interacts with the protein. In our model, the difference of this environment free energy between the folded state and the unfolded (extended) state of a protein is shown to be proportional to the area buried by nonpolar atoms upon folding. General properties of this environment free energy are derived from statistical studies on a database of 82 well-refined protein structures. This free energy is shown to be able to discriminate misfolded from correct structural models, to provide an estimate of the stabilization due to oligomerization, and to predict the stability of mutants in which hydrophobic residues have been substituted by site-directed mutagenesis, provided that no large structural modifications occur. © 1994 Wiley-Liss, Inc.     

pK(a) values for the unfolded state under native conditions explain the pH-dependent stability of PGB1

Understanding the role of electrostatics in protein stability requires knowledge of these interactions in both the folded and unfolded states. Electrostatic interactions can be probed experimentally by characterizing ionization equilibria of titrating groups, parameterized as pK_a values. However, pK_a values of the unfolded state are rarely accessible under native conditions, where the unfolded state has a very low population. Here, we report pK_a values under nondenaturing conditions for two unfolded fragments of the protein G B1 domain that mimic the unfolded state of the intact protein. pK_a values were determined for carboxyl groups by monitoring their pH-dependent ¹³C chemical shifts. Monte Carlo simulations using a Gaussian chain model provide corrections for changes in electrostatic interactions that arise from fragmentation of the protein. Most pK_a values for the unfolded state agree well with model values, but some residues show significant perturbations that can be rationalized by local electrostatic interactions. The pH-dependent stability was calculated from the experimental pK_a values of the folded and unfolded states and compared to experimental stability data. The use of experimental pK_a values for the unfolded state results in significantly improved agreement with experimental data, as compared to calculations based on model data alone.     

Density functional theory study of the potassium complexation of an unsymmetrical 1,3-alternate calix[4]-crown-5-N-azacrown-5 bearing two different crown rings

Theoretical studies of an unsymmetrical calix[4]-crown-5-N-azacrown-5 (1) in a fixed 1,3-alternate conformation and the complexes 1·K⁺(a), 1·K⁺(b), 1·K⁺(c) and 1·K⁺K⁺ were performed using density functional theory (DFT) at the B3LYP/6-31G* level. The fully optimized geometric structures of the free macroligand and its 1:1 and 1:2 complexes, as obtained from DFT calculations, were used to perform natural bond orbital (NBO) analysis. The two main types of driving force metal–ligand and cation–π interactions were investigated. NBO analysis indicated that the stabilization interaction energies (E ₂) for O…K⁺ and N…K⁺ are larger than the other intermolecular interactions in each complex. The significant increase in electron density in the RY* or LP* orbitals of K⁺ results in strong host–guest interactions. In addition, the intermolecular interaction thermal energies (ΔE, ΔH, ΔG) were calculated by frequency analysis at the B3LYP/6-31G* level. For all structures, the most pronounced changes in the geometric parameters upon interaction are observed in the calix[4]arene molecule. The results indicate that both the intermolecular electrostatic interactions and the cation–π interactions between the metal ion and π orbitals of the two pairs that face the inverted benzene rings play a significant role.     

Quantum-chemical and empirical calculations on phospholipids IV. Glycero-phosphorylethanolamine in a quasihexagonal planar lattice

The energetically most stable head group conformations of a racemic mixture of diacyl-glycero-phosphorylethanolamine in a planar quasihexagonal lattice were calculated using empirical 1-6-12 atom-atom potential functions for intra- and intermolecular interactions. The results demonstrate that the conformation of phospholipid head groups in bilayer systems is determined by intramolecular interactions as well as by intermolecular interactions with neighbouring phospholipid molecules and with solvent molecules. The most stable conformers are that with a φ2 = guache^? conformation of the phosphodiester group. All conformers with a φ2 = gauche⁽⁺⁾ or trans conformation have energies more than 15 kcal ☆ mol^?1 above that of the global minimum. The calculated torsional angles ?1 and φ1 are in very good agreement with the results of the X-ray diffraction analysis of 1,2-dilauroyl-DL-phosphatidylethanolamine (DLPE) acetic acid single crystals.     

Salt-Dependent Folding Energy Landscape of RNA Three-Way Junction

RNAs are highly negatively charged chain molecules. Metal ions play a crucial role in RNA folding stability and conformational changes. In this work, we employ the recently developed tightly bound ion (TBI) model, which accounts for the correlation between ions and the fluctuation of ion distributions, to investigate the ion-dependent free energy landscape for the three-way RNA junction in a 16S rRNA domain. The predicted electrostatic free energy landscape suggests that 1), ion-mediated electrostatic interactions cause an ensemble of unfolded conformations narrowly populated around the maximally extended structure; and 2), Mg²⁺ ion-induced correlation effects help bring the helices to the folded state. Nonelectrostatic interactions, such as noncanonical interactions within the junctions and between junctions and helix stems, might further limit the conformational diversity of the unfolded state, resulting in a more ordered unfolded state than the one predicted from the electrostatic effect. Moreover, the folded state is predominantly stabilized by the coaxial stacking force. The TBI-predicted folding stability agrees well with the experimental measurements for the different Na⁺ and Mg²⁺ ion concentrations. For Mg²⁺ solutions, the TBI model, which accounts for the Mg²⁺ ion correlation effect, gives more improved predictions than the Poisson-Boltzmann theory, which tends to underestimate the role of Mg²⁺ in stabilizing the folded structure. Detailed control tests indicate that the dominant ion correlation effect comes from the charge-charge Coulombic correlation rather than the size (excluded volume) correlation between the ions. Furthermore, the model gives quantitative predictions for the ion size effect in the folding energy landscape and folding cooperativity.     

Role of the amino-terminal extrahelical region of type I collagen in directing the 4D overlap in fibrillogenesis

The amino-terminal telopeptide of the collagen α1(I) chain has a highly conserved sequence. This sequence was analyzed by the Chou-Fasman criteria, and a folded β-sheet conformation, including a β-turn, was predicted. This folded “hairpin” region favors both ionic and hydrophobic intermolecular interactions with α1(I) chain residues 930–938 on a neighboring, end-overlapped molecule. An end-overlap interaction of this nature could direct the initial step in fibril formation. The predicted structure also places the potential crosslink-forming lysyl residue, 9^N, in a unique site at the β-turn end of the telopeptide.     

Logical analysis of the mechanism of protein folding III. Prediction of the strong long-range interactions.

It has been found that strong long-range interactions occur in regions having large β-structural potentials. As has been described previously (Nagano, 1974), interactions among regions having both helical and β-structural potentials (αβ-gaβ interactions) are also very important. Accordingly, an idea is presented in this paper that the relative stability of a protein conformation could be estimated by a relatively simple mathematical function of sequence and conformation. The function P(p,q) is called the non-energy part of pseudo-free energy, because minimization of the sum of P(p,q) and energy functions (cf. Levitt, 1974; Warme &; Scheraga, 1974) can be expected to lead to a plausible model of a protein. A merit of the function is that it can help us decide which way to go in manipulating a temporarily built model, e.g. towards a helix-rich protein or towards a β-structure-rich protein. The estimation of P(p,q) as an artificial potential does not use much computer time because only the co-ordinates of the β-carbon atoms (α-carbon atoms if the residue is Gly) are used. It is composed of terms of the long-range interactions P_L and short-range interactions P_S. The term P_L represents the relative strength of helix-helix interactions, helix-β-candidate interactions and β-candidate-β-candidate interactions. It is assumed that both helical and β-structural potentials can be measured as the differences between the predicted function for helix and β-structure, respectively, as defined previously (Nagano, 1973), and the corresponding largest values ever found. A hypothesis that two residues distantly separated in the primary sequence contribute less to the stability of the whole molecule is finally discarded because the true conformation of concanavalin A becomes very unstable compared with its false conformation folded like the main part of subtilisin. The parameters thus determined indicate that the helix-β-candidate interactions are almost as important as the β-candidate-β-candidate interactions. Both helix and loop prediction functions are combined to give the short-range interactions term, P_S, according to whether the region is really helical or not, and to whether it is really looped or not. The function P(p,q) can be used as a criterion for judging whether the predicted conformation is realistic or false, because the parameters can be adjusted to give, within limits, reasonable values of −10 kcal/residue for true conformations and higher than −5 kcal/residue for false conformations.As an application of the present theory of protein folding, the tertiary structure of bacteriophage T4 lysozyme is predicted and presented in Figure 1, prior to the X-ray structure becoming available.     

Intensity fluctuation spectroscopy and transfer RNA conformation. III. Influence of NaCl concentration on the size and shape of the initially salt-free tRNA in solution.

The influence of sodium ion concentration in solution on the initially salt-free conformation of bulk tRNA from baker's yeast has been investigated by means of photon correlation spectroscopy. From the measured values of translational (D^T) and rotational (D^R) diffusion coefficients, the semiaxes of an ellipsoid of revolution, which are hydrodynamically equivalent to the tRNA molecule, were calculated for tRNA solutions in pure H₂O as well as in 0.005, 0.1, 0.5M NaCl and 0.01M MgCl₂ solutions at pH 4.2 and 7.5. These data, combined with our previous studies, suggested a model which describes the formation of an ordered tRNA structure due to increasing NaCl concentrations. Furthermore, we have obtained information concerning intermolecular interactions between tRNA molecules in solution. In low-salt or salt-free tRNA solutions, we detected in the linewidth distribution function an extra-fast component which can be attributed as possibly due to charge fluctuations related to the reaction of ionization of organic bases. In our light-scattering linewidth measurements, we do not see fluctuations of charged and uncharged states directly as concentration fluctuations. Rather, we postulate a modulation of long-range intermolecular electrostatic interactions between the tRNA molecules due to such charge fluctuations. It is this modulation which is related to the fast component of the time correlation function at finite concentrations. A quantitative theory is needed to provide a more definitive explanation of the dynamical behavior of tRNA in salt-free or low-salt solutions.     

Monte Carlo simulation of protein folding in the presence of residue-specific binding sites

Ensemble growth Monte Carlo (EGMC) and dynamic Monte Carlo (DMC) simulations are used to study sequential folding and thermodynamic stability of hydrophobic-polar (HP) chains that fold to a compact structure. Molecularly imprinted cavities are modeled as hard walls having sites that are attractive to specific polar residues on the chain. Using EGMC simulation, we find that the folded conformation can be stabilized using a small number of carefully selected residue-specific sites while a random selection of surface-bound residues may only slightly contribute toward stabilizing the folded conformation, and in some cases may hinder the folding of the chain. DMC simulations of the surface-bound chain confirm increased stability of the folded conformation over a free chain. However, a different trend of the equilibrium population of folded chains as a function of residue-external site interactions is predicted with the two simulation methods.     

Lens Crystallins and Their Microbial Homologs: Structure,Stability, and Function

Referee: Franz Schmid, Biochemicshes Laboratorium, Universitaet Bayeuth, D-95440 Bayeuth, Germany

abg-Crystallins are the major protein components in the vertebrate eye lens — a as a molecular chaperone and b and g as structural proteins. Surprisingly, the latter two share some structural characteristics with a number of microbial stress proteins. The common denominator is not only the Greek key topology of their polypeptide chains but also their high intrinsic stability, which, in certain microbial crystallin homologs, is further enhanced by high-affinity Ca²⁺-binding. Recent studies of natural and mutant vertebrate bg-crystallins as well as spherulin 3a from Physarum polycephalum and Protein S from Myxococcus xanthus allowed the correlation of structure and stability of crystallins to be elucidated in some detail. From the thermo-dynamic point of view, stability increments come from (1) local interactions involved in the close packing of the cooperative units, (2) the all-b secondary structure of the Greek-key motif, (3) intramolecular interactions between domains, (4) intermolecular domain interactions, including 3D domain swapping and (v) excluded volume effects due to “molecular crowding” at the high cellular protein concentrations. Apart from these contributions to the Gibbs free energy of stability, significant kinetic stabilization originates from the high activation energy barrier determining the rate of unfolding from the native to the unfolded state. From the functional point of view, the high stability is responsible for the long-term transparency of the eye lens, on the one hand, and the stress resistance of the microorganisms in their dormant state on the other. Local structural perturbations due to chemical modification, wrong protein interactions, or other irreversible processes may lead to protein aggregation. A leading cataract hypothesis is that only after a-crystallin, a member of the small heat-shock protein family, is titrated out does pathological opacity occur. Understanding the structural basis of protein stability in the healthy eye lens is the route to solve the enormous medical and economical problem of cataract.     

Intermolecular interactions during ultrafiltration of pegylated proteins

Recent studies have demonstrated the feasibility of using membrane ultrafiltration for the purification of pegylated proteins; however, the separations have all been performed at relatively low protein concentrations where intermolecular interactions are unimportant. The objective of this study was to examine the behavior at higher PEG concentrations and to develop an appropriate theoretical framework to describe the effects of intermolecular interactions. Ultrafiltration experiments were performed using pegylated α‐lactalbumin as a model protein with both neutral and charged composite regenerated cellulose membranes. The transmission of the pegylated α‐lactalbumin, PEG, and α‐lactalbumin all increase with increasing PEG concentration due to the increase in the solute partition coefficient arising from unfavorable intermolecular interactions in the bulk solution. The experimental results were in good agreement with a simple model that accounts for the change in Gibbs free energy associated with these intermolecular interactions, including the effects of concentration polarization on the local solute concentrations upstream of the membrane. These intermolecular interactions are shown to cause a greater than expected loss of pegylated product in a batch ultrafiltration system, and they alter the yield and purification factor that can be achieved during a diafiltration process to remove unreacted PEG. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 29:655–663, 2013     

Introducing folding stability into the score function for computational design of RNA‐binding peptides boosts the probability of success

A computational strategy that integrates our peptide search algorithm with atomistic molecular dynamics simulation was used to design rational peptide drugs that recognize and bind to the anticodon stem and loop domain (ASL^Lys3) of human for the purpose of interrupting HIV replication. The score function of the search algorithm was improved by adding a peptide stability term weighted by an adjustable factor λ to the peptide binding free energy. The five best peptide sequences associated with five different values of λ were determined using the search algorithm and then input in atomistic simulations to examine the stability of the peptides' folded conformations and their ability to bind to ASL^Lys3. Simulation results demonstrated that setting an intermediate value of λ achieves a good balance between optimizing the peptide's binding ability and stabilizing its folded conformation during the sequence evolution process, and hence leads to optimal binding to the target ASL^Lys3. Thus, addition of a peptide stability term significantly improves the success rate for our peptide design search. Proteins 2016; 84:700–711. © 2016 Wiley Periodicals, Inc.     

Homology modeling and docking studies of FabH (β-ketoacyl-ACP synthase III) enzyme involved in type II fatty acid biosynthesis of Chlorella variabilis: a potential algal feedstock for biofuel production

The concept of using microalgae as an alternative renewable source of biofuel has gained much importance in recent years. However, its commercial feasibility is still an area of concern for researchers. Unraveling the fatty acid metabolic pathway and understanding structural features of various key enzymes regulating the process will provide valuable insights to target microalgae for augmented oil content. FabH (β-ketoacyl-acyl carrier protein synthase; KAS III) is a condensing enzyme catalyzing the initial elongation step of type II fatty acid biosynthetic process and acyl carrier protein (ACP) facilitates the shuttling of the fatty acyl intermediates to the active site of the respective enzymes in the pathway. In the present study, a reliable three-dimensional structure of FabH from Chlorella variabilis, an oleaginous green microalga was modeled and subsequently the key residues involved in substrate binding were determined by employing protein–protein docking and molecular dynamics (MD) simulation protocols. The FabH-ACP complex having the lowest docking energy score showed the binding of ACP to the electropositive FabH surface with strong hydrogen bond interactions. The MD simulation results indicated that the substrate-complexed FabH adopted a more stable conformation than the free enzyme. Further, the FabH structure retained its stability throughout the simulation although noticeable displacements were observed in the loop regions. Molecular simulation studies suggested the importance of crucial hydrogen bonding of the conserved Arg⁹¹ of FabH with Glu⁵³ and Asp⁵⁶ of ACP for exhibiting high affinity between the enzyme and substrate. The molecular modeling results are consistent with available experimental results on the flexibility of FabH and the present study provides first in silico insights into the structural and dynamical aspect of catalytic mechanism of FabH, which could be used for further site-specific mutagenic experiments to develop engineered high oil-yielding microalgal strains for biofuel production.     

The determinants of bond angle variability in protein/peptide backbones: A comprehensive statistical/quantum mechanics analysis

The elucidation of the mutual influence between peptide bond geometry and local conformation has important implications for protein structure refinement, validation, and prediction. To gain insights into the structural determinants and the energetic contributions associated with protein/peptide backbone plasticity, we here report an extensive analysis of the variability of the peptide bond angles by combining statistical analyses of protein structures and quantum mechanics calculations on small model peptide systems. Our analyses demonstrate that all the backbone bond angles strongly depend on the peptide conformation and unveil the existence of regular trends as function of ψ and/or φ. The excellent agreement of the quantum mechanics calculations with the statistical surveys of protein structures validates the computational scheme here employed and demonstrates that the valence geometry of protein/peptide backbone is primarily dictated by local interactions. Notably, for the first time we show that the position of the H^α hydrogen atom, which is an important parameter in NMR structural studies, is also dependent on the local conformation. Most of the trends observed may be satisfactorily explained by invoking steric repulsive interactions; in some specific cases the valence bond variability is also influenced by hydrogen‐bond like interactions. Moreover, we can provide a reliable estimate of the energies involved in the interplay between geometry and conformations. Proteins 2015; 83:1973–1986. © 2015 Wiley Periodicals, Inc.     

The contribution of proline residues to protein stability is associated with isomerization equilibrium in both unfolded and folded states

It has long been understood that the proline residue has lower configurational entropy than any other amino acid residue due to pyrrolidine ring hindrance. The peptide bond between proline and its preceding amino acid (Xaa-Pro) typically exists as a mixture of cis- and trans-isomers in the unfolded protein. Cis–trans isomerization of Xaa-Pro peptide bonds are infrequent, but still occur in folded proteins. Therefore, the effects of the cis–trans isomerization equilibrium in both unfolded and folded states should be taken into account when estimating the stability contribution of a specific proline residue. In order to study the stability contribution of the four proline residues to the hyperthermophilic protein Ssh10b, in this work, we expressed and purified a series of Pro→Ala mutants of Ssh10b, and performed correlative unfolding experiments in detail. We proposed a new unfolding model including proline isomerization. The model predicts that the contribution of a proline residue to protein stability is associated with the thermodynamic equilibrium between cis- and trans-isomers both in the unfolded and folded states, agreeing well with the experimental results.