pH‐dependent conformational dynamics of beta‐secretase 1: A molecular dynamics study

Beta‐secretase 1 (BACE‐1) is an aspartyl protease implicated in the overproduction of β‐amyloid fibrils responsible for Alzheimer disease. The process of β‐amyloid genesis is known to be pH dependent, with an activity peak between solution pH of 3.5 and 5.5. We have studied the pH‐dependent dynamics of BACE‐1 to better understand the pH dependent mechanism. We have implemented support for graphics processor unit (GPU) accelerated constant pH molecular dynamics within the AMBER molecular dynamics software package and employed this to determine the relative population of different aspartyl dyad protonation states in the pH range of greatest β‐amyloid production, followed by conventional molecular dynamics to explore the differences among the various aspartyl dyad protonation states. We observed a difference in dynamics between double‐protonated, mono‐protonated, and double‐deprotonated states over the known pH range of higher activity. These differences include Tyr 71‐aspartyl dyad proximity and active water lifetime. This work indicates that Tyr 71 stabilizes catalytic water in the aspartyl dyad active site, enabling BACE‐1 activity.     

UDP‐Api/UDP‐Xyl synthases affect plant development by controlling the content of UDP‐Api to regulate the RG‐II‐borate complex

Rhamnogalacturonan‐II (RG‐II) is structurally the most complex glycan in higher plants, containing 13 different sugars and 21 distinct glycosidic linkages. Two monomeric RG‐II molecules can form an RG‐II‐borate diester dimer through the two apiosyl (Api) residues of side chain A to regulate cross‐linking of pectin in the cell wall. But the relationship of Api biosynthesis and RG‐II dimer is still unclear. In this study we investigated the two homologous UDP‐D‐apiose/UDP‐D‐xylose synthases (AXSs) in Arabidopsis thaliana that synthesize UDP‐D‐apiose (UDP‐Api). Both AXSs are ubiquitously expressed, while AXS2 has higher overall expression than AXS1 in the tissues analyzed. The homozygous axs double mutant is lethal, while heterozygous axs1/+ axs2 and axs1 axs2/+ mutants display intermediate phenotypes. The axs1/+ axs2 mutant plants are unable to set seed and die. By contrast, the axs1 axs2/+ mutant plants exhibit loss of shoot and root apical dominance. UDP‐Api content in axs1 axs2/+ mutants is decreased by 83%. The cell wall of axs1 axs2/+ mutant plants is thicker and contains less RG‐II‐borate complex than wild‐type Col‐0 plants. Taken together, these results provide direct evidence of the importance of AXSs for UDP‐Api and RG‐II‐borate complex formation in plant growth and development.     

Genetic alteration of UDP‐rhamnose metabolism in Botrytis cinerea leads to the accumulation of UDP‐KDG that adversely affects development and pathogenicity

Botrytis cinerea is a model plant‐pathogenic fungus that causes grey mould and rot diseases in a wide range of agriculturally important crops. A previous study has identified two enzymes and corresponding genes (bcdh, bcer) that are involved in the biochemical transformation of uridine diphosphate (UDP)‐glucose, the major fungal wall nucleotide sugar precursor, to UDP‐rhamnose. We report here that deletion of bcdh, the first biosynthetic gene in the metabolic pathway, or of bcer, the second gene in the pathway, abolishes the production of rhamnose‐containing glycans in these mutant strains. Deletion of bcdh or double deletion of both bcdh and bcer has no apparent effect on fungal development or pathogenicity. Interestingly, deletion of the bcer gene alone adversely affects fungal development, giving rise to altered hyphal growth and morphology, as well as reduced sporulation, sclerotia production and virulence. Treatments with wall stressors suggest the alteration of cell wall integrity. Analysis of nucleotide sugars reveals the accumulation of the UDP‐rhamnose pathway intermediate UDP‐4‐keto‐6‐deoxy‐glucose (UDP‐KDG) in hyphae of the Δbcer strain. UDP‐KDG could not be detected in hyphae of the wild‐type strain, indicating fast conversion to UDP‐rhamnose by the BcEr enzyme. The correlation between high UDP‐KDG and modified cell wall and developmental defects raises the possibility that high levels of UDP‐KDG result in deleterious effects on cell wall composition, and hence on virulence. This is the first report demonstrating that the accumulation of a minor nucleotide sugar intermediate has such a profound and adverse effect on a fungus. The ability to identify molecules that inhibit Er (also known as NRS/ER) enzymes or mimic UDP‐KDG may lead to the development of new antifungal drugs.     

