Exploring the inhibition mechanism of SARS-CoV-2 main protease by ebselen: A molecular docking, molecular dynamics simulation and DFT approach.

The main protease (Mpro) of SARS-CoV-2 plays an essential role in the virus life cycle and is considered a key target for therapeutic development. This study explores the inhibition mechanism of SARS-CoV-2 Mpro by ebselen, an organoselenium drug that shows potent inhibitory activity. By using a combination of multiple computational methods including molecular docking, molecular dynamics simulations, and density functional theory calculations, the complete covalent inhibition process of ebselen is simulated for the first time. Two possible pathways with different bound conformations of ebselen are identified. The hydrolysis of the enzyme-ebselen adduct is found to be the rate-determining step. The simulation results show that the behavior of water molecules at the hydrolysis site is crucial to distinguish the two paths energetically. Our simulations, which are in agreement with existing experimental results, provide a theoretical basis for the rational design and mechanism exploration of ebselen-based inhibitors.

Theoretical understanding of the effect of specifically adsorbed halide anions on Cu-catalyzed CO2 electroreduction activity and product selectivity.

Initial CO2 electroreduction into CO and its subsequent electroreduction pathways were selected to study the effect of specifically adsorbed halide anions X- (X = F, Cl, Br, I) on CO2 electroreduction activity and product selectivity at Cu(111)/H2O interfaces via DFT calculations. The calculated results show that the presence of halide anions can exert a notable effect on the CO2 adsorption characteristics and that chemically adsorbed CO2 molecules can be formed. Furthermore, the halide-anion-modified Cu(111)/H2O interfaces could significantly enhance the initial CO2 electroreduction into CO activity, which is regarded as the rate-determining step during CO2 electroreduction at clean Cu(111)/H2O interfaces. Analysis of the initial CO2 electroreduction and Volmer reaction pathways showed that the halide-anion-modified Cu(111)/H2O interfaces could suppress the HER and thus improve the CO2 electroreduction activity and product selectivity. It is speculated that the enhanced initial CO2 electroreduction activity at the F--, Cl--, Br--, and I--modified Cu(111)/H2O interfaces may originate from the decreased work functions and anion radical ·CO2- formations. Simultaneously, we concluded that dimer OCCO formations in the presence of halide anions were more favorable than CHO during CO electroreduction according to the order of I- > Br- > Cl- > F- and could result in the production of C2 product, suggesting an improved CO2 electroreduction product selectivity. The present analyses of electronic structure may explain the more favorable OCCO formations in the order of I- > Br- > Cl- > F-. The present understanding of this effect will provide an improved scientific guideline for the control of CO2 electroreduction pathways and design of more efficient electrocatalysts.

Study on the Biomolecular Competitive Mechanism of Polybrominated Diphenyl Ethers and Their Derivatives on Thyroid Hormones.

Polybrominated diphenyl ethers (PBDEs) are widely used brominated flame retardants. PBDEs and their derivatives, hydroxylated PBDEs (OH-PBDEs), can bind to hormone receptors and impact hormone secretion, transportation, and metabolism, leading to endocrine disruption and the development of various diseases. They have particularly strong interference effects on thyroid hormones. This study used decabromodiphenyl ether (BDE-209); 2,2',4,4'-tetrabromodiphenyl ether (BDE-47); and 6-OH-BDE-47 as representative compounds of PBDEs and their derivatives, OH-PBDEs. A fluorescence probe, fluorescein-isothiocyanate-L-thyroxine (FITC-T4, F-T4), specific for binding to transthyretin (TTR), a thyroid transport protein, was prepared. The binding capacity of PBDEs and their derivatives, OH-PBDEs, to TTR was quantitatively measured using fluorescence spectroscopy. The principle of quenching the fluorescence intensity of F-T4 after binding to TTR was used to analyze the competitive interaction between the probe and BDE-209, BDE-47, and 6-OH-BDE-47, thereby evaluating the toxic effects of PBDEs and their derivatives on the thyroid system. Additionally, AutoDock molecular docking software (1.5.6) was used to further analyze the interference mechanism of OH-PBDEs on T4. The results of the study are as follows: (1) Different types of PBDEs and OH-PBDEs exhibit varying degrees of interference with T4. Both the degree of bromination and hydroxylation affect their ability to competitively bind to TTR. Higher bromination and hydroxylation degrees result in stronger competitive substitution. (2) The competitive substitution ability of the same disruptor varies at different concentrations. Higher concentrations lead to stronger substitution ability, but there is a threshold beyond which the substitution ability no longer increases. (3) When OH-PBDEs have four or more bromine atoms and exhibit the most structural similarity to T4, their binding affinity to TTR is stronger than that of T4.

