Hydrophobic regions on protein surfaces. Derivation of the solvation energy from their area distribution in crystallographic protein structures

Supra-biological performance of immobilized enzymes enabled by chaperone-like specific non-covalent interactions

Designing complex synthetic materials for enzyme immobilization could unlock the utility of biocatalysis in extreme environments. Inspired by biology, we investigate the use of random copolymer brushes as dynamic immobilization supports that enable supra-biological catalytic performance of immobilized enzymes. This is demonstrated by immobilizing Bacillus subtilis Lipase A on brushes doped with aromatic moieties, which can interact with the lipase through multiple non-covalent interactions. Incorporation of aromatic groups leads to a 50 °C increase in the optimal temperature of lipase, as well as a 50-fold enhancement in enzyme activity. Single-molecule FRET studies reveal that these supports act as biomimetic chaperones by promoting enzyme refolding and stabilizing the enzyme's folded and catalytically active state. This effect is diminished when aromatic residues are mutated out, suggesting the importance of π-stacking and π-cation interactions for stabilization. Our results underscore how unexplored enzyme-support interactions may enable uncharted opportunities for using enzymes in industrial biotransformations.

Contribution of Hydrophobic Interactions to Protein Mechanical Stability

The role of hydrophobic and polar interactions in providing thermodynamic stability to folded proteins has been intensively studied, but the relative contribution of these interactions to the mechanical stability is less explored. We used steered molecular dynamics simulations with constant-velocity pulling to generate force-extension curves of selected protein domains and monitor hydrophobic surface unravelling upon extension. Hydrophobic contribution was found to vary between one fifth and one third of the total force while the rest of the contribution is attributed primarily to hydrogen bonds. Moreover, hydrophobic force peaks were shifted towards larger protein extensions with respect to the force peaks attributed to hydrogen bonds. The higher importance of hydrogen bonds compared to hydrophobic interactions in providing mechanical resistance is in contrast with the relative importance of the hydrophobic interactions in providing thermodynamic stability of proteins. The different contributions of these interactions to the mechanical stability are explained by the steeper free energy dependence of hydrogen bonds compared to hydrophobic interactions on the relative positions of interacting atoms. Comparative analyses for several protein domains revealed that the variation of hydrophobic forces is modest, while the contribution of hydrogen bonds to the force peaks becomes increasingly important for mechanically resistant protein domains.

Understanding Design Rules for Optimizing the Interface between Immobilized Enzymes and Random Copolymer Brushes

A long-standing goal in the field of biotechnology is to develop and understand design rules for the stabilization of enzymes upon immobilization to materials. While immobilization has sometimes been successful as a strategy to stabilize enzymes, the design of synthetic materials that stabilize enzymes remains largely empirical. We sought to overcome this challenge by investigating the mechanistic basis for the stabilization of immobilized lipases on random copolymer brush surfaces comprised of poly(ethylene glycol) methacrylate (PEGMA) and sulfobetaine methacrylate (SBMA), which represent novel heterogeneous supports for immobilized enzymes. Using several related but structurally diverse lipases, including Bacillus subtilis lipase A (LipA), Rhizomucor miehei lipase, Candida rugosa lipase, and Candida antarctica lipase B (CALB), we showed that the stability of each lipase at elevated temperatures was strongly dependent on the fraction of PEGMA in the brush layer. This dependence was explained by developing and applying a new algorithm to quantify protein surface hydrophobicity, which involved using unsupervised cluster analysis to identify clusters of hydrophobic atoms. Characterization of the lipases showed that the optimal brush composition correlated with the free energy of solvation per enzyme surface area, which ranged from -17.1 kJ/mol·nm² for LipA to -11.8 kJ/mol·nm² for CALB. Additionally, using this algorithm, we found that hydrophobic patches consisting of aliphatic residues had a higher free energy than patches consisting of aromatic residues. By providing the basis for rationally tuning the interface between enzymes and materials, this understanding will transform the use of materials to reliably ruggedize enzymes under extreme conditions.

Thermodynamic profiling for fragment-based lead discovery and optimization

