Assessing the predictive capabilities of design heuristics and coarse-grain simulation toward understanding and optimizing site-specific covalent immobilization of ß-lactamase.

For industrial applications, covalent immobilization of enzymes provides minimum leakage, recoverability, reusability, and high stability. Yet, the suitability of a given site on the enzyme for immobilization remains a trial-and-error procedure. Here, we investigate the reliability of design heuristics and a coarse-grain molecular simulation in predicting the optimum sites for covalent immobilization of TEM-1 ß-lactamase. We utilized Escherichia coli-lysate-based cell-free protein synthesis (CFPS) to produce variants containing a site-specific incorporated unnatural amino acid with a unique moiety to facilitate site directed covalent immobilization. To constrain the number of potential immobilization sites, we investigated the predictive capability of several design heuristics. The suitability of immobilization sites was determined by analyzing expression yields, specific activity, immobilization efficiency, and stability of variants. These experimental findings are compared with coarse-grain simulation of TEM-1 domain stability and thermal stability and analyzed for a priori predictive capabilities. This work demonstrates that the design heuristics successfully identify a subset of locations for experimental validation. Specifically, the nucleotide following amber stop codon and domain stability correlate well with the expression yield and specific activity of the variants, respectively. Our approach highlights the advantages of combining coarse-grain simulation and high-throughput experimentation using CFPS to identify optimal enzyme immobilization sites.

Assessing site-specific PEGylation of TEM-1 ß-lactamase with cell-free protein synthesis and coarse-grained simulation.

PEGylation is a broadly used strategy to enhance the pharmacokinetic properties of therapeutic proteins. It is well established that the location and extent of PEGylation have a significant impact on protein properties. However, conventional PEGylation techniques have limited control over PEGylation sites. Emerging site-specific PEGylation technology provides control of PEG placement by conjugating PEG polymers via click chemistry reaction to genetically encoded non-canonical amino acids. Unfortunately, a method to rapidly determine the optimal PEGylation location has yet to be established. Here we seek to address this challenge. In this work, coarse-grained molecular dynamic simulations are paired with high-throughput experimental screening utilizing cell-free protein synthesis to investigate the effect of site-specific PEGylation on the two-state folder protein TEM-1 ß-lactamase. Specifically, the conjugation efficiency, thermal stability, and enzymatic activity are studied for the enzyme PEGylated at several different locations. The results of this analysis confirm that the physical properties of the PEGylated protein vary considerably with PEGylation site and that traditional design recommendations are insufficient to predict favorable PEGylation sites. In this study, the best predictor of the most favorable conjugation site is coarse-grained simulation. Thus, we propose a dual combinatorial screening approach in which coarse-grained molecular simulation informs site selection for high-throughput experimental verification.

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 p-Azidophenylalanine Incorporation on Protein Structure and Stability.

Functionalization is often needed to harness the power of proteins for beneficial use but can cause losses to stability and/or activity. State of the art methods to limit these deleterious effects accomplish this by substituting an amino acid in the wild-type molecule into an unnatural amino acid, such as p-azidophenylalanine (pAz), but selecting the residue for substitution a priori remains an elusive goal of protein engineering. The results of this work indicate that all-atom molecular dynamics simulation can be used to determine whether substituting pAz for a natural amino acid will be detrimental to experimentally determined protein stability. These results offer significant hope that local deviations from wild-type structure caused by pAz incorporation observed in simulations can be a predictive metric used to reduce the number of costly experiments that must be done to find active proteins upon substitution with pAz and subsequent functionalization.

Parameterization of Unnatural Amino Acids with Azido and Alkynyl R-Groups for Use in Molecular Simulations.

