Jitrayut Jitonnom

Quantum Mechanics/Molecular Mechanics Modeling of Substrate-Assisted Catalysis in Family 18 Chitinases: Conformational Changes and the Role of Asp142 in Catalysis in ChiB

Family 18 chitinases catalyze the hydrolysis of β-1,4-glycosidic bonds in chitin. The mechanism has been proposed to involve the formation of an oxazolinium ion intermediate via an unusual substrate-assisted mechanism, in which the substrate itself acts as an intramolecular nucleophile (instead of an enzyme residue). Here, we have modeled the first step of the chitin hydrolysis catalyzed by Serratia marcescens chitinase B for the first time using a combined quantum mechanics/molecular mechanics approach. The calculated reaction barriers based on multiple snapshots are 15.8-19.8 kcal mol(-1) [B3LYP/6-31+G(d)//AM1-CHARMM22], in good agreement with the activation free energy of 16.1 kcal mol(-1) derived from experiment. The enzyme significantly stabilizes the oxazolinium intermediate. Two stable conformations ((4)C(1)-chair and B(3,O)-boat) of the oxazolinium ion intermediate in subsite -1 were unexpectedly observed. The transition state structure has significant oxacarbenium ion-like character. The glycosyl residue in subsite -1 was found to follow a complex conformational pathway during the reaction ((1,4)B → [(4)H(5)/(4)E](++) → (4)C(1) ↔ B(3,O)), indicating complex conformational behavior in glycoside hydrolases that utilize a substrate-assisted catalytic mechanism. The D142N mutant is found to follow the same wild-type-like mechanism: the calculated barriers for reaction in this mutant (16.0-21.1 kcal mol(-1)) are higher than in the wild type, in agreement with the experiment. Asp142 is found to be important in transition state and intermediate stabilization.

QM/MM free-energy simulations of reaction in Serratia marcescens Chitinase B reveal the protonation state of Asp142 and the critical role of Tyr214.

Serratia marcescens Chitinase B (ChiB), belonging to the glycosidase family 18 (GH18), catalyzes the hydrolysis of β-1,4-glycosidic bond, with retention of configuration, via an unusual substrate-assisted mechanism, in which the substrate itself acts as an intramolecular nucleophile. Here, both elementary steps (glycosylation and deglycosylation) of the ChiB-catalyzed reaction are investigated by means of combined quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations at the SCC-DFTB/CHARMM22 level of theory. We examine the influence of the Asp142 protonation state on the reaction and the role that this residue performs in the reaction. Our simulations show that reaction with a neutral Asp142 is preferred and demonstrate that this residue provides electrostatic stabilization of the oxazolinium ion intermediate formed in the reaction. Insight into the conformational itinerary ((1,4)B↔(4)H5↔(4)C1) adopted by the substrate (bound in subsite -1) along the preferred reaction pathway is also provided by the simulations. The relative energies of the stationary points found along the reaction pathway calculated with SCC-DFTB and B3LYP were compared. The results suggest that SCC-DFTB is an accurate method for estimating the relative barriers for both steps of the reaction; however, it was found to overestimate the relative energy of an intermediate formed in the reaction when compared with the higher level of theory. Glycosylation is suggested to be a rate-determining step in the reaction with calculated overall reaction free-energy barrier of 20.5 kcal/mol, in a reasonable agreement with the 16.1 kcal/mol barrier derived from the experiment. The role of Tyr214 in catalysis was also investigated with the results, indicating that the residue plays a critical role in the deglycosylation step of the reaction. Simulations of the enzyme-product complex were also performed with an unbinding event suggested to have been observed, affording potential new mechanistic insight into the release of the product of ChiB.

Pairwise decomposition of residue interaction energies of single chain Fv with HIV-1 p17 epitope variants.

Computational assisted modeling was carried out to investigate the importance of specific residues in the binding site of scFv. In this study, scFv against HIV-1 epitope at the C-terminal on p17 (scFv anti-p17) was used as a candidate molecule for evaluating the method. The wild-type p17 and its nine natural mutants were docked with scFv anti-p17. Potential mean force (PMF) scores predicted the most favorable binding interaction, and the correlation agreed well with the corresponding activity data from the peptide based ELISA. In the interaction with solvent molecules, the 3D structures of scFv anti-p17 and selected peptide epitopes were further investigated by molecular dynamics (MDs) simulation with the AMBER 9 program. Post-processing of the snapshot at equilibrium was performed to evaluate the binding free energy and pairwise decomposition or residue-based energy calculation of complexes in solution using the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) protocol. Our results demonstrated that the specific residues located in the complementary determining regions (CDRs) of scFv anti-p17, MET100, LYS101, ASN169, HIS228, and LEU229, play a crucial role in the effective binding interaction with the absolute relative decomposed energy more than 2.00 kcal/mol in comparison to the original substrate.

