Electronic structure benchmark calculations of CO2 fixing elementary chemical steps in RuBisCO using the projector-based embedding approach.

Ribulose 1,5-bisphosphate carboxylase-oxygenase (RuBisCO) is the main enzyme involved in atmospheric carbon dioxide (CO2 ) fixation in the biosphere. This enzyme catalyzes a set of five chemical steps that take place in the same active-site within magnesium (II) coordination sphere. Here, a set of electronic structure benchmark calculations have been carried out on a reaction path proposed by Gready et al. by means of the projector-based embedding approach. Activation and reaction energies for all main steps catalyzed by RuBisCO have been calculated at the MP2, SCS-MP2, CCSD, and CCSD(T)/aug-cc-pVDZ and cc-pVDZ levels of theory. The treatment of the magnesium cation with post-HF methods is explored to determine the nature of its involvement in the mechanism. With the high-level ab initio values as a reference, we tested the performance of a set of density functional theory (DFT) exchange-correlation (xc) functionals in reproducing the reaction energetics of RuBisCO carboxylase activity on a set of model fragments. Different DFT xc-functionals show large variation in activation and reaction energies. Activation and reaction energies computed at the B3LYP level are close to the reference SCS-MP2 results for carboxylation, hydration and protonation reactions. However, for the carbon-carbon bond dissociation reaction, B3LYP and other functionals give results that differ significantly from the ab initio reference values. The results show the applicability of the projector-based embedding approach to metalloenzymes. This technique removes the uncertainty associated with the selection of different DFT xc-functionals and so can overcome some of inherent limitations of DFT calculations, complementing, and potentially adding to modeling of enzyme reaction mechanisms with DFT methods.

Taming the Reactivity of Monoterpene Synthases To Guide Regioselective Product Hydroxylation.

Monoterpenoids are industrially important natural products with applications in the flavours, fragrances, fuels and pharmaceutical industries. Most monoterpenoids are produced by plants, but recently two bacterial monoterpene synthases have been identified, including a cineole synthase (bCinS). Unlike plant cineole synthases, bCinS is capable of producing nearly pure cineole from geranyl diphosphate in a complex cyclisation cascade that is tightly controlled. Here we have used a multidisciplinary approach to show that Asn305 controls water attack on the α-terpinyl cation and subsequent cyclisation and deprotonation of the α-terpineol intermediate, key steps in the cyclisation cascade which direct product formation towards cineole. Mutation of Asn305 results in variants that no longer produce α-terpineol or cineole. Molecular dynamics simulations revealed that water coordination is disrupted in all variants tested. Quantum mechanics calculations indicate that Asn305 is most likely a (transient) proton acceptor for the final deprotonation step. Our synergistic approach gives unique insight into how a single residue, Asn305, tames the promiscuous chemistry of monoterpene synthase cyclisation cascades. It does this by tightly controlling the final steps in cineole formation catalysed by bCinS to form a single hydroxylated monoterpene product.

Projector-Based Embedding Eliminates Density Functional Dependence for QM/MM Calculations of Reactions in Enzymes and Solution.

Combined quantum mechanics/molecular mechanics (QM/MM) methods are increasingly widely utilized in studies of reactions in enzymes and other large systems. Here, we apply a range of QM/MM methods to investigate the Claisen rearrangement of chorismate to prephenate, in solution, and in the enzyme chorismate mutase. Using projector-based embedding in a QM/MM framework, we apply treatments up to the CCSD(T) level. We test a range of density functional QM/MM methods and QM region sizes. The results show that the calculated reaction energetics are significantly more sensitive to the choice of density functional than they are to the size of the QM region in these systems. Projector-based embedding of a wave function method in DFT reduced the 13 kcal/mol spread in barrier heights calculated at the DFT/MM level to a spread of just 0.3 kcal/mol, essentially eliminating dependence on the functional. Projector-based embedding of correlated ab initio methods provides a practical method for achieving high accuracy for energy profiles derived from DFT and DFT/MM calculations for reactions in condensed phases.

Visualizing protein-ligand binding with chemical energy-wise decomposition (CHEWD): application to ligand binding in the kallikrein-8 S1 Site.

Kallikrein-8, a serine protease, is a target for structure-based drug design due to its therapeutic potential in treating Alzheimer's disease and is also useful as a biomarker in ovarian cancer. We present a binding assessment of ligands to kallikrein-8 using a residue-wise decomposition of the binding energy. Binding of four putative inhibitors of kallikrein-8 is investigated through molecular dynamics simulation and ligand binding energy evaluation with two methods (MM/PBSA and WaterSwap). For visualization of the residue-wise decomposition of binding energies, chemical energy-wise decomposition or CHEWD is introduced as a plugin to UCSF Chimera and Pymol. CHEWD allows easy comparison between ligands using individual residue contributions to the binding energy. Molecular dynamics simulations indicate one ligand binds stably to the kallikrein-8 S1 binding site. Comparison with other members of the kallikrein family shows that residues responsible for binding are specific to kallikrein-8. Thus, ZINC02927490 is a promising lead for development of novel kallikrein-8 inhibitors.

Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors.

