Effective fragment potentials and the enzyme active site

Optimization of the binding conformation of a substrate in an enzyme active site using ab initio quantum chemistry methods are intractable since the active site comprises several hundred atoms. However, the active site can be decomposed into an active and spectator region where the spectator residues are represented by effective fragment potentials and reducing the number of all-electron atoms involved in the chemistry to a reasonable level. The effective fragment potentials for electrostatics and polarization are implemented in GAMESS but the repulsive and charge transfer potentials are fit to interaction energies of water with models of the residues. These repulsive/charge transfer potentials are generated for the protein residues and the EFP are then used to optimize binding of a transition state analogue to chorismate mutase (B. subtilis) and small dianions to ribonuclease A. For chorismate mutase the calculated binding conformation compares well to the comparable X-ray structure. The binding of the inhibitor to the glutamate/glutamine mutant active site is then predicted with the optimization including the glutamine residue constrained only at the C alpha atom. The binding conformations suggest important roles for tyr108 and arg63, which have not been noted earlier. The electrostatic stabilization of the transition state by the active site charge distribution has to be augmented by a specific electronic activation by glu78. In ribonuclease A, the protons are found to move to provide a clustering of the charges to bind the small dianions, phosphate, thiophosphate, and sulfate.     

Computational docking,molecular dynamics simulation and subsite structure analysis of a maltogenic amylase from Bacillus lehensis G1 provide insights into substrate and product specificity

Maltogenic amylase (MAG1) from Bacillus lehensis G1 displayed the highest hydrolysis activity on β-cyclodextrin (β-CD) to produce maltose as a main product and exhibited high transglycosylation activity on malto-oligosaccharides with polymerization degree of three and above. These substrate and product specificities of MAG1 were elucidated from structural point of view in this study. A three-dimensional structure of MAG1 was constructed using homology modeling. Docking of β-CD and malto-oligosaccharides was then performed in the MAG1 active site. An aromatic platform in the active site was identified which is responsible in substrate recognition especially in determining the enzyme’s preference toward β-CD. Molecular dynamics (MD) simulation showed MAG1 structure is most stable when docked with β-CD and least stable when docked with maltose. The docking analysis and MD simulation showed that the main subsites for substrate stabilization in the active site are −2, −1, +1 and +2. A bulky residue, Trp359 at the +2 subsite was identified to cause steric interference to the bound linear malto-oligosaccharides thus prevented it to occupy subsite +3, which can only be reached by a highly bent glucose molecule such as β-CD. The resulted modes of binding from docking simulation show a good correlation with the experimentally determined hydrolysis pattern. The subsite structure generated from this study led to a possible mode of action that revealed how maltose was mainly produced during hydrolysis. Furthermore, maltose only occupies subsite +1 and +2, therefore could not be hydrolyzed or transglycosylated by the enzyme. This important knowledge has paved the way for a novel structure-based molecular design for modulation of its catalytic activities.     

2-(4-甲氧基苯基)亚肼基乙酸与乙醛酸反应密度泛函理论研究

本文应用量子化学密度泛函理论(DFT)研究了2-(4-甲氧基苯基)亚肼基乙酸与乙醛酸的反应机理及动力学行为,并研究了不同取代基对反应的影响.在RB3LYP/6-31G*水平上,优化了各反应通道的反应物、中间体、过渡态及产物的几何构型,计算了各驻点的振动频率、零点能和电荷分布.计算结果表明,该反应有2条反应通道,分别生成(2Z)和(2E)-3(4-甲氧基苯基二氮烯基)丙烯酸.产物发生了键长的平均化和电荷的重新分布.两反应通道具有相同的反应入口,Z式产物为主要产物,2-(4-甲氧基苯基)亚肼基乙酸中苯环对位被给电子基团取代,有利于反应进行.     

A molecular dynamics study of the BACE1 conformational change from Apo to closed form induced by hydroxyethylamine derived compounds

