Postsynaptic protein assembly in three and two dimensions studied by mesoscopic simulations.

Recently, cellular biomolecular condensates formed via phase separation have received considerable attention. While they can be formed either in cytosol (denoted as 3D) or beneath the membrane (2D), the underlying difference between the two has not been well clarified. To compare the phase behaviors in 3D and 2D, postsynaptic density (PSD) serves as a model system. PSD is a protein condensate located under the postsynaptic membrane that influences the localization of glutamate receptors and thus contributes to synaptic plasticity. Recent in vitro studies have revealed the formation of droplets of various soluble PSD proteins via liquid-liquid phase separation. However, it is unclear how these protein condensates are formed beneath the membrane and how they specifically affect the localization of glutamate receptors in the membrane. In this study, focusing on the mixture of a glutamate receptor complex, AMPAR-TARP, and a ubiquitous scaffolding protein, PSD-95, we constructed a mesoscopic model of protein-domain interactions in PSD and performed comparative molecular simulations. The results showed a sharp contrast in the phase behaviors of protein assemblies in 3D and those under the membrane (2D). A mixture of a soluble variant of the AMPAR-TARP complex and PSD-95 in the 3D system resulted in a phase-separated condensate, which was consistent with the experimental results. However, with identical domain interactions, AMPAR-TARP embedded in the membrane formed clusters with PSD-95, but did not form a stable separated phase. Thus, the cluster formation behaviors of PSD proteins in the 3D and 2D systems were distinct. The current study suggests that, more generally, stable phase separation can be more difficult to achieve in and beneath the membrane than in 3D systems.

Fuzzy supertertiary interactions within PSD-95 enable ligand binding.

The scaffold protein PSD-95 links postsynaptic receptors to sites of presynaptic neurotransmitter release. Flexible linkers between folded domains in PSD-95 enable a dynamic supertertiary structure. Interdomain interactions within the PSG supramodule, formed by PDZ3, SH3, and Guanylate Kinase domains, regulate PSD-95 activity. Here we combined discrete molecular dynamics and single molecule Förster resonance energy transfer (FRET) to characterize the PSG supramodule, with time resolution spanning picoseconds to seconds. We used a FRET network to measure distances in full-length PSD-95 and model the conformational ensemble. We found that PDZ3 samples two conformational basins, which we confirmed with disulfide mapping. To understand effects on activity, we measured binding of the synaptic adhesion protein neuroligin. We found that PSD-95 bound neuroligin well at physiological pH while truncated PDZ3 bound poorly. Our hybrid structural models reveal how the supertertiary context of PDZ3 enables recognition of this critical synaptic ligand.

Human Serine Racemase Weakly Binds the Third PDZ Domain of PSD-95.

Human serine racemase (hSR) is a pyridoxal-5'-phosphate (PLP)-dependent dimer that catalyzes the formation of D-serine from L-serine, as well as the dehydration of both L- and D-serine to pyruvate and ammonia. As D-serine is a co-agonist of N-methyl-D-aspartate receptors (NMDARs), hSR is a key enzyme in glutamatergic neurotransmission. hSR activity is finely regulated by Mg2+, ATP, post-translational modifications, and the interaction with protein partners. In particular, the C-terminus of murine SR binds the third PDZ domain (PDZ3) of postsynaptic density protein 95 (PSD-95), a member of the membrane-associated guanylate kinase (MAGUK) family involved in the trafficking and localization of glutamate receptors. The structural details of the interaction and the stability of the complex have not been elucidated yet. We evaluated the binding of recombinant human PSD-95 PDZ3 to hSR by glutaraldehyde cross-linking, pull-down assays, isothermal titration calorimetry, nuclear magnetic resonance, and enzymatic assays. Overall, a weak interaction was observed, confirming the binding for the human orthologs but supporting the hypothesis that a third protein partner (i.e., stargazin) is required for the regulation of hSR activity by PSD-95 and to stabilize their interaction.

Zinc-chelating postsynaptic density-95 N-terminus impairs its palmitoyl modification.

