Predictions of structural elements for the binding of Hin recombinase with the hix site of DNA

Molecular dynamics simulations were coupled with experimental data from biochemistry and genetics to generate a theoretical structure for the binding domain of Hin recombinase complexed with the hix site of DNA. The theoretical model explains the observed sequence specificity of Hin recombinase and leads to a number of testable predictions concerning altered sequence selectivity for various mutants of protein and DNA. Combining molecular dynamics simulations with constraints based on current knowledge of protein structure leads to a theoretical structure of the binding domain of Hin recombinase with the hix site of DNA. The model offers a mechanistic explanation of the presently known characteristics of Hin and predicts the effects of specific mutations of both protein and DNA. The predictions can be tested by currently feasible experiments that should lead to refinements in and improvements on the current theoretical model. Because current experimental and theoretical methods are all limited to providing only partial information about protein-DNA interactions, we believe that this approach of basing molecular simulations on experimental knowledge and using the results of these simulations to design new, more precise experimental tests will be of general utility. These results provide additional evidence for the generality of the helix-turn-helix motif in DNA recognition and stabilization of proteins on DNA

A microscopic description of acid-base equilibrium

Acid-base reactions are ubiquitous in nature. Understanding their mechanisms is crucial in many fields, from biochemistry to industrial catalysis. Unfortunately, experiments only give limited information without much insight into the molecular behaviour. Atomistic simulations could complement experiments and shed precious light on microscopic mechanisms. The large free energy barriers connected to proton dissociation however make the use of enhanced sampling methods mandatory. Here we perform an ab initio molecular dynamics (MD) simulation and enhance sampling with the help of methadynamics. This has been made possible by the introduction of novel descriptors or collective variables (CVs) that are based on a conceptually new outlook on acid-base equilibria. We test successfully our approach on three different aqueous solutions of acetic acid, ammonia, and bicarbonate. These are representative of acid, basic, and amphoteric behaviour

Critical Role of the Secondary Binding Pocket in Modulating the Enzymatic Activity of DUSP5 toward Phosphorylated ERKs

DUSP5 is an inducible nuclear dual-specificity phosphatase that specifically interacts with and deactivates extracellular signal-regulated kinases ERK1 and ERK2, which are responsible for cell proliferation, differentiation, and survival. The phosphatase domain (PD) of DUSP5 has unique structural features absent from other nuclear DUSPs, such as the presence of a secondary anion-binding site in the proximity of the reaction center and a glutamic acid E264 positioned next to the catalytic cysteine C263, as well as a remote intramolecular disulfide linkage. The overall 400 ns molecular dynamics simulations indicate that the secondary binding site of DUSP5 PD acts as an allosteric regulator of the phosphatase activity of DUSP5. Our studies have identified E264 as a critical constituent of the dual binding pocket, which regulates the catalytic activity of DUSP5 by forming a salt bridge with arginine R269. Molecular dynamics studies showed that initial occupation of the secondary binding pocket leads to the breakage of the salt bridge, which then allows the occupation of the active site. Indeed, biochemical analysis using the pERK assay on mutant E264Q demonstrated that mutation of glutamic acid E264 leads to an increase in the DUSP5 catalytic activity. The role of the secondary binding site in assembling the DUSP5–pERK pre-reactive complex was further demonstrated by molecular dynamics simulations that showed that the remote C197–C219 disulfide linkage controls the structure of the secondary binding pocket based on its redox state (i.e., disulfide/dithiol) and, in turn, the enzymatic activity of DUSP5

Pathogenic mutations in the hydrophobic core of the human prion protein can promote structural instability and misfolding

Transmissible spongiform encephalopathies, or prion diseases, are caused by misfolding and aggregation of the prion protein PrP. These diseases can be hereditary in humans and four of the many disease-associated missense mutants of PrP are in the hydrophobic core: V180I, F198S, V203I and V210I. The T183A mutation is related to the hydrophobic core mutants as it is close to the hydrophobic core and known to cause instability. We have performed extensive molecular dynamics simulations of these five PrP mutants and compared their dynamics and conformations to wild-type PrP. The simulations highlight the changes that occur upon introduction of mutations and help to rationalize experimental findings. Changes can occur around the mutation site, but they can also be propagated over long distances. In particular, the F198S and T183A mutations lead to increased flexibility in parts of the structure that are normally stable, and the short β-sheet moves away from the rest of the protein. Mutations V180I, V210I and, to a lesser extent, V203I cause changes similar to those observed upon lowering the pH, which has been linked to misfolding. Early misfolding is observed in one V180I simulation. Overall, mutations in the hydrophobic core have a significant effect on the dynamics and stability of PrP, including the propensity to misfold, which helps to explain their role in the development of familial prion diseases