Phylogenomic analysis of UDP‐dependent glycosyltransferases provides insights into the evolutionary landscape of glycosylation in plant metabolism

Glycosylated metabolites generated by UDP‐dependent glycosyltransferases (UGTs) play critical roles in plant interactions with the environment as well as human and animal nutrition. The evolution of plant UGTs has previously been explored, but with a limited taxon sampling. In this study, 65 fully sequenced plant genomes were analyzed, and stringent criteria for selection of candidate UGTs were applied to ensure a more comprehensive taxon sampling and reliable sequence inclusion. In addition to revealing the overall evolutionary landscape of plant UGTs, the phylogenomic analysis also resolved the phylogenetic association of UGTs from free‐sporing plants and gymnosperms, and identified an additional UGT group (group R) in seed plants. Furthermore, lineage‐specific expansions and contractions of UGT groups were detected in angiosperms, with the total number of UGTs per genome remaining constant generally. The loss of group Q UGTs in Poales and Brassicales, rather than functional convergence in the group Q containing species, was supported by a gene tree of group Q UGTs sampled from many species, and further corroborated by the absence of group Q homologs on the syntenic chromosomal regions in Arabidopsis thaliana (Brassicales). Branch‐site analyses of the group Q UGT gene tree allowed for identification of branches and amino acid sites that experienced episodic positive selection. The positively selected sites are located on the surface of a representative group Q UGT (PgUGT95B2), away from the active site, suggesting their role in protein folding/stability or protein–protein interactions.     

Crystal structures of the archaeal UDP‐GlcNAc 2‐epimerase from Methanocaldococcus jannaschii reveal a conformational change induced by UDP‐GlcNAc

Uridine diphosphate N ‐ acetylglucosamine (UDP‐GlcNAc) 2‐epimerase catalyzes the interconversion of UDP‐GlcNAc to UDP‐N‐acetylmannosamine (UDP‐ManNAc), which is used in the biosynthesis of cell surface polysaccharides in bacteria. Biochemical experiments have demonstrated that mutation of this enzyme causes changes in cell morphology and the thermoresistance of the cell wall. Here, we present the crystal structures of Methanocaldococcus jannaschii UDP‐GlcNAc 2‐epimerase in open and closed conformations. A comparison of these crystal structures shows that upon UDP and UDP‐GlcNAc binding, the enzyme undergoes conformational changes involving a rigid‐body movement of the C‐terminal domain. We also present the crystal structure of Bacillus subtilis UDP‐GlcNAc 2‐epimerase in the closed conformation in the presence of UDP and UDP‐GlcNAc. Although a structural overlay of these two closed‐form structures reveals that the substrate‐binding site is evolutionarily conserved, some areas of the allosteric site are distinct between the archaeal and bacterial UDP‐GlcNAc 2‐epimerases. This is the first report on the crystal structure of archaeal UDP‐GlcNAc 2‐epimerase, and our results clearly demonstrate the changes between the open and closed conformations of this enzyme. Proteins 2014; 82:1519–1526. © 2014 Wiley Periodicals, Inc.     

Dynamic behavior of the post‐SET loop region of NSD1: Implications for histone binding and drug development