Chemical synaptic transmission represents the most sophisticated dynamic process and is highly regulated with optimized neurotransmitter balance. Imbalanced transmitters can lead to transmission impairments, for example, intracellular zinc accumulation is a hallmark of degenerating neurons. However, the underlying mechanisms remain elusive. Postsynaptic density protein-95 (PSD-95) is a primary postsynaptic membrane-associated protein and the major scaffolding component in the excitatory postsynaptic densities, which performs substantial functions in synaptic development and maturation. Its membrane association induced by palmitoylation contributes largely to its regulatory functions at postsynaptic sites. Unlike other structural domains in PSD-95, the N-terminal region (PSD-95NT) is flexible and interacts with various targets, which modulates its palmitoylation of two cysteines (C3/C5) and glutamate receptor distributions in postsynaptic densities. PSD-95NT contains a putative zinc-binding motif (C2H2) with undiscovered functions. This study is the first effort to investigate the interaction between Zn2+ and PSD-95NT. The NMR titration of 15 N-labeled PSD-95NT by ZnCl2 was performed and demonstrated Zn2+ binds to PSD-95NT with a binding affinity (Kd ) in the micromolar range. The zinc binding was confirmed by fluorescence and mutagenesis assays, indicating two cysteines and two histidines (H24, H28) are critical residues for the binding. These results suggested the concentration-dependent zinc binding is likely to influence PSD-95 palmitoylation since the binding site overlaps the palmitoylation sites, which was verified by the mimic PSD-95 palmitoyl modification and intact cell palmitoylation assays. This study reveals zinc as a novel modulator for PSD-95 postsynaptic membrane association by chelating its N-terminal region, indicative of its importance in postsynaptic signaling.

Dissecting Inter-domain Cooperativity in the Folding of a Multi Domain Protein.

Correct protein folding underlies all cellular functions. While there are detailed descriptions and a good understanding of protein folding pathways for single globular domains there is a paucity of quantitative data regarding folding of multidomain proteins. We have here investigated the folding of a three-domain supramodule from the protein PSD-95, consisting of one PDZ domain, one SH3 domain and one guanylate kinase-like (GK) domain. This supramodule has previously been shown to work as one functional unit with regard to ligand binding. We used equilibrium and kinetic folding experiments to demonstrate that the PDZ domain folds faster and independently from the SH3-GK tandem, which folds as one cooperative unit. However, concurrent folding of the PDZ domain slows down folding of SH3-GK by non-native interactions, resulting in an off-pathway folding intermediate. Our data contribute to an emerging description of multidomain protein folding in which individual domains cannot a priori be viewed as separate folding units.

A novel cell-penetrating peptide targeting calpain-cleavage of PSD-95 induced by excitotoxicity improves neurological outcome after stroke.

Postsynaptic density protein-95 (PSD-95) is a multidomain protein critical to the assembly of signaling complexes at excitatory synapses, required for neuronal survival and function. However, calpain-processing challenges PSD-95 function after overactivation of excitatory glutamate receptors (excitotoxicity) in stroke, a leading cause of death, disability and dementia in need of efficient pharmacological treatments. A promising strategy is neuroprotection of the infarct penumbra, a potentially recoverable area, by promotion of survival signaling. Interference of PSD-95 processing induced by excitotoxicity might thus be a therapeutic target for stroke and other excitotoxicity-associated pathologies. Methods: The nature and stability of PSD-95 calpain-fragments was analyzed using in vitro assays or excitotoxic conditions induced in rat primary neuronal cultures or a mouse model of stroke. We then sequenced PSD-95 cleavage-sites and rationally designed three cell-penetrating peptides (CPPs) containing these sequences. The peptides effects on PSD-95 stability and neuronal viability were investigated in the cultured neurons, subjected to acute or chronic excitotoxicity. We also analyzed the effect of one of these peptides in the mouse model of stroke by measuring infarct size and evaluating motor coordination and balance. Results: Calpain cleaves three interdomain linker regions in PSD-95 and produces stable fragments corresponding to previously described PSD-95 supramodules (PDZ1-2 and P-S-G) as well as a truncated form SH3-GK. Peptide TP95414, containing the cleavage site in the PDZ3-SH3 linker, is able to interfere PSD-95 downregulation and reduces neuronal death by excitotoxicity. Additionally, TP95414 is delivered to mice cortex and, in a severe model of permanent ischemia, significantly improves the neurological outcome after brain damage. Conclusions: Interference of excitotoxicity-induced PSD-95-processing with specific CPPs constitutes a novel and promising therapeutic approach for stroke treatment.

The Speed of Allosteric Signaling Within a Single-Domain Protein.

