Molecular dynamics simulations of the whey protein {beta}-lactoglobulin

Molecular dynamics simulations have been used to model the motionsand conformational behavior of the whey protein bovine ß-lactoglobulin.Simulations were performed for the protein by itself and complexedto a single retinol ligand located in a putative interior bindingpocket. In the absence of the retinol ligand, the backbone loopsaround the opening of this ulterior pocket shifted inward topartially close off this cavity, similar to the shifts observedin previously reported molecular dynamics simulations of theuncomplexed form of the homologous retinol binding protein.The protein complexed with retinol does not exhibit the sameconformational shifts. Conformational changes of this type couldserve as a recognition signal allowing in vivo discriminationbetween the free and retinol complexed forms of the 3-lactoglobulinmolecule. The unusual bending of the single a-helix observedin the simulations of retinol binding protein were not observedin the present calculations     

Thermostable variants of bovine {beta}-lactoglobulin

The thermal stability of bovine ß-lactoglobulin (BLG)has been enhanced by the introduction of an additional disulfidebond. Wild-type BLG has two disulfide bonds, C106–C119and C66–C160, with a free cysteine at position 121. Wehave designed, with the aid of molecular modeling calculations,two mutants of a recombinant BLG (rBLG), L104C and A132C. Moleculardynamics simulations were performed at 300K to study the effectof these alterations on the conformation of the protein. Thesemutants were then created by site-directed mutagenesis and purifiedfrom Escherichia coli carrying a tac expression vector usinga two-step renaturation method. Formation of disulfide linkagesin the correct arrangement, as designed, was confirmed by peptidemapping. In contrast to wild-type rBLG, which polymerizes attemperatures >65°C, neither of the mutant proteins polymerized.The conformational stability of the L104C and A132C mutant proteinsagainst thermal denaturation has been substantially increased(8- 10°C) as compared with wild-type rBLG. Furthermore,the A132C rBLG exhibits an enhanced stability against denaturationby guanidine hydrocnloride as compared with the wild-type orL104C rBLG     

Insight into the mechanism of the IMP-1 metallo-{beta}-lactamase by molecular dynamics simulations

Two models, a purely nonbonded model and a cationic dummy atomapproach, were examined for the modeling of the binuclear zinc-containingIMP-1 metallo-ß-lactamase in complex with a mercaptocarboxylateinhibitor. The cationic dummy atom approach had substantialadvantages as it maintained the initial, experimentally determinedgeometry of the metal-containing active site during moleculardynamics simulations in water. The method was extended to themodeling of the free enzyme and the enzyme in complex with acephalosporin substrate docked in an intermediate structure.For all three systems, the modeled complexes and the tetrahedralcoordination of the zinc ions were stable. The average zinc–zincdistance increased by     

Molecular dynamics simulations of the protein unfolding/folding reaction

All-atom molecular dynamics simulations of proteins in solvent are now able to realistically map the protein-unfolding pathway. The agreement with experiments probing both folding and unfolding suggests that these simulated unfolding events also shed light on folding. The simulations have produced detailed models of protein folding transition, intermediate, and denatured states that are in both qualitative and quantitative agreement with experiment. The various studies presented here highlight how such simulations both complement and extend experiment.     

Characterization of glycosylated variants of {beta}-lactoglobulin expressed in Pichia pastoris

Glycosylated variants of ß-lactoglobulin (BLG) wereproduced in the methylotrophic yeast Pichia pastoris to mimicthe glycosylation pattern of glycodelin, a homologue of BLGfound in humans. Glycodelin has three sites for glycosylation,corresponding to amino acids 63–65 (S1), 85–87 (S2)and 28–30 (S3) of BLG. These three sites were engineeredinto BLG to produce the variants S2, S12 and S123, which carriedone, two and three glycosylation sites, respectively. The oligosaccharideson these BLG variants ranged from (mannose)₉(N-acetylglucosamine)₂(Man₉GN₂) to Man₁₅GN₂ and were of the     

The elusive role of the N-terminal extension of {beta}A3- and {beta}Al-crystallin

