Osmotic second virial cross‐coefficient measurements for binary combination of lysozyme,ovalbumin, and α‐amylase in salt solutions

Interactions measurement is a valuable tool to predict equilibrium phase separation of a desired protein in the presence of unwanted macromolecules. In this study, cross‐interactions were measured as the osmotic second virial cross‐coefficients (B₂₃) for the three binary protein systems involving lysozyme, ovalbumin, and α‐amylase in salt solutions (sodium chloride and ammonium sulfate). They were correlated with solubility for the binary protein mixtures. The cross‐interaction behavior at different salt concentrations was interpreted by either electrostatic or hydrophobic interaction forces. At low salt concentrations, the protein surface charge dominates cross‐interaction behavior as a function of pH. With added ovalbumin, the lysozyme solubility decreased linearly at low salt concentration in sodium chloride and increased at high salt concentration in ammonium sulfate. The B₂₃ value was found to be proportional to the slope of the lysozyme solubility against ovalbumin concentration and the correlation was explained by preferential interaction theory. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 29:1203–1211, 2013     

Protein Structure Prediction with a Combined Solvation Free Energy-Molecular Mechanics Force Field

Abstract

Models of protein structure are frequently used to determine the physical characteristics of a protein when the crystal structure is not available. We developed a procedure to optimize such models, by use of a combined solvation free energy and molecular mechanics force field. Appropriately chosen atomic solvation parameters were defined using the criterion that the resulting protein model should deviate least from the crystal structure upon a forty picosecond molecular dynamics simulation carried out using the combined force field. Several tests were performed to refine the set of atomic solvation parameters which best complement the molecular mechanics forces. Four sets of parameters from the literature were tested and an empirically optimized set is proposed. The parameters are defined on a well characterized small molecule (alanyl dipeptide) and on the highly refined crystal structure of rat trypsin, and then tested on a second highly refined crystal structure of α-lytic protease. The new set of atomic solvation parameters provides a significant improvement over molecular mechanics alone in energy minimization of protein structures. This combined force field also has advantages over the use of explicit solvent as it is possible to take solvent effects into account during energetic conformational searching when modeling a homologous protein structure from a known crystal structure.     

Protein oligomerization through domain swapping: role of inter-molecular interactions and protein concentration

Domain swapping has been shown to be an important mechanism controlling multiprotein assembly and has been suggested recently as a possible mechanism underlying protein aggregation. Understanding oligomerization via domain swapping is therefore of theoretical and practical importance. By using a symmetrized structure-based (Gō) model, we demonstrate that in the free-energy landscape of domain swapping, a large free-energy barrier separates monomeric and domain-swapped dimeric configurations. We investigate the effect of finite monomer concentration, by implementing a new semi-analytical method, which involves computing the second virial coefficient, a thermodynamic indicator of inter-molecular interactions. This method, together with the symmetrized structure-based (Gō) model, minimizes the need for expensive many-protein simulations, providing a convenient framework to investigate concentration effect. Finally, we perform direct simulations of domain-swapped trimer formation, showing that this modeling approach can be used for higher-order oligomers.     

A model binding site for testing scoring functions in molecular docking

Prediction of interaction energies between ligands and their receptors remains a major challenge for structure-based inhibitor discovery. Much effort has been devoted to developing scoring schemes that can successfully rank the affinities of a diverse set of possible ligands to a binding site for which the structure is known. To test these scoring functions, well-characterized experimental systems can be very useful. Here, mutation-created binding sites in T4 lysozyme were used to investigate how the quality of atomic charges and solvation energies affects molecular docking. Atomic charges and solvation energies were calculated for 172,118 molecules in the Available Chemicals Directory using a semi-empirical quantum mechanical approach by the program AMSOL. The database was first screened against the apolar cavity site created by the mutation Leu99Ala (L99A). Compared to the electronegativity-based charges that are widely used, the new charges and desolvation energies improved ranking of known apolar ligands, and better distinguished them from more polar isosteres that are not observed to bind. To investigate whether the new charges had predictive value, the non-polar residue Met102, which forms part of the binding site, was changed to the polar residue glutamine. The structure of the resulting Leu99Ala and Met102Gln double mutant of T4 lysozyme (L99A/M102Q) was determined and the docking calculation was repeated for the new site. Seven representative polar molecules that preferentially docked to the polar versus the apolar binding site were tested experimentally. All seven bind to the polar cavity (L99A/M102Q) but do not detectably bind to the apolar cavity (L99A). Five ligand-bound structures of L99A/M102Q were determined by X-ray crystallography. Docking predictions corresponded to the crystallographic results to within 0.4A RMSD. Improved treatment of partial atomic charges and desolvation energies in database docking appears feasible and leads to better distinction of true ligands. Simple model binding sites, such as L99A and its more polar variants, may find broad use in the development and testing of docking algorithms.     

Correlation of second virial coefficient with solubility for proteins in salt solutions

