Photochemical Reaction Dynamics of the Primary Event of Vision Studied by Means of a Hybrid Molecular Simulation

The photoisomerization reaction dynamics of a retinal chromophore in the visual receptor rhodopsin was investigated by means of hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations. The photoisomerization reaction of retinal constitutes the primary step of vision and is known as one of the fastest reactions in nature. To elucidate the molecular mechanism of the high efficiency of the reaction, we carried out hybrid ab initio QM/MM MD simulations of the complete reaction process from the vertically excited state to the photoproduct via electronic transition in the entire chromophore-protein complex. An ensemble of reaction trajectories reveal that the excited-state dynamics is dynamically homogeneous and synchronous even in the presence of thermal fluctuation of the protein, giving rise to the very fast formation of the photoproduct. The synchronous nature of the reaction dynamics in rhodopsin is found to originate from weak perturbation of the protein surroundings and from dynamic regulation of volume-conserving motions of the chromophore. The simulations also provide a detailed view of time-dependent modulations of hydrogen-out-of-plane vibrations during the reaction process, and identify molecular motions underlying the experimentally observed dynamic spectral modulations.     

Ab initio models for receptor-ligand interactions in proteins. 4. Model assembly study of the catalytic mechanism of triosephosphate isomerase

The catalytic mechanism of triosephosphate isomerase (TIM) was investigated with ab initio quantum mechanical calculations. Electrostatic interactions between the quantum mechanical active site and the protein and solvent environment were modeled using the finite difference Poission-Boltzman method. The complexes of TIM with the substrate dihydroxyacetone phosphate (DHAP), five possible intermediates and the product glyceraldehyde-3-phosphate (GAP) were optimized in the active-site model at the 3-21G^(*) level and energy profile for the proton abstraction from DHAP by the active-site Glu167 was calculated at the MP2/3-21G^(*)//3-21G^(*) level. Calculated energetics of the enzyme reaction were found to be in reasonable agreement with the experimental findings. Calculations revealed that an enediol of the substrate is a probable intermediate in the enzyme reaction. It was suggested that the proton abstracted from the substrate by the active-site glutamate goes to the carbonyl oxygen of the substrate producing enediol intermediate either directly or after it is exchanged with solvent. © 1996 Wiley-Liss, Inc.     

Acetylcholinesterase: mechanisms of covalent inhibition of H447I mutant determined by computational analyses

The reaction mechanisms of two inhibitor TFK(+) and TFK(0) binding to H447I mutant mouse acetylcholinesterase (mAChE) have been investigated by using a combined ab initio quantum mechanical/molecular mechanical (QM/MM) approach and classical molecular dynamics (MD) simulations. TFK(+) binding to the H447I mutant may proceed with a different reaction mechanism from the wild-type. A water molecule takes over the role of His447 and participates in the bond breaking and forming as a "charge relayer". Unlike in the wild-type mAChE case, Glu334, a conserved residue from the catalytic triad, acts as a catalytic base in the reaction. The calculated energy barrier for this reaction is about 8kcal/mol. These predictions await experimental verification. In the case of the neutral ligand TFK(0), however, multiple MD simulations on the TFK(0)/H447I complex reveal that none of the water molecules can be retained in the active site as a "catalytic" water. Taken together our computational studies confirm that TFK(0) is almost inactive in the H447I mutant, and also provide detailed mechanistic insights into the experimental observations.     

Theoretical investigations on Azotobacter vinelandii ferredoxin I: effects of electron transfer on protein dynamics