Introduction: The enthalpic and entropic components of the ligand-protein binding free energy reflect the type and quality of the interactions and relate to the physicochemical properties of the ligands. These findings have significance in medicinal chemistry optimizations since they suggest that the thermodynamic profiling of the binding may help monitor and control the unfavorable size and hydrophobicity increase typically accompanying affinity improvements and leading to suboptimal pharmacokinetic properties.Areas covered: This review describes the ligand-protein binding event in terms of elementary steps, their associated interactions, and their enthalpic and entropic consequences. The relationships among the breaking and forming interactions, the binding thermodynamic profile, and the physicochemical properties of the ligands are also discussed.Expert opinion: Analysis of the size dependence of available affinity and favorable enthalpy highlights the limitation of the simultaneous optimization of these quantities. Indeed, moderate, rather than very high affinities can be conciliated with favorable physicochemical and pharmacokinetic profiles as it is supported by the affinity range of historical oral drugs. Although thermodynamic quantities are not suitable endpoints for medicinal chemistry optimizations owing to the complexity of the binding thermodynamics, thermodynamic profiling together with structural studies can be advantageously used to understand the details of the binding process and to optimize it.

Determination of Protein-Protein Interactions in a Mixture of Two Monoclonal Antibodies

The coformulation of monoclonal antibody (mAb) mixtures provides an attractive route to achieving therapeutic efficacy where the targeting of multiple epitopes is necessary. Controlling and predicting the behavior of such mixtures requires elucidating the molecular basis for the self- and cross-protein-protein interactions and how they depend on solution variables. While self-interactions are now beginning to be well understood, systematic studies of cross-interactions between mAbs in solution do not exist. Here, we have used static light scattering to measure the set of self- and cross-osmotic second virial coefficients in a solution containing a mixture of two mAbs, mAbA and mAbB, as a function of ionic strength and pH. mAbB exhibits strong association at a low ionic strength, which is attributed to an electrostatic attraction that is enhanced by the presence of a strong short-ranged attraction of nonelectrostatic origin. Under all solution conditions, the measured cross-interactions are intermediate self-interactions and follow similar patterns of behavior. There is a strong electrostatic attraction at higher pH values, reflecting the behavior of mAbB. Protein-protein interactions become more attractive with an increasing pH due to reducing the overall protein net charges, an effect that is attenuated with an increasing ionic strength due to the screening of electrostatic interactions. Under moderate ionic strength conditions, the reduced cross-virial coefficient, which reflects only the energetic contribution to protein-protein interactions, is given by a geometric average of the corresponding self-coefficients. We show the relationship can be rationalized using a patchy sphere model, where the interaction energy between sites i and j is given by the arithmetic mean of the i-i and j-j interactions. The geometric mean does not necessarily apply to all mAb mixtures and is expected to break down at a lower ionic strength due to the nonadditivity of electrostatic interactions.

Discovery of a novel splice variant of Fcar (CD89) unravels sequence segments necessary for efficient secretion: A story of bad signal peptides and good ones that nevertheless do not make it

The IgA receptor, Fcar (CD89) consists of 5 sequence segments: 2 segments (S1, S2) forming the potential signal peptide, 2 extracellular EC domains that include the IgA binding site, and the transmembrane and cytoplasmic tail (TM/C) region. Numerous Fcar splice variants have been reported with various combinations of the sequence segments mentioned above. Here, we report a novel splice variant termed variant APD isolated from a healthy volunteer that lacks only the IgA-binding EC1 domain. Despite possessing the complete signal peptide S1+S2, the variant APD is only found in the intracellular space whereas the wild-type variant 1 is efficiently secreted and variant 4 leaks to the extracellular space. Further mutational experiments involving signal peptide replacements, cleavage site modifications, and studies on alternative isoforms demonstrate that despite the completeness of the signal peptide motif, the presence of the EC1 domain is essential for efficient extracellular export.

Elastic properties and heterogeneous stiffness of the phi29 motor connector channel