While much is known about different allosteric regulation mechanisms, the nature of the allosteric signal and the time scale on which it propagates remains elusive. The PDZ3 domain from postsynaptic density-95 protein is a small protein domain with a terminal third α-helix, i.e., the α3-helix, which is known to be allosterically active. By cross-linking the allosteric helix with an azobenzene moiety, we obtained a photocontrollable PDZ3 variant. Photoswitching triggers its allosteric transition, resulting in a change in binding affinity of a peptide to the remote binding pocket. Using time-resolved infrared and UV/vis spectroscopy, we follow the allosteric signal transduction and reconstruct the timeline in which the allosteric signal propagates through the protein within 200 ns.

Sensing the allosteric force.

Allosteric regulation is an innate control in most metabolic and signalling cascades that enables living organisms to adapt to the changing environment by tuning the affinity and regulating the activity of target proteins. For a microscopic understanding of this process, a protein system has been designed in such a way that allosteric communication between the binding and allosteric site can be observed in both directions. To that end, an azobenzene-derived photoswitch has been linked to the α3-helix of the PDZ3 domain, arguably the smallest allosteric protein with a clearly identifiable binding and allosteric site. Photo-induced trans-to-cis isomerisation of the photoswitch increases the binding affinity of a small peptide ligand to the protein up to 120-fold, depending on temperature. At the same time, ligand binding speeds up the thermal cis-to-trans back-isomerisation rate of the photoswitch. Based on the energetics of the four states of the system (cis vs trans and ligand-bound vs free), the concept of an allosteric force is introduced, which can be used to drive chemical reactions.

Pressure Sensitivity of SynGAP/PSD-95 Condensates as a Model for Postsynaptic Densities and Its Biophysical and Neurological Ramifications.

Biomolecular condensates consisting of proteins and nucleic acids can serve critical biological functions, so that some condensates are referred as membraneless organelles. They can also be disease-causing, if their assembly is misregulated. A major physicochemical basis of the formation of biomolecular condensates is liquid-liquid phase separation (LLPS). In general, LLPS depends on environmental variables, such as temperature and hydrostatic pressure. The effects of pressure on the LLPS of a binary SynGAP/PSD-95 protein system mimicking postsynaptic densities, which are protein assemblies underneath the plasma membrane of excitatory synapses, were investigated. Quite unexpectedly, the model system LLPS is much more sensitive to pressure than the folded states of typical globular proteins. Phase-separated droplets of SynGAP/PSD-95 were found to dissolve into a homogeneous solution already at ten-to-hundred bar levels. The pressure sensitivity of SynGAP/PSD-95 is seen here as a consequence of both pressure-dependent multivalent interaction strength and void volume effects. Considering that organisms in the deep sea are under pressures up to about 1 kbar, this implies that deep-sea organisms have to devise means to counteract this high pressure sensitivity of biomolecular condensates to avoid harm. Intriguingly, these findings may shed light on the biophysical underpinning of pressure-related neurological disorders in terrestrial vertebrates.

Ligand-dependent intra- and interdomain motions in the PDZ12 tandem regulate binding interfaces in postsynaptic density protein-95.

The postsynaptic density protein-95 (PSD-95) regulates synaptic plasticity through interactions mediated by its peptide-binding PDZ domains. The two N-terminal PDZ domains of PSD-95 form an autonomous structural unit, and their interdomain orientation and dynamics depend on ligand binding. To understand the mechanistic details of the effect of ligand binding, we generated conformational ensembles using available experimentally determined nuclear Overhauser effect interatomic distances and S2 order parameters. In our approach, the fast dynamics of the two domains is treated independently. We find that intradomain structural changes induced by ligand binding modulate the probability of the occurrence of specific domain-domain orientations. Our results suggest that the ß2-ß3 loop in the PDZ domains is a key regulatory region, which influences both intradomain motions and supramodular rearrangement.

Functional interplay between protein domains in a supramodular structure involving the postsynaptic density protein PSD-95.

