Comparison of multiple crystal structures with NMR data for engrailed homeodomain

Anionic lipids induce a fold-unfold transition in the membrane-translocating Engrailed homeodomain

Homeoprotein transcription factors have the property of interacting with membranes through their DNA-binding homeodomain, which is involved in unconventional internalization and secretion. Both processes depend on membrane-translocating events but their detailed molecular mechanisms are still poorly understood. We have previously characterized the conformational properties of Engrailed 2 homeodomain (EnHD) in aqueous solution and in micelles as membrane-mimetic environments. In the present study, we used small isotropic lipid bicelles as a more relevant membrane-mimetic model to characterize the membrane-bound state of EnHD. We show that lipid bicelles, in contrast to micelles, adequately reproduce the requirement of anionic lipids in the membrane binding and conformational transition of EnHD. The fold-unfold transition of EnHD induced by anionic lipids was characterized by NMR using ¹H, ¹³C, ¹⁵N chemical shifts, nuclear Overhauser effects, residual dipolar couplings, intramolecular and intermolecular paramagnetic relaxation enhancements induced by site-directed spin-label or paramagnetic lipid probe, respectively. A global unpacking of EnHD helices is observed leading to a loss of the native fold. However, near-native propensities of EnHD backbone conformation are maintained in membrane environment, including not only the three helices but also the turn connecting helices H2 and H3. NMR and coarse-grained molecular dynamics simulations reveal that the EnHD adopts a shallow insertion in the membrane, with the three helices oriented parallel to the membrane. EnHD explores extended conformations and closed U-shaped conformations, which are stabilized by anionic lipid recruitment.

A Dynamic Hydrophobic Core and Surface Salt Bridges Thermostabilize a Designed Three-Helix Bundle

Thermostable proteins are advantageous in industrial applications, as pharmaceuticals or biosensors, and as templates for directed evolution. As protein-design methodologies improve, bioengineers are able to design proteins to perform a desired function. Although many rationally designed proteins end up being thermostable, how to intentionally design de novo, thermostable proteins is less clear. UVF is a de novo-designed protein based on the backbone structure of the Engrailed homeodomain (EnHD) and is highly thermostable (T_m > 99°C vs. 52°C for EnHD). Although most proteins generally have polar amino acids on their surfaces and hydrophobic amino acids buried in their cores, protein engineers followed this rule exactly when designing UVF. To investigate the contributions of the fully hydrophobic core versus the fully polar surface to UVF’s thermostability, we built two hybrid, chimeric proteins combining the sets of buried and surface residues from UVF and EnHD. Here, we determined a structural, dynamic, and thermodynamic explanation for UVF’s thermostability by performing 4 μs of all-atom, explicit-solvent molecular dynamics simulations at 25 and 100°C, Tanford-Kirkwood solvent accessibility Monte Carlo electrostatic calculations, and a thermodynamic analysis of 40 temperature runs by the weighted-histogram analysis method of heavy-atom, structure-based models of UVF, EnHD, and both chimeric proteins. Our models showed that UVF was highly dynamic because of its fully hydrophobic core, leading to a smaller loss of entropy upon folding. The charged residues on its surface made favorable electrostatic interactions that contributed enthalpically to its thermostability. In the chimeric proteins, both the hydrophobic core and charged surface independently imparted thermostability.

Validating Molecular Dynamics Simulations against Experimental Observables in Light of Underlying Conformational Ensembles