The DNA packaging motor of the bacteriophage ϕ29, comprising head-tail connector, ATPase, and pRNA, transports the viral DNA inside the procapsid against pressure differences of up to ∼60 atm during replication. Several models for the DNA packaging mechanism have been proposed, which attribute different roles to the connector, and require specific mechanical properties of the connector. To characterize these properties at the atomic level, and to understand how the connector withstands this large pressure, we have carried out molecular dynamics simulations of the whole connector both in equilibrium and under mechanical stress. The simulations revealed a quite heterogeneous distribution of stiff and soft regions, resembling that of typical composite materials that are also optimized to resist mechanical stress. In particular, the conserved middle α-helical region is found to be remarkably stiff, similar only to structural proteins forming viral shell, silk, or collagen. In contrast, large parts of the peripheral interface to the ϕ29 procapsid turned out to be rather soft. Force probe and umbrella sampling simulations showed that large connector deformations are remarkably reversible, and served to calculate the free energies required for these deformations. In particular, for an untwisting deformation by 12°, as postulated by the untwist-twist model, more than four times' larger energy is required than is available from hydrolysis of one ATP molecule. Combined with previous experiments, this result is incompatible with the untwist-twist model. In contrast, our simulations support the recently proposed one-way revolution model and suggest in structural terms how the connector blocks DNA leakage. In particular, conserved loops at the rim of the central channel, which are in direct contact with the DNA, are found to be rather flexible and tightly anchored to the rigid central region. These findings suggest a check-valve mechanism, with the flexible loops obstructing the channel by interacting with the viral DNA.

Empirical estimation of the energetic contribution of individual interface residues in structures of protein-protein complexes

We report a simple algorithm to scan interfaces in protein-protein complexes for identifying binding 'hot spots'. The change in side-chain solvent accessible area (DeltaASA) of interface residues has been related to change in binding energy due to mutating interface residues to Ala (DeltaDeltaG (X --> ALA)) based on two criteria-hydrogen bonding across the interface and location in the interface core-both of which are major determinants in specific, high-affinity binding. These relationships are used to predict the energetic contribution of individual interface residues. The predictions are tested against 462 experimental X --> ALA mutations from 28 interfaces with an average unsigned error of 1.04 kcal/mol. More than 80% of interface hot spots (with experimental DeltaDeltaG > or = 2 kcal/mol) could be identified as being energetically important. From the experimental values, Asp, Lys, Tyr and Trp are found to contribute most of the binding energy, burying >45 A2 on average. The method described here would be useful to understand and interfere with protein interactions by assessing the energetic importance of individual interface residues.

Preservation of protein clefts in comparative models

BACKGROUND

Comparative, or homology, modelling of protein structures is the most widely used prediction method when the target protein has homologues of known structure. Given that the quality of a model may vary greatly, several studies have been devoted to identifying the factors that influence modelling results. These studies usually consider the protein as a whole, and only a few provide a separate discussion of the behaviour of biologically relevant features of the protein. Given the value of the latter for many applications, here we extended previous work by analysing the preservation of native protein clefts in homology models. We chose to examine clefts because of their role in protein function/structure, as they are usually the locus of protein-protein interactions, host the enzymes' active site, or, in the case of protein domains, can also be the locus of domain-domain interactions that lead to the structure of the whole protein.

RESULTS

We studied how the largest cleft of a protein varies in comparative models. To this end, we analysed a set of 53507 homology models that cover the whole sequence identity range, with a special emphasis on medium and low similarities. More precisely we examined how cleft quality - measured using six complementary parameters related to both global shape and local atomic environment, depends on the sequence identity between target and template proteins. In addition to this general analysis, we also explored the impact of a number of factors on cleft quality, and found that the relationship between quality and sequence identity varies depending on cleft rank amongst the set of protein clefts (when ordered according to size), and number of aligned residues.

CONCLUSION

We have examined cleft quality in homology models at a range of seq.id. levels. Our results provide a detailed view of how quality is affected by distinct parameters and thus may help the user of comparative modelling to determine the final quality and applicability of his/her cleft models. In addition, the large variability in model quality that we observed within each sequence bin, with good models present even at low sequence identities (between 20% and 30%), indicates that properly developed identification methods could be used to recover good cleft models in this sequence range.

Conservation and relative importance of residues across protein-protein interfaces

