Structure and stability of a model three-helix-bundle protein on tailored surfaces

Coarse-grained simulation of PEGylated and tethered protein devices at all experimentally accessible surface residues on β-lactamase for stability analysis and comparison

PEGylated and surface-tethered proteins are used in a variety of biotechnological applications, but traditional methods offer little control over the placement of the functionalization sites on the protein. Fortunately, recent experimental methods functionalize the protein at any location on the amino acid sequence, so the question becomes one of selecting the site that will result in the best protein function. This work shows how molecular simulation can be used to screen potential attachment sites for surface tethering or PEGylation. Previous simulation work has shown promise in this regard for a model protein, but these studies are limited to screening only a few of the surface-accessible sites or only considered surface tethering or PEGylation separately rather than their combined effects. This work is done to overcome these limitations by screening all surface-accessible functionalization sites on a protein of industrial and therapeutic importance (TEM-1) and to evaluate the effects of tethering and PEGylation simultaneously in an effort to create a more accurate screen. The results show that functionalization site effectiveness appears to be a function of super-secondary and tertiary structures rather than the primary structure, as is often currently assumed. Moreover, sites in the middle of secondary structure elements, and not only those in loops regions, are shown to be good options for functionalization-a fact not appreciated in current practice. Taken as a whole, the results show how rigorous molecular simulation can be done to identify candidate amino acids for functionalization on a protein to facilitate the rational design of protein devices.

The effects of antigen size, binding site valency, and flexibility on fab-antigen binding near solid surfaces

Next generation antibody microarray devices have the potential to outperform current molecular detection methods and realize new applications in medicine, scientific research, and national defense. However, antibody microarrays, or arrays of antibody fragments ("fabs"), continue to evade mainstream use in part due to persistent reliability problems despite improvements to substrate design and protein immobilization strategies. Other factors could be disrupting microarray performance, including effects resulting from antigen characteristics. Target molecules embody a wide range of sizes, shapes, number of epitopes, epitope accessibility, and other physical and chemical properties. As a result, it may not be ideal for microarray designs to utilize the same substrate or immobilization strategy for all of the capture molecules. This study investigates how three antigen properties, such as size, binding site valency, and molecular flexibility, affect fab binding. The work uses an advanced, experimentally validated, coarse-grain model and umbrella sampling to calculate the free energy of ligand binding and how this energy landscape is different on the surface compared to in the bulk. The results confirm that large antigens interact differently with immobilized fabs compared to smaller antigens. Analysis of the results shows that despite these differences, tethering fabs in an upright orientation on hydrophilic surfaces is the best configuration for antibody microarrays.

The effects of tether placement on antibody stability on surfaces

Despite their potential benefits, antibody microarrays have fallen short of performing reliably and have not found widespread use outside of the research setting. Experimental techniques have been unable to determine what is occurring on the surface of an atomic level, so molecular simulation has emerged as the primary method of investigating protein/surface interactions. Simulations of small proteins have indicated that the stability of the protein is a function of the residue on the protein where a tether is placed. The purpose of this research is to see whether these findings also apply to antibodies, with their greater size and complexity. To determine this, 24 tethering locations were selected on the antibody Protein Data Bank (PDB) ID: 1IGT. Replica exchange simulations were run on two different surfaces, one hydrophobic and one hydrophilic, to determine the degree to which these tethering sites stabilize or destabilize the antibody. Results showed that antibodies tethered to hydrophobic surfaces were in general less stable than antibodies tethered to hydrophilic surfaces. Moreover, the stability of the antibody was a function of the tether location on hydrophobic surfaces but not hydrophilic surfaces.

Modeling and simulation of protein-surface interactions: achievements and challenges

