Cassette mutagenesis of Met121 in azurin from Pseudomonas aeruginosa

Cassette mutagenesis was used to exchange the suggested copperligand Met121 in azurin to all other amino acids, and a stopcodon. The mutant proteins were characterized by optical absorptionspectroscopy and EPR. At low pH, all mutants exhibit the characteristicsof a blue type 1 copper protein, indicating that methionineis not needed to create a blue copper site. At high pH, theGlu121 and the Lys121 mutants constitute a new form of protein-boundcopper that is outside the range of type1 copper.     

Molecular dynamics simulation of winter flounder antifreeze protein variants in solution: correlation between side chain spacing and ice lattice

The solution structure of the 38 amino acid C-terminal regionof the precursor for the HPLC-6 antifreeze protein from winterflounder has been investigated with molecular dynamics usingthe AMBER software. The simulation for the peptide in aqueoussolution was carried out at a constant temperature of 0°Cand at atmospheric pressure. The simulation covered 120 ps andthe results were analyzed based on data sampled upon reachinga stable equilibrium phase. Information has been obtained onthe quality of constant temperature and pressure simulations,the solution structure and dynamics, the hydrogen bonding network,the helix-stabilizing role of terminal charges and the interactionwith the surrounding water molecules. The Lys18–Glu22interactions and the terminal charged residues are found tostabilize a helical structure with the side chains of Thr2,Thr13, Thr24 and Thr35 equally spaced on one side of the helix.The spacing between oxygen atoms in the hydroxyl group of thethreonine side chains exhibits fluctuations of the order of2–3 Å during the 120 ps of simulation, but valuessimultaneously close to the repeat distance of 16.6 Åbetween oxygen atoms along the [0112] direction in ice are observed.Furthermore, two engineered variants were studied using thesame simulation protocol.     

Synthesis of a squash-type protease inhibitor by gene engineering and effects of replacements of conserved hydrophobic amino acid residues on its inhibitory activity

Cucurbita maxima trypsin inhibitor I (CMTI-I), a member of thesquash-type protease inhibitor family, is composed of 29 aminoacids and shows strong inhibition of trypsin by its compactstructure. To study the structure–function relationshipof this inhibitor using protein engineering methods, we constructedan expression system for CMTI-I as a fused protein with porcineadenylate kinase (ADK). A Met residue was introduced into thejunction of ADK and CMTI-I to cleave the fusion protein withCNBr, whereas a Met at position 8 of authentic CMTI-I was replacedby Leu. Escherichia coli JM109 transformed with the constructedplasmid expressed the fused protein as an inclusion body. Aftercleavage of the expressed protein with CNBr, fully reduced speciesof CMTI-I were purified by reversed-phase HPLC and then oxidizedwith air by shaking. For efficient refolding of CMTI-I, we used50 mM NH₄HCO₃ (pH 7.8) containing 0.1% PEG 6000 at higher proteinconcentration. Strong inhibitory activity toward trypsin wasdetected only in the first of three HPLC peaks. The inhibitorconstant of CMTI-I thus obtained, in which Met8 was replacedby Leu, was 1.4x10-¹⁰ M. The effect of replacement of Met withLeu at position 8 was shown to be small by comparison of theinhibitor constant of authentic CMTI-III bearing Lys at position9 (8.9x10-¹¹ M) with that of its mutant bearing Leu at position8 and Lys at position 9 (1.8x10-¹⁰ M). To investigate the roleof the well conserved hydrophobic residues of CMTI-I in itsinteraction with trypsin, CMTI-I mutants in which one or allof the four hydrophobic residues were replaced by Ala were prepared.The inhibitor constants of these mutants indicated that thosewith single replacements were 5–40 times less effectiveas trypsin inhibitors and that the quadruple mutant was –450times less effective, suggesting that the hydrophobic residuesin CMTI-I contribute to its tight binding with trypsin. However,each mutant was not converted to a temporary inhibitor.     

Amino acid composition of protein termini are biased in different manners