NSD1 is a SET‐domain histone methyltransferase that methylates lysine 36 of histone 3. In the crystal structure of NSD1, the post‐SET loop is in an autoinhibitory position that blocks binding of the histone peptide as well as the entrance to the lysine‐binding channel. The conformational dynamics preceding histone binding and the mechanism by which the post‐SET loop moves to accommodate the target lysine is currently unknown, although potential models have been proposed. Using molecular dynamics simulations, we have identified potential conformations of the post‐SET loop differing from those of previous studies, as well as proposed a model of peptide‐bound NSD1. Our simulations illustrate the dynamic behavior of the post‐SET loop and the presence of a few distinct conformations. In every case, the post‐SET loop remains in an autoinhibitory position blocking the peptide‐binding cleft, suggesting that another interaction is required to optimally position NSD1 in an active conformation. This finding provides initial evidence for a mechanism by which NSD1 preferentially binds nucleosomal substrates.     

Molecular dynamics studies of alanine racemase: a structural model for drug design

Alanine racemase (AlaR) is a bacterial enzyme that catalyzes the interconversion of L- and D-alanine, which is an essential constituent of the peptidoglycan layer of the bacterial cell wall and requires pyridoxal 5'-phosphate (PLP) as a cofactor. The enzyme is universal to bacteria, including mycobacteria, making it an attractive target for drug design. To investigate the effects of flexibility on the binding modes of the substrate and an inhibitor and to analyze how the active site is affected by the presence of the substrate versus inhibitor, a molecular dynamics simulation on the full AlaR dimer from Bacillus stearothermophilus (pdb code: 1SFT) with a D-alanine molecule in one active site and the noncovalent inhibitor, propionate, in the second site has been carried out. Within the time scale of the simulation, we show that the active site becomes more stabilized in the presence of substrate versus inhibitor. The results of this simulation are in agreement with the proposed mechanism of alanine racemase reaction in which the substrate carboxyl group directly participates in the catalysis by acting cooperatively with Tyr 265' and Lys 39. A structural water molecule in contact with both substrate and inhibitor (i.e., in both active sites) and bridging residues in both active sites was identified. It shows a remarkably low mobility and does not exchange with bulk water. This water molecule can be taken into account for the design of specific AlaR inhibitors by either utilizing it as a bridging group or displacing it with an inhibitor atom. The results presented here provide insights into the dynamics of the alanine racemase in the presence of substrate/inhibitor, which will be used for the rational design of novel inhibitors.     

The allosteric mechanism of yeast chorismate mutase: a dynamic analysis

The effector-regulated allosteric mechanism of yeast chorismate mutase (YCM) was studied by normal mode analysis and targeted molecular dynamics. The normal mode analysis shows that the conformational change between YCM in the R state and in the T state can be represented by a relatively small number of low-frequency modes. This suggests that the transition is coded in the structure and is likely to have a low energetic barrier. Quantitative comparisons (i.e. frequencies) between the low-frequency modes of YCM with and without effectors (modeled structures) reveal that the binding of Trp increases the global flexibility, whereas Tyr decreases global flexibility. The targeted molecular dynamics simulation of substrate analog release from the YCM active site suggests that a series of residues are critical for orienting and "recruiting" the substrate. The simulation led to the switching of a series of substrate-release-coupled salt-bridge partners in the ligand-binding domain; similar changes occur in the transition between YCM R-state and T-state crystal structures. Thus, the normal mode analysis and targeted molecular dynamics results provide evidence that the effectors regulate YCM activity by influencing the global flexibility. The change in flexibility is coupled to the binding of substrate to the T state and release of the product from the R state, respectively.     

Chiral Inhibition of Rivaroxaban Derivatives Towards UDP‐Glucuronosyltransferase (UGT) Isoforms