A core region surrounded by a rim characterizes biological interfaces. We ascertain the importance of the core by showing the sequence entropies of the residues comprising the core to be smaller than those in the rim. Such a distinction is not seen in the 2-fold-related, nonphysiological interfaces formed in crystal lattices of monomeric proteins, thereby providing a procedure for characterizing the oligomeric state from crystal structures of protein molecules. This method is better than those that rely on the comparison of the sequence entropies in the interface and the rest of the protein surface, especially in cases where the surface harbors additional binding sites. To a good approximation there is a correlation between the accessible surface area lost because of complexation and DeltaDeltaG values obtained through alanine-scanning mutagenesis (26-38 cal per A(2) of the surface buried) for residues located in the core, a relationship that is not discernable for rim residues. If, however, a residue participates in hydrogen bonding across the interface, the extent of stabilization is 52 cal/mol per 1 A(2) of the nonpolar surface area buried by the residue. As opposed to an amino acid classification used earlier, an environment-based grouping of residues yields a better discrimination in the sequence entropy between the core and the rim.

A double-deletion method to quantifying incremental binding energies in proteins from experiment: example of a destabilizing hydrogen bonding pair

The contribution of a specific hydrogen bond in apoflavodoxin to protein stability is investigated by combining theory, experiment and simulation. Although hydrogen bonds are major determinants of protein structure and function, their contribution to protein stability is still unclear and widely debated. The best method so far devised to estimate the contribution of side-chain interactions to protein stability is double mutant cycle analysis, but the interaction energies so derived are not identical to incremental binding energies (the energies quantifying net contributions of two interacting groups to protein stability). Here we introduce double-deletion analysis of 'isolated' residue pairs as a means to precisely quantify incremental binding. The method is exemplified by studying a surface-exposed hydrogen bond in a model protein (Asp96/Asn128 in apoflavodoxin). Combined substitution of these residues by alanines slightly destabilizes the protein due to a decrease in hydrophobic surface burial. Subtraction of this effect, however, clearly indicates that the hydrogen-bonded groups in fact destabilize the native conformation. In addition, molecular dynamics simulations and classic double mutant cycle analysis explain quantitatively that, due to frustration, the hydrogen bond must form in the native structure because when the two groups get approximated upon folding their binding becomes favorable. We would like to remark that 1), this is the first time the contribution of a specific hydrogen bond to protein stability has been measured by experiment; and 2), more hydrogen bonds need to be analyzed to draw general conclusions on protein hydrogen bond energetics. To that end, the double-deletion method should be of help.

A thermodynamic and kinetic analysis of the folding pathway of an SH3 domain entropically stabilised by a redesigned hydrophobic core

The folding thermodynamics and kinetics of the alpha-spectrin SH3 domain with a redesigned hydrophobic core have been studied. The introduction of five replacements, A11V, V23L, M25V, V44I and V58L, resulted in an increase of 16% in the overall volume of the side-chains forming the hydrophobic core but caused no remarkable changes to the positions of the backbone atoms. Judging by the scanning calorimetry data, the increased stability of the folded structure of the new SH3-variant is caused by entropic factors, since the changes in heat capacity and enthalpy upon the unfolding of the wild-type and mutant proteins were identical at 298 K. It appears that the design process resulted in an increase in burying both the hydrophobic and hydrophilic surfaces, which resulted in a compensatory effect upon the changes in heat capacity and enthalpy. Kinetic analysis shows that both the folding and unfolding rate constants are higher for the new variant, suggesting that its transition state becomes more stable compared to the folded and unfolded states. The phi(double dagger-U) values found for a number of side-chains are slightly lower than those of the wild-type protein, indicating that although the transition state ensemble (TSE) did not change overall, it has moved towards a more denatured conformation, in accordance with Hammond's postulate. Thus, the acceleration of the folding-unfolding reactions is caused mainly by an improvement in the specific and/or non-specific hydrophobic interactions within the TSE rather than by changes in the contact order. Experimental evidence showing that the TSE changes globally according to its hydrophobic content suggests that hydrophobicity may modulate the kinetic behaviour and also the folding pathway of a protein.