Structural, energetic, and dynamical studies of Azotobacter vinelandii ferredoxin I are presented for native and mutant forms. The protein contains two iron-sulfur clusters, one of which ([3Fe-4S]) is believed to play a central role in the electron-coupled proton transfer. Different charge sets for the [3Fe-4S] cluster in its reduced and oxidized state are calculated with broken symmetry ab initio density functional theory methods and used in molecular dynamics (MD) simulations. The validity of the ab initio calculations is assessed by comparing partially optimized structures of the [3Fe-4S] clusters with x-ray structures. Possible proton transfer pathways between the protein and the iron-sulfur cluster are examined by both MD simulations and ab initio calculations. The MD simulations identify three main-chain hydrogen atoms--HN(13), HN(14), and HN(16)--that are within H-bonding distance of the [3Fe-4S] cluster throughout the MD simulations. They could thus play a role in the proton transfer from the protein to the iron-sulfur cluster. By contrast, the HD2(15) atom of the Asp-15 is seldom close enough to the [3Fe-4S] cluster to transfer a proton. Poisson-Boltzmann calculations indicate that there is a low, but nonzero probability, that Asp-15 is protonated at pH 7; this is a requirement for it to serve as a proton donor. Ab initio calculations with a fragment model for the protein find similar behavior for the transfer of a proton from the OH of the protonated side chain and the main-chain NH of Asp-15. The existence of a stable salt bridge between Asp-15 and Lys-84 in the D15E mutant, versus its absence in the wild-type, has been suggested as the cause of the difference in the rate of proton transfer. Extensive MD simulations were done to test this idea; the results do not support the proposal. The present findings, together with the available data, serve as the basis for an alternative proposal for the mechanism of the coupled electron-proton transfer reaction in ferredoxin I.     

Atomic resolution crystal structures and quantum chemistry meet to reveal subtleties of hydroxynitrile lyase catalysis

Hydroxynitrile lyases are versatile enzymes that enantiospecifically cope with cyanohydrins, important intermediates in the production of various agrochemicals or pharmaceuticals. We determined four atomic resolution crystal structures of hydroxynitrile lyase from Hevea brasiliensis: one native and three complexes with acetone, isopropyl alcohol, and thiocyanate. We observed distinct distance changes among the active site residues related to proton shifts upon substrate binding. The combined use of crystallography and ab initio quantum chemical calculations allowed the determination of the protonation states in the enzyme active site. We show that His(235) of the catalytic triad must be protonated in order for catalysis to proceed, and we could reproduce the cyanohydrin synthesis in ab initio calculations. We also found evidence for the considerable pK(a) shifts that had been hypothesized earlier. We envision that this knowledge can be used to enhance the catalytic properties and the stability of the enzyme for industrial production of enantiomerically pure cyanohydrins.     

Kinemage of action - proposed reaction mechanism of glutamate-1-semialdehyde aminomutase at an atomic level

Glutamate-1-semialdehyde aminomutase (GSAM), a key enzyme in tetrapyrrole cofactor biosynthesis, performs a unique transamination on a single substrate. The substrate, glutamate-1-semialdehyde (GSA), undergoes a reaction that exchanges the position of an amine and a carbonyl group to produce 5-aminolevulinic acid (ALA). This transamination reaction is unique in the fact that is does not require an external cofactor to act as a nitrogen donor or acceptor in this transamination reaction. One of the other remarkable features of the catalytic mechanism is the release free in the enzyme active site of the intermediate 4,5-diaminovaleric acid (DAVA). The action of a gating loop prevents the escape of DAVA from the active site. In a MD simulation approach, using snapshots provided by X-ray crystallography and protein crystal absorption spectrometry data, the individual catalytic steps in this unique intramolecular transamination have been elucidated.     

Flexibility of prolyl oligopeptidase: molecular dynamics and molecular framework analysis of the potential substrate pathways

The flexibility of prolyl oligopeptidase has been investigated using molecular dynamics (MD) and molecular framework approaches to delineate the route of the substrate to the active site. The selectivity of the enzyme is mediated by a seven-bladed beta-propeller that in the crystal structure does not indicate the possible passage for the substrate to the catalytic center. Its open topology however, could allow the blades to move apart and let the substrate into the large central cavity. Flexibility analysis of prolyl oligopeptidase structure using the FIRST (Floppy Inclusion and Rigid Substructure Topology) approach and the atomic fluctuations derived from MD simulations demonstrated the rigidity of the propeller domain, which does not permit the substrate to approach the active site through this domain. Instead, a smaller tunnel at the inter-domain region comprising the highly flexible N-terminal segment of the peptidase domain and a facing hydrophilic loop from the propeller (residues 192-205) was identified by cross-correlation analysis and essential dynamics as the only potential pathway for the substrate. The functional importance of the flexible loop has been also verified by kinetic analysis of the enzyme with a split loop. Catalytic effect of engineered disulfide bridges was rationalized by characterizing the concerted motions of the two domains.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

A Wrench in the Works of Human Acetylcholinesterase: Soman Induced Conformational Changes Revealed by Molecular Dynamics Simulations

