PubMed Health⌕ Search

Biomedical subjects

John M Prausnitz

Publications and source records attributed to John M Prausnitz.

13 recordsLinked to original sources

Molecular simulation of protein aggregation.

Computer simulation offers unique possibilities for investigating molecular-level phenomena difficult to probe experimentally. Drawing from a wealth of studies concerning protein folding, computational studies of protein aggregation are emerging. These studies have been successful in capturing aspects of aggregation known from experiment and are being used to refine experimental methods aimed at abating aggregation. Here we review molecular-simulation studies of protein aggregation conducted in our laboratory. Specific attention is devoted to issues with implications for biotechnology.

Binding Sites↗

Effect of salt identity on the phase diagram for a globular protein in aqueous electrolyte solution.

Monte Carlo simulations are used to establish the potential of mean force between two globular proteins in an aqueous electrolyte solution. This potential includes nonelectrostatic contributions arising from dispersion forces first, between the globular proteins, and second, between ions in solution and between each ion and the globular protein. These latter contributions are missing from standard models. The potential of mean force, obtained from simulation, is fitted to an analytic equation. Using our analytic potential of mean force and Barker-Henderson perturbation theory, we obtain phase diagrams for lysozyme solutions that include stable and metastable fluid-fluid and solid-fluid phases when the electrolyte is 0.2 M NaSCN or NaI or NaCl. The nature of the electrolyte has a significant effect on the phase diagram.

Computer Simulation↗

Effect of single-point sequence alterations on the aggregation propensity of a model protein.

Sequences of contemporary proteins are believed to have evolved through a process that optimized their overall fitness, including their resistance to deleterious aggregation. Biotechnological processing may expose therapeutic proteins to conditions that are much more conducive to aggregation than those encountered in a cellular environment. An important task of protein engineering is to identify alternative sequences that would protect proteins when processed at high concentrations without altering their native structure associated with specific biological function. Our computational studies exploit parallel tempering simulations of coarse-grained model proteins to demonstrate that isolated amino acid residue substitutions can result in significant changes in the aggregation resistance of the protein in a crowded environment while retaining protein structure in isolation. A thermodynamic analysis of protein clusters subject to competing processes of folding and association shows that moderate mutations can produce effects similar to those caused by changes in system conditions, including temperature, concentration, and solvent composition, that affect the aggregation propensity. The range of conditions where a protein can resist aggregation can therefore be tuned by sequence alterations, although the protein generally may retain its generic ability for aggregation.

Amino Acid Substitution↗

Protein-folding landscapes in multichain systems.

Computational studies of proteins have significantly improved our understanding of protein folding. These studies are normally carried out by using chains in isolation. However, in many systems of practical interest, proteins fold in the presence of other molecules. To obtain insight into folding in such situations, we compare the thermodynamics of folding for a Miyazawa-Jernigan model 64-mer in isolation to results obtained in the presence of additional chains. The melting temperature falls as the chain concentration increases. In multichain systems, free-energy landscapes for folding show an increased preference for misfolded states. Misfolding is accompanied by an increase in interprotein interactions; however, near the folding temperature, the transition from folded chains to misfolded and associated chains is entropically driven. A majority of the most probable interprotein contacts are also native contacts, suggesting that native topology plays a role in early stages of aggregation.

Amino Acids↗

Interactions of lysozyme in guanidinium chloride solutions from static and dynamic light-scattering measurements.

The interactions of partially unfolded proteins provide insight into protein folding and protein aggregation. In this work, we studied partially unfolded hen egg lysozyme interactions in solutions containing up to 7 M guanidinium chloride (GdnHCl). The osmotic second virial coefficient (B(22)) of lysozyme was measured using static light scattering in GdnHCl aqueous solutions at 20 degrees C and pH 4.5. B(22) is positive in all solutions, indicating repulsive protein-protein interactions. At low GdnHCl concentrations, B(22) decreases with rising ionic strength: in the absence of GdnHCl, B(22) is 1.1 x 10(-3) mLmol/g(2), decreasing to 3.0 x 10(-5) mLmol/g(2) in the presence of 1 M GdnHCl. Lysozyme unfolds in solutions at GdnHCl concentrations higher than 3 M. Under such conditions, B(22) increases with ionic strength, reaching 8.0 x 10(-4) mLmol/g(2) at 6.5 M GdnHCl. Protein-protein hydrodynamic interactions were evaluated from concentration-dependent diffusivity measurements, obtained from dynamic light scattering. At moderate GdnHCl concentrations, lysozyme interparticle interactions are least repulsive and hydrodynamic interactions are least attractive. The lysozyme hydrodynamic radius was calculated from infinite-dilution diffusivity and did not change significantly during protein unfolding. Our results contribute toward better understanding of protein interactions of partially unfolded states in the presence of a denaturant; they may be helpful for the design of protein refolding processes that avoid protein aggregation.

