Nonequilibrium Modeling of the Elementary Step in PDZ3 Allosteric Communication (original) (raw)

Molecular Dynamics-Markov State Model of Protein Ligand Binding and Allostery in CRIB-PDZ: Conformational Selection and Induced Fit

The journal of physical chemistry. B, 2017

Conformational selection and induced fit are well-known contributors to ligand binding and allosteric effects in proteins. Molecular dynamics (MD) simulations now enable the theoretical study of protein-ligand binding in terms of ensembles of interconverting microstates and the population shifts characteristic of "dynamical allostery." Here we investigate protein-ligand binding and allostery based on a Markov state model (MSM) with states and rates obtained from all-atom MD simulations. As an exemplary case, we consider the single domain protein par-6 PDZ with and without ligand and allosteric effector. This is one of the smallest proteins in which allostery has been experimentally observed. In spite of the increased complexity intrinsic to a statistical ensemble perspective, we find that conformational selection and induced fit mechanisms can be readily identified in the analysis. In the nonallosteric pathway, MD-MSM shows that PDZ binds ligand via conformational selectio...

Long-Range Conformational Response of a PDZ Domain to Ligand Binding and Release: A Molecular Dynamics Study

Journal of chemical theory and computation, 2016

The binding of a ligand to a protein may induce long-range structural or dynamical changes in the biomacromolecule even at sites physically well separated from the binding pocket. A system for which such behavior has been widely discussed is the PDZ2 domain of human tyrosine phosphatase 1E. Here, we present results from equilibrium trajectories of the PDZ2 domain in the free and ligand-bound state, as well as nonequilibrium simulations of the relaxation of PDZ2 after removal of its peptide ligand. The study reveals changes in inter-residue contacts, backbone dihedral angles, and Cα positions upon ligand release. Our findings show a long-range conformational response of the PDZ2 domain to ligand release in the form of a collective shift of the secondary structure elements α2, β2, β3, α1-β4, and the C terminal loop relative to the rest of the protein away from the N-terminus, and a shift of the loops β2-β3 and β1-β2 in the opposite direction. The shifts lead to conformational changes ...

Real-time observation of ligand-induced allosteric transitions in a PDZ domain

Proceedings of the National Academy of Sciences, 2020

Significance Allostery is a process in which a signal sensed upon ligand-binding at a distal site is transduced to the effector site, allowing for regulation of the activity of the latter. The propagation of an allosteric signal is a nonequilibrium process, but neither the nature of the signal is known on a molecular level (e.g., whether it is structural or dynamical properties that change) nor its speed. The real-time observation of such a signal requires the design of protein systems, in which one can synchronize ligand (un)binding events. Such a design is presented here, allowing us to investigate its allosteric transition in unprecedented detail.

Long-range conformational transition of a photoswitchable allosteric protein: molecular dynamics simulation study

The journal of physical chemistry. B, 2014

A local perturbation of a protein may lead to functional changes at some distal site. An example is the PDZ2 domain of human tyrosine phosphatase 1E, which shows an allosteric transition upon binding to a peptide ligand. Recently Buchli et al. presented a time-resolved study of this transition by covalently linking an azobenzene photoswitch across the binding groove and using a femtosecond laser pulse that triggers the cis-trans photoisomerization of azobenzene. To aid the interpretation of these experiments, in this work seven microsecond runs of all-atom molecular dynamics simulations each for the wild-type PDZ2 in the ligand-bound and -free state, as well as the photoswitchable protein (PDZ2S) in the cis and trans states of the photoswitch, in explicit water were conducted. First the theoretical model is validated by recalculating the available NMR data from the simulations. By comparing the results for PDZ2 and PDZ2S, it is analyzed to what extent the photoswitch indeed mimics t...

Conserved tertiary couplings stabilize elements in the PDZ fold, leading to characteristic patterns of domain conformational flexibility

Protein Science, 2010

Single-domain allostery has been postulated to occur through intramolecular pathways of signaling within a protein structure. We had previously investigated these pathways by introducing a local thermal perturbation and analyzed the anisotropic propagation of structural changes throughout the protein. Here, we develop an improved approach, the Rotamerically Induced Perturbation (RIP), that identifies strong couplings between residues by analyzing the pathways of heat-flow resulting from thermal excitation of rotameric rotations at individual residues. To explore the nature of these couplings, we calculate the complete coupling maps of 5 different PDZ domains. Although the PDZ domain is a well conserved structural fold that serves as a scaffold in many protein-protein complexes, different PDZ domains display unique patterns of conformational flexibility in response to ligand binding: some show a significant shift in a set of a-helices, while others do not. Analysis of the coupling maps suggests a simple relationship between the computed couplings and observed conformational flexibility. In domains where the a-helices are rigid, we find couplings of the a-helices to the body of the protein, whereas in domains having ligand-responsive a-helices, no couplings are found. This leads to a model where the a-helices are intrinsically dynamic but can be damped if sidechains interact at key tertiary contacts. These tertiary contacts correlate to high covariation contacts as identified by the statistical coupling analysis method. As these dynamic modules are exploited by various allosteric mechanisms, these tertiary contacts have been conserved by evolution.