Cell scaffolding and signaling are governed by protein-protein interactions. Although a particular interaction is often defined by two specific domains binding to each other, this interaction often occurs in the context of other domains in multidomain proteins. How such adjacent domains form supertertiary structures and modulate protein-protein interactions has only recently been addressed and is incompletely understood. The postsynaptic density protein PSD-95 contains a three-domain supramodule, denoted PSG, which consists of PDZ, Src homology 3 (SH3), and guanylate kinase-like domains. The PDZ domain binds to the C terminus of its proposed natural ligand, CXXC repeat-containing interactor of PDZ3 domain (CRIPT), and results from previous experiments using only the isolated PDZ domain are consistent with the simplest scenario for a protein-protein interaction; namely, a two-state mechanism. Here we analyzed the binding kinetics of the PSG supramodule with CRIPT. We show that PSG binds CRIPT via a more complex mechanism involving two conformational states interconverting on the second timescale. Both conformational states bound a CRIPT peptide with similar affinities but with different rates, and the distribution of the two conformational states was slightly shifted upon CRIPT binding. Our results are consistent with recent structural findings of conformational changes in PSD-95 and demonstrate how conformational transitions in supertertiary structures can shape the ligand-binding energy landscape and modulate protein-protein interactions.

PtdIns4P-mediated electrostatic forces influence S-acylation of peripheral proteins at the Golgi complex.

Protein S-acylation is a reversible post-translational modification involving the addition of fatty acids to cysteines and is catalyzed by transmembrane protein acyltransferases (PATs) mainly expressed at the Golgi complex. In case of soluble proteins, S-acylation confers stable membrane attachment. Myristoylation or farnesylation of many soluble proteins constitutes the initial transient membrane adsorption step prior to S-acylation. However, some S-acylated soluble proteins, such as the neuronal growth-associated protein Growth-associated protein-43 (GAP-43), lack the hydrophobic modifications required for this initial membrane interaction. The signals for GAP-43 S-acylation are confined to the first 13 amino acids, including the S-acylatable cysteines 3 and 4 embedded in a hydrophobic region, followed by a cluster of basic amino acids. We found that mutation of critical basic amino acids drastically reduced membrane interaction and hence S-acylation of GAP-43. Interestingly, acute depletion of phosphatidylinositol 4-phosphate (PtdIns4P) at the Golgi complex reduced GAP-43 membrane binding, highlighting a new, pivotal role for this anionic lipid and supporting the idea that basic amino acid residues are involved in the electrostatic interactions between GAP-43 and membranes of the Golgi complex where they are S-acylated.

Disruption of Striatal-Enriched Protein Tyrosine Phosphatase Signaling Might Contribute to Memory Impairment in a Mouse Model of Sepsis-Associated Encephalopathy.

Sepsis-associated encephalopathy (SAE) is a potentially irreversible acute cognitive dysfunction with unclear mechanism. Striatal-enriched protein tyrosine phosphatase (STEP) is a brain-specific phosphatase which normally opposes synaptic strengthening by regulating key signaling molecules involved in synaptic plasticity and neuronal function. Thus, we hypothesized that abnormal STEP signaling pathway was involved in sepsis-induced cognitive impairment evoked by lipopolysaccharides (LPS) injection. The levels of STEP, phosphorylation of GluN2B (pGluN2B), the kinases extracellular signal-regulated kinase 1/2 (pERK), cAMP-response element binding protein (CREB), synaptophysin, brain derived neurotrophic factor (BDNF), and post-synaptic density protein 95 (PSD95) in the hippocampus, prefrontal cortex, and striatum were determined at the indicated time points. In the present study, we found that STEP levels were significantly increased in the hippocampus, prefrontal cortex, and striatum following LPS injection, which might resulted from the disruption of the ubiquitin-proteasome system. Notably, a STEP inhibitor TC-2153 treatment alleviated sepsis-induced memory impairment by increasing phosphorylation of GluN2B and ERK1/2, CREB/BDNF, and PSD95. In summary, our results support the key role of STEP in sepsis-induced memory impairment in a mouse model of SAE, whereas inhibition of STEP may provide a novel therapeutic approach for this disorder and possible other neurodegenerative diseases.

Postsynaptic density protein 95 (PSD-95) is transported by KIF5 to dendritic regions.

Postsynaptic density protein 95 (PSD-95) is a pivotal postsynaptic scaffolding protein in excitatory neurons. Although the transport and regulation of PSD-95 in synaptic regions is well understood, dendritic transport of PSD-95 before synaptic localization still remains to be clarified. To evaluate the role of KIF5, conventional kinesin, in the dendritic transport of PSD-95 protein, we expressed a transport defective form of KIF5A (ΔMD) that does not contain the N-terminal motor domain. Expression of ΔMD significantly decreased PSD-95 level in the dendrites. Consistently, KIF5 was associated with PSD-95 in in vitro and in vivo assays. This interaction was mediated by the C-terminal tail regions of KIF5A and the third PDZ domain of PSD-95. Additionally, the ADPDZ3 (the association domain of NMDA receptor and PDZ3 domain) expression significantly reduced the levels of PSD-95, glutamate receptor 1 (GluA1) in dendrites. The association between PSD-95 and KIF5A was dose-dependent on Staufen protein, suggesting that the Staufen plays a role as a regulatory role in the association. Taken together, our data suggest a new mechanism for dendritic transport of the AMPA receptor-PSD-95.