Far from the static, idealized conformations deposited into structural databases, proteins are highly dynamic molecules that undergo conformational changes on temporal and spatial scales that may span several orders of magnitude. These conformational changes, often intimately connected to the functional roles that proteins play, may be obscured by traditional biophysical techniques. Over the past 40 years, molecular dynamics (MD) simulations have complemented these techniques by providing the "hidden" atomistic details that underlie protein dynamics. However, there are limitations of the degree to which molecular simulations accurately and quantitatively describe protein motions. Here we show that although four molecular dynamics simulation packages (AMBER, GROMACS, NAMD, and ilmm) reproduced a variety of experimental observables for two different proteins (engrailed homeodomain and RNase H) equally well overall at room temperature, there were subtle differences in the underlying conformational distributions and the extent of conformational sampling obtained. This leads to ambiguity about which results are correct, as experiment cannot always provide the necessary detailed information to distinguish between the underlying conformational ensembles. However, the results with different packages diverged more when considering larger amplitude motion, for example, the thermal unfolding process and conformational states sampled, with some packages failing to allow the protein to unfold at high temperature or providing results at odds with experiment. While most differences between MD simulations performed with different packages are attributed to the force fields themselves, there are many other factors that influence the outcome, including the water model, algorithms that constrain motion, how atomic interactions are handled, and the simulation ensemble employed. Here four different MD packages were tested each using best practices as established by the developers, utilizing three different protein force fields and three different water models. Differences between the simulated protein behavior using two different packages but the same force field, as well as two different packages with different force fields but the same water models and approaches to restraining motion, show how other factors can influence the behavior, and it is incorrect to place all the blame for deviations and errors on force fields or to expect improvements in force fields alone to solve such problems.

Atomistic molecular dynamics simulations of bioactive engrailed 1 interference peptides (EN1-iPeps)

The neural-specific transcription factor Engrailed 1 - is overexpressed in basal-like breast tumours. Synthetic interference peptides - comprising a cell-penetrating peptide/nuclear localisation sequence and the Engrailed 1-specific sequence from the N-terminus have been engineered to produce a strong apoptotic response in tumour cells overexpressing EN1, with no toxicity to normal or non Engrailed 1-expressing cells. Here scaled molecular dynamics simulations were used to study the conformational dynamics of these interference peptides in aqueous solution to characterise their structure and dynamics. Transitions from disordered to α-helical conformation, stabilised by hydrogen bonds and proline-aromatic interactions, were observed throughout the simulations. The backbone of the wild-type peptide folds to a similar conformation as that found in ternary complexes of anterior Hox proteins with conserved hexapeptide motifs important for recognition of pre-B-cell leukemia Homeobox 1, indicating that the motif may possess an intrinsic preference for helical structure. The predicted NMR chemical shifts of these peptides are consistent with the Hox hexapeptides in solution and Engrailed 2 NMR data. These findings highlight the importance of aromatic residues in determining the structure of Engrailed 1 interference peptides, shedding light on the rational design strategy of molecules that could be adopted to inhibit other transcription factors overexpressed in other cancer types, potentially including other transcription factor families that require highly conserved and cooperative protein-protein partnerships for biological activity.

Solution x-ray scattering and structure formation in protein dynamics

We propose a computationally effective approach that builds on Landau mean-field theory in combination with modern nonequilibrium statistical mechanics to model and interpret protein dynamics and structure formation in small- to wide-angle x-ray scattering (S/WAXS) experiments. We develop the methodology by analyzing experimental data in the case of Engrailed homeodomain protein as an example. We demonstrate how to interpret S/WAXS data qualitatively with a good precision and over an extended temperature range. We explain experimental observations in terms of protein phase structure, and we make predictions for future experiments and for how to analyze data at different ambient temperature values. We conclude that the approach we propose has the potential to become a highly accurate, computationally effective, and predictive tool for analyzing S/WAXS data. For this, we compare our results with those obtained previously in an all-atom molecular dynamics simulation.

Protein environment affects the water–tryptophan binding mode. MD, QM/MM, and NMR studies of engrailed homeodomain mutants

Water molecules can interact with the π-face of tryptophan either forming an O–H⋯π hydrogen bond or by a lone-pair⋯π interaction. Surrounding amino acids can favor the one or the other interaction type.

Deconvoluting Protein (Un)folding Structural Ensembles Using X-Ray Scattering, Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation

The folding and unfolding of protein domains is an apparently cooperative process, but transient intermediates have been detected in some cases. Such (un)folding intermediates are challenging to investigate structurally as they are typically not long-lived and their role in the (un)folding reaction has often been questioned. One of the most well studied (un)folding pathways is that of Drosophila melanogaster Engrailed homeodomain (EnHD): this 61-residue protein forms a three helix bundle in the native state and folds via a helical intermediate. Here we used molecular dynamics simulations to derive sample conformations of EnHD in the native, intermediate, and unfolded states and selected the relevant structural clusters by comparing to small/wide angle X-ray scattering data at four different temperatures. The results are corroborated using residual dipolar couplings determined by NMR spectroscopy. Our results agree well with the previously proposed (un)folding pathway. However, they also suggest that the fully unfolded state is present at a low fraction throughout the investigated temperature interval, and that the (un)folding intermediate is highly populated at the thermal midpoint in line with the view that this intermediate can be regarded to be the denatured state under physiological conditions. Further, the combination of ensemble structural techniques with MD allows for determination of structures and populations of multiple interconverting structures in solution.

Bayesian inference of protein structure from chemical shift data

Protein chemical shifts are routinely used to augment molecular mechanics force fields in protein structure simulations, with weights of the chemical shift restraints determined empirically. These weights, however, might not be an optimal descriptor of a given protein structure and predictive model, and a bias is introduced which might result in incorrect structures. In the inferential structure determination framework, both the unknown structure and the disagreement between experimental and back-calculated data are formulated as a joint probability distribution, thus utilizing the full information content of the data. Here, we present the formulation of such a probability distribution where the error in chemical shift prediction is described by either a Gaussian or Cauchy distribution. The methodology is demonstrated and compared to a set of empirically weighted potentials through Markov chain Monte Carlo simulations of three small proteins (ENHD, Protein G and the SMN Tudor Domain) using the PROFASI force field and the chemical shift predictor CamShift. Using a clustering-criterion for identifying the best structure, together with the addition of a solvent exposure scoring term, the simulations suggests that sampling both the structure and the uncertainties in chemical shift prediction leads more accurate structures compared to conventional methods using empirical determined weights. The Cauchy distribution, using either sampled uncertainties or predetermined weights, did, however, result in overall better convergence to the native fold, suggesting that both types of distribution might be useful in different aspects of the protein structure prediction.

Investigation of homeodomain membrane translocation properties: insights from the structure determination of engrailed-2 homeodomain in aqueous and membrane-mimetic environments

In addition to their well-known DNA-binding properties, homeodomains have the ability to efficiently translocate across biological membranes through still poorly-characterized mechanisms. To date, most biophysical studies addressing the mechanisms of internalization have focused on small synthetic peptides rather than full-length globular homeodomains. In this work, we characterized the conformational properties of chicken Engrailed 2 homeodomain (En2HD) in aqueous solution and in membrane mimetic environments using circular dichroism, Trp fluorescence, and NMR spectroscopy. En2HD adopts a well-defined three-helical bundle fold in aqueous solution. The Trp-48 residue, which is critical for internalization, is fully buried in the hydrophobic core. Circular dichroism and fluorescence reveal that a conformational transition occurs in anionic lipid vesicles and in micelles. En2HD loses its native three-dimensional structure in micellar environments but, remarkably, near-native helical secondary structures are maintained. Long-range interactions could be detected using site-directed spin labels, indicating that the three helices do not adopt extended orientations. Noncovalent paramagnetic probes yielded information about helix positioning and unveiled the burial of critical aromatic and basic residues within the micelles. Our results suggest that electrostatic interactions with membranes may be determinant in inducing a conformational change enabling Trp-48 to insert into membranes.

Fast protein folding kinetics

Fast-folding proteins have been a major focus of computational and experimental study because they are accessible to both techniques: they are small and fast enough to be reasonably simulated with current computational power, but have dynamics slow enough to be observed with specially developed experimental techniques. This coupled study of fast-folding proteins has provided insight into the mechanisms, which allow some proteins to find their native conformation well <1 ms and has uncovered examples of theoretically predicted phenomena such as downhill folding. The study of fast folders also informs our understanding of even 'slow' folding processes: fast folders are small; relatively simple protein domains and the principles that govern their folding also govern the folding of more complex systems. This review summarizes the major theoretical and experimental techniques used to study fast-folding proteins and provides an overview of the major findings of fast-folding research. Finally, we examine the themes that have emerged from studying fast folders and briefly summarize their application to protein folding in general, as well as some work that is left to do.