Recent new methods to functionalize proteins at specific amino acid locations use unnatural amino acids that contain azido and alkynyl groups. This capability is unprecedented and enables the creation of site-specific protein devices. Because of the high specificity of these devices, many protein configurations are possible and in silico screens have shown promise in predicting optimal attachment site locations. Therefore, there is significant interest in improving current molecular dynamics (MD) models to include the unique chemistries of these linear moieties. This work uses the force field tool kit to obtain the bonded and nonbonded CHARMM parameters for small molecules that contain azido and alkynyl groups. Next, the reliability of these parameters is tested by running simulated MD analysis to prove that the modeled structures match those found in the literature and quantum theory. Finally, the protein MD simulation compares this parameter set with crystallographic data to give a greater understanding of unnatural amino acid influence on the protein structure.

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 Locational Impact of Site-Specific PEGylation: Streamlined Screening with Cell-Free Protein Expression and Coarse-Grain Simulation.

Although polyethylene glycol (PEG) is commonly used to improve protein stability and therapeutic efficacy, the optimal location for attaching PEG onto proteins is not well understood. Here, we present a cell-free protein synthesis-based screening platform that facilitates site-specific PEGylation and efficient evaluation of PEG attachment efficiency, thermal stability, and activity for different variants of PEGylated T4 lysozyme, including a di-PEGylated variant. We also report developing a computationally efficient coarse-grain simulation model as a potential tool to narrow experimental screening candidates. We use this simulation method as a novel tool to evaluate the locational impact of PEGylation. Using this screen, we also evaluated the predictive impact of PEGylation site solvent accessibility, conjugation site structure, PEG size, and double PEGylation. Our findings indicate that PEGylation efficiency, protein stability, and protein activity varied considerably with PEGylation site, variations that were not well predicted by common PEGylation guidelines. Overall our results suggest current guidelines are insufficiently predictive, highlighting the need for experimental and simulation screening systems such as the one presented here.

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.

Uncertainty quantification and propagation of errors of the Lennard-Jones 12-6 parameters for n-alkanes.

Molecular simulation has the ability to predict various physical properties that are difficult to obtain experimentally. For example, we implement molecular simulation to predict the critical constants (i.e., critical temperature, critical density, critical pressure, and critical compressibility factor) for large n-alkanes that thermally decompose experimentally (as large as C48). Historically, molecular simulation has been viewed as a tool that is limited to providing qualitative insight. One key reason for this perceived weakness in molecular simulation is the difficulty to quantify the uncertainty in the results. This is because molecular simulations have many sources of uncertainty that propagate and are difficult to quantify. We investigate one of the most important sources of uncertainty, namely, the intermolecular force field parameters. Specifically, we quantify the uncertainty in the Lennard-Jones (LJ) 12-6 parameters for the CH4, CH3, and CH2 united-atom interaction sites. We then demonstrate how the uncertainties in the parameters lead to uncertainties in the saturated liquid density and critical constant values obtained from Gibbs Ensemble Monte Carlo simulation. Our results suggest that the uncertainties attributed to the LJ 12-6 parameters are small enough that quantitatively useful estimates of the saturated liquid density and the critical constants can be obtained from molecular simulation.

Probing the effects of surface hydrophobicity and tether orientation on antibody-antigen binding.

Antibody microarrays have the potential to revolutionize molecular detection for many applications, but their current use is limited by poor reliability, and efforts to change this have not yielded fruitful results. One difficulty which limits the rational engineering of next-generation devices is that little is known, at the molecular level, about the antibody-antigen binding process near solid surfaces. Atomic-level structural information is scant because typical experimental techniques (X-ray crystallography and NMR) cannot be used to image proteins bound to surfaces. To overcome this limitation, this study uses molecular simulation and an advanced, experimentally validated, coarse-grain, protein-surface model to compare fab-lysozyme binding in bulk solution and when the fab is tethered to hydrophobic and hydrophilic surfaces. The results show that the tether site in the fab, as well as the surface hydrophobicity, significantly impacts the binding process and suggests that the optimal design involves tethering fabs upright on a hydrophilic surface. The results offer an unprecedented, molecular-level picture of the binding process and give hope that the rational design of protein-microarrays is possible.

