Evaluating adsorbed-phase activity coefficient models using a 2D-lattice model

Despite the wide use of the real adsorbed solution theory to predict multicomponent adsorption equilibrium, the models used for the adsorbed phase activity coefficients are usually borrowed from the gas–liquid phase equilibria. In this work, the accuracy of the Wilson and NRTL models for evaluating adsorbed phase activity coefficients is tested using a 2D-lattice model. An accurate model for adsorbed-phase activity coefficients should have no problem in fitting adsorption data obtained using this simple lattice model. The results, however, show that the commonly used Wilson and NRTL models cannot describe the adsorbed phase activity coefficients for slightly non-ideal to strong non-ideal mixtures. Therefore, until new models for adsorbed phase activity coefficients are developed, we should use existing models for liquids with care. In the second part of this work, the use of Monte Carlo simulations on a segregated 2D-lattice model, for predicting adsorption of mixtures is investigated. The segregated model assumes that the competition for adsorption occurs at isolated adsorption sites, and that the molecules from each adsorption site interact with the bulk phase independently. Two binary mixtures in two adsorbent materials were used as case studies for testing the predictions of the segregated 2D-lattice model: the binary system CO₂–N₂ in the hypothetical pure silica zeolite PCOD8200029, with isolated adsorption sites and normal preference for adsorption, and the binary system CO₂–C₃H₈ in pure silica mordenite (MOR), with isolated adsorption sites and inverse site preference. The segregated 2D-lattice model provides accurate predictions for the system CO₂–N₂ in PCOD8200029 but fails in predicting the adsorption behaviour of CO₂–C₃H₈ in pure silica MOR. The predictions of the segregated ideal adsorbed solution theory model are superior to those of the 2D-lattice model.     

Domains in folding of model proteins.

By means of Monte Carlo simulation, we investigated the equilibrium between folded and unfolded states of lattice model proteins. The amino acid sequences were designed to have pronounced energy minimum target conformations of different length and shape. For short fully compact (36-mer) proteins, the all-or-none transition from the unfolded state to the native state was observed. This was not always the case for longer proteins. Among 12 designed sequences with the native structure of a fully compact 48-mer, a simple all-or-none transition was observed in only three cases. For the other nine sequences, three states of behavior-the native, denatured, and intermediate states-were found. The contiguous part of the native structure (domain) was conserved in the intermediate state, whereas the remaining part was completely unfolded and structureless. These parts melted separately from each other.     

Monte Carlo Simulation of Binary Gas Adsorption in Zeolite Cavities

Abstract

Grand canonical Monte Carlo simulations have been performed for binary adsorption of Lennard-Jones molecules with point multipole moments in zeolite cavities of type X. Fluid-solid electrostatic interactions were taken into account. Phase diagrams and total coverage were calculated for three binaries and compared with experimental measurements. MC simulations gave good agreement with experiment for two mixtures (C₂H₄-CO₂ and CO₂-CH₄) but there were discrepancies between simulation and experiment for the system i-C₄H₁₀-C₂H₄. The dependence of excess Gibbs free energy on the composition and pressure was studied. Negative deviations from ideality are due to energetic heterogeneity and size effects. Unlike liquid-vapor equilibrium, deviations from the Lorentz-Berthelot mixing rules for the adsorbates have little effect upon the phase behavior. Density distributions show that the components compete for the high energy sites inside the cavity; depending on its relative strength of adsorption, one component may be excluded from such positions (CH₄ in CO₂-CH₄), or the two species may share sites inside the cavity (C₂H₄-CO₂).     

Crystallographic analysis of a full-length streptavidin with its C-terminal polypeptide bound in the biotin binding site

The structure of a full-length streptavidin has been determined at 1.7 A resolution and shows that the 20 residue extension at the C terminus forms a well-ordered polypeptide loop on the surface of the tetramer. Residues 150-153 of the extension are bound to the ligand-binding site, possibly competing with exogenous ligands. The binding mode of these residues is compared with that of biotin and peptidic ligands. The observed structure helps to rationalize the observations that full-length mature streptavidin binds biotinylated macromolecules with reduced affinity.     

A fast and simple method for probing the interaction of peptides and proteins with lipids and membrane-mimetics using GB1 fusion proteins and NMR spectroscopy