Understanding protein-inorganic surface interactions is central to the rational design of new tools in biomaterial sciences, nanobiotechnology and nanomedicine. Although a significant amount of experimental research on protein adsorption onto solid substrates has been reported, many aspects of the recognition and interaction mechanisms of biomolecules and inorganic surfaces are still unclear. Theoretical modeling and simulations provide complementary approaches for experimental studies, and they have been applied for exploring protein-surface binding mechanisms, the determinants of binding specificity towards different surfaces, as well as the thermodynamics and kinetics of adsorption. Although the general computational approaches employed to study the dynamics of proteins and materials are similar, the models and force-fields (FFs) used for describing the physical properties and interactions of material surfaces and biological molecules differ. In particular, FF and water models designed for use in biomolecular simulations are often not directly transferable to surface simulations and vice versa. The adsorption events span a wide range of time- and length-scales that vary from nanoseconds to days, and from nanometers to micrometers, respectively, rendering the use of multi-scale approaches unavoidable. Further, changes in the atomic structure of material surfaces that can lead to surface reconstruction, and in the structure of proteins that can result in complete denaturation of the adsorbed molecules, can create many intermediate structural and energetic states that complicate sampling. In this review, we address the challenges posed to theoretical and computational methods in achieving accurate descriptions of the physical, chemical and mechanical properties of protein-surface systems. In this context, we discuss the applicability of different modeling and simulation techniques ranging from quantum mechanics through all-atom molecular mechanics to coarse-grained approaches. We examine uses of different sampling methods, as well as free energy calculations. Furthermore, we review computational studies of protein-surface interactions and discuss the successes and limitations of current approaches.

Equilibration of complexes of DNA and H-NS proteins on charged surfaces: a coarse-grained model point of view

The Histone-like Nucleoid Structuring protein (H-NS) is a nucleoid-associated protein, which is involved in both gene regulation and DNA compaction. Although it is a key player in genome organization by forming bridges between DNA duplexes, the precise structure of complexes of DNA and H-NS proteins is still not well understood. In particular, it is not clear whether the structure of DNA/H-NS complexes in the living cell is similar to that of complexes deposited on mica surfaces, which may be observed by AFM microscopy. A coarse-grained model, which helps getting more insight into this question, is described and analyzed in the present paper. This model is able of describing both the bridging of bacterial DNA by H-NS in the bulk and the deposition and equilibration of the complex on a charged surface. Simulations performed with the model reveal that a slight attraction between DNA and the charged surface is sufficient to let DNA/H-NS complexes reorganize from 3D coils to planar plasmids bridged by H-NS proteins similar to those observed by AFM microscopy. They furthermore highlight the antagonistic effects of the interactions between DNA and the surface. Indeed, increasing these interactions slows down the equilibration of naked plasmids on the surface but, on the other hand, enables a faster equilibration of DNA/H-NS complexes. Based on the distribution of the lifetimes of H-NS bridges and the time evolution of the number of trans-binding protein dimers during equilibration of the complexes on the surface, it is argued that the decrease of the equilibration time of the complex upon increase of the interaction strength between DNA and the surface is ascribable to the associated decrease of the probability to form new bridges between DNA and the proteins.

Effect of thermal stability on protein adsorption to silica using homologous aldo-keto reductases

Gaining more insight into the mechanisms governing the behavior of proteins at solid/liquid interfaces is particularly relevant in the interaction of high-value biologics with storage and delivery device surfaces, where adsorption-induced conformational changes may dramatically affect biocompatibility. The impact of structural stability on interfacial behavior has been previously investigated by engineering nonwild-type stability mutants. Potential shortcomings of such approaches include only modest changes in thermostability, and the introduction of changes in the topology of the proteins when disulfide bonds are incorporated. Here we employ two members of the aldo-keto reductase superfamily (alcohol dehydrogenase, AdhD and human aldose reductase, hAR) to gain a new perspective on the role of naturally occurring thermostability on adsorbed protein arrangement and its subsequent impact on desorption. Unexpectedly, we find that during initial adsorption events, both proteins have similar affinity to the substrate and undergo nearly identical levels of structural perturbation. Interesting differences between AdhD and hAR occur during desorption and both proteins exhibit some level of activity loss and irreversible conformational change upon desorption. Although such surface-induced denaturation is expected for the less stable hAR, it is remarkable that the extremely thermostable AdhD is similarly affected by adsorption-induced events. These results question the role of thermal stability as a predictor of protein adsorption/desorption behavior.

Effects of tethering a multistate folding protein to a surface

