Evaluation of catalytic free energies in genetically modified proteins

A combination of the empirical valence bond method and a free energy perturbation approach is used to simulate the activity of genetically modified enzymes. The simulations reproduce in a semiquantitative way the observed effects of mutations on the activity and binding free energies of trypsin and subtilisin. This suggests that we are approaching a stage of quantitative structure-function correlation of enzymes. The analysis of the calculations points towards the electrostatic energy of the reacting system as the key factor in enzyme catalysis. The changes in the charges of the reacting system and the corresponding changes in "solvation" free energy (generalized here as the interaction between the charges and the given microenvironment) are emphasized. It is argued that a reliable evaluation of these changes might be sufficient for correlating structure and catalysis. The use of free energy perturbation methods and thermodynamic cycles for evaluation of solvation energies and reactivity is discussed, pointing out our early contributions. The apparent elaborated nature of our treatment is clarified, explaining that such a treatment is essential for consistent calculations of chemical reactions in polar environments. The problems associated with seemingly more rigorous quantum mechanical methods are discussed, emphasizing the inconsistency associated with using gas phase charge distributions. The importance of dynamic aspects is examined by evaluating the autocorrelation of the protein "reaction field" on the reacting substrate. It is found that, at least in the present case, dynamic effects are not important. The nature of the catalytic free energy is considered, arguing that the protein provides preoriented dipoles (polarized to stabilize the transition state charge distribution) and small reorganization energy, thus reducing the activation free energy. The corresponding catalytic free energy is related to the folding free energy, which is being invested in aligning the active site dipoles.     

Evaluation of protein-protein association energies by free energy perturbation calculations

The association energy upon binding of different amino acids in the specificity pocket of trypsin was evaluated by free energy perturbation calculations on complexes between bovine trypsin (BT) and bovine pancreatic trypsin inhibitor (BPTI). Three simulations of mutations of the primary binding residue (P(1)) were performed (P(1)-Ala to Gly, P(1)-Met to Gly and P(1)-Met to Ala) and the resulting differences in association energy (DeltaDeltaG(a)) are 2. 28, 5.08 and 2.93 kcal/mol for P(1)-Ala to Gly, P(1)-Met to Gly and to Ala with experimental values of 1.71, 4.62 and 2.91 kcal/mol, respectively. The calculated binding free energy differences are hence in excellent agreement with the experimental binding free energies. The binding free energies, however, were shown to be highly dependent on water molecules at the protein-protein interface and could only be quantitatively estimated if the correct number of such water molecules was included. Furthermore, the cavities that were formed when a large amino acid side-chain is perturbed to a smaller one seem to create instabilities in the systems and had to be refilled with water molecules in order to obtain reliable results. In addition, if the protein atoms that were perturbed away were not replaced by water molecules, the simulations dramatically overestimated the initial state of the free energy perturbations.     

Poisson-Boltzmann calculations of nonspecific salt effects on protein-protein binding free energies

The salt dependence of the binding free energy of five protein-protein hetero-dimers and two homo-dimers/tetramers was calculated from numerical solutions to the Poisson-Boltzmann equation. Overall, the agreement with experimental values is very good. In all cases except one involving the highly charged lactoglobulin homo-dimer, increasing the salt concentration is found both experimentally and theoretically to decrease the binding affinity. To clarify the source of salt effects, the salt-dependent free energy of binding is partitioned into screening terms and to self-energy terms that involve the interaction of the charge distribution of a monomer with its own ion atmosphere. In six of the seven complexes studied, screening makes the largest contribution but self-energy effects can also be significant. The calculated salt effects are found to be insensitive to force-field parameters and to the internal dielectric constant assigned to the monomers. Nonlinearities due to high charge densities, which are extremely important in the binding of proteins to negatively charged membrane surfaces and to nucleic acids, make much smaller contributions to the protein-protein complexes studied here, with the exception of highly charged lactoglobulin dimers. Our results indicate that the Poisson-Boltzmann equation captures much of the physical basis of the nonspecific salt dependence of protein-protein complexation.     

The orientationally controlled assembly of genetically modified enzymes in an amperometric biosensor