Role of Molecular Interactions and Protein Rearrangement in the Dissociation Kinetics of p38α MAP Kinase Type-I/II/III Inhibitors.

Understanding the governing factors of fast or slow inhibitor binding/unbinding assists in developing drugs with preferred kinetic properties. For inhibitors with the same binding affinity targeting different binding sites of the same protein, the kinetic behavior can profoundly differ. In this study, we investigated unbinding kinetics and mechanisms of fast (type-I) and slow (type-II/III) binders of p38α mitogen-activated protein kinase, where the crystal structures showed that type-I and type-II/III inhibitors bind to pockets with different conformations of the Asp-Phe-Gly (DFG) motif. The work used methods that combine conventional molecular dynamics (MD), accelerated molecular dynamics (AMD) simulations, and the newly developed pathway search guided by internal motions (PSIM) method to find dissociation pathways. The study focuses on revealing key interactions and molecular rearrangements that hinder ligand dissociation by using umbrella sampling and post-MD processing to examine changes in free energy during ligand unbinding. As anticipated, the initial dissociation steps all require breaking interactions that appeared in crystal structures of the bound complexes. Interestingly, for type-I inhibitors such as SB2, p38α keeps barrier-free conformational fluctuation in the ligand-bound complex and during ligand dissociation. In contrast, with a type-II/III inhibitor such as BIRB796, with the rearrangements of p38α in its bound state, ligand unbinding features energetically unfavorable protein-ligand concerted movement. Our results also show that the type-II/III inhibitors preferred dissociation pathways through the allosteric channel, which is consistent with an existing publication. The study suggests that the level of required protein rearrangement is one major determining factor of drug binding kinetics in p38α systems, providing useful information for development of inhibitors.

Visualizing the tunnel in tryptophan synthase with crystallography: Insights into a selective filter for accommodating indole and rejecting water.

Four new X-ray structures of tryptophan synthase (TS) crystallized with varying numbers of the amphipathic N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F6) molecule are presented. These structures show one of the F6 ligands threaded into the tunnel from the ß-site and reveal a distinct hydrophobic region. Over this expanse, the interactions between F6 and the tunnel are primarily nonpolar, while the F6 phosphoryl group fits into a polar pocket of the ß-subunit active site. Further examination of TS structures reveals that one portion of the tunnel (T1) binds clusters of water molecules, whereas waters are not observed in the nonpolar F6 binding region of the tunnel (T2). MD simulation of another TS structure with an unobstructed tunnel also indicates the T2 region of the tunnel excludes water, consistent with a dewetted state that presents a significant barrier to the transfer of water into the closed ß-site. We conclude that hydrophobic molecules can freely diffuse between the α- and ß-sites via the tunnel, while water does not. We propose that exclusion of water serves to inhibit reaction of water with the α-aminoacrylate intermediate to form ammonium ion and pyruvate, a deleterious side reaction in the αß-catalytic cycle. Finally, while most TS structures show ßPhe280 partially blocking the tunnel between the α- and ß-sites, new structures show an open tunnel, suggesting the flexibility of the ßPhe280 side chain. Flexible docking studies and MD simulations confirm that the dynamic behavior of ßPhe280 allows unhindered transfer of indole through the tunnel, therefore excluding a gating role for this residue.