A procedure for detection and quantitation of cavity volumes proteins. Application to measure the strength of the hydrophobic driving force in protein folding

Accurate identification of cavities is important in the study of protein structure, stability, design, and ligand binding. Identification and quantitation of cavities is a nontrivial problem because most cavities are connected to the protein exterior. We describe a computational procedure for quantitating cavity volumes and apply this to derive an estimate of the hydrophobic driving force in protein folding. A grid-based Monte Carlo procedure is used to position water molecules on the surface of a protein. A Voronoi procedure is used to identify and quantitate empty space within the solvated protein. Additional cavities not detected by other existing procedures can be identified. Most of these are close to surface concavities. Residue volumes for both the interior and the surface residues as well as cavity volumes are in good agreement with volumes calculated from fully hydrated protein structures obtained from molecular dynamic simulations. We show that the loss of stability because of cavity-creating mutations correlates better with cavity volumes determined by this procedure than with cavity volumes determined by other methods. Available structural and thermodynamic data for a number of cavity-containing mutants were analyzed to obtain estimates of 26.1 cal x mol(-1) x A(-3) and 18.5 cal x mol(-1) x A(-2) for the relative contributions of cavity formation and the hydrophobic effect to the observed stability changes. The present estimate for the hydrophobic driving force is at the lower end of estimates derived from model compound studies and considerably lower than previous estimates of approximately 50 cal x mol(-1) x A(-2) derived from protein mutational data. In the absence of structural rearrangement, on average, deletion of a single methylene group is expected to result in losses in stability of 0.41 and 0.70 kcal x mol(-1) resulting from decrease in hydrophobicity and packing, respectively.

Protein-protein interactions in concentrated electrolyte solutions

Protein-protein interactions were measured for ovalbumin and for lysozyme in aqueous salt solutions. Protein-protein interactions are correlated with a proposed potential of mean force equal to the free energy to desolvate the protein surface that is made inaccessible to the solvent due to the protein-protein interaction. This energy is calculated from the surface free energy of the protein that is determined from protein-salt preferential-interaction parameter measurements. In classical salting-out behavior, the protein-salt preferential interaction is unfavorable. Because addition of salt raises the surface free energy of the protein according to the surface-tension increment of the salt, protein-protein attraction increases, leading to a reduction in solubility. When the surface chemistry of proteins is altered by binding of a specific ion, salting-in is observed when the interactions between (kosmotrope) ion-protein complexes are more repulsive than those between the uncomplexed proteins. However, salting-out is observed when interactions between (chaotrope) ion-protein complexes are more attractive than those of the uncomplexed proteins.

Hydrophobic forces between protein molecules in aqueous solutions of concentrated electrolyte

Protein-protein interactions have been measured for a mutant (D101F) lysozyme and for native lysozyme in concentrated solutions of ammonium sulfate at pH 7 and sodium chloride at pH 4.5. In the mutant lysozyme, a surface aspartate residue has been replaced with a hydrophobic phenylalanine residue. The protein-protein interactions of D101F lysozyme are more attractive than those of native lysozyme for all conditions studied. The salt-induced attraction is correlated with a solvation potential of mean force given by the work required to desolvate the part of the protein surfaces that is buried by the protein-protein interaction. This work is proportional to the aqueous surface-tension increment of the salt and the fractional non-polar surface coverage of the protein. Experimental measurements of osmotic second virial coefficients validate a proposed potential of mean force that ascribes the salt-induced attraction between protein molecules to an enhancement of the hydrophobic attraction. This model provides a first approximation for predicting the protein-protein potential of mean force in concentrated aqueous electrolyte solutions; this potential is useful for determining solution conditions favorable for protein crystallization.

Unspecific hydrophobic stabilization of folding transition states