The expression of peptides and proteins as fusions to the B1 domain of streptococcal protein G (GB1) is very popular since GB1 often improves the solubility of the target protein and because the first purification step using IgG affinity chromatography is simple and efficient. However, the following protease digest is not always complete or can result in a digest of the target protein. In addition, a further purification step such as RP-HPLC has to be used to get rid of the GB1 tag and undigested fusion protein. Because the protease digest and the following purification step are not only time-consuming but generally also expensive, we tested if GB1 fusion proteins can directly be used for NMR interaction studies using lipids or membrane-mimetics. Based on NMR binding studies using only the GB1 part, this fusion tag does not significantly interact with different membrane-mimetics such as micelles, bicelles, or liposomes. Thus spectral changes observed using GB1-fusion proteins indicate lipid- and membrane interactions of the target protein. The method was initially established to probe membrane interactions of a large number of mutants of the FATC domain of the ser/thr kinase TOR. To demonstrate the usefulness of the approach, we show NMR binding data for the wild type protein and a leucine to alanine mutant.     

Specificity of helix packing in transmembrane dimer of the cell death factor BNIP3: a molecular modeling study

BNIP3 is a mitochondrial 19-kDa proapoptotic protein, a member of the Bcl-2 family. It has a single COOH-terminal transmembrane (TM) alpha-helical domain, which is required for membrane targeting, proapoptotic activity, hetero- and homo-dimerization in membrane. The role and the molecular details of association of TM helices of BNIP3 are yet to be established. Here, we present a molecular modeling study of helix interactions in its membrane domain. The approach combines Monte Carlo conformational search in an implicit hydrophobic slab followed by molecular dynamics simulations in a hydrated full-atom lipid bilayer. The former technique was used for exhaustive sampling of the peptides' conformational space and for generation of putative "native-like" structures of the dimer. The latter ones were taken as realistic starting points to assess stability and dynamic behavior of the complex in explicit lipid-water surrounding. As a result, several groups of tightly packed right-handed structures of the dimer were proposed. They have almost similar helix-helix interface, which includes the motif A(176)xxxG(180)xxxG(184) and agrees well with previous mutagenesis data and preliminary NMR analysis. Molecular dynamics simulations of these structures reveal perfect adaptation of most of them to heterogeneous membrane environment. A remarkable feature of the predicted dimeric structures is the occurrence of a cluster of H-bonded histidine 173 and serines 168 and 172 on the helix interface, near the N-terminus. Because of specific polar interactions between the monomers, this part of the dimer has no such dense packing as the C-terminal one, thus allowing penetration of water from the extramembrane side into the membrane interior. We propose that the ionization state of His(173) can mediate structural and dynamic properties of the dimer. This, in turn, may be related to pH-dependent proapoptotic activity of BNIP3, which is triggering on by acidosis appearing under hypoxic conditions.     

Brownian dynamics simulations of glycolytic enzyme subsets with F-actin

Previous Brownian dynamics (BD) simulations identified specific basic residues on fructose-1,6-bisphophate aldolase (aldolase) (I. V. Ouporov et al., Biophysical Journal, 1999, Vol. 76, pp. 17-27) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (I. V. Ouporov et al., Journal of Molecular Recognition, 2001, Vol. 14, pp. 29-41) involved in binding F-actin, and suggested that the quaternary structure of the enzymes may be important. Herein, BD simulations of F-actin binding by enzyme dimers or peptides matching particular sequences of the enzyme and the intact enzyme triose phosphate isomerase (TIM) are compared. BD confirms the experimental observation that TIM has little affinity for F-actin. For aldolase, the critical residues identified by BD are found in surface grooves, formed by subunits A/D and B/C, where they face like residues of the neighboring subunit enhancing their electrostatic potentials. BD simulations between F-actin and aldolase A/D dimers give results similar to the native tetramer. Aldolase A/B dimers form complexes involving residues that are buried in the native structure and are energetically weaker; these results support the importance of quaternary structure for aldolase. GAPDH, however, placed the critical residues on the corners of the tetramer so there is no enhancement of the electrostatic potential between the subunits. Simulations using GAPDH dimers composed of either S/H or G/H subunits show reduced binding energetics compared to the tetramer, but for both dimers, the sets of residues involved in binding are similar to those found for the native tetramer. BD simulations using either aldolase or GAPDH peptides that bind F-actin experimentally show complex formation. The GAPDH peptide bound to the same F-actin domain as did the intact tetramer; however, unlike the tetramer, the aldolase peptide lacked specificity for binding a single F-actin domain.     

Brownian dynamics of interactions between aldolase mutants and F-actin