An exhaustive statistical analysis of the amino acid sequencesat the carboxyl (C) and amino (N) termini of proteins and ofcoding nucleic acid sequences at the 5' side of the stop codonswas undertaken. At the N ends, Met and Ala residues are over-representedat the first (+1) position whereas at positions 2 and 5 Thris preferred. These peculiarities at N-termini are most probablyrelated to the mechanism of initiation of translation (for Met)and to the mechanisms governing the life-span of proteins viaregulation of their degradation (for Ala and Thr). We assumethat the C-terminal bias facilitates fixation of the C endson the protein globule by a preference for charged and Cys residues.The terminal biases, a novel feature of protein structure, haveto be taken into account when molecular evolution, three-dimensionalstructure, initiation and termination of translation, proteinfolding and life-span are concerned. In addition, the bias ofprotein termini composition is an important feature which shouldbe considered in protein engineering experiments.     

Can three-dimensional contacts in protein structures be predicted by analysis of correlated mutations?

A method has been developed to detect pairs of positions withcorrelated mutations in protein multiple sequence alignments.The method is based on reconstruction of the phylogenetic treefor a set of sequences and statistical analysis of the distributionof mutations in the branches of the tree. The database of homology-derivedprotein structures (HSSP) is used as the source of multiplesequence alignments for proteins of known three-dimensionalstructure. We analyse pairs of positions with correlated mutationsin 67 protein families and show quantitatively that the presenceof such positions is a typical feature of protein families.A significant but weak tendency is observed for correlated residuepairs to be close in the three-dimensional structure. With furtherimprovements, methods of this type may be useful for the predictionof residue-residue contacts and subsequent prediction of proteinstructure using distance geometry algorithms. In conclusion,we suggest a new experimental approach to protein structuredetermination in which selection of functional mutants afterrandom mutagenesis and analysis of correlated mutations providesufficient proximity constraints for calculation of the proteinfold     

Modulating D-amino acid oxidase substrate specificity: production of an enzyme for analytical determination of all D-amino acids by directed evolution

Recent research on the flavoenzyme D-amino acid oxidase from Rhodotorula gracilis (RgDAAO) has revealed new, intriguing properties of this catalyst and offers novel biotechnological applications. Among them, the reaction of RgDAAO has been exploited in the analytical determination of the D-amino acid content in biological samples. However, because the enzyme does not oxidize acidic D-amino acids, it cannot be used to detect the total amount of D-amino acids. We now present the results obtained using a random mutagenesis approach to produce RgDAAO mutants with a broader substrate specificity. The libraries of RgDAAO mutants were generated by error-prone PCR, expressed in BL21(DE3)pLysS Escherichia coli cells and screened for their ability to oxidize different substrates by means of an activity assay. Five random mutants that have a 'modified' substrate specificity, more useful for the analytical determination of the entire content of D-amino acids than wild-type RgDAAO, have been isolated. With the only exception of Y223 and G199, none of the effective amino acid substitutions lie in segments predicted to interact directly with the bound substrate. The substitutions appear to cluster on the protein surface: it would not have been possible to predict that these substitutions would enhance DAAO activity. We can only conclude that these substitutions synergistically generate small structural changes that affect the dynamics and/or stability of the protein in a way that enhances substrate binding or subsequently catalytic turnover.     

Fluorescence resonance energy transfer analysis of the folding pathway of Engrailed Homeodomain

The Engrailed Homeodomain folds on the microsecond time scale via an intermediate that is experimentally well characterised using structural Engrailed-Homeodomain mimics. Here, we analysed directly the changes in distance between key residues during the kinetics of unfolding and at equilibrium using fluorescence resonance energy transfer (FRET). Trp was the donor and 5-(((acetylamino)ethyl)amino) naphthalene-1-sulphate, the acceptor, substituted in positions that caused little change in stability. Distances calculated for the native state were in good agreement with those derived from the NMR structure. The distances between the N- and C-termini of Helix I and of Helix III increased, then decreased and finally increased again with increasing GdmCl concentration on equilibrium denaturation. This behaviour implied that there was a folding intermediate on the folding pathway and that this intermediate was populated at low concentrations of GdmCl concentration ( approximately 1 M). We analysed the changes in distance during temperature-jump relaxation kinetics, using a qualitative and very conservative procedure that drew conclusions only when changes in fluorescence of mutants containing either the donor or the acceptor alone would not obscure the change in the FRET signal when both donor and acceptor were present. The distance changes obtained under equilibrium and kinetic measurements were self-consistent and also consistent with the known high-resolution structures of the mimics of the folding intermediates. We showed that for analysing distances in disordered ensembles, it is important to use FRET probes with a critical distance close to the average separation in the ensemble. Otherwise, average distances could be over or underestimated.     

