Domain motions in dihydrofolate reductase: a molecular dynamics study.

Molecular dynamics simulations have been carried out on the enzyme dihydrofolate reductase from Lactobacillus casei complexed with methotrexate, NADPH and 264 crystallographic water molecules. Analysis of correlations in atomic fluctuations reveal the presence of highly correlated motion (correlation coefficient > 0.6) in the region between residues 30 to 35 and 85 to 90 leading to the identification of two domains, an "adenosine-binding domain" and a "large domain", which rotate by 3 to 4 degrees with respect to each other. The strongest correlation (> 0.6) within the large domain involves a coupling between the motions of the "teen-loop", and the spatially contiguous loops linking beta 6-beta 7 and beta 7-beta 8. Moreover, there is a significant correlation (approximately 0.5) between the adenosine fragment of NADPH and the pteridine and p-aminobenzoyl fragments of methotrexate, which are separated by approximately 17 A, and is lost on removal of "rigid-body" motion from the original trajectory. This provides support for the idea that the relative motion of the two domains is a means by which the occupation of the binding site for the adenosine end of the coenzyme can affect methotrexate binding and vice versa. Quasiharmonic vibrational analysis of the trajectory reveals that the overall dynamics of the system are governed by domain motions whose contributions are dominant at low frequencies. In addition, different low-frequency modes are responsible for separately coupling the adenosine-binding site and parts of methotrexate.

The crystal structure of the L1 metallo-beta-lactamase from Stenotrophomonas maltophilia at 1.7 A resolution.

The structure of the L1 metallo-beta-lactamase from the opportunistic pathogen Stenotrophomonas maltophilia has been determined at 1.7 A resolution by the multiwavelength anomalous dispersion (MAD) approach exploiting both the intrinsic binuclear zinc centre and incorporated selenomethionine residues. L1 is unique amongst all known beta-lactamases in that it exists as a tetramer. The protein exhibits the alphabeta/betaalpha fold found only in the metallo-beta-lactamases and displays several unique features not previously observed in these enzymes. These include a disulphide bridge and two substantially elongated loops connected to the active site of the enzyme. Two closely spaced zinc ions are bound at the active site with tetrahedral (Zn1) and trigonal bipyramidal (Zn2) co-ordination, respectively; these are bridged by a water molecule which we propose acts as the nucleophile in the hydrolytic reaction. Ligation of the second zinc ion involves both residues and geometry which have not been previously observed in the metallo-beta-lactamases. Simulated binding of the substrates ampicillin, ceftazidime and imipenem suggests that the substrate is able to bind to the enzyme in a variety of different conformations whose common features are direct interactions of the beta-lactam carbonyl oxygen and nitrogen with the zinc ions and of the beta-lactam carboxylate with Ser187. We describe a catalytic mechanism whose principal features are a nucleophilic attack of the bridging water on the beta-lactam carbonyl carbon, electrostatic stabilisation of a negatively charged tetrahedral transition state and protonation of the beta-lactam nitrogen by a second water molecule co-ordinated by Zn2. Further, we propose that direct metal:substrate interactions provide a substantial contribution to substrate binding and that this may explain the lack of specificity which is a feature of this class of enzyme.

Crystal structure of GerE, the ultimate transcriptional regulator of spore formation in Bacillus subtilis.

The small, DNA-binding protein GerE regulates gene transcription in the terminally differentiated mother-cell compartment during late stages of sporulation in Bacillus subtilis. This versatile transcription factor shares sequence homology with the LuxR/FixJ/UhpA family of activators and modulates the expression of a number of genes, in particular those encoding the components of the coat that surrounds the mature spore. GerE orchestrates the final stages of coat deposition and maturation that lead to a spore with remarkable resistance properties but that must be responsive to low levels of germination signals. As this germination process is largely passive and can occur in the absence of de novo protein synthesis, the correct assembly of germination machinery, including germinant receptors and energy storage compounds, is crucial to the survival of the cell. The crystal structure of GerE has been solved at 2.05 A resolution using multi-wavelength anomalous dispersion techniques and reveals the nature of the GerE dimer. Each monomer comprises four alpha-helices, of which the central pair forms a helix-turn-helix DNA-binding motif. Implications for DNA-binding and the structural organisation of the LuxR/FixJ/UhpA family of transcription activator domains are discussed.