Rivaroxaban is an oral direct factor Xa (FXa) inhibitor clinically used to prevent and treat thromboembolic disorders. Drug–drug interaction (DDI) exist for rivaroxaban and the inhibitors of CYP3A4/5. This study aims to investigate the inhibition of rivaroxaban and its derivatives with a chiral center towards UDP‐glucuronosyltransferases (UGTs). Chemical synthesis was performed to obtain rivaroxaban derivatives with different chiral centers. UGTs supersomes‐catalyzed 4‐methylumbelliferone (4‐MU) glucuronidation was employed to evaluate the inhibition potential towards various UGT isoforms. A significant influence of rivaroxaban derivatives towards UGT1A3 was observed. Chiral centers produce different effects towards the effect of four pairs of rivaroxaban derivatives towards UGT1A3 activity, with stronger inhibition potential of S1 than R1, but stronger inhibition capability of R2, R3, R4 than S2, S3, and S4. Competitive inhibition of R3 and R4 towards UGT1A3 was demonstrated by Dixon and Lineweaver‐Burk plots. In conclusion, the significant influence of rivaroxaban derivatives towards UGT1A3's activity was demonstrated in the present study. The chirality centers highly affected the inhibition behavior of rivaroxaban derivatives towards UGT1A3. Chirality 27:936–943, 2015. © 2015 Wiley Periodicals, Inc.     

Molecular modeling,dynamics, and an insight into the structural inhibition of cofactor independent phosphoglycerate mutase isoform 1 from Wuchereria bancrofti using cheminformatics and mutational studies

Phosphoglycerate mutase catalyzes the interconversion between 2-phosphoglycerate and 3-phosphoglycerate in the glycolytic and gluconeogenic pathways. They exist in two unrelated forms, that is either cofactor (2,3-diphosphoglycerate) dependent or cofactor-independent. These two enzymes have no similarity in amino acid sequence, tertiary structure, and in catalytic mechanism. Wuchereria bancrofti (WB) contains the cofactor-independent form, whereas other organisms can possess the dependent form or both. Since, independent phosphoglycerate mutase (iPGM) is an essential gene for the survival of nematodes, and it has no sequence or structural similarity to the cofactor-dependent phosphoglycerate mutase found in mammals, it represents an attractive drug target for the filarial nematodes. In this current study, a putative cofactor-iPGM gene was identified in the protein sequence of the WB. In the absence of crystal structure, a three-dimensional structure was determined using the homology modeling approximation, and the most stable protein conformation was identified through the molecular dynamics simulation studies, using GROMACS 4.5. Further, the functional or characteristic residues were identified through the sequence analysis, potential inhibitors were short-listed and validated, and potential inhibitors were ranked using the cheminformatics and molecular dynamics simulations studies, Prime MM-GBSA approach, respectively.     

G‐LoSA: An efficient computational tool for local structure‐centric biological studies and drug design

Molecular recognition by protein mostly occurs in a local region on the protein surface. Thus, an efficient computational method for accurate characterization of protein local structural conservation is necessary to better understand biology and drug design. We present a novel local structure alignment tool, G‐LoSA. G‐LoSA aligns protein local structures in a sequence order independent way and provides a GA‐score, a chemical feature‐based and size‐independent structure similarity score. Our benchmark validation shows the robust performance of G‐LoSA to the local structures of diverse sizes and characteristics, demonstrating its universal applicability to local structure‐centric comparative biology studies. In particular, G‐LoSA is highly effective in detecting conserved local regions on the entire surface of a given protein. In addition, the applications of G‐LoSA to identifying template ligands and predicting ligand and protein binding sites illustrate its strong potential for computer‐aided drug design. We hope that G‐LoSA can be a useful computational method for exploring interesting biological problems through large‐scale comparison of protein local structures and facilitating drug discovery research and development. G‐LoSA is freely available to academic users at http://im.compbio.ku.edu/GLoSA/ .     

Restrained molecular dynamics simulations of HIV-1 protease: the first step in validating a new target for drug design