Comparison of systematic search and database methods for constructing segments of protein structure

Two principal methods of determining the conformation of shortpieces of polypeptide backbone in proteins have been developed:using a database of known structures and systematicallygeneratingall conformations. In this paper, we compare the effectivenessof these two techniques. The completeness of the database forsegments of different lengths is examined and it is found tocontain most conformations for segments seven residues long,but to deteriorate rapidly for longer regions. When the databasesegment is to be incorporated into the rest of a structure,at least seven residues are required to build four new residues,because of the need to positionthe segment relative to the restof the structure.It is found that such positioning using flankingresidues results in large errors in the inserted region. Weconclude that the database method is currently not effectivefor comparative modeling, even for short segments. The systematicsearchprocedure is found to generate almost all structures of shortsegments found in proteinsIn contrast to the database method,low root mean square error structures are obtained for a setof trial segments embedded in the rest of a protein structure.Thus, it should be considered the method of choice.     

Protein engineering loops in aspartic proteinases: site-directed mutagenesis, biochemical characterization and X-ray analysis of chymosin with a replaced loop from rhizopuspepsin

The loop exchange mutant chymosm 155–164 rhizopuspepsinwas expressed in Trichoderma reesei and exported into the mediumto yield a correctly folded and active product. The biochemicalcharacterization and crystal structure determination at 2.5Å resolution confirm that the mutant enzyme adopts a nativefold. However, the conformation of the mutated loop is unlikethat in native rhizopuspepsin and involves the chelation ofa water molecule in the loop. Kinetic analysis using two syntheticpeptide substrates (six and 15 residues long) and the naturalsubstrate, milk, revealed a reduction in the activity of themutant enzyme with respect to the native when acting on boththe long peptide substrate and milk. This may be a consequenceof the different charge distribution of the mutated loop, itsincreased size and/or its different conformation.     

Antagonism, non-native interactions and non-two-state folding in S6 revealed by double-mutant cycle analysis

When folding to the native state N in the presence of salt, the apparent two-state folder S6 transiently forms a transient off-pathway state C with substantial secondary and tertiary structure. Fifteen double mutant cycles were analysed to compare side-chain interaction energies DeltaDeltaG(int) in C, N and TS (the transition state between N and the denatured state). The kinetic signatures of these destabilizing mutants suggest folding scenarios involving unfolding intermediates and even alternative unfolding pathways. However, restricting the kinetic data to linear parts of the chevron plot allows reliable extrapolation to zero molar denaturant of rate constants of folding, unfolding and misfolding. Side-chain interactions appear to contribute to the stability of C, but in a substantially non-native environment, as shown by changes in the sign of DeltaDeltaG(int) between C and N. Remarkably, there appear to be significant (0.7-2 kcal/mol) antagonistic interactions between the two residues Leu30 and Leu75 in N and TS, which may be linked to subtle structural changes seen in the crystal structures of the mutants. A small number of overlapping residues are involved in these kinds of antagonistic interactions in N, TS and C, suggesting that repulsive interactions are coded into the protein topology whether the protein folds or misfolds. Destabilizing double mutants indicate that apparent two-state folders can be induced to behave in more complex ways provided that the native state is suitably destabilized.     

The PASTA server for protein aggregation prediction