Protein/surface interactions are important in a variety of fields and devices, yet fundamental understanding of the relevant phenomena remains fragmented due to resolution limitations of experimental techniques. Molecular simulation has provided useful answers, but such studies have focused on proteins that fold through a two-state process. This study uses simulation to show how surfaces can affect proteins which fold through a multistate process by investigating the folding mechanism of lysozyme (PDB ID: 7LZM). The results demonstrate that in the bulk 7LZM folds through a process with four stable states: the folded state, the unfolded state, and two stable intermediates. The folding mechanism remains the same when the protein is tethered to a surface at most residues; however, in one case the folding mechanism changes in such a way as to eliminate one of the intermediates. An analysis of the molecular configurations shows that tethering at this site is advantageous for protein arrays because the active site is both presented to the bulk phase and stabilized. Taken as a whole, the results offer hope that rational design of protein arrays is possible once the behavior of the protein on the surface is ascertained.

Protein stability at a carbon nanotube interface

The interactions of proteins with solid surfaces occur in a variety of situations. Motivated by the many nanoengineering applications of protein-carbon nanotube hybrids, we investigate the conformational transitions of hen egg white lysozyme adsorbed on a carbon nanotube. Using a C(α) structure-based model and replica exchange molecular dynamics, we show how the folding/unfolding equilibrium of the adsorbed protein varies with the strength of its coupling to the surface. The stability of the native state depends on the balance between the favorable entropy and unfavorable enthalpy change on adsorption. In the case of a weakly attractive surface when the former dominates, the protein is stabilized. In this regime, the protein can fold and unfold while maintaining the same binding fraction. With increasing surface attraction, the unfavorable enthalpic effect dominates, the native state is destabilized, and the protein has to extensively unbind before changing states from unfolded to folded. At the highest surface coupling, the entropic penalty of folding vanishes, and a folding intermediate is strongly stabilized. In this intermediate state, the α-domain of lysozyme is disrupted, while the β-sheet remains fully structured. We rationalize the relative stability of the two domains on the basis of the residue contact order.

Surface orientation of magainin 2: molecular dynamics simulation and sum frequency generation vibrational spectroscopic studies

We combined molecular dynamics based free energy calculations with sum frequency generation (SFG) spectroscopy to study the orientational distribution of solvated peptides near hydrophobic surfaces. Using a simplified atomistic model of the polystyrene (PS) surface, molecular dynamics simulations have been applied to compute the orientational probability of an α-helical peptide, magainin 2, with respect to the PS/water interface. Free energy calculations revealed that the preferred (horizontal) peptide orientation was driven by the favorable interactions between the hydrophobic PS surface and the hydrophobic residues on the helix, and additional simulations examined the importance of small aggregate formation. Concentration-dependent measurements obtained via SFG vibrational spectroscopy suggest that, at very low peptide concentrations, magainin molecules tend to lie down at the PS/solution interface, which correlates well with the simulation results. When the concentration is increased, peptides exhibit behavior not captured by MD simulations using single helical peptides. A combination of simulations and experiments was shown to yield more reliable results with molecular-level insights into interaction between peptides and polymer surfaces.

End-to-surface reaction dynamics of a single surface-attached DNA or polypeptide

The dynamics of surface-attached polymers play a key role in the operation of a number of biological sensors, yet its current understanding is rather limited. Here we use computer simulations to study the dynamics of a reaction between the free end of a polymer chain and a surface, to which its other end has been attached. We consider two limiting cases, the diffusion-controlled limit, where the reaction is accomplished whenever the free chain end diffuses to within a specified distance from the surface, and the reaction-controlled limit, where slow, intrinsic reaction kinetics rather than diffusion of the chain is rate limiting. In the diffusion-controlled limit, we find that the overall rate scales as N(-b), where N is the number of monomers in the chain and b approximately = 2.2 for excluded volume chains. This value of the scaling exponent b is close to that derived from a simple approximate theory treating the dynamics of the chain end relative to the surface as one-dimensional diffusion in an effective potential. In the reaction-controlled limit, the value of the scaling exponent b is close to 1. We compare our findings with those for the related (and better studied) problem of end-to-end reactions within an unconstrained polymer chain and discuss their implications for electrochemical DNA sensors.

An entropic perspective of protein stability on surfaces