ß-Crystallins are structural lens proteins with aconserved two-domain structure and variable N- and C-terminalextensions. These extensions are assumed to be involved in quaternaryinteractions within the ß-crystallin oligomers orwith other lens proteins. Therefore, the production of ßA3-and ßAl-crystallin from the single ßA3/A1mRNA by dual translation initiation is of interest. These crystallinsare identical, except that ßAl has a much shorterN-terminal extension than ßA3. This rare mechanismhas been conserved for over 250 million years during the evolutionof the ßA3/A1 gene, suggesting that the generationof different N-terminal extensions confers a selective advantage.We therefore compared the stability and association behaviourof recombinant ßA3- and ßAl-crystallin.Both proteins are equally stable in urea- and pH-induced denaturationexperiments. Gel filtration and analytical ultracentrifugationestablished that ßA3 and ßA1 both form homodimers.In the water-soluble proteins of bovine lens, ßA3and ßA1 are present in the same molecular weight fractions,indicating that they oligomerize equally with other ß-crystallins.¹H-NMR spectroscopy showed that residues Met1 to Asn22 of theN-terminal extension of ßA3 have great flexibilityand are solvent exposed, excluding them from protein interactionsin the homodimer. These results indicate that the differentN-terminal extensions of ßA3 and ßA1 donot affect their homo- or heteromeric interactions.     

Cloning of rabbit interleukin-1{beta}: differential evolution of IL-1{alpha} and IL-1{beta} proteins

We have cloned the rabbit IL-1ß cDNA, which encodesa 268 amino acid precursor similar in length to other sequencedIL-1 precursors. Comparison of all published IL-1 and IL-1ßsequences respectively indicates that the IL-1 gene family isevolving faster than the IL-1ß family, and that thetwo genes diverged –270 million years ago. Surprisingly,there are differences in the regions preferentially conservedwithin the two families. The IL-1 family is most conserved atthe amino terminus whereas the IL-1ß family is mostconserved in the carboxy-terminal half. This is despite thefact that the carboxy-terminal half encodes the active portionof both molecules and would be expected to adopt a similar ß-sheetstructure in IL-1 as in the published X-ray structure of matureIL-1ß. These findings suggest that differences inthe function and properties of the IL-1 and IL-1ßprecursor molecules may have been conserved. These differencesmay therefore provide an explanation for the existence of twoIL-1 molecules.     

Estimating the twist of {beta}-strands embedded within a regular parallel {beta}-barrel structure

The parallel ß-barrel is a recurrent structural motiffound in a large variety of different enzymes belonging to thefamily of /ß-proteins. It has been shown previouslythat the hyperboloid can be considered as a scaffold describingthe parallel ß-barrel structure. To assess restraintson ß-strand twist imposed by a given scaffold geometry,the notion of scaffold twist, T_s, is introduced. From T_s, theß-strand twist (Tw_ß) expected for a givenscaffold geometry can be derived and it is verified that thiscomputed twist can be used to identify ß-barrels characterizedby good hydrogen bonding. It is noted that Tw_ß isonly slightly affected for ß-barrels differing inthe number (N) of ß-strands, suggesting that restraintson main-chain conformation of ß-strands are not likelyto account for the N = 8 invariability observed in natural parallelß-barrels thereby strengthening previous work rationalizingthis constancy.     

Molecular dynamics simulations of phospholipases A2

An extensive molecular dynamics study of phospholipases A₂ frompancreatic bovine and Crotalus atrox venom has shown that thewell-conserved homologous core of the phospholipases A₂, includingthe so called catalytic network, is very stable during the courseof the calculations. The fluctuations which occur are locatedin segments which have significantly different three-dimensionalconformations in the two phospholipases A₂ studied, suggestingthat a particularly stable core conformation gives rise to alarge homologous family of similar three-dimensional structure.The calcium ion, which exhibits a crucial structural role inthe monomeric phospholipases A₂, appears not to be requiredto stabilize the C.atrox dimer. Moreover, the behaviour of thedimeric structure during the dynamics raises the question ofa possible dissociation of the two subunits into functionalmonomers.     

Automating the identification and analysis of protein {beta}-barrels

ßBarrels are widespread and well-studied featuresof a great many protein structures. In this paper an unsuper-visedmethod for the detection of P-barrels is developed based ontechniques from graph theory. The hydrogen bonded connectivityof ß-sheets is derived using standard pattern recognitiontechniques and expressed as a graph. Barrels correspond to topologicalrings in these connectivity graphs and can thus be identifiedusing ring perception algorithms. Following from this, the characteristictopological structure of a barrel can be expressed using a novelform of reduced nomenclature that counts sequence separationsbetween successive members of the ring set These techniquesare tested by applying them to the detection of barrels in anon-redundant subset of the Brookhaven database. Results indicatethat topological rings do seem to correspond uniquely to ß-barrelsand that the technique, as implemented, finds the majority ofbarrels present in the dataset.     

Molecular modelling of Staphylococcal {delta}-toxin ion channels by restrained molecular dynamics