A novel, nitroreductase (NTR) containing a sequence of six cysteine amino acids, enabling strong thiolate bonds to form on a gold electrode surface without the loss of enzyme activity, was genetically engineered. The enzyme was directly immobilised at a gold electrode without the need for pre-treatment of the surface with a self-assembled monolayer or a conducting polymer. The ensemble was used to develop an amperometric biosensor for the detection of explosives containing nitroaromatic compounds. Preliminary results demonstrate detection levels down to 100 parts per trillion, signifying tremendous promise towards an in situ sensor for the detection of explosives.     

Detection of glycoalkaloids using disposable biosensors based on genetically modified enzymes

In this work we present a rapid, selective, and highly sensitive detection of α-solanine and α-chaconine using cholinesterase-based sensors. The high sensitivity of the devices is brought by the use of a genetically modified acetylcholinesterase (AChE), combined with a one-step detection method based on the measurement of inhibition slope. The selectivity was obtained by using butyrylcholinesterase (BChE), an enzyme able to detect these two toxins with differential inhibition kinetics. The enzymes were immobilized via entrapment in PVA-AWP polymer directly on the working electrode surface. The analysis of the resulting inhibition slope was performed employing linear regression function included in Matlab. The high toxicity of α-chaconine compared to α-solanine due to a better affinity to the active site was proved. The inhibition of glycoalkaloids (GAs) mixture was performed over AChE enzyme wild-type AChE and BChE biosensors resulting in the detection of synergism effect. The developed method allows the detection of (GAs) at 50 ppb in potato matrix.     

Characterization of digestive enzymes of bruchid parasitoids-initial steps for environmental risk assessment of genetically modified legumes

Genetically modified (GM) legumes expressing the α-amylase inhibitor 1 (αAI-1) from Phaseolus vulgaris L. or cysteine protease inhibitors are resistant to several bruchid pests (Coleoptera: Chrysomelidae). In addition, the combination of plant resistance factors together with hymenopteran parasitoids can substantially increase the bruchid control provided by the resistance alone. If the strategy of combining a bruchid-resistant GM legume and biological control is to be effective, the insecticidal trait must not adversely affect bruchid antagonists. The environmental risk assessment of such GM legumes includes the characterization of the targeted enzymes in the beneficial species and the assessment of the in vitro susceptibility to the resistance factor. The digestive physiology of bruchid parasitoids remain relatively unknown, and their susceptibility to αAI-1 has never been investigated. We have detected α-amylase and serine protease activities in all five bruchid parasitoid species tested. Thus, the deployment of GM legumes expressing cysteine protease inhibitors to control bruchids should be compatible with the use of parasitoids. In vitro inhibition studies showed that sensitivity of α-amylase activity to αAI-1 in the parasitoids was comparable to that in the target species. Direct feeding assays revealed that harmful effects of α-amylase inhibitors on bruchid parasitoids cannot be discounted and need further evaluation.     

Cutaneous melanoma in genetically modified animals

Cutaneous melanomas are tumors originating from skin melanocytes which are present in hair follicles, and interfollicular epidermal and dermal layers. Experimental work in model systems involves in silico, in vitro and in vivo analyses. Such models allow to mimick melanocytic aberrations characteristic of melanoma, and to potentially exploit novel therapies. Transgenic technologies can be used to modify specifically the genome of the model organism and thereby generate transgenic strains, and combinations of such strains, which may develop metastasizing melanoma. In such strains, metastasizing melanoma either arises spontaneously after a period of latency or requires additional physical or chemical induction. In this review, we summarize the work of currently available transgenic melanoma models and discuss the most recent progress in creating improved and/or inducible models reflecting the human disease.     

Analysis of genetically modified oils

Genetically modified oils with altered functional or nutritional characteristics are being introduced into the marketplace. A wide array of analytical techniques has been utilized to facilitate developing these oils. This article attempts to review the utilization of these analytical procedures for characterizing both the chemistry and some functionality of these oils. Although techniques to assess oxidative stability in frying and food applications are covered, measurement of nutritional characteristics are not.     

Ultrafast catalytic processes in enzymes

