Exploration of swapping enzymatic function between two proteins: a simulation study of chorismate mutase and isochorismate pyruvate lyase.

The enzyme chorismate mutase EcCM from Escherichia coli catalyzes one of the few pericyclic reactions in biology, the transformation of chorismate to prephenate. The isochorismate pyruvate lyase PchB from Pseudomonas aeroginosa catalyzes another pericyclic reaction, the isochorismate to salicylate transformation. Interestingly, PchB possesses weak chorismate mutase activity as well thus being able to catalyze two distinct pericyclic reactions in a single active site. EcCM and PchB possess very similar folds, despite their low sequence identity. Using molecular dynamics simulations of four combinations of the two enzymes (EcCM and PchB) with the two substrates (chorismate and isochorismate) we show that the electrostatic field due to EcCM at atoms of chorismate favors the chorismate to prephenate transition and that, analogously, the electrostatic field due to PchB at atoms of isochorismate favors the isochorismate to salicylate transition. The largest differences between EcCM and PchB in electrostatic field strengths at atoms of the substrates are found to be due to residue side chains at distances between 0.6 and 0.8 nm from particular substrate atoms. Both enzymes tend to bring their non-native substrate in the same conformation as their native substrate. EcCM and to a lower extent PchB fail in influencing the forces on and conformations of the substrate such as to favor the other chemical reaction (isochorismate pyruvate lyase activity for EcCM and chorismate mutase activity for PchB). These observations might explain the difficulty of engineering isochorismate pyruvate lyase activity in EcCM by solely mutating active site residues.

Multi-resolution simulation of biomolecular systems: a review of methodological issues.

Theoretical-computational modeling with an eye to explaining experimental observations in regard to a particular chemical phenomenon or process requires choices concerning essential degrees of freedom and types of interactions and the generation of a Boltzmann ensemble or trajectories of configurations. Depending on the degrees of freedom that are essential to the process of interest, for example, electronic or nuclear versus atomic, molecular or supra-molecular, quantum- or classical-mechanical equations of motion are to be used. In multi-resolution simulation, various levels of resolution, for example, electronic, atomic, supra-atomic or supra-molecular, are combined in one model. This allows an enhancement of the computational efficiency, while maintaining sufficient detail with respect to particular degrees of freedom. The basic challenges and choices with respect to multi-resolution modeling are reviewed and as an illustration the differential catalytic properties of two enzymes with similar folds but different substrates with respect to these substrates are explored using multi-resolution simulation at the electronic, atomic and supra-molecular levels of resolution.

Molecular dynamics simulation of the last step of a catalytic cycle: product release from the active site of the enzyme chorismate mutase from Mycobacterium tuberculosis.

The protein chorismate mutase MtCM from Mycobacterium tuberculosis catalyzes one of the few pericyclic reactions known in biology: the transformation of chorismate to prephenate. Chorismate mutases have been widely studied experimentally and computationally to elucidate the transition state of the enzyme catalyzed reaction and the origin of the high catalytic rate. However, studies about substrate entry and product exit to and from the highly occluded active site of the enzyme have to our knowledge not been performed on this enzyme. Crystallographic data suggest a possible substrate entry gate, that involves a slight opening of the enzyme for the substrate to access the active site. Using multiple molecular dynamics simulations, we investigate the natural dynamic process of the product exiting from the binding pocket of MtCM. We identify a dominant exit pathway, which is in agreement with the gate proposed from the available crystallographic data. Helices H2 and H4 move apart from each other which enables the product to exit from the active site. Interestingly, in almost all exit trajectories, two residues arginine 72 and arginine 134, which participate in the burying of the active site, are accompanying the product on its exit journey from the catalytic site.

Preferential affinity of the components of liquid mixtures at a rigid non-polar surface: enthalpic and entropic driving forces.

The modulation of the properties of lipid membranes by polyhydroxylated cosolutes such as sugars is a phenomenon of considerable biological, technological and medicinal relevance. A few years ago, we proposed the sugar-like mechanism--binding driven by the release of water molecules--as an attempt to rationalize the preferential affinity of carbohydrate molecules compared to water molecules for the surface of lipid bilayers, which is presumably related to the bioprotective action of these compounds. The goal herein is to gain a better understanding of the driving force underlying this mechanism, in terms of specific interactions or effects, as well as in terms of the energy-entropy partitioning. This is done in the simplest possible context of an apolar rigid-wall model representing the membrane, and mixtures of closely related and possibly artificial species in solution, namely monomers or dimers of Lennard-Jones particles, water with physical or reduced charges, and hydroxymethyl groups. The results indicate that although the sugar-like mechanism seems phenomenologically reasonable, the main driving force underlying this mechanism is not the entropy gain upon releasing water molecules into the bulk, as originally suggested, but rather the hydrophobic effect. Note that the latter effect is a generic concept and may in principle involve both a solvent release and an interaction component, depending on the solute considered.

Definition and testing of the GROMOS force-field versions 54A7 and 54B7.