Multimolecule test-tube simulations of protein unfolding and aggregation

Molecular dynamics simulations of protein folding or unfolding, unlike most in vitro experimental methods, are performed on a single molecule. The effects of neighboring molecules on the unfolding/folding pathway are largely ignored experimentally and simply not modeled computationally. Here, we present two all-atom, explicit solvent molecular dynamics simulations of 32 copies of the Engrailed homeodomain (EnHD), an ultrafast-folding and -unfolding protein for which the folding/unfolding pathway is well-characterized. These multimolecule simulations, in comparison with single-molecule simulations and experimental data, show that intermolecular interactions have little effect on the folding/unfolding pathway. EnHD unfolded by the same mechanism whether it was simulated in only water or also in the presence of other EnHD molecules. It populated the same native state, transition state, and folding intermediate in both simulation systems, and was in good agreement with experimental data available for each of the three states. Unfolding was slowed slightly by interactions with neighboring proteins, which were mostly hydrophobic in nature and ultimately caused the proteins to aggregate. Protein-water hydrogen bonds were also replaced with protein-protein hydrogen bonds, additionally contributing to aggregation. Despite the increase in protein-protein interactions, the protein aggregates formed in simulation did not do so at the total exclusion of water. These simulations support the use of single-molecule techniques to study protein unfolding and also provide insight into the types of interactions that occur as proteins aggregate at high temperature at an atomic level.

Promiscuous contacts and heightened dynamics increase thermostability in an engineered variant of the engrailed homeodomain

A thermostabilized variant (UVF) of the engrailed homeodomain (EnHD) was previously engineered by Mayo and co-workers. The melting temperature of the non-natural, designed protein is 50°C higher than the natural wild-type protein (>99 vs. 52°C), and the two proteins share 22% sequence identity. We have performed extensive (1 μs) all-atom, explicit solvent molecular dynamics simulations of the wild-type and engineered proteins to investigate their structural and dynamic properties at room temperature and at 100°C. Our simulations are in good agreement with nuclear magnetic resonance data available for the two proteins [nuclear Overhauser effect crosspeaks (NOEs), J-coupling constants and order parameters for EnHD; and NOEs for UVF], showing that we reproduce the backbone dynamics and side chain packing in the native state of both proteins. UVF was more dynamic at room temperature than EnHD, with respect to both its backbone and side chain motion. When the temperature was raised, the thermostable protein maintained this mobility while retaining its native conformation. EnHD, on the other hand, was unable to maintain its more rigid native structure at higher temperature and began to unfold. Heightened protein dynamics leading to promiscuous and dynamically interchangeable amino acid contacts makes UVF more tolerant to increasing temperature, providing a molecular explanation for heightened thermostability of this protein.

Refolding the engrailed homeodomain: structural basis for the accumulation of a folding intermediate

The ultrafast folding pathway of the engrailed homeodomain has been exceptionally well characterized by experiment and simulation. Helices II and III of the three-helix bundle protein form the native helix-turn-helix motif as an on-pathway intermediate within a few microseconds. The slow step is then the proper docking of the helices in approximately 15 mus. However, there is still the unexplained puzzle of why helix docking is relatively slow, which is part of the more general question as to why rearrangements of intermediates occur slowly. To address this problem, we performed 46 all-atom molecular dynamics refolding simulations in explicit water, for a total of 15 micros of simulation time. The simulations started from an intermediate state structure that was generated in an unfolding simulation at 498 K and was then quenched to folding-permissive temperatures. The protein refolded successfully in only one of the 46 simulations, and in that case the refolding pathway mirrored the unfolding pathway at high temperature. In the 45 simulations in which the protein did not fully fold, nonnative salt bridges trapped the protein, which explains why the protein folds relatively slowly from the intermediate state.