Conformational Compression and Barrier Height Heterogeneity in the N-Acetylglutamate Kinase

Accounting for the instantaneous disorder in the enzyme-substrate Michaelis complex to calculate the Gibbs free energy barrier of an enzyme reaction

Many enzyme reactions present instantaneous disorder. These dynamic fluctuations in the enzyme-substrate Michaelis complexes generate a wide range of energy barriers that cannot be experimentally observed, but that determine the measured kinetics of the reaction. These individual energy barriers can be calculated using QM/MM methods, but then the problem is how to deal with this dispersion of energy barriers to provide kinetic information. So far, the most usual procedure has implied the so-called exponential average of the energy barriers. In this paper, we discuss the foundations of this method, and we use the free energy perturbation theory to derive an alternative equation to get the Gibbs free energy barrier of the enzyme reaction. In addition, we propose a practical way to implement it. We have chosen four enzyme reactions as examples. In particular, we have studied the hydrolysis of a glycosidic bond catalyzed by the enzyme Thermus thermophilus β-glycosidase, and the mutant Y284P Ttb-gly, and the hydrogen abstraction reactions from C13 and C7 of arachidonic acid catalyzed by the enzyme rabbit 15-lipoxygenase-1.

Asymmetric dynamic coupling promotes alternative evolutionary pathways in an enzyme dimer

The importance of dynamic factors in enzyme evolution is gaining recognition. Here we study how the evolution of a new enzymatic activity exploits conformational tinkering and demonstrate that conversion of a dimeric phosphotriesterase to an arylesterase in Pseudomonas diminuta is accompanied by structural divergence between the two subunits. Deviations in loop conformations increase with promiscuity, leading to functionally distinct states, while they decrease during specialisation for the new function. We show that opposite loop movements in the two subunits are due to a dynamic coupling with the dimer interface, the importance of which is also corroborated by the co-evolution of the loop and interface residues. These results illuminate how protein dynamics promotes conformational heterogeneity in a dimeric enzyme, leading to alternative evolutionary pathways for the emergence of a new function.

A protocol to obtain multidimensional quantum tunneling corrections derived from QM(DFT)/MM calculations for an enzyme reaction

The multidimensional small-curvature tunneling (SCT) method with Electrostatic Embedding calculations is a compromise between an accessible computational cost and the attainment of an accurate enough estimation of tunneling for an enzyme reaction.

Phosphorylation Mechanism of N-Acetyl-l-glutamate Kinase, a QM/MM Study

In microorganisms and plants, N-acetyl-l-glutamate kinase (NAGK) catalyzes the second step in l-arginine synthesis, the phosphorylation of N-Acetyl-l-glutamate (NAG) to give N-acetyl-l-glutamate-5-phosphate. NAGK is only present in microorganisms and plants but absent in mammals, which makes it an attractive target for antimicrobial or biocidal development. Understanding the substrate binding mode and reaction mechanism of NAGK is crucial for targeting the kinase to develop potential therapies. Here, the substrate binding mode was studied by comparing the conformational change of NAGK in the presence and in the absence of the NAG substrate based on molecular dynamics simulations. We revealed that with substrate binding, the catalytic site of the kinase involving three loops in NAGK exhibits a closed conformation, which is predominantly controlled by an interaction between Arg98 and the α-COO^- of NAG. Lys41 is found to guide phosphate transfer through the interactions with the β-,γ-, and γ-phosphate oxygen atoms of adenosine 5'-triphosphate surrounded by two highly conserved glycine residues (Gly44 and Gly76), while Arg98 helps to position the NAG substrate in the catalytic site, which facilitates the phosphate transfer. Furthermore, we elucidated phosphate-transfer reaction mechanism using hybrid density functional theory-based quantum mechanics/molecular mechanics calculations (B97D/AMBER99) and found that the catalysis follows a dissociative mechanism.

How Many Conformations of Enzymes Should Be Sampled for DFT/MM Calculations? A Case Study of Fluoroacetate Dehalogenase

The quantum mechanics/molecular mechanics (QM/MM) method (e.g., density functional theory (DFT)/MM) is important in elucidating enzymatic mechanisms. It is indispensable to study "multiple" conformations of enzymes to get unbiased energetic and structural results. One challenging problem, however, is to determine the minimum number of conformations for DFT/MM calculations. Here, we propose two convergence criteria, namely the Boltzmann-weighted average barrier and the disproportionate effect, to tentatively address this issue. The criteria were tested by defluorination reaction catalyzed by fluoroacetate dehalogenase. The results suggest that at least 20 conformations of enzymatic residues are required for convergence using DFT/MM calculations. We also tested the correlation of energy barriers between small QM regions and big QM regions. A roughly positive correlation was found. This kind of correlation has not been reported in the literature. The correlation inspires us to propose a protocol for more efficient sampling. This saves 50% of the computational cost in our current case.

QM/MM study of the mechanism of reduction of 3-hydroxy-3-methylglutaryl coenzyme A catalyzed by human HMG-CoA reductase

Detailing with atomistic resolution the reaction mechanism of human HMG-CoA reductase (HMG-CoA-R) might provide valuable insights for the development of new cholesterol-lowering drugs.

An integrated structural proteomics approach along the druggable genome of Corynebacterium pseudotuberculosis species for putative druggable targets

Background

The bacterium Corynebacterium pseudotuberculosis (Cp) causes caseous lymphadenitis (CLA), mastitis, ulcerative lymphangitis, and oedema in a number of hosts, comprising ruminants, thereby intimidating economic and dairy industries worldwide. So far there is no effective drug or vaccine available against Cp. Previously, a pan-genomic analysis was performed for both biovar equi and biovar ovis and a Pathogenicity Islands (PAIS) analysis within the strains highlighted a large set of proteins that could be relevant therapeutic targets for controlling the onset of CLA. In the present work, a structural druggability analysis pipeline was accomplished along 15 previously sequenced Cp strains from both biovar equi and biovar ovis.

Methods and results

We computed the whole modelome of a reference strain Cp1002 (NCBI Accession: NC_017300.1) and then the homology models of proteins, of 14 different Cp strains, with high identity (≥ 85%) to the reference strain were also done. Druggability score of all proteins pockets was calculated and only those targets that have a highly druggable (HD) pocket in all strains were kept, a set of 58 proteins. Finally, this information was merged with the previous PAIS analysis giving two possible highly relevant targets to conduct drug discovery projects. Also, off-targeting information against host organisms, including Homo sapiens and a further analysis for protein essentiality provided a final set of 31 druggable, essential and non-host homologous targets, tabulated in table S4, additional file 1. Out of 31 globally druggable targets, 9 targets have already been reported in other pathogenic microorganisms, 3 of them (3-isopropylmalate dehydratase small subunit, 50S ribosomal protein L30, Chromosomal replication initiator protein DnaA) in C. pseudotuberculosis.

Conclusion

Overall we provide valuable information of possible targets against C. pseudotuberculosis where some of these targets have already been reported in other microorganisms for drug discovery projects, also discarding targets that might be physiologically relevant but are not amenable for drug binding. We propose that the constructed in silico dataset might serve as a guidance for the scientific community to have a better understanding while selecting putative therapeutic protein candidates as druggable ones as effective measures against C. pseudotuberculosis.