Two polymorphisms facilitate differences in plasticity between two chicken major histocompatibility complex class I proteins

Major histocompatibility complex class I molecules (MHC I) present peptides to cytotoxic T-cells at the surface of almost all nucleated cells. The function of MHC I molecules is to select high affinity peptides from a large intracellular pool and they are assisted in this process by co-factor molecules, notably tapasin. In contrast to mammals, MHC homozygous chickens express a single MHC I gene locus, termed BF2, which is hypothesised to have co-evolved with the highly polymorphic tapasin within stable haplotypes. The BF2 molecules of the B15 and B19 haplotypes have recently been shown to differ in their interactions with tapasin and in their peptide selection properties. This study investigated whether these observations might be explained by differences in the protein plasticity that is encoded into the MHC I structure by primary sequence polymorphisms. Furthermore, we aimed to demonstrate the utility of a complimentary modelling approach to the understanding of complex experimental data. Combining mechanistic molecular dynamics simulations and the primary sequence based technique of statistical coupling analysis, we show how two of the eight polymorphisms between BF2*15:01 and BF2*19:01 facilitate differences in plasticity. We show that BF2*15:01 is intrinsically more plastic than BF2*19:01, exploring more conformations in the absence of peptide. We identify a protein sector of contiguous residues connecting the membrane bound ?3 domain and the heavy chain peptide binding site. This sector contains two of the eight polymorphic residues. One is residue 22 in the peptide binding domain and the other 220 is in the ?3 domain, a putative tapasin binding site. These observations are in correspondence with the experimentally observed functional differences of these molecules and suggest a mechanism for how modulation of MHC I plasticity by tapasin catalyses peptide selection allosterically

Temperature Dependent Conformational Transitions and Hydrogen Bond Dynamics of the Elastin-Like Octapeptide GVG(VPGVG): a Molecular Dynamics Study

A joint experimental / theoretical investigation of the elastin-like octapeptide GVG(VPGVG) was carried out. In this paper a comprehensive molecular dynamics study of the temperature dependent folding and unfolding of the octapeptide is presented. The current study, as well as its experimental counterpart find that this peptide undergoes an "inverse temperature transition", ITT, leading to a folding at about 310-330 K. In addition, an unfolding transition is identified at unusually high temperatures approaching the boiling point of water. Due to the small size of the system two broad temperature regimes are found: the "ITT regime" (at about 280-320 K) and the "unfolding regime" at about T > 330 K, where the peptide has a maximum probability of being folded at approximately 330 K. A detailed molecular picture involving a thermodynamic order parameter, or reaction coordinate, for this process is presented along with a time-correlation function analysis of the hydrogen bond dynamics within the peptide as well as between the peptide and solvating water molecules. Correlation with experimental evidence and ramifications on the properties of elastin are discussed.Comment: 15 pages, 1 table, 8 figure

Picosecond fluctuating protein energy landscape mapped by pressure–temperature molecular dynamics simulation

Microscopic statistical pressure fluctuations can, in principle, lead to corresponding fluctuations in the shape of a protein energy landscape. To examine this, nanosecond molecular dynamics simulations of lysozyme are performed covering a range of temperatures and pressures. The well known dynamical transition with temperature is found to be pressure-independent, indicating that the effective energy barriers separating conformational substates are not significantly influenced by pressure. In contrast, vibrations within substates stiffen with pressure, due to increased curvature of the local harmonic potential in which the atoms vibrate. The application of pressure is also shown to selectively increase the damping of the anharmonic, low-frequency collective modes in the protein, leaving the more local modes relatively unaffected. The critical damping frequency, i.e., the frequency at which energy is most efficiently dissipated, increases linearly with pressure. The results suggest that an invariant description of protein energy landscapes should be subsumed by a fluctuating picture and that this may have repercussions in, for example, mechanisms of energy dissipation accompanying functional, structural, and chemical relaxation

Anharmonicity and self-similarity of the free energy landscape of protein G