An improved statistical analysis for predicting the critical temperature and critical density with Gibbs ensemble Monte Carlo simulation.

A rigorous statistical analysis is presented for Gibbs ensemble Monte Carlo simulations. This analysis reduces the uncertainty in the critical point estimate when compared with traditional methods found in the literature. Two different improvements are recommended due to the following results. First, the traditional propagation of error approach for estimating the standard deviations used in regression improperly weighs the terms in the objective function due to the inherent interdependence of the vapor and liquid densities. For this reason, an error model is developed to predict the standard deviations. Second, and most importantly, a rigorous algorithm for nonlinear regression is compared to the traditional approach of linearizing the equations and propagating the error in the slope and the intercept. The traditional regression approach can yield nonphysical confidence intervals for the critical constants. By contrast, the rigorous algorithm restricts the confidence regions to values that are physically sensible. To demonstrate the effect of these conclusions, a case study is performed to enhance the reliability of molecular simulations to resolve the n-alkane family trend for the critical temperature and critical density.

Communication: Antibody stability and behavior on surfaces.

Antibody microarrays have the potential to revolutionize molecular detection in scientific, medical, and other biosensor applications, but their current use is limited because of poor reliability. It is hypothesized that one reason for their poor performance results from strong antibody-surface interactions that destabilize the antibody structure and create steric interference for antigen recognition. Using a recently developed coarse-grain protein-surface model that has been parameterized against experimental data, antibody-surface interactions for two antibody orientations on two types of surfaces have been investigated. The results show that regardless of attachment geometry, antibodies tend to collapse onto hydrophobic surfaces and exhibit lower overall stability compared to antibodies on hydrophilic surfaces or in bulk solution. The results provide an unprecedented view into the dynamics of antibodies on surfaces and offer new insights into the poor performance exhibited by current antibody microarrays.

The effects of multiple probes on the hybridization of target DNA on surfaces.

DNA microarrays have disruptive potential in many fields including genetics and medicine, but the technology has yet to find widespread clinical use due to poor reliability. Microarrays work on the principle of hybridization and can only be as dependable as this process is reliable. As such, a significant amount of theoretical research has been done to understand hybridization on surfaces on the molecular level. Previous simulations of a target strand with a single, surface-tethered probe molecule have yielded valuable insights, but such is an ideal system and little is known about the effects of multiple probes-a situation that more closely approximates the real system. This work uses molecular simulation to determine the specific differences in duplex stability between one, three, six, and nine tethered probes on a surface. The results show that it is more difficult for a single target to hybridize to a probe as the number of probes on the surface increases due to crowding effects; however, once hybridized, the duplex is more stable than when fewer probes are present. The data also indicate that hybridization of a target to a probe on the face of a group of probes is more stable than hybridization to probes at the edge or center locations. Taken as a whole, the results offer new insights into the cause of the poor reproducibility exhibited by microarrays.

Communication: Using multiple tethers to stabilize proteins on surfaces.

Protein surface interactions are important in many applications in biotechnology including protein arrays, but these technologies have not lived up to their transformative potential because it is difficult to attach proteins to surfaces in a manner that preserves function and theoretical understanding of the relevant phenomena remains limited. Here is reported the effect of using multiple tethers to attach a protein (lysozyme) to a surface and the effects on the structure and stability of the molecule. The simulations show how using two tethers can drastically change the folding mechanism such that a protein that is initially unstable and inactive when attached using a single tether can become more stable and functional when two tethers are used. The results offer hope that the rational design of protein arrays is possible.

A coarse grain model for protein-surface interactions.

The interaction of proteins with surfaces is important in numerous applications in many fields-such as biotechnology, proteomics, sensors, and medicine--but fundamental understanding of how protein stability and structure are affected by surfaces remains incomplete. Over the last several years, molecular simulation using coarse grain models has yielded significant insights, but the formalisms used to represent the surface interactions have been rudimentary. We present a new model for protein surface interactions that incorporates the chemical specificity of both the surface and the residues comprising the protein in the context of a one-bead-per-residue, coarse grain approach that maintains computational efficiency. The model is parameterized against experimental adsorption energies for multiple model peptides on different types of surfaces. The validity of the model is established by its ability to quantitatively and qualitatively predict the free energy of adsorption and structural changes for multiple biologically-relevant proteins on different surfaces. The validation, done with proteins not used in parameterization, shows that the model produces remarkable agreement between simulation and experiment.