BACE1 is an aspartyl protease which is a therapeutic target for Alzheimer’s disease (AD) because of its participation in the rate-limiting step in the production of Aβ-peptide, the accumulation of which produces senile plaques and, in turn, the neurodegenerative effects associated with AD. The active site of this protease is composed in part by two aspartic residues (Asp93 and Asp289). Additionally, the catalytic site has been found to be covered by an antiparallel hairpin loop called the flap. The dynamics of this flap are fundamental to the catalytic function of the enzyme. When BACE1 is inactive (Apo), the flap adopts an open conformation, allowing a substrate or inhibitor to access the active site. Subsequent interaction with the ligand induces flap closure and the stabilization of the macromolecular complex. Further, the protonation state of the aspartic dyad is affected by the chemical nature of the species entering the active site, so that appropriate selection of protonation states for the ligand and the catalytic residues will permit the elucidation of the inhibitory pathway for BACE1. In the present study, comparative analysis of different combinations of protonation states for the BACE1-hydroxyethylamine (HEA) system is reported. HEAs are potent inhibitors of BACE1 with favorable pharmacological and kinetic properties, as well as oral bioavailability. The results of Molecular Dynamics (MD) simulations and population density calculations using 8 different parameters demonstrate that the LnAsp289 configuration (HEA with a neutral amine and the Asp289 residue protonated) is the only one which permits the expected conformational change in BACE1, from apo to closed form, after flap closure. Additionally, differences in their capacities to establish and maintain interactions with residues such as Asp93, Gly95, Thr133, Asp289, Gly291, and Asn294 during this step allow differentiation among the inhibitory activities of the HEAs. The results and methodology here reported will serve to elucidate the inhibitory pathway of other families of compounds that act as BACE1 inhibitors, as well as the design of better leader compounds for the treatment of AD.     

Ab initio QM/MM modelling of acetyl-CoA deprotonation in the enzyme citrate synthase