Understanding the mechanism of ice binding by type III antifreeze proteins.

Type III antifreeze proteins (AFPs) are present in the body fluids of some polar fishes where they inhibit ice growth at subzero temperatures. Previous studies of the structure of type III AFP by NMR and X-ray identified a remarkably flat surface on the protein containing amino acids that were demonstrated to be important for interaction with ice by mutational studies. It was proposed that this protein surface binds onto the (1 0 [\bar 1] 0) plane of ice with the key amino acids interacting directly with the water molecules in the ice crystal. Here, we show that the mechanism of type III AFP interaction with ice crystals is more complex than that proposed previously. We report a high-resolution X-ray structure of type III AFP refined at 1.15 A resolution with individual anisotropic temperature factors. We report the results of ice-etching experiments that show a broad surface coverage, suggesting that type III AFP binds to a set of planes that are parallel with or inclined at a small angle to the crystallographic c-axis of the ice crystal. Our modelling studies, performed with the refined structure, confirm that type III AFP can make energetically favourable interactions with several ice surfaces.

Structure of a slow processing precursor penicillin acylase from Escherichia coli reveals the linker peptide blocking the active-site cleft.

Penicillin G acylase is a periplasmic protein, cytoplasmically expressed as a precursor polypeptide comprising a signal sequence, the A and B chains of the mature enzyme (209 and 557 residues respectively) joined by a spacer peptide of 54 amino acid residues. The wild-type AB heterodimer is produced by proteolytic removal of this spacer in the periplasm. The first step in processing is believed to be autocatalytic hydrolysis of the peptide bond between the C-terminal residue of the spacer and the active-site serine residue at the N terminus of the B chain. We have determined the crystal structure of a slowly processing precursor mutant (Thr263Gly) of penicillin G acylase from Escherichia coli, which reveals that the spacer peptide blocks the entrance to the active-site cleft consistent with an autocatalytic mechanism of maturation. In this mutant precursor there is, however, an unexpected cleavage at a site four residues from the active-site serine residue. Analyses of the stereochemistry of the 260-261 bond seen to be cleaved in this precursor structure and of the 263-264 peptide bond have suggested factors that may govern the autocatalytic mechanism.

Rational design of thermally stable proteins: relevance to bionanotechnology.

Design of thermally stable proteins is spurred by their applications in bionanotechnology. There are three major issues governing this: first, the upper limit on the temperature at which proteins remain physiologically active and are available for technological applications (answers may emerge from the discovery of new, natural hyperthermophilic enzymes that are active above 125 degrees C or from the selection of mutants of hyperthermophilic enzymes that are more stable); second, the use of hyperthermophilic enzymes as molecular templates to design highly stable enzymes that have high activity at low temperatures; third, the link between rigidity and flexibility to thermostability and activity, respectively. We review progress in these areas.

Conformational change in the activation of lipase: an analysis in terms of low-frequency normal modes.

The interfacial activation of Rhizomucor miehei lipase (RmL) involves the motion of an alpha-helical region (residues 82-96) which acts as a "lid" over the active site of the enzyme, undergoing a displacement from a "closed" to an "open" conformation upon binding of substrate. Normal mode analyses performed in both low and high dielectric media reveal that low-frequency vibrational modes contribute significantly to the conformational transition between the closed and open conformations. In these modes, the lid displacement is coupled to local motions of active site loops as well as global breathing motions. Atomic fluctuations of the first hinge of the lid (residues 83-84) are substantially larger in the low dielectric medium than in the high dielectric medium. Our results also suggest that electrostatic interactions of Arg86 play an important role in terms of both the intrinsic stability of the lid and its displacement, through enhancement of hinge mobility in a high dielectric medium. Additional calculations demonstrate that the observed patterns of atomic fluctuations are an intrinsic feature of the protein structure and not dependent on the nature of specific energy minima.

