The Genetic Algorithm and the Conformational Search of Polypeptides and Proteins

Abstract

The genetic algorithm is a technique of function optimization derived from the principles of evolutionary theory. We have adapted it to perform conformational search on polypeptides and proteins. The algorithm was first tested on several small polypeptides and the 46 amino acid protein crambin under the AMBER potential energy function. The probable global minimum conformations of the polypeptides were located 90% of the time and a non-native conformation of crambin was located that was 150kcal/mol lower in potential energy than the minimized crystal structure conformation. Next, we used a knowledge-based potential function to predict the structures of melittin, pancreatic polypeptide, and crambin. A 2.31 Å ΔRMS conformation of melittin and a 5.33 Å ΔRMS conformation of pancreatic polypeptide were located by genetic algorithm-based conformational search under the knowledge-based potential function. Although the ΔRMS of pancreatic polypeptide was somewhat high, most of the secondary structure was correct. The secondary structure of crambin was predicted correctly, but the potential failed to promote packing interactions. Finally, we tested the packing aspects of our potential function by attempting to predict the tertiary structure of cytochrome b ₅₆₂ given correct secondary structure as a constraint. The final predicted conformation of cytochrome b ₅₆₂ was an almost completely extended continuous helix which indicated that the knowledge-based potential was useless for tertiary structure prediction. This work serves as a warning against testing potential functions designed for tertiary structure prediction on small proteins.     

Acid phosphatase from maize scutellum: Succinylation and some kinetic properties of the active enzyme

Acid phosphatase purified from maize scutellum did not dissociate into subunits upon acylation with succinic anhydride. The enzyme maintained its catalytic activity after succinylation of 52 free amino groups permolecule. The results also showed that free amino groups may play an important role in the maintenance of enzyme stability at pH values greater than 5.4.     

Design and evaluation of succinylated soy protein tablets as delayed drug delivery systems

The impact of succinylation on soy proteins as excipients for delayed delivery of drugs in the gastrointestinal tract was studied. Succinylation decreased protein solubility and protein charge density at pH 1.2 and increased solubility and zeta potential at pH above 4.5. Tablet erosion and swelling were decreased at pH 1.2 and increased at pH 7.5. FTIR analysis indicated polypeptide chain unfolding as a result of succinylation. Tablets of protein succinylated 50% or 100% released less than 10% of loaded riboflavin or rifampicin in 2 h at gastric pH in the presence of pepsin but released these compounds rapidly at intestinal pH. Succinylated soy protein tablets were thus gastroresistant, suggesting their use as excipients for controlled release of medicinal or nutraceutical agents.     

Effect of succinylation of antibodies on their conformation and interaction with the antigen

Using succinic anhydride, a succinylated derivative of anti-urease IgG having 49 ± 6% modification was prepared and its physicochemical and immunological properties were studied. IgG undergoes substantial changes in its native conformation on succinylation, which was mainly attributed to electrostatic destabilization of the native protein conformation. The modified IgG exhibited a decrease in its cross-reactivity with urease. This decrease is attributed to the conformational change in IgG upon succinylation and/or is due to the disruption of the lysine residues in the antigen-binding site of IgG upon succinylation, which may be involved in binding the antigen. IgG was able to bind to the specific antigen although its conformation was partially modified. Therefore, partial modification of the conformation of the antigen-binding site of IgG is permissible in order to bind to the antigen. Published in Russian in Biokhimiya, 2006, Vol. 71, No. 12, pp. 1642–1647.     

Changes in conformation and immunological activity of ovalbumin during its modification with different acid anhydrides.

