Conformational Biases of Linear Motifs

Conformational dynamics of Tetracenomycin aromatase/cyclase regulate polyketide binding and enzyme aggregation propensity

BACKGROUND

The N-terminal domain of Tetracenomycin aromatase/cyclase (TcmN), an enzyme derived from Streptomyces glaucescens, is involved in polyketide cyclization, aromatization, and folding. Polyketides are a diverse class of secondary metabolites produced by certain groups of bacteria, fungi, and plants with various pharmaceutical applications. Examples include antibiotics, such as tetracycline, and anticancer drugs, such as doxorubicin. Because TcmN is a promising enzyme for in vitro production of polyketides, it is important to identify conditions that enhance its thermal resistance and optimize its function.

METHODS

TcmN unfolding, stability, and dynamics were evaluated by fluorescence spectroscopy, circular dichroism, nuclear magnetic resonance ¹⁵N relaxation experiments, and microsecond molecular dynamics (MD) simulations.

RESULTS

TcmN thermal resistance was enhanced at low protein and high salt concentrations, was pH-dependent, and denaturation was irreversible. Conformational dynamics on the μs-ms timescale were detected for residues in the substrate-binding cavity, and two predominant conformers representing opened and closed cavity states were observed in the MD simulations.

CONCLUSION

Based on the results, a mechanism was proposed in which the thermodynamics and kinetics of the TcmN conformational equilibrium modulate enzyme function by favoring ligand binding and avoiding aggregation.

GENERAL SIGNIFICANCE

Understanding the principles underlying TcmN stability and dynamics may help in designing mutants with optimal properties for biotechnological applications.

The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.

Prediction of Binding Energy of Keap1 Interaction Motifs in the Nrf2 Antioxidant Pathway and Design of Potential High-Affinity Peptides

Nuclear factor erythroid 2-related factor 2 (Nrf2) is a transcription factor and principal regulator of the antioxidant pathway. The Kelch domain of Kelch-like ECH-associated protein 1 (Keap1) binds to motifs in the N-terminal region of Nrf2, promoting its degradation. There is interest in developing ligands that can compete with Nrf2 for binding to Kelch, thereby activating its transcriptional activities and increasing antioxidant levels. Using experimental Δ G_bind values of Kelch-binding motifs determined previously, a revised hydrophobicity-based model was developed for estimating Δ G_bind from amino acid sequence and applied to rank potential uncharacterized Kelch-binding motifs identified from interaction databases and BLAST searches. Model predictions and molecular dynamics (MD) simulations suggested that full-length MAD2A binds Kelch more favorably than a high-affinity 20-mer Nrf2 E78P peptide, but that the motif in isolation is not a particularly strong binder. Endeavoring to develop shorter peptides for activating Nrf2, new designs were created based on the E78P peptide, some of which showed considerable propensity to form binding-competent structures in MD, and were predicted to interact with Kelch more favorably than the E78P peptide. The peptides could be promising new ligands for enhancing the oxidative stress response.

Investigation of the binding mode of a novel cruzain inhibitor by docking, molecular dynamics, ab initio and MM/PBSA calculations

Chagas disease remains a major health problem in South America, and throughout the world. The two drugs clinically available for its treatment have limited efficacy and cause serious adverse effects. Cruzain is an established therapeutic target of Trypanosoma cruzi, the protozoan that causes Chagas disease. Our group recently identified a competitive cruzain inhibitor (compound 1) with an IC₅₀ = 15 µM that is also more synthetically accessible than the previously reported lead, compound 2. Prior studies, however, did not propose a binding mode for compound 1, hindering understanding of the structure-activity relationship and optimization. Here, the cruzain binding mode of compound 1 was investigated using docking, molecular dynamics (MD) simulations with ab initio derived parameters, ab initio calculations, and MM/PBSA. Two ligand protonation states and four binding poses were evaluated. A careful ligand parameterization method was employed to derive more physically meaningful parameters than those obtained by automated tools. The poses of unprotonated 1 were unstable in MD, showing large conformational changes and diffusing away from the binding site, whereas the protonated form showed higher stability and interaction with negatively charged residues Asp161 and Cys25. MM/PBSA also suggested that these two residues contribute favorably to binding of compound 1. By combining results from MD, ab initio calculations, and MM/PBSA, a binding mode of 1 is proposed. The results also provide insights for further optimization of 1, an interesting lead compound for the development of new cruzain inhibitors.