Characterization of Promiscuous Binding of Phosphor Ligands to Breast-Cancer-Gene 1 (BRCA1) C-Terminal (BRCT): Molecular Dynamics, Free Energy, Entropy and Inhibitor Design.

Inhibition of the protein-protein interaction (PPI) mediated by breast-cancer-gene 1 C-terminal (BRCT) is an attractive strategy to sensitize breast and ovarian cancers to chemotherapeutic agents that induce DNA damage. Such inhibitors could also be used for studies to understand the role of this PPI in DNA damage response. However, design of BRCT inhibitors is challenging because of the inherent flexibility associated with this domain. Several studies identified short phosphopeptides as tight BRCT binders. Here we investigated the thermodynamic properties of 18 phosphopeptides or peptide with phosphate mimic and three compounds with phosphate groups binding to BRCT to understand promiscuous molecular recognition and guide inhibitor design. We performed molecular dynamics (MD) simulations to investigate the interactions between inhibitors and BRCT and their dynamic behavior in the free and bound states. MD simulations revealed the key role of loops in altering the shape and size of the binding site to fit various ligands. The mining minima (M2) method was used for calculating binding free energy to explore the driving forces and the fine balance between configuration entropy loss and enthalpy gain. We designed a rigidified ligand, which showed unfavorable experimental binding affinity due to weakened enthalpy. This was because it lacked the ability to rearrange itself upon binding. Investigation of another phosphate group containing compound, C1, suggested that the entropy loss can be reduced by preventing significant narrowing of the energy well and introducing multiple new compound conformations in the bound states. From our computations, we designed an analog of C1 that introduced new intermolecular interactions to strengthen attractions while maintaining small entropic penalty. This study shows that flexible compounds do not always encounter larger entropy penalty, compared with other more rigid binders, and highlights a new strategy for inhibitor design.

Amino acid networks in a (ß/α)8 barrel enzyme change during catalytic turnover.

Proteins can be viewed as small-world networks of amino acid residues connected through noncovalent interactions. Nuclear magnetic resonance chemical shift covariance analyses were used to identify long-range amino acid networks in the α subunit of tryptophan synthase both for the resting state (in the absence of substrate and product) and for the working state (during catalytic turnover). The amino acid networks observed stretch from the surface of the protein into the active site and are different between the resting and working states. Modification of surface residues on the network alters the structural dynamics of active-site residues over 25 Å away and leads to changes in catalytic rates. These findings demonstrate that amino acid networks, similar to those studied here, are likely important for coordinating structural changes necessary for enzyme function and regulation.

CmCML11 interacts with CmCAMTA5 to enhance γ-aminobutyric acid (GABA) accumulation by regulating GABA shunt in fresh-cut cantaloupe.

The effect of calcium chloride (CaCl2) treatment on γ-aminobutyric acid (GABA) accumulation in fresh-cut cantaloupe and the involved mechanisms were investigated. The result showed that 1% (w/v) CaCl2 treatment increased GABA content and activities of glutamate decarboxylase (GAD) and succinate semialdehyde dehydrogenase (SSADH), while decreased glutamate (Glu) content and GABA transaminase (GABA-T) activities in fresh-cut cantaloupe. CmCML11 and CmCAMTA5 expressions of CaCl2-treated fruit increased by 187.4% and 165.6% than control fruit in the initial 6 h. Besides, expressions of GABA shunt genes, including CmGAD1, CmGAD2, CmGABA-T and CmSSADH were also up-regulated by CaCl2 treatment during early storage. Moreover, acting as a transcriptional activator, CmCAMTA5 could bind to the CG-box in promoters of CmGAD1, CmGABA-T and CmSSADH and activate their transcription. Furthermore, the interaction between CmCML11 and CmCAMTA5 could enhance the transcriptional activation on GABA shunt genes which were regulated by CmCAMTA5. Collectively, our findings revealed that CaCl2 treatment promoted GABA accumulation in fresh-cut cantaloupe via the combined effect of CmCML11 and CmCAMTA5 in the regulation of expressions of CmGAD1, CmGABA-T, and CmSSADH in GABA shunt.