Irreversible inactivation of human acetylcholinesterase (hAChE) by organophosphorous pesticides (OPs) and chemical weapon agents (CWA) has severe morbidity and mortality consequences. We present data from quantum mechanics/molecular mechanics (QM/MM) and 80 classical molecular dynamics (MD) simulations of the apo and soman-adducted forms of hAChE to investigate the effects on the dynamics and protein structure when the catalytic Serine 203 is phosphonylated. We find that the soman phosphonylation of the active site Ser203 follows a water assisted addition-elimination mechanism with the elimination of the fluoride ion being the highest energy barrier at 6.5 kcal/mole. We observe soman-dependent changes in backbone and sidechain motions compared to the apo form of the protein. These alterations restrict the soman-adducted hAChE to a structural state that is primed for the soman adduct to be cleaved and removed from the active site. The altered motions and resulting structures provide alternative pathways into and out of the hAChE active site. In the soman-adducted protein both side and back door pathways are viable for soman adduct access. Correlation analysis of the apo and soman adducted MD trajectories shows that the correlation of gorge entrance and back door motion is disrupted when hAChE is adducted. This supports the hypothesis that substrate and product can use two different pathways as entry and exit sites in the apo form of the protein. These alternative pathways have important implications for the rational design of medical countermeasures.     

Exosite interactions determine the affinity of factor X for the extrinsic Xase complex

The initiation of coagulation results from the activation of factor X by an enzyme complex (Xase) composed of the trypsin-like serine proteinase, factor VIIa, bound to tissue factor (TF) on phospholipid membranes. We have investigated the basis for the protein substrate specificity of Xase using TF reconstituted into vesicles of phosphatidylcholine, phosphatidylserine, or pure phosphatidylcholine. We show that occupation of the active site of VIIa within Xase by a reversible inhibitor or an alternate peptidyl substrate is sufficient to exclude substrate interactions at the active site but does not alter the affinity of Xase for factor X. This is evident as classical competitive inhibition of peptidyl substrate cleavage but as classical noncompetitive inhibition of factor X activation by active site-directed ligands. This implies that the productive recognition of factor X by Xase arises from a multistep reaction requiring an initial interaction at sites on the enzyme complex distinct from the active site (exosites), followed by active site interactions and bond cleavage. Exosite interactions determine protein substrate affinity, whereas the second binding step influences the maximum catalytic rate for the reaction. We also show that competitive inhibition can be achieved by interfering with exosite binding using factor X derivatives that are expected to have limited or abrogated interactions with the active site of VIIa within Xase. Thus, substrate interactions at exosites, sites removed from the active site of VIIa within the enzyme complex, determine affinity and binding specificity in the productive recognition of factor X by the VIIa-TF complex. This may represent a prevalent strategy through which distinctive protein substrate specificities are achieved by the homologous enzymes of coagulation.     

Ab initio study of the reaction mechanism of ribonuclease A with cytidyl-3',5'-adenosine. I. Geometry optimization of cytidyl-3', 5'-adenosine.

As a first part of the ab initio study of the reaction mechanism of ribonuclease A with cytidyl-3',5'-adenosine, the geometry of the cytidyl-3',5'-adenosine substrate has been optimized using the Hartree-Fock method. Eleven different starting structures of cytidyl-3',5'-adenosine have been studied. To guarantee a proper alignment with the active site of the ribonuclease A enzyme, a part of the substrate was fixed during the geometry optimization. The geometry and intramolecular interactions of the refined conformations have been evaluated and two possible prototype structures have been proposed. One of these prototypes is more in accordance with the results of a molecular dynamics simulation and is therefore presented as a model for the geometry of cytidyl-3', 5'-adenosine in the initial step of the reaction with ribonuclease A.     

A theoretical study of glucosamine synthase