Stabilizing the eIF4G1 α-helix increases its binding affinity with eIF4E: implications for peptidomimetic design strategies.

Eukaryotic initiation factor (eIF)4E is overexpressed in many types of cancer such as breast, head and neck, and lung. A consequence of increased levels of eIF4E is the preferential translation of pro-tumorigenic proteins such as c-Myc, cyclin D1, and vascular endothelial growth factor. Inhibition of eIF4E is therefore a potential therapeutic target for human cancers. A novel peptide based on the eIF4E-binding peptide eIF4G1, where the α-helix was stabilized by the inclusion of α-helix inducers as shown by CD measurements, was synthesized. The helically stabilized peptide binds with an apparent K(d) of 9.43±2.57 nM, which is ∼15.7-fold more potent than the template peptide from which it is designed. The helically stabilized peptide showed significant biological activity at a concentration of 400 µM, unlike the naturally occurring eIFG1 peptide when measured in cell-based cap-dependent translational reporter and WST-1 (4-[3-(4-iodophenyl)-2-(4-nitrophenyl)-2H-5-tetrazolio]-1,3-benzene disulfonate) assays. Fusion of the template peptide and the stabilized peptide to the cell-penetrating peptide TAT produced more active but equally potent inhibitors of cap-dependent translation in cell lines. They also equally disrupted cell metabolism as measured in a WST-1 assay. Propidium iodide staining revealed that the TAT-fused, helically stabilized peptide caused more cell death than the TAT-fused eIF4G1 template peptide with substantial decreases in the G1 and G2 cell populations. Annexin-staining experiments also indicated that the TAT-fused eIF4G1 derivative peptides caused cell death by apoptosis. The results presented should offer further insight into peptidomimetics development for eIF4E.

Updates on p53: modulation of p53 degradation as a therapeutic approach.

The p53 pathway is aberrant in most human tumours with over 50% expressing mutant p53 proteins. The pathway is critically controlled by protein degradation. Here, we discuss the latest developments in the search for small molecules that can modulate p53 pathway protein stability and restore p53 activity for cancer therapy.

Protein flexibility: its role in structure and mechanism revealed by molecular simulations.

Computer simulations at the atomic level have arrived at a stage where they provide realistic modeling of flexibility in proteins (and the mobility of their associated solvent) that is important in understanding the nature of molecular motions. This can now be extended to the molecular and atomic motions that are associated with protein mechanisms. Moreover, the derived data agree reasonably accurately with experimental measurements of several kinetic and thermodynamic parameters. Fundamental insights emerge on the roles that this intrinsic flexibility plays in the thermodynamic characteristics of macromolecules in solution; these equip the investigator to probe the consequences of cognate interactions and ligand binding on entropy and enthalpy. Thus simulations can now provide a powerful tool for investigating protein mechanisms that complements the existing and the emerging experimental techniques.

Binding of buried structural water increases the flexibility of proteins.

Water deeply buried in proteins is considered to be an integral part of the folded structure. Such structural water molecules make strong H bonds with polar groups of the surrounding protein and therefore are believed to tighten the protein matrix. Surprisingly, our computational analysis of the binding of a buried water molecule to bovine pancreatic trypsin inhibitor shows that the protein actually becomes more flexible, as revealed by an increase in the vibrational entropy. We find that this effect must be common in proteins, because the large entropic cost of immobilizing a single water molecule [-TDeltaS = 20.6 kcal/mol (1 kcal = 4.18 kJ) for the lost translational and rotational degrees of freedom] can only be partly compensated by water-protein interactions, even when they are nearly perfect, as in the case of bovine pancreatic trypsin inhibitor (DeltaE = -19.8 kcal/mol), leaving no room for a further decrease in entropy from protein tightening. This study illustrates the importance of considering changes in protein flexibility (which in this case favor binding by 3.5 kcal/mol) for the prediction of ligand binding affinities.