-Toxin is a 26-residue channel-forming peptide from Staphylococcusaureus which forms an amphipathic -helix in a membrane environment.Channel formation in planar bilayers suggests that an averageof six -toxin helices self-assemble to form transbilayer pores.Molecular models for channels formed by -toxin and by a syntheticanalogue have been generated using a simulated annealing protocolapplied via restrained molecular dynamics. These models areanalysed in terms of the predicted geometric and energetic propertiesof the transbilayer pores. Pore radius calculations of the modelsdemonstrate that rings of channel-lining residues contributea series of constrictions along the pore. Electrostatic propertiesof the pores are determined both by pore-lining charged sidechains and by the aligned helix dipoles of the parallel helixbundle. Molecular dynamics simulations (100 ps) of -toxin modelscontaining intra-pore water were performed. Analysis of theresultant dynamics trajectories further supports the proposalthat alternative conformations of pore-constricting side chainsmay be responsible for the observed conductance heterogeneityof -toxin ion channels.     

Display of {beta}-lactamase on the Escherichia coli surface: outer membrane phenotypes conferred by Lpp'-OmpA'-{beta}-lactamase fusions

Bacterial cell-surface exposure of foreign peptides and solubleproteins has been achieved recently by employing a fusion proteinmethodology. An Lpp'–OmpA(46–159)–Bla fusionprotein has been shown previously to display the normally periplasmicenzyme ß-lactamase (Bla) on the cell surface of theGram-negative bacterium Escherichia coli. Here, we have investigatedthe role of the OmpA domain of the tripartite fusion proteinin the surface display of the passenger domain (Bla) and havecharacterized the effects of the fusion proteins on the integrityand permeability of the outer membrane. We show that in additionto OmpA(46–159), a second OmpA segment, consisting ofamino acids 46–66, can also mediate the display of Blaon the cell surface. Other OmpA domains of various lengths (aminoacids 46–84, 46–109, 46–128, 46–141and 46–145) either anchored the Bla domain on the periplasmicface of the outer membrane or caused a major disruption of theouter membrane, allowing the penetration of antibodies intothe cell. Detergent and antibiotic sensitivity and periplasmicleakage assays showed that changes in the permeability of theouter membrane are an unavoidable consequence of displayinga large periplasmic protein on the surface of E.coli. This isthe first systematic report on the effects that cell surfaceengineering may have on the integrity and permeability propertiesof bacterial outer membranes.     

Nuclear magnetic resonance studies and molecular dynamics simulations of the solution conformation of a 'designed', {alpha}-helical peptide

The complete three-dimensional structure in methanol of an amphipathica-helical peptide, that has been designed by taking into accountthe three-dimensional structures of small haemolytic peptides,secondary structure prediction algorithms and the well documentedliterature on -helix stabilizing factors, has been elucidatedby two-dimensional NMR spectroscopy. Initially various two-dimensionalspectra (COSY, TOCSY, and NOESY) allowed the complete sequencespecific assignment of all signals in the ¹H spectrum. Consequentlytrial structures were generated which were then subjected tomolecular dynamics simulations using 121 NOE-derived distancesand 25 vicinal coupling constant values as structural restraintsto give a final set of calculated structures. These structuresare in complete agreement with the results of a circular dichroismstudy and reveal that the peptide adopted a highly ordered a-helicalconformation. Details of the structure which throw light onfuture peptide/protein design are discussed.     

An SS1-SS2 {beta}-barrel structure for the voltage-activated potassium channel

To examine the feasibility of a ß structure for thepore-lining region of the voltage-gated potassium channel, wehave characterized a family of 12 antiparallel ß-barrels.Each is comprised of four identical pairs of ß-strandsorganized with approximate 4-fold symmetry about a channel axis.The Cand N-termini of the ß-strand pairs are assumedto be at the extracellular end of the channel, and each pairis connected by a hairpin turn at the intracellular end of thechannel. The models differ in the residues located in the hairpinturn and in the orientation of the two strands of each pairin the barrel, i.e. whether the C-terminus of a pair is clockwise(CW) or counterclockwise (CCW) from the N-terminus when thechannel is viewed from outside the cell. Following known structureprecedents and potential energy predictions, the barrel is assumedto be right-twisting in all cases. All models have crowded layersof inward-projecting aromatic sidechains near the center ofthe channel which could regulate channel selectivity. The modelswith an odd number of amino acids in the hairpin turn have theadvantage of predicting that F433 points into the barrel, butthe disadvantage that V438 does not. Of these models, two ofthe models are most consistent with the external tetraethylanunonhim(TEA) block data, and of those, one (T439 CCW 3:5) is most consistentwith the internal TEA block data.     

Alcohol-induced changes of {beta}-lactoglobulin - retinol-binding stoichiometry