Continuing our theoretical studies of glucosamine synthase catalysis, we have carried out MNDO and ab initio calculations of the first stage of the reaction, which involves the attack of a cysteine thiol group from the enzyme active site on the side chain carboxyamide group of glutamine, producing ammonia and thioester. The reactants were modelled by methyl mercaptate and acetamide, respectively. For two considered mechanisms of the reaction the energy surfaces were evaluated. Mechanism I, proposed by Chmara et al. (1985) involves the nucleophilic attack of a deprotonated thiol group on the carbonyl carbon atom. Mechanism II, postulated in our previous work (Tempczyk et al. 1989), assumes the concerted binding of the mercaptate sulphur to the carbonyl carbon and the sulfhydryl hydrogen to the amide nitrogen with simultaneous breaking of the S-H bond. The energy surface of mechanism I shows no minimum on the approach of the mercaptide anion towards the carbonyl carbon, which is also consistent with ab initio calculations in a 4-31 G basis set. Therefore, mechanism I seems to be unlikely. The same analysis of mechanism II shows that it leads to the desired products: methyl thioacetate and ammonia. The presence of a sulfhydryl hydrogen causes apparent pyramidicity of the amido nitrogen and lengthening of the C-N bond in the transition state, making conditions for the release of the ammonia molecule. The MNDO calculated energy barrier of the reaction is 50.1 kcal/mol and the approximate 4-31 G ab initio barrier (at the MNDO geometries of the substrate complex and the transition state) is 63 kcal/mol. The biggest energy contribution to the barrier comes from the breaking of the S-H bond, which also causes a large charge separation in the transition state. The latter affect may result in the stabilisation of the transition state in a real enzymatic environment when compared to the gas phase, e.g. by the interaction of the reacting center with a pair of oppositely charged amino acid side chains such as lysinium and glutamate (aspartate), which are present in the enzyme studied. To estimate the magnitude of this effect, molecular mechanics calculations were carried out on the reaction center at the transition state in our proposed model of the enzymatic active site. The site was supplemented by ammonium and acetate ion, which were to mimic the lysinium and glutamate/aspartate side chains. A transition state stabilization energy of 20 kcal/mol was obtained and this lowers the energy barrier to about 30 kcal/mol. This value is within the thermal energy range of an average protein and indicates that our mechanism is a possible route of glucosamine synthase catalysis. Offprint requests to: E. Borowski     

Coupling interactions of distal residues enhance dihydrofolate reductase catalysis: mutational effects on hydride transfer rates

Recently, the participation in catalysis of residues spatially removed from the enzyme's active site has received considerable attention. The influence of the distal Gly-121 residue on the chemical step of hydride transfer in dihydrofolate reductase (DHFR) catalysis had been demonstrated previously [Cameron, C. E., and Benkovic, S. J. (1997) Biochemistry 36, 15792-15800]. In our continuing effort to identify catalytically important residues that are distal from the active site, we used sequence conservation information, kinetic data on site-directed mutants, dynamic motion information from NMR methods, and correlated motions from MD simulations to identify a subset of residues. Among them, the region spanning positions 41-45 is distal to the active site and was chosen as the focus for the mutagenesis and kinetic studies reported here. Specifically, the highly conserved Met-42 was selected for site-directed mutagenesis. While the reaction kinetics for the M42F mutant enzyme did not deviate from wild-type behavior, a 41-fold reduction in the forward hydride transfer rate was found for the M42W mutant. Given the established role of Gly-121 in the hydride transfer process, double mutant enzymes involving positions 42 and 121 were constructed and characterized. These double mutant enzymes generally showed little changes in substrate and cofactor binding but synergistic decreases in forward hydride transfer rates, while the decreases in reverse rates were additive. Along with supporting information from mixed quantum/classical MD simulations [Agarwal, P. K., Billeter, S. R., Rajagopalan, P. T., Benkovic, S. J., and Hammes-Schiffer, S. (2002) Proc. Natl. Acad. Sci. U.S.A. 99, 2794-2799], the data suggest that a coupled dynamic motion of these distal residues enhances crossing of the chemical reaction barrier and imply an expanded nonstatic role for the protein fold in catalysis.     

A quantum chemical study of the enzymatic deamination of benzoadenine derivatives. A theoretical model of the interactions occurring between nucleosides and the active site of adenosine deaminase