Here we present a method for determining the inference of non-native conformations in the folding of a small domain, alpha-spectrin Src homology 3 domain. This method relies on the preservation of all native interactions after Tyr/Phe exchanges in solvent-exposed, contact-free positions. Minor changes in solvent exposure and free energy of the denatured ensemble are in agreement with the reverse hydrophobic effect, as the Tyr/Phe mutations slightly change the polypeptide hydrophilic/hydrophobic balance. Interestingly, more important Gibbs energy variations are observed in the transition state ensemble (TSE). Considering the small changes induced by the H/OH replacements, the observed energy variations in the TSE are rather notable, but of a magnitude that would remain undetected under regular mutations that alter the folded structure free energy. Hydrophobic residues outside of the folding nucleus contribute to the stability of the TSE in an unspecific nonlinear manner, producing a significant acceleration of both unfolding and refolding rates, with little effect on stability. These results suggest that sectors of the protein transiently reside in non-native areas of the landscape during folding, with implications in the reading of phi values from protein engineering experiments. Contrary to previous proposals, the principle that emerges is that non-native contacts, or conformations, could be beneficial in evolution and design of some fast folding proteins.

Ab initio prediction of the solution structures and populations of a cyclic pentapeptide in DMSO based on an implicit solvation model

Using a recently developed statistical mechanics methodology, the solution structures and populations of the cyclic pentapeptide cyclo(D-Pro(1)-Ala(2)-Ala(3)-Ala(4)-Ala(5)) in DMSO are obtained ab initio, i.e., without using experimental restraints. An important ingredient of this methodology is a novel optimization of implicit solvation parameters, which in our previous publication [Baysal, C.; Meirovitch, H. J Am Chem Soc 1998, 120, 800-812] has been applied to a cyclic hexapeptide in DMSO. The molecule has been described by the simplified energy function E(tot) = E(GRO) + summation operator(k) sigma(k)A(k), where E(GRO) is the GROMOS force-field energy, sigma(k) and A(k) are the atomic solvation parameter (ASP) and the solvent accessible surface area of atom k. This methodology, which relies on an extensive conformational search, Monte Carlo simulations, and free energy calculations, is applied here with E(tot) based on the ASPs derived in our previous work, and for comparison also with E(GRO) alone. For both models, entropy effects are found to be significant. For E(tot), the theoretical values of proton-proton distances and (3)J coupling constants agree very well with the NMR results [Mierke, D. F.; Kurz, M.; Kessler, H. J Am Chem Soc 1994, 116, 1042-1049], while the results for E(GRO) are significantly worse. This suggests that our ASPs might be transferrable to other cyclic peptides in DMSO as well, making our methodology a reliable tool for an ab initio structure prediction; obviously, if necessary, parts of this methodology can also be incorporated in a best-fit analysis where experimental restraints are used.

The effect of vicinal polar and charged groups on hydrophobic hydration

The use of a linear relationship between free energy of hydrophobic hydration and solvent-accessible apolar surface area has been helpful in interpreting the thermodynamics of biological macromolecules. However, a recent study (Y.-K. Cheng, P. J. Rossky, Nature 1998, Vol. 392, pp. 696-699) has established a substantial enthalpic dependence on biomolecular surface topography, originating from solvent hydrogen-bonding loss in a restrictive geometry. In this study, we use molecular dynamics simulations of 2-Zn insulin in water solvent to explore the further effect of vicinal polar or charged groups on hydrophobic hydration at a biomolecular surface. In contrast to the case for solvent more isolated from such polar solute influences, the binding energies of the water that is proximal to the hydrophobic dimeric interface of insulin and vicinal to polar and charged groups are comparable to the bulk solvent value, a result of compensating interaction primarily with the solute counterions. The results suggest a special importance for such polar/charged groups in biological processes involving hydrophobic surface regions of restricted geometry and also suggest a general route for tuning the hydrophobicity of interfaces.

Free energy based populations of interconverting microstates of a cyclic peptide lead to the experimental NMR data