Many different proteins aggregate into amyloid fibrils characterized by cross-beta structure. beta-strands contributed by distinct protein molecules are generally found in a parallel in-register alignment. Here, we describe the web server for a novel algorithm, prediction of amyloid structure aggregation (PASTA), to predict the most aggregation-prone portions and the corresponding beta-strand inter-molecular pairing for a given input sequence. PASTA was previously shown to yield results in excellent agreement with available experimental observations, when tested on both natively unfolded and structured proteins. The web server and downloadable source code are freely accessible from the URL: http://protein.cribi.unipd.it/pasta/.     

Secondary structure prediction for modelling by homology

An improved method of secondary structure prediction has beendeveloped to aid the modelling of proteins by homology. Selecteddata from four published algorithms are scaled and combinedas a weighted mean to produce consensus algorithms. Each consensusalgorithm is used to predict the secondary structure of a proteinhomologous to the target protein and of known structure. Bycomparison of the predictions to the known structure, accuracyvalues are calculated and a consensus algorithm chosen as theoptimum combination of the composite data for prediction ofthe homologous protein. This customized algorithm is then usedto predict the secondary structure of the unknown protein. Inthis manner the secondary structure prediction is initiallytuned to the required protein family before prediction of thetarget protein. The method improves statistical secondary structureprediction and can be incorporated into more comprehensive systemssuch as those involving consensus prediction from multiple sequencealignments. Thirty one proteins from five families were usedto compare the new method to that of Garnier, Osguthorpe andRobson (GOR) and sequence alignment. The improvement over GORis naturally dependent on the similarity of the homologous protein,varying from a mean of 3% to 7% with increasing alignment significancescore.     

Evolutionary divergence and conservation of trypsin

The trypsin sequences currently available in the data bankshave been collected and aligned using first the amino acid sequencehomology and, subsequently, the superposed crystal structuresof trypsins from the cow, the bacterium Streptomyces griseusand the fungus Fusarium oxysporum. The phylogenetic tree constructedaccording to this multiple alignment is consistent with a continuousevolutionary divergence of trypsin from a common ancestor ofboth prokaryotes and eukaryotes. Comparison of crystal structuresreveals a strict conservation of secondary structure. Similarly,in the alignment of all the sequences, insertions and deletionsoccur only in regions corresponding to loops between the secondarystructure elements in the known crystal structures. The conservedresidues cluster around the active site. Almost all conservedresidues can be associated with one of the basic functionalfeatures of the protein: zymogen activation, catalysis and substratespecificity. In contrast, the residues of the hydrophobic coreof the protein and the calcium ion binding sites are generallynot conserved. The conserved features of trypsin and the natureof the conservation are discussed In detail     

Spectral tuning of photoproteins by partnering site-directed mutagenesis strategies with the incorporation of chromophore analogs

Aequorin and obelin are photoproteins whose calcium controlled bioluminescent light emission is used for labeling in assays, for the determination of calcium concentrations in vivo, and as a reporter in cellular imaging. Both of these photoproteins emit blue light from a 2-hydroperoxycoelenterazine chromophore, which is non-covalently bound in the hydrophobic core of the proteins. In an effort to produce aequorin and obelin variants with improved analytical properties, such as alternative emission colors and altered decay kinetics, seven mutants of aequorin and obelin were prepared and combined with 10 different coelenterazine analogs. These semi-synthetic photoprotein mutants exhibited shifts in bioluminescent properties when compared with wild-type proteins. The bioluminescent parameters determined for these semi-synthetic photoprotein mutants included specific activity, emission spectra and decay half-life time. This spectral tuning strategy resulted in semi-synthetic photoprotein mutants that had significantly altered bioluminescent properties. The largest emission maxima shift obtained was 44 nm, and the largest decay half-life difference was 23.91 s.     

Study of cysteine residues in the {alpha} subunit of Escherichia coli tryptophan synthase. 1. Role in conformational stability

