Molecular Dynamics Simulation of Phosphorylated KID Post-Translational Modification

BACKGROUND:Kinase-inducible domain (KID) as transcriptional activator can stimulate target gene expression in signal transduction by associating with KID interacting domain (KIX). NMR spectra suggest that apo-KID is an unstructured protein. After post-translational modification by phosphorylation, KID undergoes a transition from disordered to well folded protein upon binding to KIX. However, the mechanism of folding coupled to binding is poorly understood. METHODOLOGY:To get an insight into the mechanism, we have performed ten trajectories of explicit-solvent molecular dynamics (MD) for both bound and apo phosphorylated KID (pKID). Ten MD simulations are sufficient to capture the average properties in the protein folding and unfolding. CONCLUSIONS:Room-temperature MD simulations suggest that pKID becomes more rigid and stable upon the KIX-binding. Kinetic analysis of high-temperature MD simulations shows that bound pKID and apo-pKID unfold via a three-state and a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound pKID folds in the order of KIX access, initiation of pKID tertiary folding, folding of helix alpha(B), folding of helix alpha(A), completion of pKID tertiary folding, and finalization of pKID-KIX binding. Our data show that the folding pathways of apo-pKID are different from the bound state: the foldings of helices alpha(A) and alpha(B) are swapped. Here we also show that Asn139, Asp140 and Leu141 with large Phi-values are key residues in the folding of bound pKID. Our results are in good agreement with NMR experimental observations and provide significant insight into the general mechanisms of binding induced protein folding and other conformational adjustment in post-translational modification

Multiple Routes and Milestones in the Folding of HIV–1 Protease Monomer

Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV–1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV–1 protease

Insights into Protein Aggregation by NMR Characterization of Insoluble SH3 Mutants Solubilized in Salt-Free Water

Protein aggregation in vivo has been extensively associated with a large spectrum of human diseases. On the other hand, mechanistic insights into protein aggregation in vitro were incomplete due to the inability in solubilizing insoluble proteins for high-resolution biophysical investigations. However, a new avenue may be opened up by our recent discovery that previously-thought insoluble proteins can in fact be solubilized in salt-free water. Here we use this approach to study the NMR structural and dynamic properties of an insoluble SH3 mutant with a naturally-occurring insertion of Val22 at the tip of the diverging turn. The obtained results reveal: 1) regardless of whether the residue is Val, Ala, Asp or Arg, the insertion will render the first hNck2 SH3 domain to be insoluble in buffers. Nevertheless, all four mutants could be solubilized in salt-free water and appear to be largely unfolded as evident from their CD and NMR HSQC spectra. 2) Comparison of the chemical shift deviations reveals that while in V22-SH3 the second helical region is similarly populated as in the wild-type SH3 at pH 2.0, the first helical region is largely unformed. 3) In V22-SH3, many non-native medium-range NOEs manifest to define non-native helical conformations. In the meanwhile a small group of native-like long-range NOEs still persists, indicating the existence of a rudimentary native-like tertiary topology. 4) Although overall, V22-SH3 has significantly increased backbone motions on the ps-ns time scale, some regions still own restricted backbone motions as revealed by analyzing 15N relaxation data. Our study not only leads to the establishment of the first high-resolution structural and dynamic picture for an insoluble protein, but also shed more light on the molecular events for the nonhierarchical folding mechanism. Furthermore, a general mechanism is also proposed for in vivo protein aggregation triggered by the genetic mutation and posttranslational modification

Correlations of Gene Expression with Blood Lead Levels in Children with Autism Compared to Typically Developing Controls

The objective of this study was to examine the correlation between gene expression and lead (Pb) levels in blood in children with autism (AU, n = 37) compared to typically developing controls (TD, n = 15). We postulated that, though lead levels did not differ between the groups, AU children might metabolize lead differently compared to TD children. RNA was isolated from blood and processed on Affymetrix microarrays. Separate analyses of covariance (ANCOVA) corrected for age and gender were performed for TD, AU, and all subjects (AU + TD). To reduce false positives, only genes that overlapped these three ANCOVAs were considered. Thus, 48 probe sets correlated with lead levels in both AU and TD subjects and were significantly different between the groups (p(Diagnosis × log2 Pb) < 0.05). These genes were related mainly to immune and inflammatory processes, including MHC Class II family members and CD74. A large number (n = 791) of probe sets correlated (P ≤ 0.05) with lead levels in TD but not in AU subjects; and many probe sets (n = 162) correlated (P ≤ 0.05) with lead levels in AU but not in TD subjects. Only 30 probe sets correlated (P ≤ 0.05) with lead levels in a similar manner in the AU and TD groups. These data show that AU and TD children display different associations between transcript levels and low levels of lead. We postulate that this may relate to the underlying genetic differences between the two groups, though other explanations cannot be excluded

Rapid Sampling of Molecular Motions with Prior Information Constraints

Proteins are active, flexible machines that perform a range of different functions. Innovative experimental approaches may now provide limited partial information about conformational changes along motion pathways of proteins. There is therefore a need for computational approaches that can efficiently incorporate prior information into motion prediction schemes. In this paper, we present PathRover, a general setup designed for the integration of prior information into the motion planning algorithm of rapidly exploring random trees (RRT). Each suggested motion pathway comprises a sequence of low-energy clash-free conformations that satisfy an arbitrary number of prior information constraints. These constraints can be derived from experimental data or from expert intuition about the motion. The incorporation of prior information is very straightforward and significantly narrows down the vast search in the typically high-dimensional conformational space, leading to dramatic reduction in running time. To allow the use of state-of-the-art energy functions and conformational sampling, we have integrated this framework into Rosetta, an accurate protocol for diverse types of structural modeling. The suggested framework can serve as an effective complementary tool for molecular dynamics, Normal Mode Analysis, and other prevalent techniques for predicting motion in proteins. We applied our framework to three different model systems. We show that a limited set of experimentally motivated constraints may effectively bias the simulations toward diverse predicates in an outright fashion, from distance constraints to enforcement of loop closure. In particular, our analysis sheds light on mechanisms of protein domain swapping and on the role of different residues in the motion