A theoretical study is presented, where, by using both ab initio and semi-empirical methodologies, the properties of benzoadenine derivatives as substrates of adenosine deaminase are discussed. The results suggest that lin-benzoadenine and lin-benzoadenosine can be recognized with an affinity similar to that of adenosine, but only if they are introduced about 0.12 nm deeper inside the active site of the enzyme than the natural substrate adenosine. This fact implies the existence of non-linear hydrogen bonds inside the active site of adenosine deaminase. Ab initio molecular electrostatic potential values suggest that these hydrogen bonds can exist, and have stability similar to that of linear hydrogen bonds. Finally, the great rate of deamination of lin-benzoadenine, comparable with that of adenosine despite the absence of the ribose, is explained in the context of the hypothesis that the protonation at the N1 atom is the rate-determining step of the whole deamination reaction.     

How does activation loop phosphorylation modulate catalytic activity in the cAMP-dependent protein kinase: a theoretical study

Phosphorylation mediates the function of many proteins and enzymes. In the catalytic subunit of cAMP-dependent protein kinase, phosphorylation of Thr 197 in the activation loop strongly influences its catalytic activity. In order to provide theoretical understanding about this important regulatory process, classical molecular dynamics simulations and ab initio QM/MM calculations have been carried out on the wild-type PKA-Mg(2) ATP-substrate complex and its dephosphorylated mutant, T197A. It was found that pThr 197 not only facilitates the phosphoryl transfer reaction by stabilizing the transition state through electrostatic interactions but also strongly affects its essential protein dynamics as well as the active site conformation.     

Proton transfer from histidine 244 may facilitate the 1,2 rearrangement reaction in coenzyme B(12)-dependent methylmalonyl-CoA mutase.

Methylmalonyl-CoA mutase is an adenosylcobalamin-dependent enzyme that catalyzes the 1,2 rearrangement of methylmalonyl-CoA to succinyl-CoA. This reaction results in the interchange of a carbonyl-CoA group and a hydrogen atom on vicinal carbons. The crystal structure of the enzyme reveals the presence of an aromatic cluster of residues in the active site that includes His-244, Tyr-243, and Tyr-89 in the large subunit. Of these, His-244 is within hydrogen bonding distance to the carbonyl oxygen of the carbonyl-CoA moiety of the substrate. The location of these aromatic residues suggests a possible role for them in catalysis either in radical stabilization and/or by direct participation in one or more steps in the reaction. The mechanism by which the initially formed substrate radical isomerizes to the product radical during the rearrangement of methylmalonyl-CoA to succinyl-CoA is unknown. Ab initio molecular orbital theory calculations predict that partial proton transfer can contribute significantly to the lowering of the barrier for the rearrangement reaction. In this study, we report the kinetic characterization of the H244G mutant, which results in an acute sensitivity of the enzyme to oxygen, indicating the important role of this residue in radical stabilization. Mutation of His-244 leads to an approximately 300-fold lowering in the catalytic efficiency of the enzyme and loss of one of the two titratable pK(a) values that govern the activity of the wild type enzyme. These data suggest that protonation of His-244 increases the reaction rate in wild type enzyme and provides experimental support for ab initio molecular orbital theory calculations that predict rate enhancement of the rearrangement reaction by the interaction of the migrating group with a general acid. However, the magnitude of the rate enhancement is significantly lower than that predicted by the theoretical studies.     

Structure and dynamics of pin1 during catalysis by NMR

The link between internal enzyme motions and catalysis is poorly understood. Correlated motions in the microsecond-to-millisecond timescale may be critical for enzyme function. We have characterized the backbone dynamics of the peptidylprolyl isomerase (Pin1) catalytic domain in the free state and during catalysis. Pin1 is a prolyl isomerase of the parvulin family and specifically catalyzes the isomerization of phosphorylated Ser/Thr-Pro peptide bonds. Pin1 has been shown to be essential for cell-cycle progression and to interact with the neuronal tau protein inhibiting its aggregation into fibrillar tangles as found in Alzheimer's disease. (15)N relaxation dispersion measurements performed on Pin1 during catalysis reveal conformational exchange processes in the microsecond timescale. A subset of active site residues undergo kinetically similar exchange processes even in the absence of a substrate, suggesting that this area is already "primed" for catalysis. Furthermore, structural data of the turning-over enzyme were obtained through inter- and intramolecular nuclear Overhauser enhancements. This analysis together with a characterization of the substrate concentration dependence of the conformational exchange allowed the distinguishing of regions of the enzyme active site that are affected primarily by substrate binding versus substrate isomerization. Together these data suggest a model for the reaction trajectory of Pin1 catalysis.     