Phosphatidylinositol-specific phospholipase C-associated phospholipid metabolism mediates DcGLRs channel to promote calcium influx under CaCl2 treatment in shredded carrots during storage.

The rapid activation of phosphatidylinositol-specific phospholipase C (PI-PLC) occurs early after the stimulation of biotic and abiotic stress in plants, which directly associated with the calcium channel-induced calcium ion (Ca2+) influx. Exogenous calcium chloride (CaCl2) mediates the calcium signaling transduction to promote the γ-aminobutyric acid accumulation and nutritional quality in shredded carrots whereas the generation mechanism remains uncertain. Therefore, the involvement of PI-PLC-associated phospholipid metabolism was investigated in present study. Our result revealed that CaCl2 treatment promoted the expression and activity of PI-PLC and increased the inositol 1,4,5-trisphosphate and hexakisphosphate content in shredded carrots. The transcripts of multi-glutamate receptor-like channels (DcGLRs), the glutamate and γ-aminobutyric acid (GABA) content, and Ca2+ influx were induced by CaCl2 treatment in shredded carrots during storage. However, PI-PLC inhibitor (U73122) treatment inhibited the activation of PI-PLC, the increase of many DcGLRs family genes expression levels, and Ca2+ influx. Moreover, the identification of DcPI-PLC4/6 and DcGLRs proteins, along with the analysis of characteristic domains such as PLCXc, PLCYc, C2 domain, transmembranous regions, and ligand binding domain, suggests their involvement in phospholipid catalysis and calcium transport in carrots. Furthermore, DcPI-PLC4/6 overexpression in tobacco leaves induced the Ca2+ influx by activating the expressions of NtGLRs and the accumulation of glutamate and GABA. These findings collectively indicate that CaCl2 treatment-induced PI-PLC activation influences DcGLRs expression levels to mediate cytosolic Ca2+ influx, thus, highlighting the "PI-PLC-GLRs-Ca2+" pathway in calcium signaling generation and GABA biosynthesis in shredded carrots.

Theoretical analysis of naproxen reaction with sulfate and hydroxyl radicals in the aqueous phase: Investigating reactive sites and reaction kinetics.

In this study, we have utilized theoretical calculations to predict the reaction active sites of naproxen when reacting with radicals and to further study the thermodynamics and kinetics of the reactions with ·OH and SO4-·. The evidence, derived from the average local ionization energy and electrostatic potential, points to the naphthalene ring as the preferred site of attack, especially for the C2, C6, C9, and C10 sites. The changes in Gibbs free energy and enthalpy of the reactions initiated by ·OH and SO4-· ranged between -19.6 kcal/mol - 26.3 kcal/mol and -22.3 kcal/mol -18.5 kcal/mol, respectively. More in-depth investigation revealed that RA2 pathway for ·OH exhibited the lowest free energy of activation, suggesting this reaction is more inclined to proceed. The second-order rate constant results indicate the ·OH attacking reaction is faster than reactions initiated by SO4·-, yet controlled by diffusion. The consistency between theoretical findings and experimental data underscores the validity of this computational method for our study.

γ-aminobutyric acid treatment inhibits browning and promotes storage quality by regulating reactive oxygen species and membrane lipid metabolism in fresh-cut stem lettuce.

The effects of γ-aminobutyric (GABA) on enzymatic browning, storage quality, membrane and reactive oxygen species (ROS) metabolism in fresh-cut stem lettuce were investigated. The results illustrated that GABA treatment delayed browning degree, polyphenol oxidase (PPO) activity and the expression of LsPPO. Meanwhile, higher chlorophyll and ascorbic acid contents were exhibited in GABA-treated stem lettuce, as well as the slower microbial propagation. Further investigation revealed that exogenous GABA application declined malondialdehyde content, electrolyte leakage and the enzyme activities of membrane metabolism, and the expression levels of related genes were also downregulated. In addition, GABA treatment scavenged ROS and strengthened the enzyme activities of ROS metabolism, as well as the expression levels of corresponding genes. Taken together, these findings implied that the repressed enzymatic browning and microbial propagation in GABA-treated stem lettuce were due to the inhibition of ROS accumulation, enhancement of membrane stability and increased resistance to oxidation.