To test the anticorrelated relationship that was recently displayed in conventional molecular dynamics (MD) simulations, several different restrained MD simulations on a wild type and on the V82F/I84V drug-resistant mutant of HIV-1 protease were performed. This anticorrelated relationship refers to the observation that compression of the peripheral ear-to-cheek region of HIV protease (i.e., the elbow of the flap to the fulcrum and the cantilever) occurred as the active site flaps were opening, and, conversely, expansion of that ear-to-cheek region occurred as both flaps were closing. An additional examination of this anticorrelated relationship was necessary to determine whether it can be harnessed in a useful manner. Consequently, six different MD experiments were performed that incorporated pairwise distance restraints in that ear-to-cheek region (i.e., the distance between the alpha-carbons of Gly40 and Gln61 was restrained to either 7.7 or 10.5 A, in both monomers). Pushing the backbones of the ear and the cheek regions away from each other slightly did force the flaps that guard the active site to remain closed in both the wild type and the mutant systems-even though there were no ligands in the active sites. Thus, these restrained MD simulations provided evidence that the anticorrelated relationship can be exploited to affect the dynamic behavior of the flaps that guard the active site of HIV-1 protease. These simulations supported our hypothesis of the mechanism governing flap motion, and they are the first step towards validating that peripheral surface as a new target for drug design.     

Binding hotspots on K‐ras: Consensus ligand binding sites and other reactive regions from probe‐based molecular dynamics analysis

We have used probe‐based molecular dynamics (pMD) simulations to search for interaction hotspots on the surface of the therapeutically highly relevant oncogenic K‐Ras G12D. Combining the probe‐based query with an ensemble‐based pocket identification scheme and an analysis of existing Ras‐ligand complexes, we show that (i) pMD is a robust and cost‐effective strategy for binding site identification, (ii) all four of the previously reported ligand binding sites are suitable for structure‐based ligand design, and (iii) in some cases probe binding and expanded sampling of configurational space enable pocket expansion and increase the likelihood of site identification. Furthermore, by comparing the distribution of hotspots in nonpocket‐like regions with known protein‐ and membrane‐interacting interfaces, we propose that pMD has the potential to predict surface patches responsible for protein‐biomolecule interactions. These observations have important implications for future drug design efforts and will facilitate the search for potential interfaces responsible for the proposed transient oligomerization or interaction of Ras with other biomolecules in the cellular milieu. Proteins 2015; 83:898–909. © 2015 Wiley Periodicals, Inc.     

Functional characterization of UDP‐rhamnose‐dependent rhamnosyltransferase involved in anthocyanin modification,a key enzyme determining blue coloration in Lobelia erinus

Because structural modifications of flavonoids are closely related to their properties, such as stability, solubility, flavor and coloration, characterizing the enzymes that catalyze the modification reactions can be useful for engineering agriculturally beneficial traits of flavonoids. In this work, we examined the enzymes involved in the modification pathway of highly glycosylated and acylated anthocyanins that accumulate in Lobelia erinus. Cultivar Aqua Blue (AB) of L. erinus is blue‐flowered and accumulates delphinidin 3‐O‐p‐coumaroylrutinoside‐5‐O‐malonylglucoside‐3′5′‐O‐dihydroxycinnamoylglucoside (lobelinins) in its petals. Cultivar Aqua Lavender (AL) is mauve‐flowered, and LC‐MS analyses showed that AL accumulated delphinidin 3‐O‐glucoside (Dp3G), which was not further modified toward lobelinins. A crude protein assay showed that modification processes of lobelinin were carried out in a specific order, and there was no difference between AB and AL in modification reactions after rhamnosylation of Dp3G, indicating that the lack of highly modified anthocyanins in AL resulted from a single mutation of rhamnosyltransferase catalyzing the rhamnosylation of Dp3G. We cloned rhamnosyltransferase genes (RTs) from AB and confirmed their UDP‐rhamnose‐dependent rhamnosyltransferase activities on Dp3G using recombinant proteins. In contrast, the RT gene in AL had a 5‐bp nucleotide deletion, resulting in a truncated polypeptide without the plant secondary product glycosyltransferase box. In a complementation test, AL that was transformed with the RT gene from AB produced blue flowers. These results suggest that rhamnosylation is an essential process for lobelinin synthesis, and thus the expression of RT has a great impact on the flower color and is necessary for the blue color of Lobelia flowers.     

Enantioselective Inhibition of Carprofen Towards UDP‐glucuronosyltransferase (UGT) 2B7