The study of biocatalysis and biotransformation in the transition-state region has been challenging and difficult, but recent advances on two important photoenzymes in nature, DNA photolyase and protochlorophyllide oxidoreductase, have enabled the investigation of their catalytic processes in real time. By following the entire evolution of substrate transformation, the functional dynamics constituting a series of elementary reactions have been mapped out. The five fundamental reactions in the enzymes, namely electron transfer, bond breaking and making, proton and hydride transfer, all occur ultrafast within subnanosecond. The direct clocking of catalytic transition states probes central, unmasked chemical processes and provides mechanistic insights into the role of the dynamics in enzyme function, which not only facilitates the formation of the enzyme-substrate complex in the transition-state configurations, but also modulates the subsequent catalytic reactions for maximum biotransformation efficiency.     

Lifespan extension in genetically modified mice

Major advances in aging research have been made by studying the effect of genetic modifications on the lifespan of organisms, such as yeast, invertebrates (worms and flies) and mice. Data from yeast and invertebrates have been the most plentiful because of the ease in which genetic manipulations can be made and the rapidity by which lifespan experiments can be performed. With the ultimate focus on advancing human health, testing genetic interventions in mammals is crucial, and the mouse has proven to be the mammal most amenable to this task. Lifespan studies in mice are resource intensive, requiring up to 4 years to complete. Therefore, it is critical that a set of scientifically-based criteria be followed to assure reliable results and establish statistically significant findings so other laboratories can replicate and build on the data. Only then will it be possible to confidently determine that the genetic modification extends lifespan and alters aging.     

Gene flow in genetically modified wheat

Understanding gene flow in genetically modified (GM) crops is critical to answering questions regarding risk-assessment and the coexistence of GM and non-GM crops. In two field experiments, we tested whether rates of cross-pollination differed between GM and non-GM lines of the predominantly self-pollinating wheat Triticum aestivum. In the first experiment, outcrossing was studied within the field by planting "phytometers" of one line into stands of another line. In the second experiment, outcrossing was studied over distances of 0.5-2.5 m from a central patch of pollen donors to adjacent patches of pollen recipients. Cross-pollination and outcrossing was detected when offspring of a pollen recipient without a particular transgene contained this transgene in heterozygous condition. The GM lines had been produced from the varieties Bobwhite or Frisal and contained Pm3b or chitinase/glucanase transgenes, respectively, in homozygous condition. These transgenes increase plant resistance against pathogenic fungi. Although the overall outcrossing rate in the first experiment was only 3.4%, Bobwhite GM lines containing the Pm3b transgene were six times more likely than non-GM control lines to produce outcrossed offspring. There was additional variation in outcrossing rate among the four GM-lines, presumably due to the different transgene insertion events. Among the pollen donors, the Frisal GM line expressing a chitinase transgene caused more outcrossing than the GM line expressing both a chitinase and a glucanase transgene. In the second experiment, outcrossing after cross-pollination declined from 0.7-0.03% over the test distances of 0.5-2.5 m. Our results suggest that pollen-mediated gene flow between GM and non-GM wheat might only be a concern if it occurs within fields, e.g. due to seed contamination. Methodologically our study demonstrates that outcrossing rates between transgenic and other lines within crops can be assessed using a phytometer approach and that gene-flow distances can be efficiently estimated with population-level PCR analyses.     

Binding free energies and free energy components from molecular dynamics and Poisson-Boltzmann calculations. Application to amino acid recognition by aspartyl-tRNA synthetase

