Molecular Dynamics Study of the Structure, Flexibility and Dynamics of Thermostable L1 Lipase at High Temperatures

Genome-level therapeutic targets identification and chimeric Vaccine designing against the Blastomyces dermatitidis

Blastomyces dermatitidis is a thermally dimorphic fungus that can cause serious and sometimes fatal infections, including blastomycosis. After spore inhalation, a pulmonary infection develops, which can be asymptomatic and have lethal effects, such as acute respiratory distress syndrome. Its most common extra-pulmonary sites are the central nervous system, bones, skin, and genito-urinary systems. Currently, no vaccine has been approved by the FDA to prevent this infection. In the study, a peptide-based vaccine was developed against blastomycosis by using subtractive proteomics and reverse vaccinology approaches. It focuses on mining the whole genome of B. dermatitidis, identifying potential therapeutic targets, and pinpointing potential epitopes for both B- and T-cells that are immunogenic, non-allergenic, non-toxic, and highly antigenic. Multi-epitope constructs were generated by incorporating appropriate linker sequences. A linker (EAAAK) was also added to incorporate an adjuvant sequence to increase immunological potential. The addition of adjuvants and linkers ultimately resulted in the formation of a vaccine construct in which the number of amino acids was 243 and the molecular weight was 26.18 kDa. The designed antigenic and non-allergenic vaccine constructs showed suitable physicochemical properties. The vaccine's structures were predicted, and further analysis verified their interactions with the human TLR-4 receptor through protein-protein docking. Additionally, MD simulation showed a potent interaction between prioritized vaccine-receptor complexes. Immune simulation predicted that the final vaccine injections resulted in significant immune responses for the T- and B-cell immune responses. Moreover, in silico cloning ensured a high expression possibility of the lead vaccine in the E. coli (K12) vector. This study offers an initiative for the development of effective vaccines against B. dermatitidis; however, it is necessary to validate the designed vaccine's immunogenicity experimentally.

Exploring the Antarctic aminopeptidase P from Pseudomonas sp. strain AMS3 through structural analysis and molecular dynamics simulation

Aminopeptidase P (APPro) is a crucial metalloaminopeptidase involved in amino acid cleavage from peptide N-termini, playing essential roles as versatile biocatalysts with applications ranging from pharmaceuticals to industrial processes. Despite acknowledging its potential for catalysis in lower temperatures, detailed molecular basis and biotechnological implications in cold environments are yet to be explored. Therefore, this research aims to investigate the molecular mechanisms underlying the cold-adapted characteristics of APPro from Pseudomonas sp. strain AMS3 (AMS3-APPro) through a detailed analysis of its structure and dynamics. In this study, structure analysis and molecular dynamics (MD) simulation of a predicted model of AMS3-APPro has been performed at different temperatures to assess structural flexibility and thermostability across a temperature range of 0-60 °C over 100 ns. The MD simulation results revealed that the structure were able to remain stable at low temperatures. Increased temperatures present a potential threat to the overall stability of AMS3-APPro by disrupting the intricate hydrogen bond networks crucial for maintaining structural integrity, thereby increasing the likelihood of protein unfolding. While the metal binding site at the catalytic core exhibits resilience at higher temperatures, highlighting its local structural integrity, the overall enzyme structure undergoes fluctuations and potential denaturation. This extensive structural instability surpasses the localized stability observed at the metal binding site. Consequently, these assessments offer in-depth understanding of the cold-adapted characteristics of AMS3-APPro, highlighting its capability to uphold its native conformation and stability in low-temperature environments. In summary, this research provides valuable insights into the cold-adapted features of AMS3-APPro, suggesting its efficient operation in low thermal conditions, particularly relevant for potential biotechnological applications in cold environments.Communicated by Ramaswamy H. Sarma.

Deciphering insights into the binding mechanism and plasticity of Telacebec with M. tuberculosis cytochrome bcc-aa3 supercomplex through an unbiased molecular dynamics simulation, free-energy analysis, and DFT study