The near-native free energy landscape of protein G is investigated through 0.4 microseconds-long atomistic molecular dynamics simulations in explicit solvent. A theoretical and computational framework is used to assess the time-dependence of salient thermodynamical features. While the quasi-harmonic character of the free energy is found to degrade in a few ns, the slow modes display a very mild dependence on the trajectory duration. This property originates from a striking self-similarity of the free energy landscape embodied by the consistency of the principal directions of the local minima, where the system dwells for several ns, and of the virtual jumps connecting them.Comment: revtex, 6 pages, 5 figure

Reconstructing the free-energy landscape of Met-enkephalin using dihedral Principal Component Analysis and Well-tempered Metadynamics

Well-Tempered Metadynamics (WTmetaD) is an efficient method to enhance the reconstruction of the free-energy surface of proteins. WTmetaD guarantees a faster convergence in the long time limit in comparison with the standard metadynamics. It still suffers however from the same limitation, i.e. the non trivial choice of pertinent collective variables (CVs). To circumvent this problem, we couple WTmetaD with a set of CVs generated from a dihedral Principal Component Analysis (dPCA) on the Ramachadran dihedral angles describing the backbone structure of the protein. The dPCA provides a generic method to extract relevant CVs built from internal coordinates. We illustrate the robustness of this method in the case of the small and very diffusive Metenkephalin pentapeptide, and highlight a criterion to limit the number of CVs necessary to biased the metadynamics simulation. The free-energy landscape (FEL) of Met-enkephalin built on CVs generated from dPCA is found rugged compared with the FEL built on CVs extracted from PCA of the Cartesian coordinates of the atoms.Comment: 17 pages, 9 figures (4 in color

Eigenvector Centrality Distribution for Characterization of Protein Allosteric Pathways

Determining the principal energy pathways for allosteric communication in biomolecules, that occur as a result of thermal motion, remains challenging due to the intrinsic complexity of the systems involved. Graph theory provides an approach for making sense of such complexity, where allosteric proteins can be represented as networks of amino acids. In this work, we establish the eigenvector centrality metric in terms of the mutual information, as a mean of elucidating the allosteric mechanism that regulates the enzymatic activity of proteins. Moreover, we propose a strategy to characterize the range of the physical interactions that underlie the allosteric process. In particular, the well known enzyme, imidazol glycerol phosphate synthase (IGPS), is utilized to test the proposed methodology. The eigenvector centrality measurement successfully describes the allosteric pathways of IGPS, and allows to pinpoint key amino acids in terms of their relevance in the momentum transfer process. The resulting insight can be utilized for refining the control of IGPS activity, widening the scope for its engineering. Furthermore, we propose a new centrality metric quantifying the relevance of the surroundings of each residue. In addition, the proposed technique is validated against experimental solution NMR measurements yielding fully consistent results. Overall, the methodologies proposed in the present work constitute a powerful and cost effective strategy to gain insight on the allosteric mechanism of proteins

On the conservation of the slow conformational dynamics within the amino acid kinase family: NAGK the paradigm

N-Acetyl-L-Glutamate Kinase (NAGK) is the structural paradigm for examining the catalytic mechanisms and dynamics of amino acid kinase family members. Given that the slow conformational dynamics of the NAGK (at the microseconds time scale or slower) may be rate-limiting, it is of importance to assess the mechanisms of the most cooperative modes of motion intrinsically accessible to this enzyme. Here, we present the results from normal mode analysis using an elastic network model representation, which shows that the conformational mechanisms for substrate binding by NAGK strongly correlate with the intrinsic dynamics of the enzyme in the unbound form. We further analyzed the potential mechanisms of allosteric signalling within NAGK using a Markov model for network communication. Comparative analysis of the dynamics of family members strongly suggests that the low-frequency modes of motion and the associated intramolecular couplings that establish signal transduction are highly conserved among family members, in support of the paradigm sequence→structure→dynamics→function © 2010 Marcos et al