Previous Brownian dynamics (BD) simulations (Ouporov IG, Knull HR and Thomasson KA 1999. Biophys. J. 76: 17-27) of complex formation between rabbit aldolase and F-actin have identified three lysine residues (K288, K293 and K341) on aldolase and acidic residues (DEDE) at the N-terminus of actin as important to binding. BD simulations of computer models of aldolase mutants with any of these lysine residues replaced by alanine show reduced binding energy; the greatest effect of a single substitution is for K341A, and replacement of all three lysines greatly reduces binding. BD simulations of wild-type rabbit aldolase vs altered F-actin show that binding is decreased if any one of the four N-terminal acidic residues is replaced by alanine and binding is greatly reduced if three or more of the N-terminal acidic residues are replaced; none of the four actin residues appear more critical for binding than the others.     

Assembly of protein structure from sparse experimental data: An efficient Monte Carlo model

A new, efficient method for the assembly of protein tertiary structure from known, loosely encoded secondary structure restraints and sparse information about exact side chain contacts is proposed and evaluated. The method is based on a new, very simple method for the reduced modeling of protein structure and dynamics, where the protein is described as a lattice chain connecting side chain centers of mass rather than Cαs. The model has implicit built-in multibody correlations that simulate short- and long-range packing preferences, hydrogen bonding cooperativity and a mean force potential describing hydrophobic interactions. Due to the simplicity of the protein representation and definition of the model force field, the Monte Carlo algorithm is at least an order of magnitude faster than previously published Monte Carlo algorithms for structure assembly. In contrast to existing algorithms, the new method requires a smaller number of tertiary restraints for successful fold assembly; on average, one for every seven residues as compared to one for every four residues. For example, for smaller proteins such as the B domain of protein G, the resulting structures have a coordinate root mean square deviation (cRMSD), which is about 3 Å from the experimental structure; for myoglobin, structures whose backbone cRMSD is 4.3 Å are produced, and for a 247-residue TIM barrel, the cRMSD of the resulting folds is about 6 Å. As would be expected, increasing the number of tertiary restraints improves the accuracy of the assembled structures. The reliability and robustness of the new method should enable its routine application in model building protocols based on various (very sparse) experimentally derived structural restraints. Proteins 32:475–494, 1998. © 1998 Wiley-Liss, Inc.     

A comparison of SCOP and CATH with respect to domain-domain interactions

The analysis and prediction of protein-protein interaction sites from structural data are restricted by the limited availability of structural complexes that represent the complete protein-protein interaction space. The domain classification schemes CATH and SCOP are normally used independently in the analysis and prediction of protein domain-domain interactions. In this article, the effect of different domain classification schemes on the number and type of domain-domain interactions observed in structural data is systematically evaluated for the SCOP and CATH hierarchies. Although there is a large overlap in domain assignments between SCOP and CATH, 23.6% of CATH interfaces had no SCOP equivalent and 37.3% of SCOP interfaces had no CATH equivalent in a nonredundant set. Therefore, combining both classifications gives an increase of between 23.6 and 37.3% in domain-domain interfaces. It is suggested that if possible, both domain classification schemes should be used together, but if only one is selected, SCOP provides better coverage than CATH. Employing both SCOP and CATH reduces the false negative rate of predictive methods, which employ homology matching to structural data to predict protein-protein interaction by an estimated 6.5%.     

Lattice simulations of cotranslational folding of single domain proteins

We use lattice protein models and Monte Carlo simulations to study cotranslational folding of small single domain proteins. We show that the assembly of native structure begins during late extrusion stages, but final formation of native state occurs during de novo folding, when all residues are extruded. There are three main results in our study. First, for the sequences displaying two-state refolding mechanism de novo cotranslational folding pathway differs from that sampled in in vitro refolding. The change in folding pathways is due to partial assembly of native interactions during extrusion that results in different starting conditions for in vitro refolding and for de novo cotranslational folding. For small single domain proteins cotranslational folding is slower than in vitro refolding, but is generally fast enough to be completed before the release from a ribosome. Second, we found that until final stages of biosynthesis cotranslational folding is essentially equilibrium. This observation is explained by low stability of structured states for partially extruded chains. Finally, our data suggest that the proteins, which refold in vitro slowly via intermediates, complete their de novo folding after the release from a ribosome. Comparison of our lattice cotranslational simulations with recent experimental and computational studies is discussed.     

The elixir phase of chain molecules