The cytochrome bcc-aa3 supercomplex, a key component in the electron transport chain pathway involved in bacterial energy production and homeostasis, is a clinically validated target for tuberculosis (TB), leading to Telacebec (Q203). Telacebec is a potent candidate drug under Phase II clinical development for the treatment of drug-sensitive and drug-resistant TB. Recently, the cryo-electron microscopy structure of this supercomplex from Mycobacterium tuberculosis (Mtb) complexed with Q203 was resolved at 6.9 Å resolution (PDB ID: 7E1W). To understand the binding site (QP site) flexibility and Q203's stability at the QP site of the Mtb cytochrome bcc complex, we conducted molecular dynamics (MD) simulation and free energy analysis on this complex in an explicit hydrated lipid bilayer environment for 500 ns. Through this study, the persistence of a range of direct and indirect interactions was observed over the course of the simulation. The significance of the interactions with His375, Tyr161, Ala178, Ala179, Ile183, His355, Leu356, and Thr313 is underlined. Electrostatic energy was the primary source of the net binding free energy, regardless of the important interacting residues. The overall binding free energy for Q203 was -112.84 ± 7.73 kcal/mol, of which the electrostatic and lipophilic energy contributions were -116.31 ± 1.14 and -21.32 ± 2.35 kcal/mol, respectively. Meanwhile, DFT calculations were utilized to elucidate Q203's molecular properties. Overall, this study deciphers key insights into the cytochrome bcc-aa3 supercomplex with Q203 on the ground of molecular mechanics and quantum mechanics that may facilitate structure-based drug design and optimization for the discovery of the next-generation antitubercular drug(s).Communicated by Ramaswamy H. Sarma.

Multiple approaches of loop region modification for thermostability improvement of 4,6-α-glucanotransferase from Limosilactobacillus fermentum NCC 3057

4,6-α-glucanotransferase (4,6-α-GT), as a member of the glycoside hydrolase 70 (GH70) family, converts starch/maltooligosaccharides into α,1-6 bond-containing α-glucan and possesses potential applications in food, medical and related industries but does not satisfy the high-temperature requirement due to its poor thermostability. In this study, a 4,6-α-GT (ΔGtfB) from Limosilactobacillus fermentum NCC 3057 was used as a model enzyme to improve its thermostability. The loops of ΔGtfB as the target region were optimized using directed evolution, sequence alignment, and computer-aided design. A total of 11 positive mutants were obtained and iteratively combined to obtain a combined mutant CM9, with high resistance to temperature (50 °C). The activity of mutant CM9 was 2.08-fold the activity of the wild type, accompanied by a 5 °C higher optimal temperature, a 5.76 °C higher melting point (T_m, 59.46 °C), and an 11.95-fold longer half-life time (t_1/2). The results showed that most of the polar residues in the loop region of ΔGtfB are mutated into rigid proline residues. Molecular dynamics simulation demonstrated that the root mean square fluctuation of CM9 significantly decreased by "Breathing" movement reduction of the loop region. This study provides a new strategy for improving the thermostability of 4,6-α-GT through rational loop region modification.

Employing an immunoinformatics approach revealed potent multi-epitope based subunit vaccine for lymphocytic choriomeningitis virus

BACKGROUND

Lymphocytic choriomeningitis virus (LCMV) infects many individuals worldwide and causes severe infection in the immunosuppressant recipient, spontaneous abortion, and congenital disabilities in infants.

OBJECTIVES

There is no specific vaccine or therapeutics available to protect against LCMV infection; thus, there is a need to design a potential vaccine to combat the virus by developing immunity in the population. Herein, we attempted to design a potent multi-epitope vaccine for LCMV using immunoinformatics methods.

METHODS

The whole proteome of the virus was screened and mapped to extract immunodominant B-cell and T-cell epitopes which were fused with appropriate linkers (EAAAK, GGGS, AAY, GPGPG, and AAY), PADRE sequence (13aa) and an adjuvant (50 S ribosomal protein L7/L12) to formulate a multi-epitope vaccine ensemble. Codon adaptation and in silico cloning of the constructed vaccine were carried out using bioinformatics tools. The secondary and tertiary structure of the vaccine construct was predicted and refined. The physicochemical profile of the designed vaccine was analyzed, and the multi-epitope vaccine's potential to bind Toll-like receptors (TLR2 and TLR4) was evaluated through molecular docking and molecular dynamics simulations. Computational immune simulation of the designed vaccine antigen was performed using the C-ImmSim server.

RESULTS

The designed multi-epitope-based vaccine (613 aa) comprised 26 immunodominant (six B-cell, nine cytotoxic T lymphocytes, and 11 helper T lymphocytes) epitopes and is predicted antigenic, non-toxic, non-allergen, soluble, and topographically accessible with a suitable physicochemical profile. The designed vaccine is expected to cover a broad worldwide population (96.35 %) and stimulate a robust adaptive immune response against the virus upon administration. In silico cloning of the constructed vaccine in PET28a (+) vector ensured its optimal expression in the Escherichia coli system. Molecular docking, molecular dynamics simulation, and binding free energy estimation collectively support the stability and energetically favourable interaction of the modeled vaccine-TLR2/4 complexes.