In this work, osmotic second virial coefficients (B(22)) were determined and correlated with the measured solubilities for the proteins, α-amylase, ovalbumin, and lysozyme. The B(22) values and solubilities were determined in similar solution conditions using two salts, sodium chloride and ammonium sulfate in an acidic pH range. An overall decrease in the solubility of the proteins (salting out) was observed at high concentrations of ammonium sulfate and sodium chloride solutions. However, for α-amylase, salting-in behavior was also observed in low concentration sodium chloride solutions. In ammonium sulfate solutions, the B(22) are small and close to zero below 2.4 M. As the ammonium sulfate concentrations were further increased, B(22) values decreased for all systems studied. The effect of sodium chloride on B(22) varies with concentration, solution pH, and the type of protein studied. Theoretical models show a reasonable fit to the experimental derived data of B(22) and solubility. B(22) is also directly proportional to the logarithm of the solubility values for individual proteins in salt solutions, so the log-linear empirical models developed in this work can also be used to rapidly predict solubility and B(22) values for given protein-salt systems.     

Simple theoretical models for biochemical systems,with applications to DNA

A set oflogically connected models, to study chemical systems ofbiological interest, is presented. The sequence in the set is dictated by a progressive reduction of details with a corresponding enlargement of the field of application. The exposition starts with models suitable for interactions among a finite number of molecules, passes then to models considering also solvent effects and ends with models specialized for DNA containing systems.     

Temperature-induced complementarity as a mechanism for biomolecular assembly

Recent advances in the measurement and theory of “hydration” interactions between biomolecules provide a basis on which to formulate mechanisms of biomolecular recognition. In this paper we have developed a mathematical formalism for analyzing specificity encoded in dynamic distributions of surface polar groups, a formalism that incorporates newly recognized properties of directly measured “hydration” forces. As expected, attraction between surfaces requires complementary patterns of surface polar groups. In contrast to usual expectations, thermal motion can create these complementary surface configurations. We have demonstrated that assembly can occur with an increase in conformational entropy of polar residues. Elevated temperature then facilitates recognition rather than hinders it. This mechanism might underlie some temperature-favored assembly reactions common in biological systems that are usually associated with the “hydrophobic effect” only. © 1994 Wiley-Liss, Inc.     

Interaction of Cardiotoxin A5 with Membrane: Role of Conformational Heterogeneity and Hydrophobic Properties

The hypothesis that local conformational differences of the snake venom cardiotoxins (cytotoxins, CT) may play a significant role in their interaction with membrane was tested by molecular modeling of the behavior of the CT A5 from the venom of Naja atra in water and at the water–membrane interface. Two models of the CT A5 spatial structure are known: the first was obtained by X-ray analysis and the second, by NMR studies in solution. A molecular dynamics (MD) analysis demonstrated that loop II of the toxin has a fixed -like shape in water, which does not depend on its initial structure. An interaction of the experimentally derived (X-ray and NMR) conformations and MD simulated conformations of CT A5 with the lipid bilayer was studied by the Monte Carlo method using the previously developed model of the implicit membrane. It is found that: (1) unlike the previously studied CT2 from the venom of cobra Naja oxiana, CT A5 has only loops I and II bound to the membrane with the involvement of a lesser number of hydrophobic residues. (2) A long hydrophobic area is formed on the surface of CT A5 due to the -like shape of loop II and the arrangement of loop I in proximity to loop II. This hydrophobic area favors the toxin embedding into the lipid bilayer. (3) The toxin retains its conformation upon interaction with the membrane. (4) The CT A5 molecule has close values of the potential energy in the membrane and in aqueous environment, which suggests dynamic character of the binding. The results of the molecular modeling indicate a definite configuration of loops I and II and, consequently, a specific character of distribution of polar and apolar properties on the toxin surface, which turns out to be the most energetically favorable.     

Molecular interactions between a plant virus movement protein and RNA: force spectroscopy investigation

RNA-protein interactions are fundamental for different aspects of molecular biology such as gene expression, assembly of biomolecular complexes or macromolecular transport. The 3a movement protein (MP) of a plant virus, Cucumber mosaic virus (CMV), forms ribonucleoprotein (RNP) complexes with viral RNA, capable of trafficking from cell-to-cell throughout the infected plant only in the presence of the CMV capsid protein (CP). However, deletion of the C-terminal 33 amino acid residues of the CMV MP (in the mutant designated 3aDeltaC33 MP) resulted in CP-independent cell-to-cell movement. The biological differences in the behaviour of CMV wild type (wt) 3a MP and 3aDeltaC33 MP could have been a consequence of differences in the RNA-binding properties of the two MPs detected previously using biochemical assays on ensembles of molecules. To investigate the physical mechanisms of MP-RNA interactions at a single molecule level, we applied atomic force microscopy to measure for the first time unbinding forces between these individual binding partners. Minimal unbinding forces determined for individual interaction of the CMV RNA molecule with the CMV wt or truncated MPs were estimated to be approximately 45 pN and approximately 90 pN, respectively, suggesting that the distinct differences in the strength of MP-RNA interactions for the wt MP and truncated MP are attributable to the molecular binding mechanism. We also demonstrated that molecules of both CMV 3a MP and 3aDeltaC33 MP were capable of self-interaction with minimal unbinding forces of approximately 50 pN and approximately 70 pN, respectively, providing a physical basis for the cooperative mechanism of the RNA binding. The significance of intermolecular force measurements for understanding the structural and functional aspects of viral RNP formation and trafficking is discussed.