Determinants of reactivity and selectivity in soluble epoxide hydrolase from quantum mechanics/molecular mechanics modeling

Soluble epoxide hydrolase (sEH) is an enzyme involved in drug metabolism that catalyzes the hydrolysis of epoxides to form their corresponding diols. sEH has a broad substrate range and shows high regio- and enantioselectivity for nucleophilic ring opening by Asp333. Epoxide hydrolases therefore have potential synthetic applications. We have used combined quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations (at the AM1/CHARMM22 level) and high-level ab initio (SCS-MP2) QM/MM calculations to analyze the reactions, and determinants of selectivity, for two substrates: trans-stilbene oxide (t-SO) and trans-diphenylpropene oxide (t-DPPO). The calculated free energy barriers from the QM/MM (AM1/CHARMM22) umbrella sampling MD simulations show a lower barrier for phenyl attack in t-DPPO, compared with that for benzylic attack, in agreement with experiment. Activation barriers in agreement with experimental rate constants are obtained only with the highest level of QM theory (SCS-MP2) used. Our results show that the selectivity of the ring-opening reaction is influenced by several factors, including proximity to the nucleophile, electronic stabilization of the transition state, and hydrogen bonding to two active site tyrosine residues. The protonation state of His523 during nucleophilic attack has also been investigated, and our results show that the protonated form is most consistent with experimental findings. The work presented here illustrates how determinants of selectivity can be identified from QM/MM simulations. These insights may also provide useful information for the design of novel catalysts for use in the synthesis of enantiopure compounds.     

Characterization of the active site of DNA polymerase beta by molecular dynamics and quantum chemical calculation

It is well established that the fully formed polymerase active site of the DNA repair enzyme, polymerase beta (pol beta), including two bound Mg2+ cations and the nucleoside triphosphate (dNTP) substrate, exists at only one point in the catalytic cycle just prior to the chemical nucleotidyl transfer step. The structure of the active conformation has been the subject of much interest as it relates to the mechanism of the chemical step and also to the question of fidelity assurance. Although crystal structures of ternary pol beta-(primer-template) DNA-dNTP complexes have provided the main structural features of the active site, they are necessarily incomplete due to intentional alterations (e.g., removal of the 3'OH groups from primer and substrate) needed to obtain a structure from midcycle. Working from the crystal structure closest to the fully formed active site [Protein Data Bank (PDB) code: 1bpy], two molecular dynamics (MD) simulations of the solvated ternary complex were performed: one with the missing 3'OHs restored, via modeling, to the primer and substrate, and the other without restoration of the 3'OHs. The results of the simulations, together with ab initio optimizations on simplified active-site models, indicate that the missing primer 3'OH in the crystal structure is responsible for a significant perturbation in the coordination sphere of the catalytic cation and allow us to suggest several corrections and additions to the active-site structure as observed by crystallography. In addition, the calculations help to resolve questions raised regarding the protonation states of coordinating ligands.     

Domain motions of hyaluronan lyase underlying processive hyaluronan translocation

Hyaluronan lyase (Hyal) is a surface enzyme occurring in many bacterial organisms including members of Streptococcus species. Streptococcal Hyal primarily degrades hyaluronan‐substrate (HA) of the extracellular matrix. This degradation appears to facilitate the spread of this bacterium throughout host tissues. Unlike purely endolytic degradation of its other substrates, unsulfated chondroitin or some chondroitin sulfates, the degradation of HA by Hyal proceeds by processive exolytic cleavage of one disaccharide at a time following an initial endolytic cut. Molecular dynamics (MD) studies of Hyal from Streptococcus pneumoniae are presented that address the enzyme's molecular mechanism of action and the role of domain motions for processive functionality. The analysis of extensive sub‐microsecond MD simulations of this enzyme action on HA‐substrates of different lengths and the connection between the domain dynamics of Hyal and the translocation of the HA‐substrate reveals that opening/closing and twisting domain motions of the Hyal are intimately linked to processive HA degradation. Enforced simulations confirmed this finding as the domain motions in SpnHyal were found to be induced by enforced substrate translocation. These results establish the dynamic interplay between Hyal flexibility and substrate translocation and provide insight into the processive mechanism of Hyal. Proteins 2009. © 2008 Wiley‐Liss, Inc.