CONCLUSION

The designed multi-epitope vaccine in the present study could serve as a potential vaccine candidate to protect against LMCV infection; however, the experimental validation and safety testing of the vaccine is warranted to validate the study's outcomes.

Contriving multi-epitope vaccine ensemble for monkeypox disease using an immunoinformatics approach

The current global outbreak of monkeypox (MPX) disease, caused by Monkeypox virus (MPXV), has resulted in 16 thousand infection cases, five deaths, and has been declared a global health emergency of international concern by the World Health Organization. Given current challenges in the safety of existing vaccines, a vaccine to prevent MPX infection and/or onset of symptoms would significantly advance disease management. In this context, a multi-epitope-based vaccine could be a well-suited approach. Herein, we searched a publicly accessible database (Virus Pathogen Database and Analysis Resource) for MPXV immune epitopes from various antigens. We prioritized a group of epitopes (10 CD8+ T cells and four B-cell epitopes) using a computer-aided technique based on desirable immunological and physicochemical properties, sequence conservation criteria, and non-human homology. Three multi-epitope vaccines were constructed (MPXV-1–3) by fusing finalized epitopes with the aid of appropriate linkers and adjuvant (beta-defensin 3, 50S ribosomal protein L7/L12, and Heparin-binding hemagglutinin). Codon optimization and in silico cloning in the pET28a (+) expression vector ensure the optimal expression of each construct in the Escherichia Coli system. Two and three-dimensional structures of the constructed vaccines were predicted and refined. The optimal binding mode of the construct with immune receptors [Toll-like receptors (TLR2, TLR3, and TLR4)] was explored by molecular docking, which revealed high docking energies of MPXV-1–TLR3 (–99.09 kcal/mol), MPXV-2–TLR3 (–98.68 kcal/mol), and MPXV-3–TLR2 (–85.22 kcal/mol). Conformational stability and energetically favourable binding of the vaccine-TLR2/3 complexes were assessed by performing molecular dynamics simulations and free energy calculations (Molecular Mechanics/Generalized Born Surface Area method). In silico immune simulation suggested that innate, adaptive, and humoral responses will be elicited upon administration of such potent multi-epitope vaccine constructs. The vaccine constructs are antigenic, non-allergen, non-toxic, soluble, topographically exposed, and possess favourable physicochemical characteristics. These results may help experimental vaccinologists design a potent MPX vaccine.

Exploring whole proteome to contrive multi-epitope-based vaccine for NeoCoV: An immunoinformtics and in-silico approach

Neo-Coronavirus (NeoCoV) is a novel Betacoronavirus (β-CoVs or Beta-CoVs) discovered in bat specimens in South Africa during 2011. The viral sequence is highly similar to Middle East Respiratory Syndrome, particularly that of structural proteins. Thus, scientists have emphasized the threat posed by NeoCoV associated with human angiotensin-converting enzyme 2 (ACE2) usage, which could lead to a high death rate and faster transmission rate in humans. The development of a NeoCoV vaccine could provide a promising option for the future control of the virus in case of human infection. In silico predictions can decrease the number of experiments required, making the immunoinformatics approaches cost-effective and convenient. Herein, with the aid of immunoinformatics and reverse vaccinology, we aimed to formulate a multi-epitope vaccine that may be used to prevent and treat NeoCoV infection. Based on the NeoCoV proteins, B-cell, cytotoxic T lymphocyte (CTL), and helper T lymphocyte (HTL) epitopes were shortlisted. Four vaccines (Neo-1-4) were devised by fusing shortlisted epitopes with appropriate adjuvants and linkers. The secondary and three-dimensional structures of final vaccines were then predicted. The binding interactions of these potential vaccines with toll-like immune receptors (TLR-2, TLR-3, and TLR-4) and major histocompatibility complex molecules (MHC-I and II) reveal that they properly fit into the receptors' binding domains. Besides, Neo-1 and Neo-4 vaccines exhibited better docking energies of -101.08 kcal/mol and -114.47 kcal/mol, respectively, with TLR-3 as compared to other vaccine constructs. The constructed vaccines are highly antigenic, non-allergenic, soluble, non-toxic, and topologically assessable with good physiochemical characteristics. Codon optimization and in-silico cloning confirmed efficient expression of the designed vaccines in Escherichia coli strain K12. In-silico immune simulation indicated that Neo-1 and Neo-4 vaccines could induce a strong immune response against NeoCoV. Lastly, the binding stability and strong binding affinity of Neo-1 and Neo-4 with TLR-3 receptor were validated using molecular dynamics simulations and free energy calculations (Molecular Mechanics/Generalized Born Surface Area method). The final vaccines require experimental validation to establish their safety and effectiveness in preventing NeoCoV infections.