The first step of the reaction catalysed by the enzyme citrate synthase is studied here with high level combined quantum mechanical/molecular mechanical (QM/MM) methods (up to the MP2/6-31+G(d)//6-31G(d)/CHARMM level). In the first step of the reaction, acetyl-CoA is deprotonated by Asp375, producing an intermediate, which is the nucleophile for attack on the second substrate, oxaloacetate, prior to hydrolysis of the thioester bond of acetyl-CoA and release of the products. A central question has been whether the nucleophilic intermediate is the enolate of acetyl-CoA, the enol, or an 'enolic' intermediate stabilized by a 'low-barrier' hydrogen bond with His274 at the active site. The imidazole sidechain of His274 is neutral, and donates a hydrogen bond to the carbonyl oxygen of acetyl-CoA in substrate complexes. We have investigated the identity of the nucleophilic intermediate by QM/MM calculations on the substrate (keto), enolate, enol and enolic forms of acetyl-CoA at the active site of citrate synthase. The transition states for proton abstraction from acetyl-CoA by Asp375, and for transfer of the hydrogen bonded proton between His274 and acetyl-CoA have been modelled approximately. The effects of electron correlation are included by MP2/6-31G(d) and MP2/6-31+G(d) calculations on active site geometries produced by QM/MM energy minimization. The results do not support the hypothesis that a low-barrier hydrogen bond is involved in catalysis in citrate synthase, in agreement with earlier calculations. The acetyl-CoA enolate is identified as the only intermediate consistent with the experimental barrier for condensation, stabilized by conventional hydrogen bonds from His274 and a water molecule.     

Mechanistic insights into the catalytic reaction of plant allene oxide synthase (pAOS) via QM and QM/MM calculations

QM cluster and QM/MM protein models have been employed to understand aspects of the reaction mechanism of plant allene oxide synthase (pAOS). In this study we have investigated two reaction mechanisms for pAOS. The standard pAOS mechanism was contrasted with an alternative involving an additional active site molecule which has been shown to facilitate proton coupled electron transfer (PCET) in related systems. Firstly, we found that the results from QM/MM protein model are comparable with those from the QM cluster model, presumably due to the large active site used. Furthermore, the results from the QM cluster model show that the Fe^III and Fe^IV pathways for the standard mechanism have similar energetic and structural properties, indicating that the reaction mechanism may well proceed via both pathways. However, while the PCET process is facilitated by an additional active site bound water in other related families, in pAOS it is not, suggesting this type of process is not general to all closely related family members.     

The three-dimensional structure of Arabidopsis thaliana O-methyltransferase predicted by homology-based modelling

O-methylation of flavonoid compounds is an important enzymatic reaction since it not only reduces the chemical reactivity of their phenolic hydroxyl groups but also increases their lipophilicity and, hence, their intracellular compartmentation. Several genes encoding flavonoid O-methyltransferases (OMTs) have been isolated and characterized both at the molecular and biochemical levels. In contrast with mammalian enzymes, plant OMTs exhibit narrow substrate specificities as well as position-specific activities, so that the homology comparison, derived using programs such as BLAST can not provide sufficient information on the enzyme function or its substrate preference. In order to study these characteristics, therefore, another approach, homology-based modelling is being carried out. We report here the determination of the 3-D structure of Arabidopsis thaliana O-methyltransferase, AtOMT1 as well as its dynamics when complexed with its substrate. The predicted structure obtained by homology-based modelling is conserved during molecular dynamics simulations. AtOMT1 exhibits a structure similar to that of caffeic acid O-methyltransferase, COMT when the latter was used as a template. Whereas COMT includes 20 alpha-helices and nine beta-sheets, AtOMT1 has 16 and 9, respectively. Although the homology between both proteins is higher than 77% and all amino acids surrounding the active sites, except one residue, are similar in their primary sequences, the two proteins exhibit different substrate preferences. The differences in substrate specificity may be explained on the basis of the predicted structures of the protein and its complex with the substrate. In addition, docking the substrate into the active site of the protein allowed the study of the structural change of the active site on the dihedral angle distribution of the residues surrounding the active site.     

External and internal electrostatic potentials of cholinesterase models

The electrostatic potentials for the three-dimensional structures of cholinesterases from various species were calculated, using the Delphi algorithm, on the basis of the Poisson–Boltzmann equation. We used structures for Torpedo californica and mouse acetylcholinesterase, and built homology models of the human, Bungarus fasciatus, and Drosophila melanogaster acetylcholinesterases and human butyrylcholinesterase. All these structures reveal a negative external surface potential, in the area around the entrance to the active-site gorge, that becomes more negative as the rim of the gorge is approached. Moreover, in all cases, the potential becomes increasingly more negative along the central axis running down the gorge, and is largest at the base of the gorge, near the active site. Ten key acidic residues conserved in the sequence alignments of AChE from various species, both in the surface area near the entrance of the active-site gorge and at its base, appear to be primarily responsible for these potentials. The potentials are highly correlated among the structures examined, down to sequence identities as low as 35%. This indicates that they are a conserved property of the cholinesterase family, could serve to attract the positively charged substrate into and down the gorge to the active site, and may play other roles important for cholinesterase function.     

Combined quantum mechanical and molecular mechanical reaction pathway calculation for aromatic hydroxylation by p-hydroxybenzoate-3-hydroxylase

The reaction pathway for the aromatic 3-hydroxylation of p-hydroxybenzoate by the reactive C4a-hydroperoxyflavin cofactor intermediate in p-hydroxybenzoate hydroxylase (PHBH) has been investigated by a combined quantum mechanical and molecular mechanical (QM/MM) method. A structural model for the C4a-hydroperoxyflavin intermediate in the PHBH reaction cycle was built on the basis of the crystal structure coordinates of the enzyme-substrate complex. A reaction pathway for the subsequent hydroxylation step was calculated by imposing a reaction coordinate that involves cleavage of the peroxide oxygen-oxygen bond and formation of the carbon-oxygen bond between the C3 atom of the substrate and the distal oxygen of the peroxide moiety of the cofactor. The geometric changes and the Mulliken charge distributions along the calculated reaction pathway are in line with an electrophilic aromatic substitution type of mechanism. The energy barrier of the calculated reaction is considerably lower when the substrate hydroxyl moiety is deprotonated, in comparison with the barrier found with a protonated hydroxyl moiety. This effect of the protonation state of the substrate on the calculated energy barrier supports experimental observations that deprotonation is required for hydroxylation of the substrate. A notable event in the calculated reaction pathway is a lengthening of the peroxide oxygen-oxygen bond at an intermediate stage. Further analysis of the reaction pathway indicates that this oxygen-oxygen bond elongation is accompanied by an increase in electrophilic reactivity on the distal oxygen of the peroxide moiety, which may assist the C-O bond formation in the reaction of the C4a-hydroperoxyflavin intermediate with the substrate. Analysis of the effect of individual active site residues on the reaction reveals a specific transition state stabilization by the backbone carbonyl moiety of Pro293. The crystal water 717 appears to drive the hydroxylation step through a stabilizing hydrogen bond interaction to the proximal oxygen of the C4a-hydroperoxyflavin intermediate, which increases in strength as the hydroperoxyflavin cofactor converts to the anionic (deprotonated) hydroxyflavin.     

Effects of amine organic groups as lattice in ZSM-5 on the hydrolysis of dimethyl ether

The effects of doping amine to ZSM-5 on its catalytic activity for hydrolysis of dimethyl ether (DME) have been studied theoretically using Density Functional Theory with the embedded cluster ONIOM(M06/6-31G(d,p):UFF) model. Doping by amine to ZSM-5 yields two new active centers, namely the protonated Z[NH₂] and non-protonated Z[NH] amine sites in addition to the normal Brønsted acid Z[OH] site. The reaction has two possible stepwise and concerted channels. The stepwise channel consists of two elementary steps; (i) the demethylation followed by (ii) the hydrolysis while the concerted channel involves in the demethylation and hydrolysis in a single step. We found that the reaction favors to proceed via the concerted channel at all three active centers. The results predict that the Z[OH] shows the best catalytic performance for the studied reaction. The Z[NH₂] is not catalytically active due to the activation barriers are extremely high for both stepwise and concerted pathways. The demethylation step is energetically favorable over the Z[NH] site, however, the product methylamonium surface intermediate is too stable to be further converted to methanol.     

Explaining the autoinhibition of the SMYD enzyme family: A theoretical study

The SMYD enzymes (SMYD1-5) are lysine methyltransferases that have diverse biological functions including gene expression and regulation of skeletal and cardiac muscle development and function. Recently, they have gained more attention as potential drug targets because of their involvement in cardiovascular diseases and in the progression of different cancer types. Their activity has been suggested to be regulated by a posttranslational mechanism and by autoinhibition. The later relies on a hinge-like movement of the N- and C-lobes to adopt an open or closed conformation, consequently, determining the accessibility of the active site and substrate specificity.In this study we aim to investigate and explain the possibility of the regulatory autoinhibition process of the SMYD enzymes by a thorough computational exploration of their dynamic, energetic, and structural changes by using extended molecular dynamics simulations; normal mode analysis (NMA); and energy correlations. Three SMYD models (SMYD1-3) were used in this study. Our results showed an obvious hinge-like motion between the N- and C-lobes. Also, we identified interaction energy pathways within the 3D structures of the proteins, and hot spots on their surfaces that could be of particular importance for the regulation of their activities via allosteric means. These results can help in a better understanding of the nature of these promising drug targets; and in designing selective drugs that can interfere with (inhibit) the function of a specific SMYD member by disrupting its dynamical and conformational behaviour without disrupting the function of the entire SMYD proteins.     

Flap opening mechanism of HIV-1 protease

The active site of aspartic proteases, such as HIV-1 protease (PR), is covered by one or more flaps, which restrict access to the active site. For HIV-1 PR, X-ray diffraction studies suggested that in the free enzyme the two flaps are packed onto each other loosely in a semi-open conformation, while molecular dynamics (MD) studies observed that the flaps can also separate into open conformations. In this study, the mechanism of flap opening and the structure and dynamics of HIV-1 PR with semi-open and open flap conformations were investigated using molecular dynamics simulations. The flaps showed complex dynamic behavior as two distinct mechanisms of flap opening and various stable flap conformations (semi-open, open and curled) were observed during the simulations. A network of weakly polar interactions between the flaps were proposed to be responsible for stabilizing the semi-open flap conformation. It is hypothesized that such interactions could be responsible for making flap opening a highly sensitive gating mechanism which control access to the active site.     

Mutational analysis of microbial hydroxycinnamoyl-CoA hydratase-lyase (HCHL) towards enhancement of binding affinity: A computational approach

Improving the industrial enzyme for better yield of the product is important and a challenging task. One of such important industrial enzymes is microbial Hydroxycinnamoyl-CoA hydratase-lyase (HCHL). It converts feruloyl-CoA to vanillin. We place our efforts towards the improvement of its catalytic activity with comprehensive computational investigation. Catalytic core of the HCHL was explored with molecular modeling and docking approaches. Site-directed mutations were introduced in the catalytic site of HCHL in a sequential manner to generate different mutants of HCHL. Basis of mutation is to increase the interaction between HCHL and substrate feruloyl-CoA through interatomic forces and hydrogen bond formation. A rigorous molecular dynamics (MD) simulation was performed to check the stability of mutant’s structure. Root mean square deviation (RMSD), root mean square fluctuation (RMSF), dynamic cross correlation (DCCM) and principal component analysis (PCA) were also performed to analyze flexibility and stability of structures. Docking studies were carried out between different mutants of HCHL and feruloyl-CoA. Investigation of the different binding sites and the interactions with mutant HCHLs and substrate allowed us to highlight the improved performance of mutants than wild type HCHL. This was further validated with MD simulation of complex consisting of different mutants and substrate. It further confirms all the structures are stable. However, mutant-2 showed better affinity towards substrate by forming hydrogen bond between active site and feruloyl-CoA. We propose that increase in hydrogen bond formation might facilitate in dissociation of vanillin from feruloyl-CoA. The current work may be useful for the future development of ‘tailor-made’ enzymes for better yield of vanillin.     

基于时空网络的地铁进出站客流量预测

客流量预测是城市智能交通系统的重要组成部分,对人们出行和交通管控有着重要的指导意义。针对地铁客流量数据具有时间维度和空间维度属性的特点,提出一种可以同时捕获数据时空特征的预测模型。该模型基于编码器解码器架构设计,其中解码器和编码器均由时空预测模块组成,在该模块中利用图卷积学习地铁站的空间拓朴结构、门控循环单元来捕获数据的时间特征。此外,模型将单位时间间隔内进站和出站客流量分别构成的两个时间序列,即进出站双时间序列作为输入,最终协同预测各站点的进站与出站人数。在上海地铁一卡通数据集上进行对比实验,实验结果表明,所提出的模型在进站与出站客流量预测上均取得了更好的效果,这表明考虑空间依赖能够有效地提高模型预测精度。     

浙江省缙云地区民居出入口形制研究

门是空间转化的媒介,门是出入口重要组成部分,由于历史,地理,社会等原因,造就了不同的建筑出入口类型,浙江地区的出入口就是一个典型的例子。通过对岩下村出入口形式的考察研究,从而了解其形成原因并进行文字记录与整理。     

Indistinguishability and identifiability of kinetic models for the MurC reaction in peptidoglycan biosynthesis

An important question in Systems Biology is the design of experiments that enable discrimination between two (or more) competing chemical pathway models or biological mechanisms. In this paper analysis is performed between two different models describing the kinetic mechanism of a three-substrate three-product reaction, namely the MurC reaction in the cytoplasmic phase of peptidoglycan biosynthesis. One model involves ordered substrate binding and ordered release of the three products; the competing model also assumes ordered substrate binding, but with fast release of the three products. The two versions are shown to be distinguishable; however, if standard quasi-steady-state assumptions are made distinguishability cannot be determined. Once model structure uniqueness is ensured the experimenter must determine if it is possible to successfully recover rate constant values given the experiment observations, a process known as structural identifiability. Structural identifiability analysis is carried out for both models to determine which of the unknown reaction parameters can be determined uniquely, or otherwise, from the ideal system outputs. This structural analysis forms an integrated step towards the modelling of the full pathway of the cytoplasmic phase of peptidoglycan biosynthesis.     

Dynamical analysis of the conformation of the active site of porcine pancreatic elastase in native and Michaelis complex states. Molecular dynamics simulations

Mechanistic details of events preceeding and occurring during enzymatic catalysis are extremely difficult to obtain experimentally. While molecular dynamics methods permit the simulation of discrete steps along the catalytic pathway, they also can overwhelm the investigator with a mountain of structural details, necessitating the development of specific analytic tools. Based on recent crystallographic data, several molecular dynamics simulations have been performed using the AMBER 3.0 program package. Two different simulations for the hydrated native enzyme (100 and 50 ps) and two simulations (30 and 70 ps) of the Michaelis complex of this enzyme with the substrate (Thr-Pro-nVal-Leu-Tyr-Thr) have been performed. Dynamical properties of the active site (especially the catalytic tetrad: Ser-195, His-57, Asp-102 and Ser-214; chymotrypsinogen numbering system) have been examined using the program MD_ANALYSIS_1. It, together with the program MDKINO, facilitates analysis of dynamical changes of conformation (especially the hydrogen bond network) of the active site. Hydrogen bonding among Asp-102, Ser-214 and His-57 was quite stable, but the catalytic Ser-195 sidechain was flexible. Therefore the catalytically crucial H-bond between HO_γ (Ser-195) and N_ϵ (His-57) is relatively labile. In the native case (i.e. without substrate) this H-bond was never formed due to competition for the acceptor atom (N_ϵ) between water molecules and the HO_γ group. In the Michaelis complex the H-bond is more readily formed, although the sidechain (Ser-195) may sometimes change its conformation do to the influence of the carbonyl group of Ser-214. Due to the dynamical motion of the enzyme there were six different short periods in which the distance between both heavy atoms in this crucial H-bond was less than 2.6 Å, which may facilitate proton transfer from Ser-195 to His-57 (the first step during the proper catalysis). These results suggest mechanistic details about the precise clockwork of a functioning enzyme.     

Gate and common pathway detection in crowd scenes and anomaly detection using motion units and LSTM predictive models

In this paper, we propose two approaches to analyze the crowd scenes. The first one is motion units and meta-tracking based approach (MUDAM Approach). In this approach, the scene is divided into a number of dynamic divisions with coherent motion dynamics called the motion units (MUs). By analyzing the relationships between these MUs using a proposed continuation likelihood, the scene entrance and exit gates are retrieved. A meta-tracking procedure is then applied and the scene dominant motion pathways are retrieved. To overcome the limitations of the MUDAM approach, and detect some of the anomalies, that may happen in these scenes, we proposed another new LSTM based approach. In this approach, the scene is divided into a number of static overlapped spatial regions named super regions (SRs), which cover the whole scene. Long Short Term Memory (LSTM) is used in defining a predictive model for each of the scene SRs. Each LSTM predictive model uses its SR tracklets in the training, such that, it can capture the whole motion dynamics of that SR. Using apriori known scene entrance segments, the proposed LSTM predictive models are applied and the scene dominant motion pathways are retrieved. an anomaly metric is formulated to be used with the LSTM predictive models to detect the scene anomalies. Prototypes of our proposed approaches were developed and evaluated on the challenging New York Grand Central station scene, in addition to four other crowded scenes. Four types of anomalies that may happen in the crowded scenes were defined in the context, and our proposed LSTM based approach was used in detecting these anomalies. Experimental results on anomalies detection were applied too on a number of data sets. Ov erall, the proposed approaches managed to outperform the state of the art methods in retrieving the scene gates and common pathways, in addition to detecting motion anomalies.

     

Molecular dynamics simulations of mutated Mycobacterium tuberculosis l-alanine dehydrogenase to illuminate the role of key residues

l-Alanine dehydrogenase from Mycobacterium tuberculosis (l-MtAlaDH) catalyzes the NADH-dependent interconversion of l-alanine and pyruvate, and it is considered to be a potential target for the treatment of tuberculosis. The experiment has verified that amino acid replacement of the conserved active-site residues which have strong stability and no great changes in biological evolutionary process, such as His96 and Asp270, could lead to inactive mutants [Ågren et al., J. Mol. Biol. 377 (2008) 1161–1173]. However, the role of these conserved residues in catalytic reaction still remains unclear. Based on the crystal structures, a series of mutant structures were constructed to investigate the role of the conserved residues in enzymatic reaction by using molecular dynamics simulations. The results show that whatever the conserved residues were mutated, the protein can still convert its conformation from open state to closed state as long as NADH is present in active site. Asp270 maintains the stability of nicotinamide ring and ribose of NADH through hydrogen bond interactions, and His96 is helpful to convert the protein conformation by interactions with Gln271, whereas, they would lead to the structural rearrangement in active site and lose the catalytic activity when they were mutated. Additionally, we deduce that Met301 plays a major role in catalytic reaction due to fixing the nicotinamide ring of NADH to prevent its rotation, and we propose that Met301 would be mutated to the hydrophobic residue with large steric hindrance in side chain to test the activity of the protein in future experiment.