Modulating the DNA Affinity of Elk-1 with Computationally Selected Mutations

Molecular dynamics simulations and in silico peptide ligand screening of the Elk-1 ETS domain

BACKGROUND

The Elk-1 transcription factor is a member of a group of proteins called ternary complex factors, which serve as a paradigm for gene regulation in response to extracellular signals. Its deregulation has been linked to multiple human diseases including the development of tumours. The work herein aims to inform the design of potential peptidomimetic compounds that can inhibit the formation of the Elk-1 dimer, which is key to Elk-1 stability. We have conducted molecular dynamics simulations of the Elk-1 ETS domain followed by virtual screening.

RESULTS

We show the ETS dimerisation site undergoes conformational reorganisation at the α1β1 loop. Through exhaustive screening of di- and tri-peptide libraries against a collection of ETS domain conformations representing the dynamics of the loop, we identified a series of potential binders for the Elk-1 dimer interface. The di-peptides showed no particular preference toward the binding site; however, the tri-peptides made specific interactions with residues: Glu17, Gln18 and Arg49 that are pivotal to the dimer interface.

CONCLUSIONS

We have shown molecular dynamics simulations can be combined with virtual peptide screening to obtain an exhaustive docking protocol that incorporates dynamic fluctuations in a receptor. Based on our findings, we suggest experimental binding studies to be performed on the 12 SILE ranked tri-peptides as possible compounds for the design of inhibitors of Elk-1 dimerisation. It would also be reasonable to consider the score-ranked tri-peptides as a comparative test to establish whether peptide size is a determinant factor of binding to the ETS domain.

The intrinsic protein flexibility of endogenous protease inhibitor TIMP-1 controls its binding interface and affects its function

Protein flexibility is thought to play key roles in numerous biological processes, including antibody affinity maturation, signal transduction, and enzyme catalysis, yet only limited information is available regarding the molecular details linking protein dynamics with function. A single point mutation at the distal site of the endogenous tissue inhibitor of metalloproteinase 1 (TIMP-1) enables this clinical target protein to tightly bind and inhibit membrane type 1 matrix metalloproteinase (MT1-MMP) by increasing only the association constant. The high-resolution X-ray structure of this complex determined at 2 A could not explain the mechanism of enhanced binding and pointed to a role for protein conformational dynamics. Molecular dynamics (MD) simulations reveal that the high-affinity TIMP-1 mutants exhibit significantly reduced binding interface flexibility and more stable hydrogen bond networks. This was accompanied by a redistribution of the ensemble of substrates to favorable binding conformations that fit the enzyme catalytic site. Apparently, the decrease in backbone flexibility led to a lower entropy cost upon formation of the complex. This work quantifies the effect of a single point mutation on the protein conformational dynamics and function of TIMP-1. Here we argue that controlling the intrinsic protein dynamics of MMP endogenous inhibitors may be utilized for rationalizing the design of selective novel protein inhibitors for this class of enzymes.

Molecular evolution of affinity and flexibility in the immune system

The immune system responds to the introduction of foreign antigens by rapidly evolving antibodies with increasing affinity for the antigen (i.e., maturation). To investigate the factors that control this process at the molecular level, we have assessed the changes in flexibility that accompany ligand binding at four stages of maturation in the 4-4-20 antibody. Our studies, based on molecular dynamics, indicate that increased affinity for the target ligand is associated with a decreased entropic cost to binding. The entropy of binding is unfavorable, opposing favorable enthalpic contributions that arise during complex formation. Computed binding free energies for the various antibody-ligand complexes qualitatively reproduce the trends observed in the experimentally derived values, although the absolute magnitude of free-energy differences is overestimated. Our results support the existence of a correlation between high-affinity interactions and decreased protein flexibility in this series of antibody molecules. This observation is likely to be a general feature of molecular association processes and key to the molecular evolution of antibody responses.

Design of Functional Ferritin-Like Proteins with Hydrophobic Cavities

Ferritin four-helix bundle subunits self-assemble to create a stable multimer with a large central hydrophilic cavity where metal ions bind. To explore the versatility of this reaction vessel, computational design was used to generate cavities with increasingly apolar surface areas inside a dodecameric ferritin-like protein, Dps. Cavity mutants, in which as many as 120 surface accessible hydrophilic residues were replaced with hydrophobic amino acids, were shown to still assemble properly using size-exclusion chromatography and dynamic light scattering measurements. Wild-type Dps exhibited highly cooperative subunit folding and assembly, which was monitored by changes in Trp fluorescence and UV circular dichroism. The hydrophobic cavity mutants showed distinctly less cooperative unfolding behavior, with one mutant forming a partially assembled intermediate upon guanidine denaturation. Although the stability of Dps to such denaturation decreased with increasing apolar surface area, all proteins exhibited high melting temperatures, T(m) = 74-90 degrees C. Despite the large number of mutations, near-native ability to mineralize iron was maintained. This work illustrates the versatility of the ferritin scaffold for engineering large protein cavities with novel properties.