Personalized Peptide Arrays for Detection of HLA Alloantibodies in Organ Transplantation

In organ transplantation, the function and longevity of the graft critically rely on the success of controlling immunological rejection reactivity against human leukocyte antigens (HLA). Histocompatibility guidelines are based on laboratory tests of anti-HLA immunity, which presents either as pre-existing or de novo generated HLA antibodies that constitute a major transplantation barrier. Current tests are built on a single-antigen beads (SAB) platform using a fixed set of ~100 preselected recombinant HLA antigens to probe transplant sera. However, in humans there exist a far greater variety of HLA types, with no two individuals other than identical twins who can share the same combination of HLA sequences. While advanced technologies for HLA typing and direct sequencing can precisely capture any mismatches in DNA sequence between a donor's and recipient's HLA, the SAB assay, due to its limited variety in sequence representation, is unable to precisely detect alloantibodies specifically against the donor HLA mismatches. We sought to develop a complementary method using a different technology to detect and characterize anti-donor HLA antibodies on a personalized basis. The screening tool is a custom peptide array of donor HLA-derived sequences for probing post-transplant sera of the organ recipient to assess the risk for antibody-mediated rejection. On a single array for one donor-recipient pair, up to 600 unique peptides are made based on the donor's HLA protein sequences, each peptide carrying at least one mismatched residue in a 15-amino acid sequence. In our pilot experiments to compare antigen patterns for pre- and post-transplant sera on these arrays, we were able to detect anti-HLA signals with the resolution that also allowed us to pinpoint the immune epitopes involved. These personalized antigen arrays allow high-resolution detection of donor-specific HLA epitopes in organ transplantation.

Aggregation and Prion-Like Properties of Misfolded Tumor Suppressors: Is Cancer a Prion Disease?

Prion diseases are disorders that share several characteristics that are typical of many neurodegenerative diseases. Recently, several studies have extended the prion concept to pathological aggregation in malignant tumors involving misfolded p53, a tumor-suppressor protein. The aggregation of p53 and its coaggregation with p53 family members, p63 and p73, have been shown. Certain p53 mutants exert a dominant-negative regulatory effect on wild-type (WT) p53. The basis for this dominant-negative effect is that amyloid-like mutant p53 converts WT p53 into an aggregated species, leading to a gain-of-function (GoF) phenotype and the loss of its tumor-suppressor function. Recently, it was shown that p53 aggregates can be internalized by cells and can coaggregate with endogenous p53, corroborating the prion-like properties of p53 aggregates. The prion-like behavior of oncogenic p53 mutants provides an explanation for its dominant-negative and GoF properties, including the high metastatic potential of cancer cells carrying p53 mutations. The inhibition of p53 aggregation appears to represent a promising target for therapeutic intervention in patients with malignant tumors.

Aggregation tendencies in the p53 family are modulated by backbone hydrogen bonds

The p53 family of proteins is comprised of p53, p63 and p73. Because the p53 DNA binding domain (DBD) is naturally unstable and possesses an amyloidogenic sequence, it is prone to form amyloid fibrils, causing loss of functions. To develop p53 therapies, it is necessary to understand the molecular basis of p53 instability and aggregation. Light scattering, thioflavin T (ThT) and high hydrostatic pressure (HHP) assays showed that p53 DBD aggregates faster and to a greater extent than p63 and p73 DBDs, and was more susceptible to denaturation. The aggregation tendencies of p53, p63, and p73 DBDs were strongly correlated with their thermal stabilities. Molecular Dynamics (MD) simulations indicated specific regions of structural heterogeneity unique to p53, which may be promoted by elevated incidence of exposed backbone hydrogen bonds (BHBs). The results indicate regions of structural vulnerability in the p53 DBD, suggesting new targetable sites for modulating p53 stability and aggregation, a potential approach to cancer therapy.