Animals↗

Thermodynamics of folding and association of lattice-model proteins.

Closely related to the "protein folding problem" is the issue of protein misfolding and aggregation. Protein aggregation has been associated with the pathologies of nearly 20 human diseases and presents serious difficulties during the manufacture of pharmaceutical proteins. Computational studies of multiprotein systems have recently emerged as a powerful complement to experimental efforts aimed at understanding the mechanisms of protein aggregation. We describe the thermodynamics of systems containing two lattice-model 64-mers. A parallel tempering algorithm abates problems associated with glassy systems and the weighted histogram analysis method improves statistical quality. The presence of a second chain has a substantial effect on single-chain conformational preferences. The melting temperature is substantially reduced, and the increase in the population of unfolded states is correlated with an increase in interactions between chains. The transition from two native chains to a non-native aggregate is entropically favorable. Non-native aggregates receive approximately 25% of their stabilizing energy from intraprotein contacts not found in the lowest-energy structure. Contact maps show that for non-native dimers, nearly 50% of the most probable interprotein contacts involve pairs of residues that form native contacts, suggesting that a domain-swapping mechanism is involved in self-association.

Binding Sites↗

Steady-state diffusion of water through soft-contact-lens materials.

Water transport through soft contact lenses (SCL) is important for acceptable performance on the human eye. Chemical-potential gradient-driven diffusion rates of water through SCL materials are measured with an evaporation-cell technique. Water is evaporated from the bottom surface of a lens membrane by impinging air at controlled flow rate and humidity. The resulting weight loss of a water reservoir covering the top surface of the contact-lens material is recorded as a function of time. New results are reported for a conventional hydrogel material (SofLens One Day, hilafilcon A, water content at saturation w10 = 70 weight %) and a silicone hydrogel material (PureVision, balafilcon A, w10 = 36%), with and without surface oxygen plasma treatment. Also, previously reported data for a conventional 2-hydroxyethyl methacrylate (HEMA)-SCL (w10 = 38%) hydrogel are reexamined and compared with those for SofLens One Day and PureVision hydrogels. Measured steady-state water fluxes are largest for SofLens One Day, followed by PureVision and HEMA. In some cases, the measured steady-state water fluxes increase with rising relative air humidity. This increase, due to an apparent mass-transfer resistance at the surface (trapping skinning), is associated with formation of a glassy skin at the air/membrane interface when the relative humidity is below 55-75%. Steady-state water fluxes are interpreted through an extended Maxwell-Stefan diffusion model for a mixture of species starkly different in size. Thermodynamic nonideality is considered through Flory-Rehner polymer-solution theory. Shrinking/swelling is self-consistently modeled by conservation of the total polymer mass. Fitted Maxwell-Stefan diffusivities increase significantly with water concentration in the contact lens.

Biocompatible Materials↗

The competition between protein folding and aggregation: off-lattice minimalist model studies.

Protein aggregation has been associated with a number of human diseases, and is a serious problem in the manufacture of recombinant proteins. Of particular interest to the biotechnology industry is deleterious aggregation that occurs during the refolding of proteins from inclusion bodies. As a complement to experimental efforts, computer simulations of multi-chain systems have emerged as a powerful tool to investigate the competition between folding and aggregation. Here we report results from Langevin dynamics simulations of minimalist model proteins. Order parameters are developed to follow both folding and aggregation. By mapping natural units to real units, the simulations are shown to be carried out under experimentally relevant conditions. Data pertaining to the contacts formed during the association process show that multiple mechanisms for aggregation exist, but certain pathways are statistically preferred. Kinetic data show that there are multiple time scales for aggregation, although most association events take place at times much shorter than those required for folding. Last, we discuss results presented here as a basis for future work aimed at rational design of mutations to reduce aggregation propensity, as well as for development of small-molecular weight refolding enhancers.

Amino Acid Sequence↗

Effect of alcohols on aqueous lysozyme-lysozyme interactions from static light-scattering measurements.