The Cytochrome P450 family of heme-containing proteins plays a major role in catalyzing phase I metabolic reactions, and the CYP3A4 subtype is responsible for the metabolism of many currently marketed drugs. Additionally, CYP3A4 has an inherent affinity for a broad spectrum of structurally diverse chemical entities, often leading to drug-drug interactions mediated by the inhibition or induction of the metabolic enzyme. The current study explores the binding of selected highly efficient CYP3A4 inhibitors by docking and molecular dynamics (MD) simulation protocols and their binding free energy calculated using the WaterSwap method. The results indicate the importance of binding pocket residues including Phe57, Arg105, Arg106, Ser119, Arg212, Phe213, Thr309, Ser312, Ala370, Arg372, Glu374, Gly481 and Leu483 for interaction with CYP3A4 inhibitors. The residue-wise decomposition of the binding free energy from the WaterSwap method revealed the importance of binding site residues Arg106 and Arg372 in the stabilization of all the selected CYP3A4-inhibitor complexes. The WaterSwap binding energies were further complemented with the MM(GB/PB)SA results and it was observed that the binding energies calculated by both methods do not differ significantly. Overall, our results could guide towards the use of multiple computational approaches to achieve a better understanding of CYP3A4 inhibition, subsequently leading to the design of highly specific and efficient new chemical entities with suitable ADMETox properties and reduced side effects.

Towards a peptide-based vaccine against Shigella sonnei: A subtractive reverse vaccinology based approach.

Shigella sonnei is one of the major causes of shigellosis in technically advanced countries and reports of its unprecedented increase are published from the Middle East, Latin America, and Asia. The pathogen exhibits resistance against first and second line antibiotics which highlights the need for the development of an effective broad-spectrum vaccine. A computational based approach comprising subtractive reverse vaccinology was used for the identification of potential peptide-based vaccine candidates in the proteome of S. sonnei reference strain (53G). The protocol revealed three essential, host non-homologous, highly virulent, antigenic, conserved and adhesive vaccine proteins: TolC, PhoE, and outer membrane porin protein. The cellular interactome of these proteins supports their direct and indirect involvement in biologically significant pathways, essential for pathogen survival. Epitope mapping of these candidates reveals the presence of surface exposed 9-mer B-cell-derived T-cell epitopes of an antigenic, virulent, non-allergen nature and have broad-spectrum potency. In addition, molecular docking studies demonstrated the deep binding of the epitopes in the binding groove and the stability of the complex with the most common binding allele in the human population, DRB1*0101. Future characterization of the screened epitopes in order to further investigate the immune protection efficacy in animal models is highly desirable.

Protein dynamics and enzyme catalysis: insights from simulations.

The role of protein dynamics in enzyme catalysis is one of the most active and controversial areas in enzymology today. Some researchers claim that protein dynamics are at the heart of enzyme catalytic efficiency, while others state that dynamics make no significant contribution to catalysis. What is the biochemist - or student - to make of the ferocious arguments in this area? Protein dynamics are complex and fascinating, as molecular dynamics simulations and experiments have shown. The essential question is: do these complex motions have functional significance? In particular, how do they affect or relate to chemical reactions within enzymes, and how are chemical and conformational changes coupled together? Biomolecular simulations can analyse enzyme reactions and dynamics in atomic detail, beyond that achievable in experiments: accurate atomistic modelling has an essential part to play in clarifying these issues. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Analysis of chorismate mutase catalysis by QM/MM modelling of enzyme-catalysed and uncatalysed reactions.