In order to probe the cause and nature of conformational changes induced by the chemical modification of amino groups in proteins, five acylated derivatives of ovalbumin namely 21% acetylated, 32% succinylated, 90% butyrated 92% succinylated, and 95% acetylated ovalbumins were prepared and their molecular and immunological properties were systematically investigated. As evidenced by the ultraviolet difference spectral, solvent perturbation, gel filtration, and viscosity data, acylation of the amino groups produced a definite conformational change in native ovalbumin whose extent was higher for higher degrees of chemical modification. The solvent pertubation data showed an exposure of 0.5 tryptophan and 3 tyrosine residues in native ovalbumin; the exposure increased to 1 tryptophan and about 5 tyrosine residues in the maximally modified proteins (i.e. 90% butyrated, 92% succinylated, and 95% acetylated ovalbumins). The Stokes radius (2.7 nm) and intrinsic viscosity (3.9 ml/g) of ovalbumin increased, respectively, to about 3.4 nm and 7.7 ml/g upon acylation of its 18 lysine residues; the intrinsic viscosity of 95% acetylated ovalbumin was 7.2 ml/g. The reduced viscosity of ovalbumin (4.2 ml/g) which remained unaltered on raising the pH to pH 11.2, increased to 7.9 ml/g on succinylation of 18 lysine residues. On raising the ionic strength from 0.15 to 1.0, the value decreased from 7.9 to 6.2 ml/g. These observations taken together with the fact that the intrinsic viscosities of 92% succinylated and 90% butyrated ovalbumins are identical, argue against the presently prevalent proposal that electrostatic effects alone are responsible for the disruption of native protein conformation during chemical modification. The immunological activity of ovalbumin towards rabbit anti-ovalbumin expectedly decreased with acylation of its amino groups but the three maximally modified ovalbumins retained 40% immunological activity. This taken along with the spectral and viscosity data showed substantial native structure (format) in the three maximally acylated derivatives. The rabbit antiserum against 95% acetylated ovalbumin did not cross-react with acetylated lysozyme and reacted poorly with the native and 92% succinylated ovalbumins suggesting that the antigenic make-up of the three maximally modified ovalbumins is different.     

Effect of lysine succinylation on the regulation of 2-oxoglutarate dehydrogenase inhibitor,OdhI, involved in glutamate production in Corynebacterium glutamicum

In Corynebacterium glutamicum, the activity of the 2-oxoglutarate dehydrogenase (ODH) complex is negatively regulated by the unphosphorylated form of OdhI protein, which is critical for L-glutamate overproduction. We examined the potential impact of protein acylation at lysine (K)-132 of OdhI in C. glutamicum ATCC13032. The K132E succinylation-mimic mutation reduced the ability of OdhI to bind OdhA, the catalytic subunit of the ODH complex, which reduced the inhibition of ODH activity. In vitro succinylation of OdhI protein also reduced the ability to inhibit ODH, and the K132R mutation blocked the effect. These results suggest that succinylation at K132 may attenuate the OdhI function. Consistent with these results, the C. glutamicum mutant strain with OdhI-K132E showed decreased L-glutamate production. Our results indicated that not only phosphorylation but also succinylation of OdhI protein may regulate L-glutamate production in C. glutamicum.     

Understanding the role of internal lysine residues of serum albumins in conformational stability and bilirubin binding

The role of internal lysine residues of different serum albumins, viz. from human, rabbit, goat, sheep and buffalo (HSA, RbSA, GSA, SSA and BuSA), in conformational stability and bilirubin binding was investigated after blocking them using acetylation, succinylation and guanidination reactions. No significant change in the secondary structure was noticed whereas the tertiary structure of these proteins was slightly altered upon acetylation or succinylation as revealed by circular dichroism (CD), fluorescence and gel filtration results. Guanidination did not affect the native protein conformation to a measurable extent. Scatchard analysis, CD and absorption spectroscopic results showed marked reductions (5-21-fold decrease in K(a) and approximately 50% decrease in the CD Cotton effect intensity) in the affinity of albumins for bilirubin upon acetylation or succinylation whereas guanidination produced a small change. Interestingly, monosignate CD spectra of bilirubin complexed with GSA, SSA and BuSA were transformed to bisignate CD spectra upon acetylation or succinylation of internal lysine residues whereas spectra remained bisignate in the case of bilirubin bound to acetylated or succinylated derivatives of HSA and RbSA. When probed by CD spectroscopy, bilirubin bound to acetylated or succinylated derivatives of GSA and SSA rapidly switched over to native albumins and not vice versa. These results suggested that salt linkage(s) contributed by internal lysine residue(s) play an important role in the high-affinity binding of bilirubin to albumin and provide stability to the native three-dimensional conformation of the bound pigment. Chloroform severely decreased the intensity of both positive and negative CD Cotton effects of bilirubin complexed with acetylated or succinylated derivatives of all albumins which otherwise increased significantly in the case of bilirubin complexed with native and guanidinated albumin derivatives, except the bilirubin-RbSA complex which showed a small decrease in intensity. These results suggest that the presence of salt linkage(s) in bilirubin-albumin complexation is(are) crucial to bring about effective and efficient stereochemical changes in the bound pigment by co-binding of chloroform which seems to have at least one conserved binding site on these albumins that is shared with bilirubin.     