Probing Protein Allostery as a Residue-Specific Concept via Residue Response Maps.

Allosteric regulation is a well-established phenomenon defined as a distal conformational or dynamical change of the protein upon allosteric effector binding. Here, we developed a novel approach to delineate allosteric effects in proteins. In this approach, we applied robust machine learning methods, including deep neural network and random forest, on extensive molecular dynamics (MD) simulations to distinguish otherwise similar allosteric states of proteins. Using the PDZ3 domain of PDS-95 as a model protein, we demonstrated that the allosteric effects could be represented as residue-specific properties through two-dimensional property-residue maps, which we refer to as "residue response maps". These maps were constructed through two machine learning methods and could accurately describe how different properties of various residues are affected upon allosteric perturbation on protein. Based on the "residue response maps", we propose allostery as a residue-specific concept, suggesting that all residues could be considered as allosteric residues because each residue "senses" the allosteric events through changing its single or multiple attributes in a quantitatively unique way. The "residue response maps" could be used to fingerprint a protein based on the unique patterns of residue responses upon binding events, providing a novel way to systematically describe the protein allosteric effects of each residue upon perturbation.

Structural determinants governing ß-arrestin2 interaction with PDZ proteins and recruitment to CRFR1.

ß-Arrestins are multifunctional adaptor proteins best know for their vital role in regulating G protein coupled receptor (GPCR) trafficking and signaling. ß-arrestin2 recruitment and receptor internalization of corticotropin-releasing factor receptor 1 (CRFR1), a GPCR whose antagonists have been shown to demonstrate both anxiolytic- and antidepressant-like effects, have previously been shown to be modulated by PDZ proteins. Thus, a structural characterization of the interaction between ß-arrestins and PDZ proteins can delineate potential mechanism of PDZ-dependent regulation of GPCR trafficking. Here, we find that the PDZ proteins PSD-95, MAGI1, and PDZK1 interact with ß-arrestin2 in a PDZ domain-dependent manner. Further investigation of such interaction using mutational analyses revealed that mutating the alanine residue at 175 residue of ß-arrestin2 to phenylalanine impairs interaction with PSD-95. Additionally, A175F mutant of ß-arrestin2 shows decreased CRF-stimulated recruitment to CRFR1 and reduced receptor internalization. Thus, our findings show that the interaction between ß-arrestins and PDZ proteins is key for CRFR1 trafficking and may be targeted to mitigate impaired CRFR1 signaling in mental and psychiatric disorders.

Conformational changes in the third PDZ domain of the neuronal postsynaptic density protein 95.

PDZ domains are protein-protein recognition modules that interact with other proteins through short sequences at the carboxyl terminus. These domains are structurally characterized by a conserved fold composed of six ß-strands and two α-helices. The third PDZ domain of the neuronal postsynaptic density protein 95 has an additional α-helix (α3), the role of which is not well known. In previous structures, a succinimide was identified in the ß2-ß3 loop instead of Asp332. The presence of this modified residue results in conformational changes in α3. In this work, crystallographic structures of the following have been solved: a truncated form of the third PDZ domain of the neuronal postsynaptic density protein 95 from which α3 has been removed, D332P and D332G variants of the protein, and a new crystal form of this domain showing the binding of Asp332 to the carboxylate-binding site of a symmetry-related molecule. Crystals of the wild type and variants were obtained in different space groups, which reflects the conformational plasticity of the domain. Indeed, the overall analysis of these structures suggests that the conformation of the ß2-ß3 loop is correlated with the fold acquired by α3. The alternate conformation of the ß2-ß3 loop affects the electrostatics of the carboxylate-binding site and might modulate the binding of different PDZ-binding motifs.

Intramolecular domain dynamics regulate synaptic MAGUK protein interactions.