Milk proteins interact with goat Binder of SPerm (BSP) proteins and decrease their binding to sperm

Seminal plasma Binder of SPerm (BSP) proteins bind to sperm at ejaculation and promote capacitation. When in excess, however, BSP proteins damage the sperm membrane. It has been suggested that milk components of semen extenders associate with BSP proteins, potentially protecting sperm. Thus, this study was conducted to investigate if milk proteins interact with BSP proteins and reduce BSP binding to goat sperm. Using gel filtration chromatography, milk was incubated with goat seminal plasma proteins and loaded onto columns with and without calcium. Milk was also fractionated into parts containing mostly whey proteins or mostly caseins, incubated with seminal plasma proteins and subjected to gel filtration. Eluted fractions were evaluated by immunoblot using anti-goat BSP antibodies, confirming milk protein-BSP protein interactions. As determined by ELISA, milk proteins coated on polystyrene wells bound to increasing of goat BSP proteins. Far-western dot blots confirmed that BSP proteins bound to caseins and β-lactoglobulin in a concentration-dependent manner. Then, cauda epididymal sperm from five goats was incubated with seminal plasma; seminal plasma followed by milk; and milk followed by seminal plasma. Sperm membrane proteins were extracted and evaluated by immunoblotting. The pattern of BSP binding to sperm membrane proteins was reduced by 59.3 % when epididymal sperm were incubated with seminal plasma and then with skimmed milk (p < 0.05). When epididymal sperm were treated with milk followed by seminal plasma, coating of sperm with BSP proteins was not significantly reduced (57.6 %; p > 0.05). In conclusion, goat BSP proteins have an affinity for caseins and whey proteins. Milk reduces BSP binding to goat sperm, depending whether or not sperm had been previously exposed to seminal plasma. Such events may explain the protective effect of milk during goat sperm preservation.

Binding of Disordered Peptides to Kelch: Insights from Enhanced Sampling Simulations

Keap1 protein plays an essential role in regulating cellular oxidative stress response and is a crucial binding hub for multiple proteins, several of which are intrinsically disordered proteins (IDP). Among Kelch's IDP binding partners, NRF2 and PTMA are the two most interesting cases. They share a highly similar binding motif; however, NRF2 binds to Kelch with a binding affinity of approximately 100-fold higher than that of PTMA. In this study, we perform an exhaustive sampling composed of 6 μs well-tempered metadynamics and 2 μs unbiased molecular dynamics (MD) simulations aiming at characterizing the binding mechanisms and structural properties of these two peptides. Our results agree with previous experimental observations that PTMA is remarkably more disordered than NRF2 in both the free and bound states. This explains PTMA's lower binding affinity. Our extensive sampling also provides valuable insights into the vast conformational ensembles of both NRF2 and PTMA, supports the hypothesis of coupled folding-binding, and confirms the essential role of linear motifs in IDP binding.

Accelerating the Conformational Sampling of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are a class of proteins lacking a well-defined secondary structure. Instead, they are able to attain multiple conformations, bind to multiple targets, and respond to changes in their surroundings. Functionally, IDPs have been associated with molecular recognition, cell regulation, and signal transduction. The dynamic conformational ensemble of IDPs is highly environmental and binding partner dependent, rendering the characterization of IDPs extremely challenging. Here, we compare the sampling efficiencies of conventional molecular dynamics (MD), well-tempered metadynamics (WT-META), and bias-exchange metadynamics (BE-META). The total simulation time was over 10 μs, and a 20-mer peptide derived from the Neh2 domain of the Nuclear factor erythroid 2-related factor 2 (Nrf2) protein was simulated. BE-META, with a neutral replica and seven biased replicas employing a set of seven relevant collective variables (CVs), provided the most reliable and efficient sampling. Finally, we propose a free-energy reconstruction method based on the probability distribution of the secondary structure contents. This postprocessing analysis confirms the presence of not only the β-hairpin conformation of the free Neh2 peptide but also its rare bound-state-like conformation, both of that have been experimentally observed. In addition, our simulations also predict other possible conformations to be verified with future experiments.