Zeta potential of yeast cells: application in cell immobilization

The zeta potential of Saccharomyces cerevisiae cells has been studied. Zeta potential was measured by microelectrophoresis with a Laser Zee meter model 500. Haploid (a and α) and diploid cells were tested and their zeta potential was found to be comparable. The influence of agents such as CaCl₂, NaCl, carboxymethylcellulose, Al₂(SO₄)₃ and cationic starch, on the zeta potential of cells has been studied. Moreover, zeta potential measurement has been used for improvement of cell immobilization by adhesion on sawdust.     

Preparation of succinylated yeast protein: Composition and solubility

A method for the succinylation of yeast Proteins during their isolation is described. Some factor affecting the extent of succinylation are described, e.g., heating of protein (80°C for 5 min) prior to dervatization reduced maximum succinylation of available ε-amino groups from 88 to 54%. Succinylation of amino groups progressively increased as the concentration of succinic anhydride was increased and there was a concomitant decrease in nucleic acid content of the precipitated protein concentrate from 12.5 to 5.3 for 0 and 88% succinylation, respectively, succinylation enhanced the solubility of the protein concentrates particularly below pH 7 and it decreased the isoelectric point pH. 4.5 to 4.0.     

Effect of environment on allosteric properties of acetylcholinesterase]

In the presence of organophosphorus inhibitors (OPI) AChE inhibition is initiated at a lower concentration of ACh; the plot reaction rate versus substrate concentration shows two maxima with a distinct minimum between them. It was shown that extremely mild conditions (short-term heating up to 50 degrees C; acidic or alkaline pH shift by 0.5 units; high concentrations of bivalent cations; erythrocyte storage) which do not affect substrate inhibition, remove this effect. The data obtained suggest that OPI effect is not directed to the site of AChE responsible for enzyme inhibition by ACh excess ("substrate inhibition site"), but to some other area. This results in a change in the conformation of the substrate inhibition site and a pronounced inhibition of the AChE activity takes place at lower substrate concentration.     

The effect of various acylating agents on the catalytic properties and structure of aspartate aminotransferase]

Changes of quaternary structure and conformation of molecule concomitant with inactivation were observed in the course of aspartate transaminase acylation by maleic, citraconic, dimethylmaleic and succinic anhydrides. It was established that acylation of 10-12 xi-amino groups of lysine did not induce the dissociation of transaminase into subunits. Further acylation of amino groups (2 groups if dimethylmaleic anhydrade was used as acylating agent) induced dissociation of transaminase dimer into subunits. These data were obtained by sedimentation analysis. The dissociation was accompanied with a sharp decrease of correlation time (from 18 nsec to 9 nsec) of the paramagnetic label covalently bound to the protein. The obtained results allow us to distinguish three types of xi-aminogroups of aspartate transaminase: exposed (about 12 residues), "contact" (2 residues) located in the vicinity to complementary surfaces of subunits and buried (about 6 residues). The stepwise inactivation occurred during the acylation as a result of conformational changes or appearance of sterical hindrances in the cataytic site of the enzyme. The thiol groups were not modified in transaminase molecule under experimental conditions used. Aspartate transaminase treated with citraconic or dimethylmaleic anhydride may be deacylated under mild conditions. After reacylation the quaternary structure was reconstituted and catalytic activity was almost fully restored.     

Identification and application of RAPD markers for anthracnose resistance in water yam (Dioscorea alata)

Anthracnose, caused by Colletotrichum gloeosporioides, is the most severe foliar disease of water yam (Dioscorea alata) worldwide. The tetraploid breeding line, TDa 95/00328, is a source of dominant genetic resistance to the moderately virulent fast growing salmon (FGS) strain of C. gloeosporioides. Bulked segregant analysis was used to search for random amplified polymorphic DNA (RAPD) markers linked to anthracnose resistance in F₁ progeny derived from a cross between TDa 95/00328 and the susceptible male parent, TDa 95–310. Two hundred and eighty decamer primers were screened using bulks obtained from pooled DNA of individuals comprising each extreme of the disease phenotype distribution. A single locus that contributes to anthracnose resistance in TDa 95/00328 was identified and tentatively named Dcg‐1. We found two RAPD markers closely linked in coupling phase with Dcg‐1, named OPI7₁₇₀₀ and OPE6₉₅₀, both of which were mapped on the same linkage group. OPI7₁₇₀₀ appeared tightly linked to the Dcg‐1 locus; it was present in all the 58 resistant F₁ individuals and absent in all but one of the 13 susceptible genotypes (genetic distance of 2.3 cM). OPE6₉₅₀ was present in 56 of the 58 resistant progeny and only one susceptible F₁ plant showed this marker (6.8 cM). Both markers successfully identified Dcg‐1 in resistant D. alata genotypes among 34 breeding lines, indicating their potential for use in marker‐assisted selection. OPI7₁₇₀₀ and OPE6₉₅₀ are the first DNA markers for yam anthracnose resistance. The use of molecular markers presents a valuable strategy for selection and pyramiding of anthracnose resistance genes in yam improvement.     