Genome-wide identification and role of HSFs in antioxidant response of hot water treated zucchini fruit during cold storage.

Zucchini squashes are cold-sensitive and vulnerable to chilling injury (CI) resulting from reactive oxygen species (ROS) and hot water (HW) immersing effectively reduce CI symptoms during cold storage. However, mechanism involved in reduced ROS due to HW treatment has not been characterized well. In this study, tender green zucchini fruit were treated with HW for 15 min at 45 ± 1 °C and stored for 15 d at 4 ± 1 °C and above 90 % relative humidity. Results showed substantial reduction in CI index, electrolyte leakage, malonaldehyde (MDA) contents and ROS accumulation along with increased activity of ROS-scavenging enzymes due to HW treatment. To gain insight into the molecular mechanism involved in antioxidant defense system, transcriptomic analysis revealed that heat shock factors (HSF) accumulated due to HW treatment regulated the ROS pathway during cold stress. CpHSFA4a was one of the highly expressed transcription factors (TF) due to HW treatment that regulated the transcription of ROS enzymes related genes. CpHSFA4a bind actively with heat shock element (HSE) in promoter regions of CpSOD, CpCAT, CpAPX1, CpAPX2, and CpAPX3, activated and increased the expression of these genes. In conclusion, HW treatment alleviated the CI by maintaining ROS homeostasis through CpHSFA4a mediated ROS pathway in zucchini squashes during cold storage.

Exogenous CaCl2 reduces the oxidative cleavage of carotenoids in shredded carrots by targeting CAMTA4-mediated transcriptional repression of carotenoid degradation pathway.

Carotenoid oxidative cleavage is a significant factor contributing to the color changes of shredded carrots and treatment with calcium chloride (CaCl2, 1% w/v) has been observed to alleviate the whitening symptom and color loss. However, the specific mechanism by which CaCl2 treatment suppresses carotenoid degradation remains unclear. In this study, the effect of CaCl2 and EGTA (calcium ion chelating agent) treatment on carotenoid biosynthesis and degradation in shredded carrots and the mechanism involved was investigated. CaCl2 treatment promoted the expression and activity of carotenoid biosynthetic enzyme (phytoene synthase, PSY), but inhibited the increases of the degradative enzyme activity of carotenoid cleavage dioxygenase (CCD) and down-regulated the corresponding transcripts, thus delayed the degradation of total carotenoid and maintaining higher levels of major carotenoid compounds including ß-carotene, α-carotene, lycopene, and lutein in shredded carrots during storage. However, EGTA treatment promoted the gene expression and enzyme activity of CCD and increased the degradation of carotenoid compounds in shredded carrots during storage. Furthermore, the CaCl2 treatment induced DcCAMTA4, identified as a calcium decoder in shredded carrots, which, in turn, suppressed the expressions of DcCCD1 and DcCCD4 by interacting with their promoters. The transient overexpression of DcCAMTA4 in tobacco leaves led to reduced expression of NtCCD1 and NtCCD4, maintaining a higher content of carotenoids. Thus, CaCl2 alleviated the oxidative cleavage of carotenoids in shredded carrots through the DcCAMTA4-mediated carotenoid degradation pathway.

Glutamate application maintains quality and antioxidant capacity of fresh-cut carrots by modulating GABA shunt, phenylpropanoid and ROS metabolism.