It has been demonstrated using CD that ethanol induces importantsecondary structure changes of ß-lactoglobulin. CDspectra indicate that ß-lactoglobulin secondary structure,which is mainly composed of ß-strands, becomes mostly-helical under the influence of the solvent polarity changes.The midpoint of ß-strand/-helix transition in ß-lactoglobulinis observed at dielectric constant {small tilde}60 (35% ethanol;v/v). According to CD measurements, the ethanol-dependent secondarystructure changes are reversible. The alkylation of lysines-NH₂ in ß-lactoglobulin weakens the central ß-barrelstructure, since the ß-strand/-helix transition midpointof alkylated ß-lactoglobulin is shifted to lower ethanolconcentration (25% ethanol; v/v). ß-Lactoglobulinstructural changes are triggering the dissociation of the ß-lactoglobulin- retinol complex as judged from complete quenching of its fluorescencein ethanol concentration >30% (v/v). However, in 20% ethanol(v/v), ß-lactoglobulin still retains most of its nativesecondary structure as shown by CD and, in this condition, oneß-lactoglobulin molecule binds an additional secondretinol molecule. This suggests that the highly populated speciesobserved around 20% ethanol (v/v) might represent an intermediatestate able to bind two molecules of retinol.     

Secondary structure characterization of {beta}-lactamase inclusion bodies

The secondary structure of proteins in E.coli inclusion bodieswas investigated via Raman spectroscopy. Inclusion bodies werepurified from cells expressing different forms of RTEM ß-lactamaseand grown at either 37 or 42° C. All of the solid phaseinclusion body samples examined gave amide I band spectra thatwere perturbed from that of the native, purified protein inboth solution and powder forms; secondary structure estimatesindicated significant decreases in a-helix and increases inß-sheet contents in the inclusion body samples. The structureestimates for inclusion bodies isolated from 37°C cultureswere similar, regardless of aggregate localization in the E.colicytoplasmic or periplasmic spaces or ß-lactamase precursorcontent. Inclusion bodies obtained from 42°C cells exhibiteda further reduction of ß-helix and augmentation of ß-sheetcontents relative to those from 37°C cultures. These resultsare consistent with the paradigm for inclusion body formationvia the self-association of intra-cellular folding intermediateshaving extensive secondary structure content. Further, the overallsecondary structure content of inclusion bodies is not significantlyaffected by subcellular compartmentalization, but may be alteredat increased temperatures     

Length preferences and periodicity in {beta}-strands. Antiparallel edge {beta}-sheets are more likely to finish in non-hydrogen bonded rings

We analysed the length distributions of different types of ß-strandin a high resolution, non-homologous set of 500 protein structures,finding differences in their mean lengths. Antiparallel edgestrands in strand–turn–strand motifs show a preferencefor an even number of residues. This propensity is enhancedif the length is corrected for ß-bulges, which insertan extra residue into the strand. Residues in antiparallel edgeß-strands alternate between being in hydrogen bondedand non-hydrogen bonded rings. Antiparallel edges with an evennumber of residues are more likely to have their final ßresidue in a non-hydrogen bonded ring. This suggests that non-hydrogenbonded rings are intrinsically more stable than hydrogen bondedrings, perhaps because its side chain packing is closer. Therefore,we suggest that a simple way to increase ß-hairpinstability, or the stability of an antiparallel edge strand,is to have a non-hydrogen bonded ring at the end of the strand. Received June 19, 2003; revised October 25, 2003; accepted November 7, 2003     

Altered specificities of genetically engineered {alpha}1 antitrypsin variants

Seven active site variants of human ₁-antitrypsin (₁AT) wereproduced in Escherichia coli following site-specific mutagenesisof the ₁AT complementary DNA. ₁AT (Ala ³⁵⁸), ₁AT (Ile³⁵⁸ and₁AT (Val³⁵⁸), were efficient inhibitors of both neutrophil andpancreatic elastases, but not of cathepsin G. ₁AT (Ala³⁵⁸, Val³⁵⁸)and ₁AT (Phe³⁵⁸ specifically inhibited pancreatic elastase andcathepsin G respectively. The most potent inhibitor of neutrophilelastase was ₁AT (Leu³⁵⁸), which also proved to be effectiveagainst cathepsin G. The ₁AT (Arg³⁵⁸) variant inactivated thrombinwith kinetics similar to antithrombin III in the presence ofheparin. Electrophoretic analysis showed that SDS-stable highmol. wt complexes were formed between the mutant inhibitorsand the cognate proteases in each case. These data indicatethat effective inhibition occurs when the ₁AT P1 residue (position358) corresponds to the primary specificity of the target protease.Moreover, alteration of the P3 residue (position 356) can furthermodify the reactivity of the inhibitor. Two of the variantshave therapeutic potential: ₁AT (Leu³⁵⁸ may be more useful thanplasma ₁AT in the treatment of destructive lung disorders and₁ (Arg³⁵⁸ could be effective in the control of thrombosis.