Effect of chain length on physicochemical properties and cytotoxicity of cationic vesicles composed of phosphatidylcholines and dialkyldimethylammonium bromides

This study investigated the physicochemical characteristics of cationic vesicles that were prepared from two phosphatidylcholines and three dialkyldimethylammonium bromides (DXDAB) with differing in dialkyl chain lengths, ranging from 2-C₁₄ to 2-C₁₈, by measuring particle size and zeta potential. The dependence of particle size, zeta potential and short-storage stability of mixed phosphatidylcholine/DXDAB vesicles on the chain length and composition were also elucidated. Transmission electron microscopy analysis verified that vesicles were formed as a phosphatidylcholine film to which DXDAB was added in a phosphate buffer saline (PBS, pH 7.4). Furthermore, the toxicity to the human keratinocytes (HaCaT) and squamous cell carcinomas (SCC25) cells that were incubated with these vesicles, evaluated by a cell viability assay, increased with the percentage of DXDAB that was incorporated and was inversely proportional to the chain length of DXDAB. The morphological features (round shape, chromatin condensation and apoptosis bodies) and results of flow cytometry analysis (increased sub-G₁ fraction) confirmed the induction of apoptosis in HaCaT and SCC25 cells by cationic vesicles. Apoptosis caused by cationic vesicles without the addition of any drugs was observed for the first time in HaCaT and SCC25 cells. The results of this investigation suggest that cytotoxicity is related to the zeta potential of the cationic vesicles.     

Variable sequences outside the SAM-binding core critically influence the conformational dynamics of the SAM-III/SMK box riboswitch

The S_MK box (SAM-III) translational riboswitches were identified in S-adenosyl-l-methionine (SAM) synthetase metK genes in members of Lactobacillales. This riboswitch switches between two alternative conformations in response to intracellular SAM concentration and controls metK expression at the level of translation initiation. We previously reported the crystal structure of the SAM-bound S_MK box riboswitch. In this study, we combined selective 2′-hydroxyl acylation analyzed by primer extension chemical probing with mutagenesis to probe the ligand-induced conformational switching mechanism. We revealed that while the majority of the apo S_MK box RNA molecules exist in an alternatively base-paired (ON) conformation, a subset of them pre-organize into a SAM-bound-like (READY) conformation, which, upon SAM exposure, is selectively stabilized into the SAM-bound (OFF) conformation through an induced-fit mechanism. Mutagenesis showed that the ON state is only slightly more stable than the READY state, as several single-nucleotide substitutions in a hypervariable region outside the SAM-binding core can alter the folding landscape to favor the READY state. Such S_MK variants display a “constitutively OFF” behavior both in vitro and in vivo. Time-resolved and temperature-dependent selective 2′-hydroxyl acylation analyzed by primer extension analyses revealed adaptation of the S_MK box RNA to its mesothermal working environment. The latter analysis revealed that the SAM-bound S_MK box RNA follows a two-step folding/unfolding process.     

Neocarzinostatin: effect of modification of side chain amino and carboxyl groups on chemical and biological properties.

The antitumor protein neocarzinostatin (NCS), isolated from Streptomyces carzinostaticus, is a single chain polypeptide with 109 amino acid residues. Complete acylation of the amino groups (alanine-1 and lysine-20) was observed when NCS was allowed to react with 3-(4-hydroxyphenyl)-propionic acid N-hydroxysuccinimide ester at pH 8.5. Since the ensuing bis[(alanine-1, lysine-20)-3-(4-hydroxyphenyl)]-propionamide NCS was fully active in antibacterial potency and in the inhibition of growth of leukemic (CCRF-CEM) cells in vitro, it appears that the two amino groups in the protein are not essential for biological activity. Radiolabeled NCS was prepared by using a tritiated or 125I-labeled acylating agent. Since the CD spectra of native and bis(alanine-1, lysine-20)-amino modified NCS were indistinguishable, there is presumably no change in the native conformation of the protein due to acylation. Reaction of NCS with ammonium chloride in the presence of 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide at pH 4.75 converted all the 10 carboxyl groups into carboxamides and produced a protein derivative of basic character. This modification caused a change in the native conformation of the protein accompanied by a loss in biological inhibitory activities.     