The effects of exogenous glutamate treatment on the quality attributes, γ-aminobutyric acid (GABA) shunt, phenylpropanoid pathway, and antioxidant capacity of fresh-cut carrots were investigated. Results showed that glutamate treatment suppressed the increases in lightness and whiteness values, inhibited the degradation of total carotenoids and maintained better flavor and taste in fresh-cut carrots. Moreover, glutamate treatment rapidly promoted the activities of glutamate decarboxylase and GABA transaminase, thus improving the GABA content. It also significantly enhanced the activities of phenylalanine ammonia-lyase, cinnamate-4-hydroxylase, and 4-coumarate coenzyme A ligase and promoted the accumulation of total phenolics as well as the main individual phenolic compounds, including chlorogenic and caffeic acid. In addition, glutamate application activated the reactive oxygen system-related enzyme including peroxidase, superoxide dismutase, ascorbate peroxidase, and catalase activities to maintain higher antioxidant capacity in fresh-cut carrots. These results demonstrated that exogenous glutamate treatment maintained better nutritional quality and alleviated color deterioration by accelerating the accumulation of GABA and phenolics and enhancing the antioxidant capacity in fresh-cut carrots.

Exogenous chlorogenic acid inhibits quality deterioration in fresh-cut potato slices.

Fresh-cut potatoes are prone to surface browning and physiological degradation. Chlorogenic acid (CGA), a natural phenolic antioxidant, has demonstrated preservative properties in various postharvest products. However, the underlying mechanisms of its application on maintaining quality remain unclear. Therefore, the effect of exogenous CGA treatment on quality deterioration of potato slices and the mechanisms involved were investigated. Results revealed CGA treatment retarded the browning coloration, suppressed microbial growth and inhibited the declines in starch, and ascorbic acid contents in potato slices. Meanwhile, the treatment activated the phenylpropanoid pathway but decreased the activities of phenolic decomposition-related enzymes such as polyphenol oxidase (PPO) and tyrosinase and downregulated StPPO expression. Moreover, the treated slices exhibited reduced accumulation of reactive oxygen species and increased activity of antioxidant enzymes. Additionally, they displayed enhanced 2,2-diphenyl-1-picrylhydrazyl radicals scavenging capacity and higher ATP levels. Therefore, these findings indicated that CGA treatment was effective for quality maintenance and antioxidant capacity enhancement in fresh-cut potatoes, thereby providing potential strategies for the preservation and processing of fresh-cut produce.

Modified stainless steel wires with superwettability for highly efficient in-tube solid-phase microextraction.

Construction of different surface wettability is meaningful for the interaction between the sorbent surface and target components. In the current study, four kinds of stainless-steel wires (SSWs) with different hydrophobic/hydrophilic property were prepared and used as the absorbents to enrich the target compounds with different polarity. Comparative extraction of six non-polar polycyclic aromatic hydrocarbons (PAHs) and six polar estrogens was carried out by in-tube solid phase microextraction (IT-SPME). The results showed that two SSWs with the superhydrophobic surfaces exhibited high extraction capacity to the non-polar PAHs with the superior enrichment factor (EF) in the range of 29-672 and 57-744, respectively. In contrast, the superhydrophilic SSWs demonstrated higher enrichment efficiency for the polar estrogens than other hydrophobic SSWs. On the basis of optimized conditions, a validated analysis method was established using six PAHs as model analytes for IT-SPME-HPLC. Acceptable linear ranges (0.5-10 µg L-1) and low detection limits (0.0056-0.32 µg L-1) were achieved using the superhydrophobic wire modified by perfluorooctyl trichlorosilane (FOTS). The relative recoveries spiked at 2, 5 and 10 µg L-1 in the lake water samples were in the range of 81.5%-113.7%. The relative standard deviation (RSD) of intraday (≤0.8%, n = 3) and interday (≤5.3%, n = 3) tests demonstrated the good extraction repeatability for the same extraction tube. Satisfactory repeatability for the preparation of extraction tubes (n = 3) was also obtained with the RSD values in the range of 3.6%-8.0%.

Potential Mean Force from Umbrella Sampling Simulations: What Can We Learn and What Is Missed?