Chorismate mutase is at the centre of current controversy about fundamental features of biological catalysts. Some recent studies have proposed that catalysis in this enzyme does not involve transition state (TS) stabilization but instead is due largely to the formation of a reactive conformation of the substrate. To understand the origins of catalysis, it is necessary to compare equivalent reactions in different environments. The pericyclic conversion of chorismate to prephenate catalysed by chorismate mutase also occurs (much more slowly) in aqueous solution. In this study we analyse the origins of catalysis by comparison of multiple quantum mechanics/molecular mechanics (QM/MM) reaction pathways at a reliable, well tested level of theory (B3LYP/6-31G(d)/CHARMM27) for the reaction (i) in Bacillus subtilis chorismate mutase (BsCM) and (ii) in aqueous solvent. The average calculated reaction (potential energy) barriers are 11.3 kcal mol(-1) in the enzyme and 17.4 kcal mol(-1) in water, both of which are in good agreement with experiment. Comparison of the two sets of reaction pathways shows that the reaction follows a slightly different reaction pathway in the enzyme than in it does in solution, because of a destabilization, or strain, of the substrate in the enzyme. The substrate strain energy within the enzyme remains constant throughout the reaction. There is no unique reactive conformation of the substrate common to both environments, and the transition state structures are also different in the enzyme and in water. Analysis of the barrier heights in each environment shows a clear correlation between TS stabilization and the barrier height. The average differential TS stabilization is 7.3 kcal mol(-1) in the enzyme. This is significantly higher than the small amount of TS stabilization in water (on average only 1.0 kcal mol(-1) relative to the substrate). The TS is stabilized mainly by electrostatic interactions with active site residues in the enzyme, with Arg90, Arg7 and Glu78 generally the most important. Conformational effects (e.g. strain of the substrate in the enzyme) do not contribute significantly to the lower barrier observed in the enzyme. The results show that catalysis is mainly due to better TS stabilization by the enzyme.

Exploration of the structural requirements of Aurora Kinase B inhibitors by a combined QSAR, modelling and molecular simulation approach.

Aurora kinase B plays an important role in the cell cycle to orchestrate the mitotic process. The amplification and overexpression of this kinase have been implicated in several human malignancies. Therefore, Aurora kinase B is a potential drug target for anticancer therapies. Here, we combine atom-based 3D-QSAR analysis and pharmacophore model generation to identify the principal structural features of acylureidoindolin derivatives that could potentially be responsible for the inhibition of Aurora kinase B. The selected CoMFA and CoMSIA model showed significant results with cross-validation values (q2) of 0.68, 0.641 and linear regression values (r2) of 0.971, 0.933 respectively. These values support the statistical reliability of our model. A pharmacophore model was also generated, incorporating features of reported crystal complex structures of Aurora kinase B. The pharmacophore model was used to screen commercial databases to retrieve potential lead candidates. The resulting hits were analyzed at each stage for diversity based on the pharmacophore model, followed by molecular docking and filtering based on their interaction with active site residues and 3D-QSAR predictions. Subsequently, MD simulations and binding free energy calculations were performed to test the predictions and to characterize interactions at the molecular level. The results suggested that the identified compounds retained the interactions with binding residues. Binding energy decomposition identified residues Glu155, Trp156 and Ala157 of site B and Leu83 and Leu207 of site C as major contributors to binding affinity, complementary to 3D-QSAR results. To best of our knowledge, this is the first comparison of WaterSwap field and 3D-QSAR maps. Overall, this integrated strategy provides a basis for the development of new and potential AK-B inhibitors and is applicable to other protein targets.

Combating tigecycline resistant Acinetobacter baumannii: A leap forward towards multi-epitope based vaccine discovery.

Global emergence of Tigecycline resistant Acinetobacter baumannii (TRAB) is on the horizon and poses a very serious threat to human health. There is a pressing demand for suitable therapeutics against this pathogen, particularly a vaccine to protect against TRAB infections. We present a comprehensive investigation of the complete proteome of a TRAB AB031 strain to predict promiscuous antigenic, non-allergenic, virulent B-cell derived T-cell epitopes and formulate a multi-epitope vaccine against the pathogen. We identified epitopes from three proteins: outer membrane protein assembly factor (BamA), fimbrial biogenesis outer membrane usher protein (FimD) and type IV secretion protein (Rhs) that are appropriate for vaccine design. These proteins constitute the core proteome of the pathogen, are essential, localized at the pathogen surface, non-homologous to humans, mice and to the beneficial probiotic bacteria that reside the human gut. Moreover, these proteins are ideal candidates for experimental investigation as they have favorable physicochemical properties and have strong cellular interacting networks. The predicted epitopes: FPLNDKPGD (BamA), FVHAEEAAA (FimD) and YVVAGTAAA (Rhs) have exo-membrane topology for efficient recognition of the host immune system and high affinity for the most prevalent allele in human populations, the DRB*0101. These epitopes were linked and attached to an adjuvant to enhance its antigenicity. The multi-epitope vaccine-construct was docked with the TLR4 receptor to assess its affinity for the protein and thus its presentation to the host immune system. Docking results were validated through molecular dynamics simulations and binding free energies were calculated using the molecular mechanics/generalized Born (MM-GBSA) method. In conclusion, we expect the designed chimeric vaccine is highly likely to be effective against infections caused by TRAB.