Time-resolved observation of protein allosteric communication

Identifying Conformational-Selection and Induced-Fit Aspects in the Binding-Induced Folding of PMI from Markov State Modeling of Atomistic Simulations

Unstructured proteins and peptides typically fold during binding to ligand proteins. A challenging problem is to identify the mechanism and kinetics of these binding-induced folding processes in experiments and atomistic simulations. In this Article, we present a detailed picture for the folding of the inhibitor peptide PMI into a helix during binding to the oncoprotein fragment ^25-109Mdm2 obtained from atomistic, explicit-water simulations and Markov state modeling. We find that binding-induced folding of PMI is highly parallel and can occur along a multitude of pathways. Some pathways are induced-fit-like with binding occurring prior to PMI helix formation, while other pathways are conformational-selection-like with binding after helix formation. On the majority of pathways, however, binding is intricately coupled to folding, without clear temporal ordering. A central feature of these pathways is PMI motion on the Mdm2 surface, along the binding groove of Mdm2 or over the rim of this groove. The native binding groove of Mdm2 thus appears as an asymmetric funnel for PMI binding. Overall, binding-induced folding of PMI does not fit into the classical picture of induced fit or conformational selection that implies a clear temporal ordering of binding and folding events. We argue that this holds in general for binding-induced folding processes because binding and folding events in these processes likely occur on similar time scales and do exhibit the time-scale separation required for temporal ordering.

Multiple Pathways and Time Scales for Conformational Transitions in apo-Adenylate Kinase

The open/close transition in adenylate kinase (AK) is regarded as a representative example for large-scale conformational transition in proteins, yet its mechanism remains unclear despite numerous experimental and computational studies. Using extensive (∼50 μs) explicit solvent atomistic simulations and Markov state analysis, we shed new lights on the mechanism of this transition in the apo form of AK. The closed basin of apo AK features an open NMP domain while the LID domain closes and rotates toward it. Therefore, although the computed structural properties of the closed ensemble are consistent with previously reported FRET and PRE measurements, our simulations suggest that NMP closure is likely to follow AMP binding, in contrast to the previous interpretation of FRET and PRE data that the apo state was able to sample the fully closed conformation for "ligand selection". The closed state ensemble is found to be kinetically heterogeneous; multiple pathways and time scales are associated with the open/close transition, providing new clues to the disparate time scales observed in different experiments. Besides interdomain interactions, a novel mutual information analysis identifies specific intradomain interactions that correlate strongly to transition kinetics, supporting observations from previous chimera experiments. While our results underscore the role of internal domain properties in determining the kinetics of open/close transition in apo AK, no evidence is observed for any significant degree of local unfolding during the transition. These observations about AK have general implications to our view of conformational states, transition pathways, and time scales of conformational changes in proteins. The key features and time scales of observed transition pathways are robust and similar from simulations using two popular fixed charge force fields.

Hinge-Shift Mechanism Modulates Allosteric Regulations in Human Pin1

Allostery, which is regulation from distant sites, plays a major role in biology. While traditional allostery is described in terms of conformational change upon ligand binding as an underlying principle, it is possible to have allosteric regulations without significant conformational change through modulating the conformational dynamics by altering the local effective elastic modulus of the protein upon ligand binding. Pin1 utilizes this dynamic allostery to regulate its function. It is a modular protein containing a WW domain and a larger peptidyl prolyl isomerase domain (PPIase) that isomerizes phosphoserine/threonine-proline (pS/TP) motifs. The WW domain serves as a docking module, whereas catalysis solely takes place within the PPIase domain. Here, we analyze the change in the dynamic flexibility profile of the PPIase domain upon ligand binding to the WW domain. Substrate binding to the WW domain induces the formation of a new rigid hinge site around the interface of the two domains and loosens the flexibility of a rigid site existing in the Apo form around the catalytic site. This hinge-shift mechanism enhances the dynamic coupling of the catalytic positions with the PPIase domain, where the rest of the domain can cooperatively respond to the local conformational changes around the catalytic site, leading to an increase in catalytic efficiency.

Metadynamics Enhanced Markov Modeling of Protein Dynamics

Enhanced sampling techniques represent a versatile approach to account for rare conformational transitions in biomolecules. A particularly promising strategy is to combine massive parallel computing of short molecular dynamics (MD) trajectories (to sample the free energy landscape of the system) with Markov state modeling (to rebuild the kinetics from the sampled data). To obtain well-distributed initial structures for the short trajectories, it is proposed to employ metadynamics MD, which quickly sweeps through the entire free energy landscape of interest. Being only used to generate initial conformations, the implementation of metadynamics can be simple and fast. The conformational dynamics of helical peptide Aib₉ is adopted to discuss various technical issues of the approach, including metadynamics settings, minimal number and length of short MD trajectories, and the validation of the resulting Markov models. Using metadynamics to launch some thousands of nanosecond trajectories, several Markov state models are constructed that reveal that previous unbiased MD simulations of in total 16 μs length cannot provide correct equilibrium populations or qualitative features of the pathway distribution of the short peptide.

High-Definition Mapping of Four Spatially Distinct Neutralizing Epitope Clusters on RiVax, a Candidate Ricin Toxin Subunit Vaccine

RiVax is a promising recombinant ricin toxin A subunit (RTA) vaccine antigen that has been shown to be safe and immunogenic in humans and effective at protecting rhesus macaques against lethal-dose aerosolized toxin exposure. We previously used a panel of RTA-specific monoclonal antibodies (MAbs) to demonstrate, by competition enzyme-linked immunosorbent assay (ELISA), that RiVax elicits similar serum antibody profiles in humans and macaques. However, the MAb binding sites on RiVax have yet to be defined. In this study, we employed hydrogen exchange-mass spectrometry (HX-MS) to localize the epitopes on RiVax recognized by nine toxin-neutralizing MAbs and one nonneutralizing MAb. Based on strong protection from hydrogen exchange, the nine MAbs grouped into four spatially distinct epitope clusters (namely, clusters I to IV). Cluster I MAbs protected RiVax's α-helix B (residues 94 to 107), a protruding immunodominant secondary structure element known to be a target of potent toxin-neutralizing antibodies. Cluster II consisted of two subclusters located on the "back side" (relative to the active site pocket) of RiVax. One subcluster involved α-helix A (residues 14 to 24) and α-helices F-G (residues 184 to 207); the other encompassed β-strand d (residues 62 to 69) and parts of α-helices D-E (154 to 164) and the intervening loop. Cluster III involved α-helices C and G on the front side of RiVax, while cluster IV formed a sash from the front to back of RiVax, spanning strands b, c, and d (residues 35 to 59). Having a high-resolution B cell epitope map of RiVax will enable the development and optimization of competitive serum profiling assays to examine vaccine-induced antibody responses across species.

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.