Changes in free energy provide valuable information for molecular recognition, including both ligand-receptor binding thermodynamics and kinetics. Umbrella sampling (US), a widely used free energy calculation method, has long been used to explore the dissociation process of ligand-receptor systems and compute binding free energy. In existing publications, the binding free energy computed from the potential of mean force (PMF) with US simulation mostly yielded "ball park" values with experimental data. However, the computed PMF values are highly influenced by factors such as initial conformations and/or trajectories provided, the reaction coordinate, and sampling of conformational space in each US window. These critical factors have rarely been carefully studied. Here we used US to study the guest aspirin and 1-butanol dissociation processes of ß-cyclodextrin (ß-CD) and an inhibitor SB2 dissociation from a p38α mitogen-activated protein kinase (MAPK) complex. For ß-CD, we used three different ß-CD conformations to generate the dissociation path with US windows. For p38α, we generated the dissociation pathway by using accelerated molecular dynamics followed by conformational relaxing with short conventional MD, steered MD, and manual pulling. We found that, even for small ß-CD complexes, different ß-CD conformations altered the height of the PMF, but the pattern of PMF was not affected if the MD sampling in each US window was well-converged. Because changing the macrocyclic ring conformation needs to rotate dihedral angles in the ring, a bound ligand largely restrains the motion of cyclodextrin. Therefore, once a guest is in the binding site, cyclodextrin cannot freely change its initial conformation, resulting in different absolute heights of the PMF, which cannot be overcome by running excessively long MD simulations for each US window. Moreover, if the US simulations were not converged, the important barrier and minimum were missed. For ligand-protein systems, our studies also suggest that the dissociation trajectories modeled by an enhanced sampling method must maintain a natural molecular movement to avoid biased PMF plots when using US simulations.

Investigation of Structural Dynamics of Enzymes and Protonation States of Substrates Using Computational Tools.

This review discusses the use of molecular modeling tools, together with existing experimental findings, to provide a complete atomic-level description of enzyme dynamics and function. We focus on functionally relevant conformational dynamics of enzymes and the protonation states of substrates. The conformational fluctuations of enzymes usually play a crucial role in substrate recognition and catalysis. Protein dynamics can be altered by a tiny change in a molecular system such as different protonation states of various intermediates or by a significant perturbation such as a ligand association. Here we review recent advances in applying atomistic molecular dynamics (MD) simulations to investigate allosteric and network regulation of tryptophan synthase (TRPS) and protonation states of its intermediates and catalysis. In addition, we review studies using quantum mechanics/molecular mechanics (QM/MM) methods to investigate the protonation states of catalytic residues of ß-Ketoacyl ACP synthase I (KasA). We also discuss modeling of large-scale protein motions for HIV-1 protease with coarse-grained Brownian dynamics (BD) simulations.

Protonation states and catalysis: Molecular dynamics studies of intermediates in tryptophan synthase.

The importance of protonation states and proton transfer in pyridoxal 5'-phosphate (PLP)-chemistry can hardly be overstated. Although experimental approaches to investigate pKa values can provide general guidance for assigning proton locations, only static pictures of the chemical species are available. To obtain the overall protein dynamics for the interpretation of detailed enzyme catalysis in this study, guided by information from solid-state NMR, we performed molecular dynamics (MD) simulations for the PLP-dependent enzyme tryptophan synthase (TRPS), whose catalytic mechanism features multiple quasi-stable intermediates. The primary objective of this work is to elucidate how the position of a single proton on the reacting substrate affects local and global protein dynamics during the catalytic cycle. In general, proteins create a chemical environment and an ensemble of conformational motions to recognize different substrates with different protonations. The study of these interactions in TRPS shows that functional groups on the reacting substrate, such as the phosphoryl group, pyridine nitrogen, phenolic oxygen and carboxyl group, of each PLP-bound intermediate play a crucial role in constructing an appropriate molecular interface with TRPS. In particular, the protonation states of the ionizable groups on the PLP cofactor may enhance or weaken the attractions between the enzyme and substrate. In addition, remodulation of the charge distribution for the intermediates may help generate a suitable environment for chemical reactions. The results of our study enhance knowledge of protonation states for several PLP intermediates and help to elucidate their effects on protein dynamics in the function of TRPS and other PLP-dependent enzymes.