To eludicate the role in conformational stability of Cys residuesburied in the interior of a protein, the thermodynamic propertiesof denaturation of mutant subunit of Escherichia coli tryptophansynthase, in which Ser, Ala, Val or Gly was substituted foreach of the three Cys residues, were analyzed using calorimetry.The mutants were less stable than the wild type, indicatingthat Cys residues contribute greatly to the stability of the subunit. In most cases, a large decrease in denaturation enthalpywas observed, compensated for by the denaturation entropy toa major extent. The extent of changes in the denaturation Gibbsenergy and denaturation enthalpy varied greatly depending onboth substituting residues and positions. Of all the mutantproteins, the Cys154Ser mutant showed the greatest decreasein denaturation enthalpy; its denaturation enthalpy was halfthat of the wild type, and was considerably repaired by addinga ligand of the subunit. Because the enthalpy of ligand bindingto Cys154Ser in the native state did not change. it seems thatthe decrease in the denaturation enthalpy of Cys154Ser and itsrecovery by ligand binding are caused by conformational changesin the denatured state due to the mutation.     

Protein engineering of the high-alkaline serine protease PB92 from Bacillus alcalophilus: functional and structural consequences of mutation at the S4 substrate binding pocket

Serine endoproteases such as trypsins and subtilisins are knownto have an extended substrate binding region that interactswith residues P6 to P3' of a substrate. In order to investigatethe structural and functional effects of replacing residuesat the S4 substrate binding pocket, the serine protease fromthe alkalophilic Bacillus strain PB92, which shows homologywith the subtilisins, was mutated at positions 102 and 126–128.Substitution of Val102 by Trp results in a 12–fold increasein activity towards succinyl-L-Ala-L-Ala-L-Pro-L-Phe-p-nitroanilide(sAAPFpNA). An X-ray structure analysis of the V102W mutantshows that the Trp side chain occupies a hydrophobic pocketat the surface of the molecule leaving a narrow crevice forthe P4 residue of a substrate. Better binding of sAAPFpNA bythe mutant compared with the wild type protein as indicatedby the kinetic data might be due to the hydrophobic interactionof Ala P4 of the substrate with the introduced Trp102 side chain.The observed difference in binding of sAAPFpNA by protease PB92and thermitase, both of which possess a Trp at position 102,is probably related to the amino acid substitutions at positions105 and 126 (in the protease PB92 numbering).Kinetic data forthe variants obtained by random mutation of residues Serl26,Prol27 and Serl28 reveal that the activity towards sAAPFpNAincreases when a hydrophobic residue is introduced at position126. An X-ray diffraction analysis was carried out for the threeprotease PB92 mutants which have residues Serl26-Prol27-Serl28replaced by Met-Ala-Gly(‘MAG’ mutant), Phe-Gln-Ser(‘FQS’ mutant) and Asn-Ser-Ala (‘NSA’mutant). Met 126 and Phel26 in the crystal structures of thecorresponding mutants are fixed in the same hydrophobic environmentas Trp102 in the V102W mutant.In contrast, Asnl26 in the ‘NSA’mutant is completely disordered in both crystal forms for whichthe structure has been determined. According to our kineticmeasurements none of the mutants with Met, Phe, Leu or Val atposition 126 binds sAAPFpNA better than the wild type enzyme.Resultsof the site-directed mutagenesis at position 127 imply thatpossible interaction of this residue with a substrate has almostno effect on activity towards sAAPFpNA and casein.     

Increasing the synthetic performance of penicillin acylase PAS2 by structure-inspired semi-random mutagenesis

A semi-random mutagenesis approach was followed to increase the performance of penicillin acylase PAS2 in the kinetically controlled synthesis of ampicillin from 6-aminopenicillanic acid (6-APA) and activated D-phenylglycine derivatives. We directed changes in amino acid residues to positions close to the active site that are expected to affect the catalytic performance of penicillin acylase: alpha R160, alpha F161 and beta F24. From the resulting triple mutant gene bank, six improved PAS2 mutants were recovered by screening only 700 active mutants with an HPLC-based screening method. A detailed kinetic analysis of the three most promising mutants, T23, TM33 and TM38, is presented. These mutants allowed the accumulation of ampicillin at 4-5 times higher concentrations than the wild-type enzyme, using D-phenylglycine methyl ester as the acyl donor. At the same time, the loss of activated acyl donor due to the competitive hydrolytic side reactions could be reduced to <20% with the mutant enzymes compared >80% wild-type PAS2. Although catalytic activity dropped by a factor of 5-10, the enhanced synthetic performance of the recovered penicillin acylase variants makes them interesting biocatalysts for the production of beta-lactam antibiotics.     