Alcohols have been widely used as protein denaturants, precipitants and crystallization reagents. We have studied the effect of alcohols on aqueous hen-egg lysozyme self-interactions by measuring the osmotic second virial coefficient (B22) using static light scattering. Addition of alcohols increases B22, indicating stronger protein-protein repulsion or weaker attraction. For the monohydric alcohols used in this study (methanol, ethanol, 1-propanol, n-butanol, iso-butanol and trifluoroethanol), B22 for lysozyme reaches a common plateau at approximately 5% (v/v) alcohol, while glycerol increases B22 more than monohydric alcohols. For a 0.05 M NaCl hen-egg lysozyme solution at pH 7, B22 increases from 2.4 x 10(-4) to 4.7 x 10(-4) ml mol/g2 upon addition of monohydric alcohols and to 5.8 x 10(-4) ml mol/g2 upon addition of glycerol. We describe the alcohol effect using a simple model that supplements the DLVO theory with an additional alcohol-dependent term representing orientation-averaged hydrophobic interactions. In this model, the increased lysozyme repulsive forces in the presence of monohydric alcohols are interpreted in terms of adsorption of alcohol molecules on hydrophobic sites on the protein surface. This adsorption reduces attractive hydrophobic protein-protein interactions. A thicker lysozyme hydration layer in aqueous glycerol solution can explain the glycerol-increased lysozyme-lysozyme repulsion.

Alcohols↗

Solvent-amino acid interaction energies in three-dimensional-lattice Monte Carlo simulations of a model 27-mer protein: Folding thermodynamics and kinetics.

Amino acid residue-solvent interactions are required for lattice Monte Carlo simulations of model proteins in water. In this study, we propose an interaction-energy scale that is based on the interaction scale by Miyazawa and Jernigan. It permits systematic variation of the amino acid-solvent interactions by introducing a contrast parameter for the hydrophobicity, C(s), and a mean attraction parameter for the amino acids, omega. Changes in the interaction energies strongly affect many protein properties. We present an optimized energy parameter set for best representing realistic behavior typical for many proteins (fast folding and high cooperativity for single chains). Our optimal parameters feature a much weaker hydrophobicity contrast and mean attraction than does the original interaction scale. The proposed interaction scale is designed for calculating the behavior of proteins in bulk and at interfaces as a function of solvent characteristics, as well as protein size and sequence.

Algorithms↗

Generalizations for the potential of mean force between two isolated colloidal particles from Monte Carlo simulations.

A substantial amount of experimental and numerical evidence has shown that the Derjaguin-Landau-Verwey-Overbeek theory is not suitable for describing those colloidal solutions that contain multivalent counterions. Toward improved understanding of such solutions, the authors report Monte Carlo calculations wherein, following Rouzina and Bloomfield, they postulate that, in the absence of van der Waals forces, the overall force between two isolated charged colloidal particles in electrolyte solutions is determined by a dimensionless parameter Gamma=z(2)l(B)/a, which measures the electrostatic repulsion between counterions adsorbed on the macroion surface, where z = counterion valence, l(B)=Bjerrum length, and a = average separation between counterions on the macroion surface calculated as if the macroion were fully neutralized. The authors find, first, that the maximum repulsion between like-charged macroions occurs at Gamma approximately 0.5 and, second, that onset of attraction occurs at Gamma approximately 1.8, essentially independent of the valence and concentration of the surrounding electrolyte. These observations might provide new understanding of interactions between electrostatic double layers and perhaps offer explanations for some electrostatic phenomena related to interactions between DNA molecules or proteins.

Journal Article↗

Amyloid fibril formation by peptide LYS (11-36) in aqueous trifluoroethanol.

Peptide LYS (11-36), derived from the beta-sheet region of T4 lysozyme, forms an amyloid fibril in aqueous trifluoroethanol (TFE) at elevated temperature. The peptide has a moderate alpha-helix content in 20 and 50% (v/v) TFE solution; large quantities of fibrils were formed after incubation at 55 degrees C for 2 weeks as monitored by a thioflavin T fluorescence assay. No fibrils were observed when the peptide initially existed predominantly as a random coil or as a complete alpha helix. Our results suggest that a moderate amount of alpha helix and random coil present in the peptide initially facilitates the fibril-formation process, but a high alpha-helix content inhibits fibril formation. Transmission electron microscopy revealed several types of fibril morphologies at different TFE concentrations. The fibrils were highly twisted and consisted of interleaved protofilaments in 50% TFE, while smooth and flat ribbonlike fibrils were found in 20% TFE. In 50% TFE, the fibril growth rate of LYS (11-36) was found to depend strongly on peptide concentration and seeding but was insensitive to solution pH and ionic strength.

Amino Acid Sequence↗