Analysis of nuclear Overhauser enhancement (NOE) intensities data of interconverting microstates of a peptide is a difficult problem in nmr. A new statistical mechanics methodology has been proposed recently, consisting of several steps: (1) potential energy wells on the energy surface of the molecule are identified (the corresponding regions are called wide microstates); (2) each wide microstate is then spanned by a Monte Carlo (MC) or molecular dynamics simulation starting from a representative structure, and the corresponding relative populations are obtained from the free energy calculated with the local states method; and (3) the overall NOEs and 3J coupling constants are obtained as averages over the corresponding contributions of the samples, weighted by the populations. Extending this methodology to cyclic peptides, we are treating here the hexapeptide cyclo(D-Pro1-Phe2-Ala3-Ser4-Phe5-Phe6) in DMSO, which was studied by Kessler et al. using nmr (Journal of the American Chemical Society, 1992, Vol. 114, pp. 4805-4818). They found that at least two structures are required to explain their NOE data, a conclusion also corroborated by our analysis (Journal of the American Chemical Society, 1998, Vol. 120, pp. 800-812) and led to a novel derivation of atomic solvation parameters (ASPs) for DMSO. Thus, the overall interactions within the peptide-solvent system are described approximately by Etot = EGRO + summation operator sigmaiAi, where EGRO is the energy of the GROMOS force field, Ai is the solvent-accessible surface area of atom i, and sigmai is the ASP. In the present paper the validity of these ASPs within the framework of the entire methodology is verified. This requires taking into account 23 microstates. A very good agreement is obtained between experimental and calculated NOEs and 3J coupling constants. The free energy based populations lead to the best results, which means that entropic effects should not be ignored. We have also studied the behavior of the internal angular fluctuations of the proton-proton vectors and discovered that they have a negligible effect on the calculated NOEs; this is due to the relatively concentrated wide microstates spanned by the MC simulations. The applicability of our ASPs to other cyclic peptides in DMSO is being studied in another work and preliminary results are discussed.

Binding of proteins to copolymers of varying hydrophobicity

Hydrophobic interactions between proteins and amphiphilic polyelectrolytes were studied by frontal analysis continuous capillary electrophoresis (Gao et al., Analytical Chemistry, 1997, Vol. 69, pp. 2945-2951). Binding isotherms were obtained for beta-lactoglobulin and for bovine serum albumin interacting with a series of alternating copolymers of maleic acid and alkyl-vinyl ethers of varying hydrophobicity. Although binding between proteins and copolymers increases with increasing alkyl chain length, a minimum alkyl chain length of 3-4 methylenes is required for significant hydrophobic interactions to occur. These copolymers, like other polyamphiphiles, can form intrapolymer micelles, and the extent of such micellization decreases with increasing degree of carboxylate ionization. Binding results obtained at different pHs suggest that competition exists between intrapolymer micelle formation and protein-polymer hydrophobic interactions.

Rigid-body oscillations of α-helices: implications for protein thermal stability

A quasi-continuity model protein consisting of two alpha-helices undergoing rigid-body torsional oscillations demonstrates that factors stabilizing the model protein, such as increased helix rigidity and hydrophobicity, are the same factors that stabilize thermophilic proteins relative to their mesophilic analogs. The model predicts oscillatory motions with frequencies in the microwave (10(10) Hz) range. These oscillations decrease in frequency with increasing helix rigidity because of compensating increases in the force constant and moment of inertia, thus explaining the retention of activity in the more rigid thermophilic enzymes. Implications for protein design, based on the predictions of the model, are discussed.

Influence of protein structure databases on the predictive power of statistical pair potentials

A long standing goal in protein structure studies is the development of reliable energy functions that can be used both to verify protein models derived from experimental constraints as well as for theoretical protein folding and inverse folding computer experiments. In that respect, knowledge-based statistical pair potentials have attracted considerable interests recently mainly because they include the essential features of protein structures as well as solvent effects at a low computing cost. However, the basis on which statistical potentials are derived have been questioned. In this paper, we investigate statistical pair potentials derived from protein three-dimensional structures, addressing in particular questions related to the form of these potentials, as well as to the content of the database from which they are derived. We have shown that statistical pair potentials depend on the size of the proteins included in the database, and that this dependence can be reduced by considering only pairs of residue close in space (i.e., with a cutoff of 8 A). We have shown also that statistical potentials carry a memory of the quality of the database in terms of the amount and diversity of secondary structure it contains. We find, for example, that potentials derived from a database containing alpha-proteins will only perform best on alpha-proteins in fold recognition computer experiments. We believe that this is an overall weakness of these potentials, which must be kept in mind when constructing a database.