Conformational dynamics of threonine 195 and the S1 subsite in functional trypsin variants

Dimerization misalignment in human glutamate-oxaloacetate transaminase variants is the primary factor for PLP release

The active form of vitamin B6, pyridoxal 5’-phosphate (PLP), plays an essential role in the catalytic mechanism of various proteins, including human glutamate-oxaloacetate transaminase (hGOT1), an important enzyme in amino acid metabolism. A recent molecular and genetic study showed that the E266K, R267H, and P300L substitutions in aspartate aminotransferase, the Arabidopsis analog of hGOT1, genetically suppress a developmentally arrested Arabidopsis RUS mutant. Furthermore, CD analyses suggested that the variants exist as apo proteins and implicated a possible role of PLP in the regulation of PLP homeostasis and metabolic pathways. In this work, we assessed the stability of PLP bound to hGOT1 for the three variant and wildtype (WT) proteins using a combined 6 μs of molecular dynamics (MD) simulation. For the variants and WT in the holo form, the MD simulations reproduced the “closed-open” transition needed for substrate binding. This conformational transition was associated with the rearrangement of the P15-R32 small domain loop providing substrate access to the R387/R293 binding motif. We also showed that formation of the dimer interface is essential for PLP affinity to the active site. The position of PLP in the WT binding site was stabilized by a unique hydrogen bond network of the phosphate binding cup, which placed the cofactor for formation of the covalent Schiff base linkage with K259 for catalysis. The amino acid substitutions at positions 266, 267, and 300 reduced the structural correlation between PLP and the protein active site and/or integrity of the dimer interface. Principal component analysis and energy decomposition clearly suggested dimer misalignment and dissociation for the three variants tested in our work. The low affinity of PLP in the hGOT1 variants observed in our computational work provided structural rationale for the possible role of vitamin B6 in regulating metabolic pathways.

Engineering trypsin for inhibitor resistance

The development of effective protease therapeutics requires that the proteases be more resistant to naturally occurring inhibitors while maintaining catalytic activity. A key step in developing inhibitor resistance is the identification of key residues in protease-inhibitor interaction. Given that majority of the protease therapeutics currently in use are trypsin-fold, trypsin itself serves as an ideal model for studying protease-inhibitor interaction. To test the importance of several trypsin-inhibitor interactions on the prime-side binding interface, we created four trypsin single variants Y39A, Y39F, K60A, and K60V and report biochemical sensitivity against bovine pancreatic trypsin inhibitor (BPTI) and M84R ecotin. All variants retained catalytic activity against small, commercially available peptide substrates [kcat /KM  = (1.2 ± 0.3) × 10(7) M(-1 ) s(-1) . Compared with wild-type, the K60A and K60V variants showed increased sensitivity to BPTI but less sensitivity to ecotin. The Y39A variant was less sensitive to BPTI and ecotin while the Y39F variant was more sensitive to both. The relative binding free energies between BPTI complexes with WT, Y39F, and Y39A were calculated based on 3.5 µs combined explicit solvent molecular dynamics simulations. The BPTI:Y39F complex resulted in the lowest binding energy, while BPTI:Y39A resulted in the highest. Simulations of Y39F revealed increased conformational rearrangement of F39, which allowed formation of a new hydrogen bond between BPTI R17 and H40 of the variant. All together, these data suggest that positions 39 and 60 are key for inhibitor binding to trypsin, and likely more trypsin-fold proteases.

Why Ser and not Thr brokers catalysis in the trypsin fold

Although Thr is equally represented as Ser in the human genome and as a nucleophile is as good as Ser, it is never found in the active site of the large family of trypsin-like proteases that utilize the Asp/His/Ser triad. The molecular basis of the preference of Ser over Thr in the trypsin fold was investigated with X-ray structures of the thrombin mutant S195T free and bound to an irreversible active site inhibitor. In the free form, the methyl group of T195 is oriented toward the incoming substrate in a conformation seemingly incompatible with productive binding. In the bound form, the side chain of T195 is reoriented for efficient substrate acylation without causing steric clash within the active site. Rapid kinetics prove that this change is due to selection of an active conformation from a preexisting ensemble of reactive and unreactive rotamers whose relative distribution determines the level of activity of the protease. Consistent with these observations, the S195T substitution is associated with a weak yet finite activity that allows identification of an unanticipated important role for S195 as the end point of allosteric transduction in the trypsin fold. The S195T mutation abrogates the Na(+)-dependent enhancement of catalytic activity in thrombin, activated protein C, and factor Xa and significantly weakens the physiologically important allosteric effects of thrombomodulin on thrombin and of cofactor Va on factor Xa. The evolutionary selection of Ser over Thr in trypsin-like proteases was therefore driven by the need for high catalytic activity and efficient allosteric regulation.

Evaluation of benzoic acid derivatives as sirtuin inhibitors

Employing a genetically modified yeast strain as a screening tool, 4-dimethylaminobenzoic acid (5) was isolated from the marine sediment-derived Streptomyces sp. CP27-53 as a weak yeast sirtuin (Sir2p) inhibitor. Using this compound as a scaffold, a series of disubstituted benzene derivatives were evaluated to elucidate the structure activity relationships for Sir2p inhibition. The results suggested that 4-alkyl or 4-alkylaminobenzoic acid is the key structure motif for Sir2p inhibitory activity. The most potent Sir2p inhibitor, 4-tert-butylbenzoic acid (20), among the tested compounds in this study turned out to be a weak but selective SIRT1 inhibitor. The calculated binding free energies between the selected compounds and the catalytic domain of SIRT1 were well correlated to their measured SIRT1 inhibitory activities.

Protein dynamics via computational microscope

The purpose of this overview is to provide a concise introduction to the methodology and current advances in molecular dynamics (MD) simulations. MD simulations emerged as a powerful and popular tool to study dynamic behavior of proteins and macromolecule complexes at the atomic resolution. This approach can extend static structural data, such as X-ray crystallography, into dynamic domains with realistic timescales (up to millisecond) and high precision, therefore becoming a veritable computational microscope. This perspective covers current advances and methodology in the simulation of protein folding and drug design as illustrated by several important published examples. Overall, recent progress in the simulation field points to the direction that MD will have significant impact on molecular biology and pharmaceutical science.