PSD-95 MAGUK family scaffold proteins are multi-domain organisers of synaptic transmission that contain three PDZ domains followed by an SH3-GK domain tandem. This domain architecture allows coordinated assembly of protein complexes composed of neurotransmitter receptors, synaptic adhesion molecules and downstream signalling effectors. Here we show that binding of monomeric CRIPT-derived PDZ3 ligands to the third PDZ domain of PSD-95 induces functional changes in the intramolecular SH3-GK domain assembly that influence subsequent homotypic and heterotypic complex formation. We identify PSD-95 interactors that differentially bind to the SH3-GK domain tandem depending on its conformational state. Among these interactors, we further establish the heterotrimeric G protein subunit Gnb5 as a PSD-95 complex partner at dendritic spines of rat hippocampal neurons. The PSD-95 GK domain binds to Gnb5, and this interaction is triggered by CRIPT-derived PDZ3 ligands binding to the third PDZ domain of PSD-95, unraveling a hierarchical binding mechanism of PSD-95 complex formation.

Molecular Dynamics Simulations of Wild Type and Mutants of SAPAP in Complexed with Shank3.

Specific interactions between scaffold protein SH3 and multiple ankyrin repeat domains protein 3 (Shank3) and synapse-associated protein 90/postsynaptic density-95⁻associated protein (SAPAP) are essential for excitatory synapse development and plasticity. In a bunch of human neurological diseases, mutations on Shank3 or SAPAP are detected. To investigate the dynamical and thermodynamic properties of the specific binding between the N-terminal extended PDZ (Post-synaptic density-95/Discs large/Zonaoccludens-1) domain (N-PDZ) of Shank3 and the extended PDZ binding motif (E-PBM) of SAPAP, molecular dynamics simulation approaches were used to study the complex of N-PDZ with wild type and mutated E-PBM peptides. To compare with the experimental data, 974QTRL977 and 966IEIYI970 of E-PBM peptide were mutated to prolines to obtain the M4P and M5P system, respectively. Conformational analysis shows that the canonical PDZ domain is stable while the ßN extension presents high flexibility in all systems, especially for M5P. The high flexibility of ßN extension seems to set up a barrier for the non-specific binding in this area and provide the basis for specific molecular recognition between Shank3 and SAPAP. The wild type E-PBM tightly binds to N-PDZ during the simulation while loss of binding is observed in different segments of the mutated E-PBM peptides. Energy decomposition and hydrogen bonds analysis show that M4P mutations only disrupt the interactions with canonical PDZ domain, but the interactions with ßN1' remain. In M5P system, although the interactions with ßN1' are abolished, the binding between peptide and the canonical PDZ domain is not affected. The results indicate that the interactions in the two-binding site, the canonical PDZ domain and the ßN1' extension, contribute to the binding between E-PBM and N-PDZ independently. The binding free energies calculated by MM/GBSA (Molecular Mechanics/Generalized Born Surface Area) are in agreement with the experimental binding affinities. Most of the residues on E-PBM contribute considerably favorable energies to the binding except A963 and D964 in the N-terminal. The study provides information to understand the molecular basis of specific binding between Shank3 and SAPAP, as well as clues for design of peptide inhibitors.

In silico identification and design of potent peptide inhibitors against PDZ-3 domain of Postsynaptic Density Protein (PSD-95).

Unique intrinsic properties of peptides like low toxicity, high biological activity, and specificity make them attractive therapeutic agents. PDZ-binding peptide inhibitors have been demonstrated for curing of Alzheimer, Parkinson, Dementia, and other central nervous system ailments. In this article, we report the successful use of an integrated computational protocol to analyze the structural basis of how peptides bind to the shallow groove of the third PDZ domain (PDZ-3) from the postsynaptic density (PSD-95) protein. This protocol employs careful and precise computational techniques for design of new strategy for predicting novel and potent peptides against PDZ protein. We attempted to generate a pharmacophore model using crystal structure of peptide inhibitor bound to the PDZ-3. A highly specific and sensitive generated pharmacophore model was used for screening virtual database generated using different combination of amino acid substitutions as well as decoy peptide database for its sensitivity and specificity. Identified hit peptides were further analyzed by docking studies, and their stability analyzed using solvated molecular dynamics. Quantum Mechanics/Molecular Mechanics (QM/MM) interaction energy and GMX-PBSA scoring schemes were used for ranking of stable peptides. Computational approach applied here generated encouraging results for identifying peptides against PDZ interaction model. The workflow can be further exercised as a virtual screening technique for reducing the search space for candidate target peptides against PDZ domains.