Structural consequences of the B5 histidine --> tyrosine mutation in human insulin characterized by X-ray crystallography and conformational analysis.

The addition of phenols to hexameric insulin solutions produces a particularly stable hexamer, resulting from a rearrangement in which residues B1-B8 change from an extended conformation (T-state) to form an alpha-helix (R-state). The R-state is, in part, stabilized by nonpolar interactions between the phenolic molecule and residue B5 His at the dimer-dimer interface. The B5 His --> Tyr mutant human insulin was constructed to see if the tyrosine side chain would mimic the effect of phenol binding in the hexamer and induce the R-state. In partial support of this hypothesis, the molecule crystallized as a half-helical hexamer (T(3)R(3)) in conditions that conventionally promote the fully nonhelical (T6) form. As expected, in the presence of phenol or resorcinol, the B5 Tyr hexamers adopt the fully helical (R6) conformation. Molecular modeling calculations were performed to investigate the conformational preference of the T-state B5 Tyr side chain in the T(3)R(3) form, this side chain being associated with structural perturbations of the A7-A10 loop in an adjacent hexamer. For an isolated dimer, several different orientations of the side chain were found, which were close in energy and readily interconvertible. In the crystal environment only one of these conformations remains low in energy; this conformation corresponds to that observed in the crystal structure. This suggests that packing constraints around residue B5 Tyr result in the observed structural rearrangements. Thus, rather than promoting the R-state in a manner analogous to phenol, the mutation appears to destabilize the T-state. These studies highlight the role of B5 His in determining hexamer conformation and in mediating crystal packing interactions, properties that are likely be important in vivo.

Probing the mechanism of ATP hydrolysis and substrate translocation in the AAA protease FtsH by modelling and mutagenesis.

We have built a homology model of the AAA domain of the ATP-dependent protease FtsH of Escherichia coli based on the crystal structure of the hexamerization domain of N-ethylmaleimide-sensitive fusion protein. The resulting model of the hexameric ring of the ATP-bound form of the AAA ATPase suggests a plausible mechanism of ATP binding and hydrolysis, in which invariant residues of Walker motifs A and B and the second region of homology, characteristic of the AAA ATPases, play key roles. The importance of these invariant residues was confirmed by site-directed mutagenesis. Further modelling suggested a mechanism by which ATP hydrolysis alters the conformation of the loop forming the central hole of the hexameric ring. It is proposed that unfolded polypeptides are translocated through the central hole into the protease chamber upon cycles of ATP hydrolysis. Degradation of polypeptides by FtsH is tightly coupled to ATP hydrolysis, whereas ATP binding alone is sufficient to support the degradation of short peptides. Furthermore, comparative structural analysis of FtsH and a related ATPase, HslU, reveals interesting similarities and differences in mechanism.

Ionization characteristics and chemical influences of aspartic acid residue 158 of papain and caricain determined by structure-related kinetic and computational techniques: multiple electrostatic modulators of active-centre chemistry.

The pK(a) of (Asp(158))-CO(2)H of papain (EC 3.4.22.2) was determined as 2.8 by using 4-chloro-7-nitrobenzofurazan (Nbf-Cl) as a reactivity probe targeted on the thiolate anion component of the Cys(25)/His(159) nucleophilic-acid/base motif of the catalytic site. The possibility of using Nbf-Cl for this purpose was established by modelling the papain-Nbf-Cl Meisenheimer intermediate by using QUANTA/CHARMM and performing molecular orbital calculations with MOPAC interfaced with Cerius 2. A pH-dependent stopped-flow kinetic study of the reaction of papain with Nbf-Cl established that the striking rate maximum at pH 3 results from reaction in a minor ionization state comprising (Cys(25))-S(-)/(His(159))-Im(+)H (in which Im represents imidazole) produced by protonic dissociation of (Cys(25))-SH/(His(159))-Im(+)H with pK(a) 3.3 and (Asp(158))-CO(2)H. Although the analogous intermediate in the reaction of caricain (EC 3.4.22.30) with Nbf-Cl has similar geometry, the pH-k profile (k being the second-order rate constant) lacks a rate maximum under acidic conditions. This precludes the experimental determination of the pK(a) value of (Asp(158))-CO(2)H of caricain, which was calculated to be 2.0 by solving the linearized Poisson-Boltzmann equation with the program UHBD ('University of Houston Brownian dynamics'). A value lower than 2.8 had been predicted by consideration of the hydrogen-bonded networks involving Asp(158) and its microenvironments in both enzymes. The difference between these pK(a) values (values not previously detected in reactions of either enzyme) accounts for the lack of the rate maximum in the caricain reaction and for the differences in the electronic absorption spectra of the two S-Nbf-enzymes under acidic conditions. The concept of control of cysteine proteinase activity by multiple electrostatic modulators, including (Asp(158))-CO(2)(-), which modifies traditional mechanistic views, is discussed.