Immunoinformatics design of a novel epitope-based vaccine candidate against dengue virus

Dengue poses a global health threat, which will persist without therapeutic intervention. Immunity induced by exposure to one serotype does not confer long-term protection against secondary infection with other serotypes and is potentially capable of enhancing this infection. Although vaccination is believed to induce durable and protective responses against all the dengue virus (DENV) serotypes in order to reduce the burden posed by this virus, the development of a safe and efficacious vaccine remains a challenge. Immunoinformatics and computational vaccinology have been utilized in studies of infectious diseases to provide insight into the host-pathogen interactions thus justifying their use in vaccine development. Since vaccination is the best bet to reduce the burden posed by DENV, this study is aimed at developing a multi-epitope based vaccines for dengue control. Combined approaches of reverse vaccinology and immunoinformatics were utilized to design multi-epitope based vaccine from the sequence of DENV. Specifically, BCPreds and IEDB servers were used to predict the B-cell and T-cell epitopes, respectively. Molecular docking was carried out using Schrödinger, PATCHDOCK and FIREDOCK. Codon optimization and in silico cloning were done using JCAT and SnapGene respectively. Finally, the efficiency and stability of the designed vaccines were assessed by an in silico immune simulation and molecular dynamic simulation, respectively. The predicted epitopes were prioritized using in-house criteria. Four candidate vaccines (DV-1-4) were designed using suitable adjuvant and linkers in addition to the shortlisted epitopes. The binding interactions of these vaccines against the receptors TLR-2, TLR-4, MHC-1 and MHC-2 show that these candidate vaccines perfectly fit into the binding domains of the receptors. In addition, DV-1 has a better binding energies of - 60.07, - 63.40, - 69.89 kcal/mol against MHC-1, TLR-2, and TLR-4, with respect to the other vaccines. All the designed vaccines were highly antigenic, soluble, non-allergenic, non-toxic, flexible, and topologically assessable. The immune simulation analysis showed that DV-1 may elicit specific immune response against dengue virus. Moreover, codon optimization and in silico cloning validated the expressions of all the designed vaccines in E. coli. Finally, the molecular dynamic study shows that DV-1 is stable with minimum RMSF against TLR4. Immunoinformatics tools are now applied to screen genomes of interest for possible vaccine target. The designed vaccine candidates may be further experimentally investigated as potential vaccines capable of providing definitive preventive measure against dengue virus infection.

In situ and real-time insight into Rhizopus chinensis lipase under high pressure and temperature: Conformational traits and biobehavioural analysis

An in situ and real-time investigation was performed using an optical cell system and in-silico analysis to reveal the impacts of pressure and temperature on the conformational state and behaviours of Rhizopus chinensis lipase (RCL). The fluorescence intensity (FI) of RCL increased remarkably under high pressure, and part of this increase was recovered after depressurization. This result displayed the partially reversible conformational change of RCL, which may be associated with the local change of Trp224 near the catalytic centre. High temperature caused a significant loss of secondary structure, whereas the α-helical segments including the lid were preserved by high pressure even at temperatures over 60 °C. The parameters of enzymatic reaction monitored by UV showed that the hydrolysis rate was remarkably enhanced by the pressure of 200 MPa. In the pressure range of 0.1-200 MPa, the active volume measured by the in situ system decreased from -2.85 to -6.73 mL/mol with the temperature increasing from 20 °C to 40 °C. The high catalytic capacity of the lipase under high pressure and high temperature was primarily attributed to pressure protection on RCL.

Exploring the substrate specificity and catalytic mechanism of imidazolonepropionase (HutI) from Bacillus subtilis

QM/MM calculations reveal that imidazolonepropionase (HutI) specifically catalyzes the hydrolytic cleavage of (S)-enantiomer of substrate to yieldl-formiminoglutamic acid.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.

Catalytic mechanism of acetolactate decarboxylase from Brevibacillus brevis towards both enantiomers of α-acetolactate

QM/MM calculations detail the conversion of both enantiomers of α-acetolactate to (R)-acetoin. (S)-α-acetolactate can be directly converted, whereas (R)-α-acetolactate has to firstly rearrange to the natural (S)-enantiomer through a carboxylate migration.