Aminoacylation and conformational properties of yeast mitochondrial tRNA mutants with respiratory deficiency

We report the identification and characterization of eight yeast mitochondrial tRNA mutants, located in mitochondrial tRNA(Gln), tRNA(Arg2), tRNA(Ile), tRNA(His), and tRNA(Cys), the respiratory phenotypes of which exhibit various degrees of deficiency. The mutations consist in single-base substitutions, insertions, or deletions, and are distributed all over the tRNA sequence and structure. To identify the features responsible for the defective phenotypes, we analyzed the effect of the different mutations on the electrophoretic mobility and efficiency of acylation of the mutated tRNAs in comparison with the respective wild-type molecules. Five of the studied mutations determine both conformational changes and defective acylation, while two have neither or limited effect. However, variations in structure and acylation are not necessarily correlated; the remaining mutation affects the tRNA conformation, but not its acylation properties. Analysis of tRNA structures and of mitochondrial and cytoplasmic yeast tRNA sequences allowed us to propose explanations for the observed defects, which can be ascribed to either the loss of identity nucleotides or, more often, of specific secondary and/or tertiary interactions that are largely conserved in native mitochondrial and cytoplasmic tRNAs.     

Enzymatic acylation of polar dipeptides: Influence of reaction media and molecular environment of functional groups

The enzymatic acylation of polar dipeptides was investigated. First, the Novozym435^®-catalyzed acylation of Lys-Ser, HCl exhibiting three potential acylable sites was carried out in organic media (2-methyl-2-butanol, oleic acid) and in an ionic liquid ([Bmim]⁺[PF₆]^?). In these reactions, the chemo-selectivity of the acylation was exclusively in favour of the N?-lysine acylation and the efficiency (substrate conversion) was demonstrated to be under control of the peptide solubility. The use of [Bmim]⁺[PF₆]^? permitted to significantly improve the dipeptide solubility, and to enhance both substrates conversion and initial rates of acylation reaction. In the three reaction media used, the O-acylated derivative of the dipeptide was never detected suggesting a weak reactivity of the serine hydroxyl group due to its molecular environment and particularly to the presence of a free carboxylic group known for its electroattractor property.Last, the acylation of a natural dipeptide (carnosine), exhibiting a very low solubility in organic solvents, was also performed. Carnosine was successfully N-acylated in 2-methyl-2-butanol, and a yield of 39% was obtained when improving the substrate solubility: a better dispersibility was obtained by application of a high pressure on the reaction medium just before starting the reaction.     

On the modulation of photoinduced fluorescence enhancement and conformational stability of albumin-bound bilirubin: effect of epsilon-NH(2) groups blocking and chloroform binding

Photoinduced fluorescence enhancement of bilirubin bound to primary binding site on human serum albumin (HSA) was completely ceased when epsilon-NH(2) groups of its internal lysine residues were covalently blocked by acetylation or succinylation though the pigment bound to these derivatives in a folded conformation akin to that bound to HSA. These photoinduced fluorescence modulations cannot be ascribed to the binding of bilirubin to secondary low affinity sites as the CD spectrum of bilirubin bound to these derivatives showed complete inversion upon addition of chloroform which binds to subdomain IIA in HSA where high affinity bilirubin binding site is located. Presence of chloroform reconciled the photoinduced alterations in the CD spectrum observed in its absence, suggesting that chloroform stabilized the bound ligand against light but the fluorescence properties of bilirubin complexed with acetylated or succinylated derivatives remained unchanged. Guanidination of internal epsilon-NH(2) groups in HSA by O-methylisourea did not alter the spectral properties of the bound ligand. These results suggest that salt linkage(s) existing between epsilon-NH(2) groups of lysine residues in HSA and carboxyl groups of bilirubin, act(s) as a potential barrier during conformational rotation of the bound ligand assisted by photoactivation and their abolishment can alter its dynamics and stereoselectivity, a hitherto unnoticed implication of salt linkage(s) in BR-HSA complex.