Theoretical Investigation on the Electronic and Optical Properties of Poly(fluorenevinylene) Derivatives as Light-Emitting Materials

Density functional theory (DFT) and time-dependent DFT (TDDFT) were employed to study ground-state properties, HOMO-LUMO gaps , excitation energies , ionization potentials (IPs), and electron affinities (EA) for PFV-alt-PDONV and PFV-alt-PDIH-PPV having different alternating groups. Excited-state properties were investigated using configuration interaction singles (CISs) while fluorescence energies were calculated using TDDFT. The results show that PFV-alt-PDONV exhibits blue-shifted energies for both HOMO-LUMO gaps and excitation energies compared with PFV-alt-PDIH-PPV. The predicted IP and EA clearly indicate that PFV-alt-PDIH-PPV has both easier hole creation and electron injection than that of PFV-alt-PDONV. The maximal absorption wavelengths of all polymers are strongly assigned to transition. The predicted radiative lifetimes of PFV-alt-PDONV and PFV-alt-PDIH-PPV for B3LYP/6-31G(d) are 0.36 and 0.61 ns, respectively, indicating that PFV-alt-PDIH-PPV should have a better performance for long-time emission than that of PFV-alt-PDONV.

A DFT study of the unusual substrate-assisted mechanism of Serratia marcescens chitinase B reveals the role of solvent and mutational effect on catalysis

Serratia marcescens chitinase B (SmChiB) catalyzes the hydrolysis of β-1,4-glycosidic bond, via an unusual substrate-assisted mechanism, in which the substrate itself acts as an intramolecular nucleophile. In this paper, the catalytic mechanism of SmChiB has been investigated by using density functional theory. The details of two consecutive steps (glycosylation and deglycosylation), the structures and energetics along the whole catalytic reaction, and the roles of solvent molecules as well as some conserved SmChiB residues (Asp142, Tyr214, Asp215, and Arg294) during catalysis are highlighted. Our calculations show that the formation of the oxazolinium cation intermediate in the glycosylation step was found to be a rate-determining step (with a barrier of 23 kcal/mol), in line with our previous computational studies (Jitonnom et al., 2011, 2014). The solvent water molecules have a significant effect on a catalytic efficiency in the degycosylation step: the catalytic water is essentially placed in a perfect position for nucleophic attack by hydrogen bond network, lowering the barrier height of this step from 11.3 kcal/mol to 2.9 kcal/mol when more water molecules were introduced. Upon the in silico mutations of the four conserved residues, their mutational effects on the relative stability of the reaction intermediates and the computed energetics can be obtained by comparing with the wild-type results. Mutations of Tyr214 to Phe or Ala have shown a profound effect on the relative stability of the oxazolinium intermediate, emphasizing a direct role of this residue in destabilizing the intermediate. In line with the experiment that the D142A mutation leads to almost complete loss of SmChiB activity, this mutation greatly decreases the stability of the intermediate, resulting in a very large increase in the activation barrier up to 50 kcal/mol. The salt-bridges residues (Asp215 and Arg294) were also found to play a role in stabilizing the oxazolinium intermediate.

Computational design of peptide inhibitor based on modifications of proregion from Plutella xylostella midgut trypsin.

Many proteases are produced as zymogens bearing the N‐terminal proregions acting both as intramolecular chaperones and as protease inhibitors. The latter role of the proregions as potent and specific inhibitors of their associated protease has been demonstrated in various peptidases and therefore has been targeted for alternative pest control. Here, we isolated amino acid sequence of Plutella xylostella midgut trypsin from the larvae of diamondback moth and tested in silico for its inhibitory activity toward propeptide models using computational modeling and docking. The propeptide models (AAAPGHR, AAAPGRR, AAAPGKR, AAPGHRI, APGHRIV, PGHRIVG, AAAAPGH, and AAAAAPG) were designed based on histidine‐mutated and frame‐shifted modifications of the 7‐amino‐acid proregion (AAAPGHR) of the Plutella xylostella trypsin. Among the eight peptides, AAAPGRR was found to give the best docking scores, showing a strong binding to the cognate enzyme. In addition, the obtained structure of trypsin–AAAPGRR complex was found to share a similar binding mode with a crystal structure of plant protease inhibitor complex. Our results may guide the experiment for the design of future peptide inhibitor with specificity and selectivity for the target enzyme.

Cationic ring-opening polymerization of cyclic carbonates and lactones by group 4 metallocenes: A theoretical study on mechanism and ring-strain effects