Thermodynamic Basis for Promiscuity and Selectivity in Protein−Protein Interactions: PDZ Domains, a Case Study

Journal of the American Chemical Society, 2006

Like other protein-protein interaction domains, PDZ domains are involved in many key cellular processes. These processes often require that specific multi-protein complexes be assembled, a task that PDZ domains accomplish by binding to specific peptide motifs in target proteins. However, a growing number of experimental studies show that PDZ domains (like other protein-protein interaction domains) can engage in a variety of interactions and bind distinct peptide motifs. Such promiscuity in ligand recognition raises intriguing questions about the molecular and thermodynamic mechanisms that can sustain it. To identify possible sources of promiscuity and selectivity underlying PDZ domain interactions, we performed molecular dynamics simulations of 20 to 25 ns on a set of 12 different PDZ domain complexes (for the proteins PSD-95, Syntenin, Erbin, GRIP, NHERF, Inad, Dishevelled and Shank). The electrostatic, non-polar and configurational entropy binding contributions were evaluated using the MM/PBSA method combined with a quasi-harmonic analysis. The results revealed that PDZ domain interactions are characterized by overwhelmingly favorable non-polar contributions and almost negligible electrostatic components, a mix that may readily sustain promiscuity. In addition, despite the structural similarity in fold and in recognition modes, the entropic and other dynamical aspects of binding were remarkably variable not only across PDZ domains but also for the same PDZ domain bound to distinct ligands. This variability suggests that entropic and dynamical components can play a role in determining selectivity either of PDZ domain interactions with peptide ligands or of PDZ domains complexes with downstream effectors.

Functional Dynamics of PDZ Binding Domains: A Normal-Mode Analysis

Biophysical Journal, 2005

Postsynaptic density-95/disks large/zonula occludens-1 (PDZ) domains are relatively small (80-120 residues) protein binding modules central in the organization of receptor clusters and in the association of cellular proteins. Their main function is to bind C-terminals of selected proteins that are recognized through specific amino acids in their carboxyl end. Binding is associated with a deformation of the PDZ native structure and is responsible for dynamical changes in regions not in direct contact with the target. We investigate how this deformation is related to the harmonic dynamics of the PDZ structure and show that one low-frequency collective normal mode, characterized by the concerted movements of different secondary structures, is involved in the binding process. Our results suggest that even minimal structural changes are responsible for communication between distant regions of the protein, in agreement with recent NMR experiments. Thus, PDZ domains are a very clear example of how collective normal modes are able to characterize the relation between function and dynamics of proteins, and to provide indications on the precursors of binding/unbinding events.

Diffusion-Limited Unbinding of Small Peptides from PDZ Domains

The Journal of Physical Chemistry B, 2007

PDZ domains are typical examples of binding motifs mediating the formation of protein-protein assemblies in many different cells. A quantitative characterization of the mechanisms intertwining structure, chemistry and dynamics with the PDZ function represent a challenge in molecular biology. Here we investigated the influence of native state topology on the thermodynamics and the dissociation kinetics for a complex PDZ-peptide via Molecular Dynamics simulations based on a coarse-grained description of PDZ domains. Our native-centric approach neglects chemical details but incorporates the basic structural information to reproduce the protein functional dynamics as it couples to the binding. We found that at physiological temperatures the unbinding of a peptide from the PDZ domain becomes increasingly diffusive rather than thermally activated, as a consequence of the significant reduction of the free energy barrier with temperature. In turn, this results in a significant slowing down of the process of two orders of magnitude with respect to the conventional Arrhenius extrapolation from low temperature calculations. Finally, a detailed analysis of a typical unbinding event based on the rupture times of single peptide-PDZ contacts allows to shed further light on the dissociation mechanism and to elaborate a coherent picture of the relation between function and dynamics in PDZ domains.

The Conformational Plasticity Vista of PDZ Domains

Life

The PDZ domain (PSD95-Discs large-ZO1) is a widespread modular domain present in the living organisms. A prevalent function in the PDZ family is to serve as scaffolding and adaptor proteins connecting multiple partners in signaling pathways. An explanation of the flexible functionality in this domain family, based just on a static perspective of the structure–activity relationship, might fall short. More dynamic and conformational aspects in the protein fold can be the reasons for such functionality. Folding studies indeed showed an ample and malleable folding landscape for PDZ domains where multiple intermediate states were experimentally detected. Allosteric phenomena that resemble energetic coupling between residues have also been found in PDZ domains. Additionally, several PDZ domains are modulated by post-translational modifications, which introduce conformational switches that affect binding. Altogether, the ability to connect diverse partners might arise from the intrinsic pl...