Structural origins of the interfacial activation in Thermomyces (Humicola) lanuginosa lipase.

The already known X-ray structures of lipases provide little evidence about initial, discrete structural steps occurring in the first phases of their activation in the presence of lipids (process referred to as interfacial activation). To address this problem, five new Thermomyces (formerly Humicola) lanuginosa lipase (TlL) crystal structures have been solved and compared with four previously reported structures of this enzyme. The bias coming from different crystallization media has been minimized by the growth of all crystals under the same crystallization conditions, in the presence of detergent/lipid analogues, with low or high ionic strength as the only main variable. Resulting structures and their characteristic features allowed the identification of three structurally distinct species of this enzyme: low activity form (LA), activated form (A), and fully Active (FA) form. The isomerization of the Cys268-Cys22 disulfide, synchronized with the formation of a new, short alpha(0) helix and flipping of the Arg84 (Arginine switch) located in the lid's proximal hinge, have been postulated as the key, structural factors of the initial transitions between LA and A forms. The experimental results were supplemented by theoretical calculations. The magnitude of the activation barrier between LA (ground state) and A (end state) forms of TlL (10.6 kcal/mol) is comparable to the enthalpic barriers typical for ring flips and disulfide isomerizations at ambient temperatures. This suggests that the sequence of the structural changes, as exemplified in various TlL crystal structures, mirror those that may occur during interfacial activation.

Variation in aspects of cysteine proteinase catalytic mechanism deduced by spectroscopic observation of dithioester intermediates, kinetic analysis and molecular dynamics simulations.

The possibility of a slow post-acylation conformational change during catalysis by cysteine proteinases was investigated by using a new chromogenic substrate, N-acetyl-Phe-Gly methyl thionoester, four natural variants (papain, caricain, actinidin and ficin), and stopped-flow spectral analysis to monitor the pre-steady state formation of the dithioacylenzyme intermediates and their steady state hydrolysis. The predicted reversibility of acylation was demonstrated kinetically for actinidin and ficin, but not for papain or caricain. This difference between actinidin and papain was investigated by modelling using QUANTA and CHARMM. The weaker binding of hydrophobic substrates, including the new thionoester, by actinidin than by papain may not be due to the well-known difference in their S2-subsites, whereby that of actinidin in the free enzyme is shorter due to the presence of Met211. Molecular dynamics simulation suggests that during substrate binding the sidechain of Met211 moves to allow full access of a Phe sidechain to the S2-subsite. The highly anionic surface of actinidin may contribute to the specificity difference between papain and actinidin. During subsequent molecular dynamics simulations the P1 product, methanol, diffuses rapidly (over<8 ps) out of papain and caricain but 'lingers' around the active centre of actinidin. Uniquely in actinidin, an Asp142-Lys145 salt bridge allows formation of a cavity which appears to constrain diffusion of the methanol away from the catalytic site. The cavity then undergoes large scale movements (over 4.8 A) in a highly correlated manner, thus controlling the motions of the methanol molecule. The changes in this cavity that release the methanol might be those deduced kinetically.