Specific amino acid binding by aminoacyl-tRNA synthetases (aaRS) is necessary for correct translation of the genetic code. Engineering a modified specificity into aminoacyl-tRNA synthetases has been proposed as a means to incorporate artificial amino acid residues into proteins in vivo. In a previous paper, the binding to aspartyl-tRNA synthetase of the substrate Asp and the analogue Asn were compared by molecular dynamics free energy simulations. Molecular dynamics combined with Poisson-Boltzmann free energy calculations represent a less expensive approach, suitable for examining multiple active site mutations in an engineering effort. Here, Poisson-Boltzmann free energy calculations for aspartyl-tRNA synthetase are first validated by their ability to reproduce selected molecular dynamics binding free energy differences, then used to examine the possibility of Asn binding to native and mutant aspartyl-tRNA synthetase. A component analysis of the Poisson-Boltzmann free energies is employed to identify specific interactions that determine the binding affinities. The combined use of molecular dynamics free energy simulations to study one binding process thoroughly, followed by molecular dynamics and Poisson-Boltzmann free energy calculations to study a series of related ligands or mutations is proposed as a paradigm for protein or ligand design.The binding of Asn in an alternate, "head-to-tail" orientation observed in the homologous asparagine synthetase is analyzed, and found to be more stable than the "Asp-like" orientation studied earlier. The new orientation is probably unsuitable for catalysis. A conserved active site lysine (Lys198 in Escherichia coli) that recognizes the Asp side-chain is changed to a leucine residue, found at the corresponding position in asparaginyl-tRNA synthetase. It is interesting that the binding of Asp is calculated to increase slightly (rather than to decrease), while that of Asn is calculated, as expected, to increase strongly, to the same level as Asp binding. Insight into the origin of these changes is provided by the component analyses. The double mutation (K198L,D233E) has a similar effect, while the triple mutation (K198L,Q199E,D233E) reduces Asp binding strongly. No binding measurements are available, but the three mutants are known to have no ability to adenylate Asn, despite the "Asp-like" binding affinities calculated here. In molecular dynamics simulations of all three mutants, the Asn ligand backbone shifts by 1-2 A compared to the experimental Asp:AspRS complex, and significant side-chain rearrangements occur around the pocket. These could reduce the ATP binding constant and/or the adenylation reaction rate, explaining the lack of catalytic activity in these complexes. Finally, Asn binding to AspRS with neutral K198 or charged H449 is considered, and shown to be less favorable than with the charged K198 and neutral H449 used in the analysis.     

Detection of genetically modified organisms in foods

Legislation enacted worldwide to regulate the presence of genetically modified organisms (GMOs) in crops, foods and ingredients, necessitated the development of reliable and sensitive methods for GMO detection. In this article, protein- and DNA-based methods employing western blots, enzyme-linked immunosorbant assay, lateral flow strips, Southern blots, qualitative-, quantitative-, real-time- and limiting dilution-PCR methods, are discussed. Where information on modified gene sequences is not available, new approaches, such as near-infrared spectrometry, might tackle the problem of detection of non-approved genetically modified (GM) foods. The efficiency of screening, identification and confirmation strategies should be examined with respect to false-positive rates, disappearance of marker genes, increased use of specific regulator sequences and the increasing number of GM foods.     

Surface free energies of oral streptococci

Abstract The surface free energies ( γ _b) of a variety of oral streptococci were determined from contact angle measurements on bacterial deposits, using the concept of dispersion and polar components. At least four strains of each species were tested. Strains of Streptococcus mutans, S. sanguis and S. salivarius possessed relatively high surface free energies (103 ± 12 mJ · m⁻²) and at the species level no significant difference was found. In contrast, the strains of S. mitis had remarkably low surface free energies (45 ± 14 mJ · m⁻²). S. milleri appeared to be a heterogeneous species, showing surface free energies over a range of 32–119 mJ · m⁻². No significant differences were observed between laboratory strains and strains freshly isolated from the oral cavity.     

Prediction of protein-protein binding free energies

We present an energy function for predicting binding free energies of protein-protein complexes, using the three-dimensional structures of the complex and unbound proteins as input. Our function is a linear combination of nine terms and achieves a correlation coefficient of 0.63 with experimental measurements when tested on a benchmark of 144 complexes using leave-one-out cross validation. Although we systematically tested both atomic and residue-based scoring functions, the selected function is dominated by residue-based terms. Our function is stable for subsets of the benchmark stratified by experimental pH and extent of conformational change upon complex formation, with correlation coefficients ranging from 0.61 to 0.66.     

Theoretical binding energies of inhibitors to enzymes

Modified quantum mechanical calculations predict the binding enthalpies of saccharide inhibitors of lysozyme to within a few kilocalories of experimental measurements.     

Advances in the calculation of binding free energies

1. Download : Download high-res image (147KB)
2. Download : Download full-size image