A phase of matter is a familiar notion for inanimate physical matter. The nature of a phase of matter transcends the microscopic material properties. For example, materials in the liquid phase have certain common properties independent of the chemistry of the constituents: liquids take the shape of the container; they flow; and they can be poured—alcohol, oil, and water as well as a Lennard-Jones computer model exhibit similar behavior when poised in the liquid phase. Here, we identify a hitherto unstudied “phase” of matter, the elixir phase, in a simple model of a polymeric chain whose backbone has the correct local cylindrical symmetry induced by the tangent to the chain. The elixir phase appears on breaking the cylindrical symmetry by adding side spheres along the negative normal direction, as in proteins. This phase, nestled between other phases, has multiple ground states made up of building blocks of helices and almost planar sheets akin to protein native folds. We discuss the similarities of this “phase” of a finite size system to the liquid crystal and spin glass phases. Our findings are relevant for understanding proteins; the creation of novel bioinspired nanomachines; and also may have implications for life elsewhere in the cosmos.     

Influence of substrate pattern on the adsorption of HP lattice proteins

With the highly simplified hydrophobic-polar model representation of a protein, we can study essential qualitative physics without an unnecessarily large computational overhead. Using Wang-Landau sampling in conjunction with a set of efficient Monte Carlo trial moves, we studied the adsorption of short HP lattice proteins on various simple patterned substrates and in particular for checkered patterned surfaces. A set of single-site mutated HP proteins is used to investigate the role of hydrophobicity of a protein chain and surface pattern for substrates with various pattern cell sizes relative to the protein’s native configuration. For most cases, we found that the adsorption transition occurs at a lower temperature, while the hydrophobic core formation is less affected. The flattening procedure after the HP protein is adsorbed is more sensitive to the change in surface patterns and single-site mutations. These observations stay valid for both strongly and weakly attractive surfaces.     

Integration of cell-free protein coexpression with an enzyme-linked immunosorbent assay enables rapid analysis of protein-protein interactions directly from DNA

Assays that integrate detection of binding with cell-free protein expression directly from DNA can dramatically increase the pace at which protein-protein interactions (PPIs) can be analyzed by mutagenesis. In this study, we present a method that combines in vitro protein production with an enzyme-linked immunosorbent assay (ELISA) to measure PPIs. This method uses readily available commodity instrumentation and generic antibody-affinity tag interactions. It is straightforward and rapid to execute, enabling many interactions to be assessed in parallel. In traditional ELISAs, reporter complexes are assembled stepwise with one layer at a time. In the method presented here, all the members of the reporter complex are present and assembled together. The signal strength is dependent on all the intercomponent interaction affinities and concentrations. Although this assay is straightforward to execute, establishing proper conditions and analysis of the results require a thorough understanding of the processes that determine the signal strength. The formation of the fully assembled reporter sandwich can be modeled as a competition between Langmuir adsorption isotherms for the immobilized components and binding equilibria of the solution components. We have shown that modeling this process provides semiquantitative understanding of the effects of affinity and concentration and can guide strategies for the development of experimental protocols. We tested the method experimentally using the interaction between a synthetic ankyrin repeat protein (Off7) and maltose-binding protein. Measurements obtained for a collection of alanine mutations in the interface between these two proteins demonstrate that a range of affinities can be analyzed.     

Protein unfolding at interfaces: slow dynamics of alpha-helix to beta-sheet transition

A two-phase sequential dynamic change in the secondary structure of hen egg lysozyme (Lys) adsorbed on solid substrates was observed. The first phase involved fast conversion of alpha-helix to random/turns (within the first minute or at very low coverage or high substrate wettability) with no perceptible change in beta-sheet content. The second phase (1-1200 min), however, involved a relatively slow conversion from alpha-helix to beta-sheet without a noticeable change in random/turns. An important finding of this work is that the concentration of lysozyme in the adsorbed state has a substantial effect on the fractional content of secondary structures. Attenuated total reflection Fourier transform infrared (ATR/FTIR) spectroscopy, along with a newly-developed optimization algorithm for predicting the content of secondary structure motifs, was used to correlate the secondary structure and the amount of adsorbed lysozyme with the surface wettability of six different flat nanoporous substrates. Although three independent variables, surface wettability, solution concentration and time for adsorption, were used to follow the fractional structural changes of lysozyme, the results were all normalized onto a single plot with the amount adsorbed as the universal independent variable. Consequently, lateral interactions among proteins likely drive the transition process. Direct intermolecular force adhesion measurements between lysozyme and different functionalized self-assembled alkanethiol monolayers confirm that hydrophobic surfaces interact strongly with proteins. The lysozyme-unfolding pathway during early adsorption appears to be similar to that predicted by published molecular modeling results.     

Effect of temperature on adsorption of mixtures in porous materials