A comparative molecular dynamics study of thermophilic and mesophilic β-fructosidase enzymes

Arabidopsis thaliana cell wall invertase 1 (AtcwINV1) and Thermotoga maritima β-fructosidase (BfrA) are among the best structurally studied members of the glycoside hydrolase family 32. Both enzymes hydrolyze sucrose as the main substrate but differ strongly in their thermal stability. Mesophilic AtcwINV1 and thermophilic BfrA have divergent sequence similarities in the N-terminal five bladed β-propeller catalytic domain (31 %) and the C-terminal β-sandwich domain (15 %) of unknown function. The two enzymes were subjected to 200 ns molecular dynamics simulations at 300 K (27 °C) and 353 K (80 °C). Regular secondary structure regions, but not loops, in AtcwINV1 and BfrA showed no significant fluctuation differences at both temperatures. BfrA was more rigid than AtcwINV1 at 300 K. The simulation at 353 K did not alter the structural stability of BfrA, but did increase the overall flexibility of AtcwINV1 exhibiting the most fluctuating regions in the β-propeller domain. The simulated heat treatment also increased the gyration radius and hydrophobic solvent accessible surface area of the plant enzyme, consistent with the initial steps of an unfolding process. The preservation of the conformational rigidity of BfrA at 353 K is linked to the shorter size of the protein loops. Shortening of BfrA loops appears to be a key mechanism for thermostability.

Lipase adsorption on different nanomaterials: a multi-scale simulation study

Adsorption orientations of lipase on different nanomaterials with different surface chemistry.

A QM/MM study of the catalytic mechanism of nicotinamidase

Nicotinamidase (Pnc1) is a member of Zn-dependent amidohydrolases that hydrolyzes nicotinamide (NAM) to nicotinic acid (NA), which is a key step in the salvage pathway of NAD(+) biosynthesis. In this paper, the catalytic mechanism of Pnc1 has been investigated by using a combined quantum-mechanical/molecular-mechanical (QM/MM) approach based on the recently obtained crystal structure of Pnc1. The reaction pathway, the detail of each elementary step, the energetics of the whole catalytic cycle, and the roles of key residues and Zn-binding site are illuminated. Our calculation results indicate that the catalytic water molecule comes from the bulk solvent, which is then deprotonated by residue D8. D8 functions as a proton transfer station between C167 and NAM, while the activated C167 serves as the nucleophile. The residue K122 only plays a role in stabilizing intermediates and transition states. The oxyanion hole formed by the amide backbone nitrogen atoms of A163 and C167 has the function to stabilize the hydroxyl anion of nicotinamide. The Zn-binding site rather than a single Zn(2+) ion acts as a Lewis acid to influence the reaction. Two elementary steps, the activation of C167 in the deamination process and the decomposition of catalytic water in the hydrolysis process, correspond to the large energy barriers of 25.7 and 28.1 kcal mol(-1), respectively, meaning that both of them contribute a lot to the overall reaction barrier. Our results may provide useful information for the design of novel and efficient Pnc1 inhibitors and related biocatalytic applications.

Biochemical characterization of a carboxylesterase from the archaeon Pyrobaculum sp. 1860 and a rational explanation of its substrate specificity and thermostability

In this work, genome mining was used to identify esterase/lipase genes in the archaeon Pyrobaculum sp. 1860. A gene was cloned and functionally expressed in Escherichia coli as His-tagged protein. The recombinant enzyme (rP186_1588) was verified by western blotting and peptide mass fingerprinting. Biochemical characterization revealed that rP186_1588 exhibited optimum activity at pH 9.0 and 80 °C towards p-nitrophenyl acetate (K(m): 0.35 mM, k(cat): 11.65 s⁻¹). Interestingly, the purified rP186_1588 exhibited high thermostability retaining 70% relative activity after incubation at 90 °C for 6 h. Circular dichroism results indicated that rP186_1588 showed slight structure alteration from 60 to 90 °C. Structural modeling showed P186_1588 possessed a typical α/β hydrolase's fold with the catalytic triad consisting of Ser97, Asp147 and His172, and was further confirmed by site-directed mutagenesis. Comparative molecular simulations at different temperatures (300, 353, 373 and 473 K) revealed that its thermostability was associated with its conformational rigidity. The binding free energy analysis by MM-PBSA method revealed that the van der Waals interaction played a major role in p-NP ester binding for P186_1588. Our data provide insights into the molecular structures of this archaeal esterase, and may help to its further protein engineering for industrial applications.

