Stochastic boundary molecular dynamics simulation of L-ribose in the active site of Actinoplanes missouriensis xylose isomerase and its Val135Asn mutant with improved reaction rate.

We used molecular dynamics simulations to study how a non-natural substrate, L-ribose, interacts with the active site of Actinoplanes missouriensis xylose isomerase. The simulations showed that L-ribose does not stay liganded in the active site in the same way as D-xylose, in which the oxygens O2 and O4 are liganded to the metal M1. The oxygen O4 of L-ribose moved away from the metal M1 to an upside down position. Furthermore, the distances of the carbons C1 and C2 of L-ribose to the catalytic metal M2 were higher than in the case of D-xylose. These findings explain the extremely low reaction rate of xylose isomerase with L-ribose. The mutation V135N close to the C5-OH of the substrate increased the reaction efficiency 2- to 4-fold with L-ribose. V135N did not affect the reaction with D-xylose and L-arabinose, whereas the reaction with D-glucose was impaired, probably due to a hydrogen bond between Asn-135 and the substrate. When L-ribose was the substrate, Asn-135 formed a hydrogen bond to Glu-181. As a consequence, O4 of L-ribose stayed liganded to the metal M1 in the V135N mutant in molecular dynamics simulations. This explains the decreased K(m) of the V135N mutant with L-ribose.

Engineering the substrate specificity of xylose isomerase.

Xylose isomerase (XI) catalyzes the isomerization and epimerization of hexoses, pentoses and tetroses. In order to clarify the reasons for the low reaction efficiency of a pentose sugar, L-arabinose, we determined the crystal structure of Streptomyces rubiginosus XI complexed with L-arabinose. The crystal structure revealed that, when compared with D-xylose and D-glucose, L-arabinose binds to the active site in a partially different position, in which the ligand has difficulties in binding the catalytic metal M2. Lys183 has been thought to stabilize the open substrate conformation by hydrogen bonding to oxygen O1. Our results with L-arabinose showed that the substrate stays in a linear form even without a hydrogen bond between Lys183 and oxygen O1. We engineered mutations to the active site of Actinoplanes missouriensis XI to improve the reaction efficiency with L-arabinose. The mutation F26W was intended to shift the position of oxygen O1 of L-arabinose closer to the catalytic metal M2. This effect of F26W was modeled by free energy perturbation simulations. In line with this, F26W increased 2-fold the catalytic efficiency of XI with L-arabinose; the increase was seen mainly in kcat. The mutation Q256D was outside the sphere of the catalytic residues and probably modified the electrostatic properties of the active site. It improved 3-fold the catalytic efficiency of XI with L-arabinose; this increase was seen in both Km and kcat. This study showed that it is possible to engineer the substrate specificity of XI.

Correlative motions and memory effects in molecular dynamics simulations of molecules: principal components and rescaled range analysis suggest that the motions of native BPTI are more correlated than those of its mutants.

In this work MD simulations of the native bovine pancreatic trypsin inhibitor (BPTI) and 16 mutants were done in vacuum in order to study memory effects in the mutants using principal component analysis (PCA) and the rescaled range analysis (Hurst exponents). Both PCA and the rescaled range analysis support our previous proposition, based on PCA of lysozyme, that the motions of a native protein are more correlated than those of mutants. The methods are compared, the nature and applications of the rule and the role of the long-range correlations in MD time series (i.e. memory) are discussed in the context of collective motions.

Stability and amino acid preferences of type VIII reverse turn: the most common turn in peptides?

Free energies of the alpha(r)beta and betabeta conformations of 14 tetrapeptides, based on the sequence SALN and protein X-ray structures, were calculated using molecular dynamics simulations and MM-PBSA calculations. The alphaalpha conformations of five of the tetrapeptides were also studied. SALN has been earlier shown by molecular dynamics simulations and NMR spectroscopy to have a tendency to form an alpha(r)beta turn. The gas-phase energy of the molecular mechanical force field (CHARMM), the electrostatic and non-polar solvation free energies and solute entropies were used to explain the free energy differences of the alphaalpha, betabeta and alpha(r)beta conformations of the peptides. The alpha(r)beta conformation of SALN and SATN was predicted to be slightly more stable than the extended conformation (betabeta), in agreement with experimental results. The SALN mutants SAIN, SAVN, SATN, SSIN and MSHV, were also predicted to be potential alpha(r)beta turn-forming peptides. We report also revised positional potentials for the type VIII turn, based on a non-homologous set of protein structures. This protein databank analysis confirms the main results of the earlier analyses and reveals several new amino acid residues with a significant positional preference. The results of this work led us to suggest that the alpha(r)beta turn may be the most common turn type in peptides. Such turns may be readily formed in aqueous solution and thereby play important roles in the protein folding process by serving as an initiation point for structure formation.

Conformational and biochemical analysis of the cyclic peptides which modulate serine protease activity.

Prostate-specific antigen (PSA), a member of the kallikrein sub-group of the trypsin serine protease family, is a widely used marker for prostate cancer. Several sequences with specific binding to PSA have been identified by using phage display peptide libraries. The GST-fusion proteins of the characterized sequences have been shown to increase the enzyme activity of PSA to a synthetic substrate. The corresponding three cyclic synthetic analogues CVFTSNYAFC (A-1), CVFAHNYNYLVC (B-2) and CVAYCIEHHCWTC (C-4) have similar PSA promoting activity. Despite differences in the amino acid sequences, all three peptides bind to the same region of PSA. The conformation of the peptides was investigated by proton NMR spectroscopy. In addition, alanine replacement was used to characterize the prerequisites for binding. It is proposed that interactions with PSA are based on the aromatic and hydrophobic features of the amino acid side chains. Furthermore, it is suggested that peptides form beta-turn structures forced by cysteine bridges directing important aromatic side chains to the same side of the turn-structure.