UDP‐glucuronosyltransferases (UGTs)‐catalyzed glucuronidation conjugation reaction plays an important role in the elimination of many important clinical drugs and endogenous substances. The present study aims to investigate the enantioselective inhibition of carprofen towards UGT isoforms. In vitro a recombinant UGT isoforms‐catalyzed 4‐methylumbelliferone (4‐MU) glucuronidation incubation mixture was used to screen the inhibition potential of (R)‐carprofen and (S)‐carprofen towards multiple UGT isoforms. The results showed that (S)‐carprofen exhibited stronger inhibition potential than (R)‐carprofen towards UGT2B7. However, no significant difference was observed for the inhibition of (R)‐carprofen and (S)‐carprofen towards other UGT isoforms. Furthermore, the inhibition kinetic behavior was compared for the inhibition of (S)‐carprofen and (R)‐carprofen towards UGT2B7. A Lineweaver–Burk plot showed that both (S)‐carprofen and (R)‐carprofen exhibited competitive inhibition towards UGT2B7‐catalyzed 4‐MU glucuronidation. The inhibition kinetic parameter (K_i) was calculated to be 7.0 μM and 31.1 μM for (S)‐carprofen and (R)‐carprofen, respectively. Based on the standard for drug–drug interaction, the threshold for (S)‐carprofen and (R)‐carprofen to induce a drug–drug interaction is 0.7 μM and 3.1 μM, respectively. In conclusion, enantioselective inhibition of carprofen towards UDP‐glucuronosyltransferase (UGT) 2B7 was demonstrated in the present study. Using the in vitro inhibition kinetic parameter, the concentration threshold of (S)‐carprofen and (R)‐carprofen to possibly induce the drug–drug interaction was obtained. Therefore, clinical monitoring of the plasma concentration of (S)‐carprofen is more important than (R)‐carprofen to avoid a possible drug–drug interaction between carprofen and the drugs mainly undergoing UGT2B7‐catalyzed metabolism. Chirality 27:189–193, 2015. © 2014 Wiley Periodicals, Inc.     

Chirality Influence of Zaltoprofen Towards UDP‐Glucuronosyltransferases (UGTs) Inhibition Potential

Zaltoprofen (ZLT) is a nonsteroidal antiinflammation drug, and has been clinically employed to treat rheumatoid arthritis, osteoarthritis, and other chronic inflammatory pain conditions. The present study aims to investigate the chirality influence of zaltoprofen towards the inhibition potential towards UDP‐glucuronosyltransferases (UGTs) isoforms. In vitro a recombinant UGT isoforms‐catalyzed 4‐methylumbelliferone (4‐MU) glucuronidation incubation system was employed to investigate the inhibition of (R)‐zaltoprofen and (S)‐zaltoprofen towards UGT isoforms. The inhibition difference capability was observed for the inhibition of (R)‐zaltoprofen and (S)‐zaltoprofen towards UGT1A8 and UGT2B7, but not for other tested UGT isoforms. (R)‐zaltoprofen exhibited noncompetitive inhibition towards UGT1A8 and competitive inhibition towards UGT2B7. The inhibition kinetic parameters were calculated to be 35.3 μM and 19.2 μM for UGT1A8 and UGT2B7. (R)‐zaltoprofen and (S)‐zaltoprofen exhibited a different inhibition type towards UGT1A7. Based on the reported maximum plasma concentration of (R)‐zaltoprofen in vivo, a high drug–drug interaction between (R)‐zaltoprofen and the drugs mainly undergoing UGT1A7, UGT1A8, and UGT2B7‐catalyzed glucuronidation was indicated. Chirality 27:359–363, 2015. © 2015 Wiley Periodicals, Inc.     

Roles of conserved tryptophans in trimerization of HIV‐1 membrane‐proximal external regions: Implications for virucidal design via alchemical free‐energy molecular simulations