Structural basis for the extreme thermostability of D-glyceraldehyde-3-phosphate dehydrogenase from Thermotoga maritima: analysis based on homology modelling

D-Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from a hyperthermophiliceubacterium, Thermotoga maritima, is remarkably heat stable(T_m = 109°C). In this work, we have applied homology modellingto predict the 3-D structure of Th.maritima GAPDH to revealthe structural basis of thermostability. Three known GAPDH structureswere used as reference proteins. First, the rough model of onesubunit was constructed using the identified structurally conservedand variable regions of the reference proteins. The holoenzymewas assembled from four subunits and the NAD molecules. Thestructure was refined by energy minimization and molecular dynamicssimulated annealing. No errors were detected in the refinedmodel using the 3-D profile method. The model was compared withthe structure of Bacillus stearothermophilus GAPDH to identifystructural details underlying the increased thermostability.In all, 12 extra ion pairs per subunit were found at the proteinsurface. This seems to be the most important factor responsiblefor thermostability. Differences in the non-specific interactions,including hydration effects, were also found. Minor changeswere detected in the secondary structure. The model predictsthat a slight increase in a-helical propensities and helix-dipoleinteractions also contribute to increased stability, but toa lesser degree.     

Structural and kinetic studies on ligand binding in wild-type and active-site mutants of penicillin acylase

Penicillin acylase catalyses the condensation of Calpha-substituted phenylacetic acids with beta-lactam nucleophiles, producing semi-synthetic beta-lactam antibiotics. For efficient synthesis a low affinity for phenylacetic acid and a high affinity for Calpha-substituted phenylacetic acid derivatives is desirable. We made three active site mutants, alphaF146Y, betaF24A and alphaF146Y/betaF24A, which all had a 2- to 10-fold higher affinity for Calpha-substituted compounds than wild-type enzyme. In addition, betaF24A had a 20-fold reduced affinity for phenylacetic acid. The molecular basis of the improved properties was investigated by X-ray crystallography. These studies showed that the higher affinity of alphaF146Y for (R)-alpha-methylphenylacetic acid can be explained by van der Waals interactions between alphaY146:OH and the Calpha-substituent. The betaF24A mutation causes an opening of the phenylacetic acid binding site. Only (R)-alpha-methylphenylacetic acid, but not phenylacetic acid, induces a conformation with the ligand tightly bound, explaining the weak binding of phenylacetic acid. A comparison of the betaF24A structure with other open conformations of penicillin acylase showed that betaF24 has a fixed position, whereas alphaF146 acts as a flexible lid on the binding site and reorients its position to achieve optimal substrate binding.     

L3 loop-mediated mechanisms of pore closing in porin: a molecular dynamics perturbation approach

L3 loop-mediated mechanisms for pore closing in porin are investigatedwith molecular dynamics simulation, using an approach that canbe related to the phenomenon of voltage gating. Voltage gatingis seen as a perturbation of the electrostatic screening insidethe porin pore where, by the influence of the potential gradient,water and counter-ion distribution can be slightly displacedfrom their equilibrium distribution. This is simulated by perturbingthe screening electrostatics of ionizable groups inside thepore. Under these conditions, a localized conformational changetakes place, involving 12 (Ile102–Alall3) out of the 44residues of the loop. The pore is reduced to a sixth of itsopen state size. The conformational change can be achieved witha small perturbation and it is reversible once the perturbationis switched off (relaxation process). Other types of behaviourpredominating at higher simulation temperatures are found forthe loop, involving an extra conformational change in the Thr92–Asp96loop segment. This conformational change completely closes thepore, but is not reversible under the simulation conditions.Both zones involved in the conformational changes contain oroverlap the zones which were described previously, using othertechniques, to be the most flexible zones of the loop.