The effect of temperature on the adsorption of a simple mixture (Ar/Kr) in disordered porous materials is investigated by means of molecular simulation. In the larger mesopores of porous silica glasses, capillary condensation occurs upon decreasing the temperature. At temperatures above the capillary condensation temperature, Kr is preferentially adsorbed at the pore surface and Ar adsorption occurs in regions of low Kr density. For temperatures below the capillary condensation temperature, Ar density surprisingly increases as temperature increases, the behaviour that is consistent with an over-solubility effect. In contrast, in the disordered sub-nanoporous carbon, filling of the pores occurs in a reversible and continuous way upon decreasing the temperature, owing to the small size and amorphous shape of the pores. These results show that the crossover between capillary condensation and continuous reversible filling observed for pure fluids in pores also exists for mixtures. We also show that the Kr selectivity exhibits a minimum in the disordered porous silica that is located at the capillary condensation temperature. In contrast, in the disordered porous carbon where no capillary condensation occurs, the selectivity decreases monotonically with increasing the temperature. These results shed light on low-temperature adsorption of mixtures confined in porous materials and provide a guide to design efficient phase separation processes.     

The toxin-antitoxin proteins relBE2Spn of Streptococcus pneumoniae: characterization and association to their DNA target

The chromosome of the pathogenic Gram-positive bacterium Streptococcus pneumoniae contains between six to 10 operons encoding toxin-antitoxin systems (TAS). TAS are widespread and redundant in bacteria and archaea and their role, albeit still obscure, may be related to important aspects of bacteria lifestyle like response to stress. One of the most abundant TAS is the relBE family, being present in the chromosome of many bacteria and archaea. Because of the high rates of morbility and mortality caused by S. pneumoniae, it has been interesting to gain knowledge on the pneumococcal TAS, among them the RelBE2Spn proteins. Here, we have analyzed the DNA binding capacity of the RelB2Spn antitoxin and the RelB2Spn-RelE2Spn proteins by band-shift assays. Thus, a DNA region encompassing the operator region of the proteins was identified. In addition, we have used analytical ultracentrifugation and native mass spectrometry to measure the oligomerization state of the antitoxin alone and the RelBE2Spn complex in solution bound or unbound to its DNA substrate. Using native mass spectrometry allowed us to unambiguously determine the stoichiometry of the RelB2Spn and of the RelBE2Spn complex alone or associated to its DNA target.     

Benchmarking of dimeric threading and structure refinement

The understanding of protein-protein interactions is a major goal in the postgenomic era. The prediction of interaction from sequence and the subsequent generation of full-length dimeric models is therefore of great interest especially because the number of structurally characterized protein-protein complexes is sparse. A quality assessment of a benchmark comprised of 170 weakly homologous dimeric target-template pairs is presented. They are predicted in a two-step method, similar to the previously described MULTIPROSPECTOR algorithm: each target sequence is assigned to a monomeric template structure by threading; then, those templates that belong to the same physically interacting dimer template are selected. Additionally we use structural alignments as the "gold standard" to assess the percentage of correctly assigned monomer and dimer templates and to evaluate the threading results with a focus on the quality of the alignments in the interfacial region. This work aims to give a quantitative picture of the quality of dimeric threading. Except for one, all monomer templates are identified correctly, but approximately 40% of the dimer templates are still problematic or incorrect. Preliminary results for three full-length dimeric models generated with the TASSER method show on average a significant improvement of the final model over the initial template.     

Intermolecular versus intramolecular interactions of the vinculin binding site 33 of talin

The cytoskeletal proteins talin and vinculin are localized at cell‐matrix junctions and are key regulators of cell signaling, adhesion, and migration. Talin couples integrins via its FERM domain to F‐actin and is an important regulator of integrin activation and clustering. The 220 kDa talin rod domain comprises several four‐ and five‐helix bundles that harbor amphipathic α‐helical vinculin binding sites (VBSs). In its inactive state, the hydrophobic VBS residues involved in binding to vinculin are buried within these helix bundles, and the mechanical force emanating from bound integrin receptors is thought necessary for their release and binding to vinculin. The crystal structure of a four‐helix bundle of talin that harbors one of these VBSs, coined VBS33, was recently determined. Here we report the crystal structure of VBS33 in complex with vinculin at 2 Å resolution. Notably, comparison of the apo and vinculin bound structures shows that intermolecular interactions of the VBS33 α‐helix with vinculin are more extensive than the intramolecular interactions of the VBS33 within the talin four‐helix bundle.