Exploring the mechanisms of DNA hybridization on a surface.

DNA microarrays are a potentially disruptive technology in the medical field, but their use in such settings is limited by poor reliability. Microarrays work on the principle of hybridization and can only be as reliable as this process is robust, yet little is known at the molecular level about how the surface affects the hybridization process. This work uses advanced molecular simulation techniques and an experimentally parameterized coarse-grain model to determine the mechanism by which hybridization occurs on surfaces. The results show that hybridization proceeds through a mechanism where the untethered (target) strand often flips orientation. For evenly lengthed strands, the surface stabilizes hybridization (compared to the bulk system) by reducing the barriers involved in the flipping event. For unevenly lengthed strands, the surface destabilizes hybridization compared to the bulk, but the degree of destabilization is dependent on the location of the matching sequence. Taken as a whole, the results offer an unprecedented view into the hybridization process on surfaces and provide some insights as to the poor reproducibility exhibited by microarrays.

Multiple molecule effects on the cooperativity of protein folding transitions in simulations.

Though molecular simulation of proteins has made notable contributions to the study of protein folding and kinetics, disagreement between simulation and experiment still exists. One of the criticisms levied against simulation is its failure to reproduce cooperative protein folding transitions. This weakness has been attributed to many factors such as a lack of polarizability and adequate capturing of solvent effects. This work, however, investigates how increasing the number of proteins simulated simultaneously can affect the cooperativity of folding transitions--a topic that has received little attention previously. Two proteins are studied in this work: phage T4 lysozyme (Protein Data Bank (PDB) ID: 7LZM) and phage 434 repressor (PDB ID: 1R69). The results show that increasing the number of proteins molecules simulated simultaneously leads to an increase in the macroscopic cooperativity for transitions that are inherently cooperative on the molecular level but has little effect on the cooperativity of other transitions. Taken as a whole, the results identify one area of consideration to improving simulations of protein folding.

Thermodynamics of DNA hybridization on surfaces.

Hybridization of single-stranded DNA (ssDNA) targets to surface-tethered ssDNA probes was simulated using an advanced coarse-grain model to identify key factors that influence the accuracy of DNA microarrays. Comparing behavior in the bulk and on the surface showed, contrary to previous assumptions, that hybridization on surfaces is more thermodynamically favorable than in the bulk. In addition, the effects of stretching or compressing the probe strand were investigated as a model system to test the hypothesis that improving surface hybridization will improve microarray performance. The results in this regard indicate that selectivity can be increased by reducing overall sensitivity by a small degree. Taken as a whole, the results suggest that current methods to enhance microarray performance by seeking to improve hybridization on the surface may not yield the desired outcomes.

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.

Predicting stability of alpha-helical, orthogonal-bundle proteins on surfaces.

The interaction of proteins with surfaces is a key phenomenon in many applications, but current understanding of the biophysics involved is lacking. At present, rational design of such emerging technologies is difficult as no methods or theories exist that correctly predict how surfaces influence protein behavior. Using molecular simulation and a coarse-grain model, this study illustrates for the first time that stability of proteins on surfaces can be correlated with tertiary structural elements for alpha-helical, orthogonal-bundle proteins. Results show that several factors contribute to stability on surfaces including the nature of the loop region where the tether is placed and the ability of the protein to freely rotate on the surface. A thermodynamic analysis demonstrates that surfaces stabilize proteins entropically and that any destabilization is an enthalpic effect. Moreover, the entropic effects are concentrated on the unfolded state of the protein while the ethalpic effects are focused on the folded state.