New parameter sets of the GROMOS biomolecular force field, 54A7 and 54B7, are introduced. These parameter sets summarise some previously published force field modifications: The 53A6 helical propensities are corrected through new φ/ψ torsional angle terms and a modification of the N-H, C=O repulsion, a new atom type for a charged -CH(3) in the choline moiety is added, the Na(+) and Cl(-) ions are modified to reproduce the free energy of hydration, and additional improper torsional angle types for free energy calculations involving a chirality change are introduced. The new helical propensity modification is tested using the benchmark proteins hen egg-white lysozyme, fox1 RNA binding domain, chorismate mutase and the GCN4-p1 peptide. The stability of the proteins is improved in comparison with the 53A6 force field, and good agreement with a range of primary experimental data is obtained.

Membrane protein dynamics in different environments: simulation study of the outer membrane protein X in a lipid bilayer and in a micelle.

The bacterial outer membrane protein OmpX from Escherichia coli has been investigated by molecular dynamics simulations when embedded in a phospholipid bilayer and as a protein-micelle aggregate. The resulting simulation trajectories were analysed in terms of structural and dynamic properties of the membrane protein. In agreement with experimental observations, highest relative stability was found for the ß-barrel region that is embedded in the lipophilic phase, whereas an extracellular protruding ß-sheet, which is a unique structural feature of OmpX that supposedly plays an important role in cell adhesion and invasion, shows larger structure fluctuations. Additionally, we investigated water permeation into the core of the ß-barrel protein, which contains a tight salt-bridge and hydrogen-bond network, so that extensive water flux is unlikely. Differences between the bilayer and the micellar system were observed in the length of the barrel and its position inside the lipid environment, and in the protein interactions with the hydrophilic part of the lipids near the lipid/water interface. Those variations suggest that micelles and other detergent environments might not offer a wholly membrane-like milieu to promote adoption of the physiological conformational state by OmpX.

Exact distances and internal dynamics of perdeuterated ubiquitin from NOE buildups.

It is proposed to convert nuclear Overhauser effects (NOEs) into relatively precise distances for detailed structural studies of proteins. To this purpose, it is demonstrated that the measurement of NOE buildups between amide protons in perdeuterated human ubiquitin using a designed (15)N-resolved HMQC-NOESY experiment enables the determination of (1)H(N)-(1)H(N) distances up to 5 A with high accuracy and precision. These NOE-derived distances have an experimental random error of approximately 0.07 A, which is smaller than the pairwise rmsd (root-mean-square deviation) of 0.24 A obtained with corresponding distances extracted from either an NMR or an X-ray structure (pdb codes: 1D3Z and 1UBQ), and also smaller than the pairwise rmsd between distances from X-ray and NMR structures (0.15 A). Because the NOE contains both structural and dynamical information, a comparison between the 3D structures and NOE-derived distances may also give insights into through-space dynamics. It appears that the extraction of motional information from NOEs by comparison to the X-ray structure or the NMR structure is challenging because the motion may be masked by the quality of the structures. Nonetheless, a detailed analysis thereof suggests motions between beta-strands and large complex motions in the alpha-helix of ubiquitin. The NOE-derived motions are, however, of smaller amplitude and possibly of a different character than those present in a 20 ns molecular dynamic simulation of ubiquitin in water using the GROMOS force field. Furthermore, a recently published set of structures representing the conformational distribution over time scales up to milliseconds (pdb: 2K39) does not satisfy the NOEs better than the single X-ray structure. Hence, the measurement of possibly thousands of exact NOEs throughout the protein may serve as an excellent probe toward a correct representation of both structure and dynamics of proteins.

Relative tolerance of an enzymatic molten globule and its thermostable counterpart to point mutation.

Enzyme structures reflect the complex interplay between the free energy of unfolding (DeltaG) and catalytic efficiency. Consequently, the effects of point mutations on structure, stability, and function are difficult to predict. It has been proposed that the mutational robustness of homologous enzymes correlates with a higher initial DeltaG. To examine this issue, we compared the tolerance of a natural thermostable chorismate mutase and an engineered molten globular variant to targeted mutation. These mutases possess similar sequence, structure, and catalytic efficiency but dramatically different DeltaG values. We find that analogous point mutations can have widely divergent effects on catalytic activity in these scaffolds. In a set of five rationally designed single-amino acid changes, the thermostable scaffold suffers activity losses ranging from 50-fold smaller, for an aspartate-to-glycine substitution at the active site, to 2-fold greater, for a phenylalanine-to-tryptophan substitution in the hydrophobic core, versus that of the molten globular scaffold. However, biophysical characterization indicates that the variations in catalytic efficiency are not caused by losses of either secondary structural integrity or thermodynamic stability. Rather, the activity differences between variant pairs are very much context-dependent and likely stem from subtle changes in the fine structure of the active site. Thus, in many cases, it may be more productive to focus on changes in local conformation than on global stability when attempting to understand and predict how enzymes respond to point mutations.