Group 4 metallocene-mediated cationic ring-opening polymerizations of a series of lactones and cyclic carbonates, with different ring sizes (n=4–8) have been theoretically studied. Using the “naked cation” approach in combination with density functional theory, the activated chain-end mechanism and the influence of transition metals, solvent and monomer ring size on the polymerizability were explored in detail. The results showed that the cationic metallocene–monomer complex, [catalyst][monomer]+, is formed, generating cationic (carbocation ion) species responsible for polymer chain growth. We found that poor polymerizability of five-membered lactone and six-membered ring carbonate depends not only on the nature of the monomer ring size but also the relative stability of the complex, which was found to correlate well with the ring strain. Subsequently, several propagation steps take place through an SN2 reaction which involves ring opening of an active monomer, via alkyl–oxygen bond cleavage. Based on the...

Comparative study on activation mechanism of carboxypeptidase A1, A2 and B: first insights from steered molecular dynamics simulations.

Different forms of procarboxypeptidases (PCPs) zymogens are observed experimentally to show different rates and modes of activation: PCPA1 shows a slow, biphasic, activation pathway compared to PCPA2 and PCPB which have a faster, monotonic activation behavior. Detailed mechanisms involved in activating these zymogen forms to the active enzymes are not well understood yet. In this work, three PCP zymogens (subtypes A1, A2 and B) were in silico converted into the primary cleavage state of zymogens using available X-ray structures. Based on these cleaved forms of zymogen, we are able to investigate their spontaneous dissociation process of the prosegment from its associated enzyme domain using steered molecular dynamics simulation. The simulations revealed the highest rupture force in PCPB followed by PCPA2 and PCPA1. We also found that the cleavage substantially destabilizes most of the hydrogen bonds at the prosegment-enzyme interface in each zymogen structure. The mechanisms of the prosegment unbinding seem to be similar in both PCPA1 and PCPB but different in PCPA2: PCPA1 and PCPB show first rupture in the connecting segment followed by the globular domain, while PCPA2 conversely shows first rupture in the globular domain and then in the connecting segment. Our simulations have included the dynamic and long range conformational effects taking place after the first proteolytic cleavage in PCPs, providing first insights into the activation of carboxypeptidase A1, A2 and B.

Influence of metal cofactors and water on the catalytic mechanism of creatininase-creatinine in aqueous solution from molecular dynamics simulation and quantum study.

The reaction mechanism of creatinine-creatininase binding to form creatine as a final product has been investigated by using a combined ab initio quantum mechanical/molecular mechanical approach and classical molecular dynamics (MD) simulations. In MD simulations, an X-ray crystal structure of the creatininase/creatinine was modified for creatininase/creatinine complexes and the MD simulations were run for free creatininase and creatinine in water. MD results reveal that two X-ray water molecules can be retained in the active site as catalytic water. The binding free energy from Molecular Mechanics Poisson-Boltzmann Surface Area calculation predicted the strong binding of creatinine with Zn2+, Asp45 and Glu183. Two step mechanisms via Mn2+/Zn2+ (as in X-ray structure) and Zn2+/Zn2+ were proposed for water adding step and ring opening step with two catalytic waters. The pathway using synchronous transit methods with local density approximations with PWC functional for the fragment in the active region were obtained. Preferable pathway Zn2+/Zn2+ was observed due to lower activation energy in water adding step. The calculated energy in the second step for both systems were comparable with the barrier of 26.03 and 24.44xa0kcal/mol for Mn2+/Zn2+ and Zn2+/Zn2+, respectively.

A computational H5N1 neuraminidase model and its binding to commercial drugs

In order to understand the mechanisms of ligand binding and interaction between two commercial drugs (ligands), zanamivir and oseltamivir and H5N1 Influenza Virus Neuraminidase subtype N1, a three-dimensional model of N1-ligand (GenBank accession no. AAS654617) was initially generated by homology modeling using the 13 high-resolution X-ray structures of neuraminidase N2 and N9 as the template. With the aid of the molecular mechanics and molecular dynamics methods, the final implicit solvent refined model was obtained. It was, then, assessed by PROCHECK, PROSA and VERIFY3D. With this model, a flexible docking study was performed. The results show strong hydrogen bond interactions between the glycerol side chains of zanamivir and Arg29 of the N1. Common hydrogen bonds between the carboxyl groups and Arg279 were found for both drugs. It was also found that the Glu30, Asp62, Arg63, Arg204, Trp310, Tyr313, Glu336, Ile338, Trp348, Ala349 were observed to facilitate the enzyme-ligand non-bonding interactions as they are located within the radius of 5 Å from all atoms of both drugs. Charge distribution was evaluated using the semi-empirical AM1 method. The results show that the total net charges of the –NH side chain of zanamivir is less negative than that of oseltamivir. This is in contrast to what is observed for the amide and alkyl (ether/glycerol) side chains. In comparison of the binding free energies between the X-ray N2-ligand and N9-ligand complexes, N1-ligand binding is found to be less potent than N2 and N9 subtypes, while N2-ligand and N9-ligand are roughly comparable. In addition, it is interesting to observe that the binding free energies for all three subtypes of the zanamivir complexes are lower than those of oseltamivir.