The interaction of proteins with surfaces regulates numerous processes in nature, science, and technology. In many applications, it is desirable to place proteins on surfaces in an active state, and tethering represents one manner in which to accomplish this. However, a clear understanding of how tether placement and design affects protein activity is lacking. Available theoretical models predict that proteins will be stabilized when tethered to substrates. Such models suggest that the surface reduces the number of states accessible to the unfolded state of the protein, thereby reducing the entropic cost of folding on the surface compared to the bulk case. Recent studies, however, have shown that this stabilization is not always seen. The purpose of this article is to determine the validity of the theory with a thorough thermodynamic analysis of the folding of peptides attached to surfaces. Configuration-temperature-density-of-states Monte Carlo simulations are used to examine the behavior of four different peptides of different secondary and tertiary structure. It is found that the surface does reduce the entropic cost of folding for tethered peptides, as the theory suggests. This effect, however, does not always translate into improved stability because the surface may also have a destabilizing enthalpic effect. The theory neglects this effect and assumes that the enthalpy of folding is the same on and off the surface. Both the enthalpic and entropic contributions to the stability are found to be topology- and tether-placement-specific; we show that stability cannot be predicted a priori. A detailed analysis of the folding of protein A shows how the same protein can be both stabilized and destabilized on a surface depending upon how the tethering enhances or hinders the ability of the peptide to form correct tertiary structures.

Peptide conformations for a microarray surface-tethered epitope of the tumor suppressor p53

Peptides or proteins near surfaces exhibit different structural properties from those present in a homogeneous solution, and these differences give rise to varied biological activity. Therefore, understanding the detailed molecular structure of these molecules tethered to a surface is important for interpreting the performance of the various microarrays based on the activities of the immobilized peptides or proteins. We performed molecular dynamics simulations of a pentapeptide, RHSVV, an epitope of the tumor suppressor protein p53, tethered via a spacer on a functionalized silica surface and free in solution, to study their structural and conformational differences. These calculations allowed analyses of the peptide-surface interactions, the sequence orientations, and the translational motions of the peptide on the surface to be performed. Conformational similarities are found among dominant structures of the tethered and free peptide. In the peptide microarray simulations, the peptide fluctuates between a parallel and tilted orientation driven in part by the hydrophobic interactions between the nonpolar peptide residues and the methyl-terminated silica surface. The perpendicular movement of the peptide relative to the surface is also restricted due to the hydrophobic nature of the microarray surface. With regard to structures available for recognition and binding, we find that similar conformations to those found in solution are available to the peptide tethered to the surface, but with a shifted equilibrium constant. Comparisons with experimental results show important implications of this for peptide microarray design and assays.

Effects of surface tethering on protein folding mechanisms

The folding mechanisms of proteins are increasingly being probed through single-molecule experiments in which the protein is immobilized on a surface. Nevertheless, a clear understanding of how the surface might affect folding, and whether or not it changes folding from its bulk behavior, is lacking. In this work, we use molecular dynamics simulations of a model beta-barrel protein tethered to a surface to systematically investigate how the surface impacts folding. In the bulk, this protein folds in a three-state manner through a compact intermediate state, and its transition state (TS) has a well formed hydrophobic core. Upon tethering, we find that folding rates and stability are impacted differently by the surface, with dependencies on both the length and location of the tether. Significant changes in folding times are observed for tether points that do not alter the folding temperature. Tethering also locally enhances the formation of structure for residues proximal to the tether point. We find that neither the folding mechanism nor the TS of this protein are altered if the tether is in a fully structured or completely unstructured region of the TS. By contrast, tethering in a partially structured region of the TS leads to dramatic changes. For one such tether point, the intermediate present in bulk folding is eliminated, leading to a two-state folding process with a heterogeneous, highly unstructured TS ensemble. These results have implications for both the design of single-molecule experiments and biotechnological applications of tethered proteins.

Confinement effects on the thermodynamics of protein folding: Monte Carlo simulations

The effects of chaperonin-like cage-induced confinement on protein stability have been studied for molecules of varying sizes and topologies. Minimalist models based on Gō-like interactions are employed for the proteins, and density-of-states-based Monte Carlo simulations are performed to accurately characterize the thermodynamic transitions. This method permits efficient sampling of conformational space and yields precise estimates of free energy and entropic changes associated with protein folding. We find that confinement-driven stabilization is not only dependent on protein size and cage radius, but also on the specific topology. The choice of the confining potential is also shown to have an effect on the observed stabilization and the scaling behavior of the stabilization with respect to the cage size.