The Dual‐Action Virolytic Entry Inhibitors, or “DAVEI's,” are a class of recombinant fusions of a lectin, a linker polypeptide, and a 15‐residue fragment from the membrane‐proximal external region (MPER) of HIV‐1 gp41. DAVEI's trigger rupture of HIV‐1 virions, and the interaction site between DAVEI MPER and HIV‐1 lies in the gp41 component of the envelope glycoprotein Env. Here, we explore the hypothesis that DAVEI MPER engages Env gp41 in a mode structurally similar to a crystallographic MPER trimer. We used alchemical free‐energy perturbation to assess the thermodynamic roles of each of the four conserved tryptophan residues on each protomer of MPER₃. We found that a W666A mutation had a large positive for all three protomers, while W672A had a large positive for only two of the three protomers, with the other tryptophans remaining unimportant contributors to MPER₃ stability. The protomer for which W672 is not important is unique in the placement of its W666 sidechain between the other two protomers. We show that the unique orientation of this W666 sidechain azimuthally rotates its protomer away from the orientation it would have if the trimer were symmetric, resulting in the diminished interaction of this W672 with the rest of MPER₃. Our findings are consistent with our previous experimental study of W‐to‐A mutants of DAVEI. This suggests that DAVEI MPER may engage HIV‐1 Env to form a mixed trimer state in which one DAVEI MPER forms a trimer by displacing a more weakly interacting protomer of the endogenous Env MPER trimer.     

Comparison of SARS and NL63 papain-like protease binding sites and binding site dynamics: inhibitor design implications

The human severe acute respiratory syndrome coronavirus (SARS-CoV) and the NL63 coronaviruses are human respiratory pathogens for which no effective antiviral treatment exists. The papain-like cysteine proteases encoded by the coronavirus (SARS-CoV: PLpro; NL63: PLP1 and PLP2) represent potential targets for antiviral drug development. Three recent inhibitor-bound PLpro structures highlight the role of an extremely flexible six-residue loop in inhibitor binding. The high binding site plasticity is a major challenge in computational drug discovery/design efforts. From conventional molecular dynamics and accelerated molecular dynamics (aMD) simulations, we find that with conventional molecular dynamics simulation, PLpro translationally samples the open and closed conformation of BL2 loop on a picosecond-nanosecond timescale but does not reproduce the peptide bond inversion between loop residues Tyr269 and Gln270 that is observed on inhibitor GRL0617 binding. Only aMD simulation, starting from the closed loop conformation, reproduced the 180° ?-ψ dihedral rotation back to the open loop state. The Tyr-Gln peptide bond inversion appears to involve a progressive conformational change of the full loop, starting at one side, and progressing to the other. We used the SARS-CoV apo X-ray structure to develop a model of the NL63-PLP2 catalytic site. Superimposition of the PLP2 model on the PLpro X-ray structure identifies binding site residues in PLP2 that contribute to the distinct substrate cleavage site specificities between the two proteases. The topological and electrostatic differences between the two protease binding sites also help explain the selectivity of non-covalent PLpro inhibitors.     

Stimulus‐dependent dynamics of p53 in single cells

Many biological networks respond to various inputs through a common signaling molecule that triggers distinct cellular outcomes. One potential mechanism for achieving specific input–output relationships is to trigger distinct dynamical patterns in response to different stimuli. Here we focused on the dynamics of p53, a tumor suppressor activated in response to cellular stress. We quantified the dynamics of p53 in individual cells in response to UV and observed a single pulse that increases in amplitude and duration in proportion to the UV dose. This graded response contrasts with the previously described series of fixed pulses in response to γ‐radiation. We further found that while γ‐triggered p53 pulses are excitable, the p53 response to UV is not excitable and depends on continuous signaling from the input‐sensing kinases. Using mathematical modeling and experiments, we identified feedback loops that contribute to specific features of the stimulus‐dependent dynamics of p53, including excitability and input‐duration dependency. Our study shows that different stresses elicit different temporal profiles of p53, suggesting that modulation of p53 dynamics might be used to achieve specificity in this network.