A proton-shuttle reaction mechanism for histone deacetylase 8 and the catalytic role of metal ions

Zinc-dependent histone deacetylase 8 (HDAC8) catalyzes the removal of acetyl moieties from histone tails, and is critically involved in regulating chromatin structure and gene expression. The detailed knowledge of its catalytic process is of high importance since it has been established as a most promising target for the development of new antitumor drugs. By employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations and umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have provided further evidence against the originally proposed general acid-base catalytic pair mechanism for Zinc-dependent histone deacetylases. Instead, our results indicated that HDAC8 employs a proton-shuttle catalytic mechanism, in which a neutral His143 first serves as the general base to accept a proton from the zinc-bound water molecule in the initial rate-determining nucleophilic attack step, and then shuttles it to the amide nitrogen atom to facilitate the cleavage of the amide bond. During the deacetylation process, the Zn(2+) ion changes its coordination mode and plays multiple catalytic roles. For the K(+) ion, which is located about 7 A from the catalytic Zn(2+) ion and conserved in class I and II HDACs, our simulations indicated that its removal would lead to the different transition state structure and a higher free energy reaction barrier for the rate-determining step. It is found that the existence of this conserved K(+) ion would enhance the substrate binding, increase the basicity of His143, strengthen the catalytic role of zinc ion, and improve the transition state stabilization by the enzyme environment.

Structural study of carboxylesterase from hyperthermophilic bacteria Geobacillus stearothermophilus by molecular dynamics simulation

Carboxylesterases are ubiquitous enzymes with important physiological, industrial and medical applications such as synthesis and hydrolysis of stereo specific compounds, including the metabolic processing of drugs, and antimicrobial agents. Here, we have performed molecular dynamics simulations of carboxylesterase from hyperthermophilic bacterium Geobacillus stearothermophilus (GsEst) for 10ns each at five different temperatures namely at 300K, 343K, 373K, 473K and 500K. Profiles of root mean square fluctuation (RMSF) identify thermostable and thermosensitive regions of GsEst. Unfolding of GsEst initiates at the thermosensitive alpha-helices and proceeds to the thermostable beta-sheets. Five ion-pairs have been identified as critical ion-pairs for thermostability and are maintained stably throughout the higher temperature simulations. A detailed investigation of the active site residues of this enzyme suggests that the geometry of this site is well preserved up to 373K. Furthermore, the hydrogen bonds between Asp188 and His218 of the active site are stably maintained at higher temperatures imparting stability of this site. Radial distribution functions (RDFs) show similar pattern of solvent ordering and water penetration around active site residues up to 373K. Principal component analysis suggests that the motion of the entire protein as well as the active site is similar at 300K, 343K and 373K. Our study may help to identify the factors responsible for thermostability of GsEst that may endeavor to design enzymes with enhanced thermostability.

Flexibility of Catalytic Zinc Coordination in Thermolysin and HDAC8: A Born-Oppenheimer ab initio QM/MM Molecular Dynamics Study

The different coordination modes and fast ligand exchange of zinc coordination has been suggested to be one key catalytic feature of the zinc ion which makes it an invaluable metal in biological catalysis. However, partly due to the well known difficulties for zinc to be characterized by spectroscopy methods, evidence for dynamic nature of the catalytic zinc coordination has so far mainly been indirect. In this work, Born-Oppenheimer ab initio QM/MM molecular dynamics simulation has been employed, which allows for a first-principle description of the dynamics of the metal active site while properly including effects of the heterogeneous and fluctuating protein environment. Our simulations have provided direct evidence regarding inherent flexibility of the catalytic zinc coordination shell in Thermolysin (TLN) and Histone Deacetylase 8 (HDAC8). We have observed different coordination modes and fast ligand exchange during the picosecond's time-scale. For TLN, the coordination of the carboxylate group of Glu166 to Zinc is found to continuously change between monodentate and bidentate manner dynamically; while for HDAC8, the flexibility mainly comes from the coordination to a non-amino-acid ligand. Such distinct dynamics in the zinc coordination shell between two enzymes suggests that the catalytic role of Zinc in TLN and HDAC8 is likely to be different in spite of the fact that both catalyze the hydrolysis of amide bond. Meanwhile, considering that such Born-Oppenheimer ab initio QM/MM MD simulations are very much desired but are widely considered to be too computationally expensive to be feasible, our current study demonstrates the viability and powerfulness of this state-of-the-art approach in simulating metalloenzymes.