Identification of specificity and promiscuity of PDZ domain interactions through their dynamic behavior

The PDZ domain of the E protein in SARS-CoV induces carcinogenesis and poor prognosis in LUAD

BACKGROUND

In both lung adenocarcinoma (LUAD) and severe acute respiratory syndrome (SARS), uncontrolled inflammation can be detected in lung tissue. The PDZ-binding motif (PBM) in the SARS-CoV-1 E protein has been demonstrated to be a virulence factor that induces a cytokine storm.

METHODS

To identify gene expression fluctuations induced by PBM, microarray sequencing data of lung tissue infected with wild-type (SARS-CoV-1-E-wt) or recombinant virus (SARS-CoV-1-E-mutPBM) were analyzed, followed by functional enrichment analysis. To understand the role of the screened genes in LUAD, overall survival and immune correlation were calculated.

RESULTS

A total of 12 genes might participate in the initial and developmental stages of LUAD through expression variation and mutation. Moreover, dysregulation of a total of 12 genes could lead to a poorer prognosis. In addition, the downregulation of MAMDC2 and ITGA8 by PBM could also affect patient prognosis. Although the conserved PBM (-D-L-L-V-) can be found at the end of the carboxyl terminus in multiple E proteins of coronaviruses, the specific function of each protein depends on the entire amino acid sequence.

CONCLUSIONS

In summary, PBM containing the SARS-CoV-1 E protein promoted the carcinogenesis of LUAD by dysregulating important gene expression profiles and subsequently influencing the immune response and overall prognosis.

E4orf1: The triple agent of adenovirus - Unraveling its roles in oncogenesis, infectious obesity and immune responses in virus replication and vector therapy

Human Adenoviruses (HAdV) are nearly ubiquitous pathogens comprising numerous sub-types that infect various tissues and organs. Among many encoded proteins that facilitate viral replication and subversion of host cellular processes, the viral E4orf1 protein has emerged as an intriguing yet under-investigated player in the complex interplay between the virus and its host. E4orf1 has gained attention as a metabolism activator and oncogenic agent, while recent research is showing that E4orf1 may play a more important role in modulating cellular pathways such as PI3K-Akt-mTOR, Ras, the immune response and further HAdV replication stages than previously anticipated. In this review, we aim to explore the structure, molecular mechanisms, and biological functions of E4orf1, shedding light on its potentially multifaceted roles during HAdV infection, including metabolic diseases and oncogenesis. Furthermore, we discuss the role of functional E4orf1 in biotechnological applications such as Adenovirus (AdV) vaccine vectors and oncolytic AdV. By dissecting the intricate relationships between HAdV types and E4orf1 proteins, this review provides valuable insights into viral pathogenesis and points to promising areas of future research.

Functional protein dynamics in a crystal

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

Functional Protein Dynamics in a Crystal

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

Allosteric Signaling in PDZ Energetic Networks: Embedding Error Analysis

Allosteric signaling in proteins has been known for some half a century, yet how the signal traverses the protein remains an active area of research. Recently, the importance of electrostatics to achieve long-range signaling has become increasingly appreciated. Our laboratory has been working on developing network approaches to capture such interactions. In this study, we turn our attention to the well-studied allosteric model protein, PDZ. We study the allosteric dynamics on a per-residue basis in key constructs involving the PDZ domain, its allosteric effector, and its peptide ligand. We utilize molecular dynamics trajectories to create the networks for the constructs to explore the allosteric effect by plotting the heat kernel results onto axes defined by principal components. We introduce a new metric to quantitate the volume sampled by a residue in the latent space. We relate our findings to PDZ and the greater field of allostery.

Network Analysis of Molecular Dynamics Sectors in the p53 Protein

Design of allosteric regulators is an emergent field in the area of drug discovery holding promise for currently untreated diseases. Allosteric regulators bind to a protein in one location and affect a distant site. The ubiquitous presence of allosteric effectors in biology and the success of serendipitously identified allosteric compounds point to the potential they hold. Although the mechanism of transmission of an allosteric signal is not unequivocally determined, one hypothesis suggests that groups of evolutionarily covarying residues within a protein, termed sectors, are conduits. A long-term goal of our lab is to allosterically modulate the activity of proteins by binding small molecules at points of allosteric control. However, methods to consistently identify such points remain unclear. Sector residues on the surfaces of proteins are a promising source of allosteric targets. Recently, we introduced molecular dynamics (MD)-based sectors; MD sectors capitalize on covariance of motion, in place of evolutionary covariance. By focusing on motional covariance, MD sectors tap into the framework of statistical mechanics afforded by the Boltzmann ensemble of structural conformations comprising the underlying data set. We hypothesized that the method of MD sectors can be used to identify a cohesive network of motionally covarying residues capable of transmitting an allosteric signal in a protein. While our initial qualitative results showed promise for the method to predict sectors, that a network of cohesively covarying residues had been produced remained an untested assumption. In this work, we apply network theory to rigorously analyze MD sectors, allowing us to quantitatively assess the biologically relevant property of network cohesiveness of sectors in the context of the tumor suppressor protein, p53. We revised the methodology for assessing and improving MD sectors. Specifically, we introduce a metric to calculate the cohesive properties of the network. Our new approach separates residues into two categories: sector residues and non-sector residues. The relatedness within each respective group is computed with a distance metric. Cohesive sector networks are identified as those that have high relatedness among the sector residues which exceeds the relatedness of the residues to the non-sector residues in terms of the correlation of motions. Our major finding was that the revised means of obtaining sectors was more efficacious than previous iterations, as evidenced by the greater cohesion of the networks. These results are discussed in the context of the development of allosteric regulators of p53 in particular and the expected applicability of the method to the drug design field in general.

Virtual screening and molecular dynamics simulation study of plant protease inhibitors against SARS-CoV-2 envelope protein

Due to the outbreak of a new strain of pandemic coronavirus, there is a huge loss of economy and health. In 2021, some vaccines are recommended as emergency licensed vaccines to protect against the virus, and efforts are continuously ongoing to evaluate the vaccine safety measures for licensed vaccines. Recently, there was an increase in the cases of a new variant of coronavirus (omicron). Envelope protein plays an important role in virus packaging and assembly. If viral assembly is blocked, there is less chance of spreading the infection to another cell.In the present study, the plant protease inhibitors (PPIs) were screened against the envelope protein of SARS CoV 2. The structures were downloaded from the protein data bank. The plant protease inhibitors cystatin-I, Eravatmin, squash, Kunitz, Bowman-Birk, Alpha-amylase inhibitors, and potato serine protease inhibitors were screened and out of them Kunitz, alpha-amylase, and squash protease inhibitors have shown maximum binding energy. The molecular dynamics simulation was performed for docked complexes showing the lowest binding energy by NMA (normal mode analysis) to visualize the motion and stability of complexes. These plant-based protease inhibitors are a good target to fight against the new emerging strain of coronavirus because plant extracted compounds are natural and there is fewer side effect than synthetic compounds.

Using azobenzene photocontrol to set proteins in motion

Controlling the activity of proteins with azobenzene photoswitches is a potent tool for manipulating their biological function. With the help of light, it is possible to change binding affinities, control allostery or manipulate complex biological processes, for example. Additionally, owing to their intrinsically fast photoisomerization, azobenzene photoswitches can serve as triggers that initiate out-of-equilibrium processes. Such switching of the activity initiates a cascade of conformational events that can be accessed with time-resolved methods. In this Review, we show how the potency of azobenzene photoswitching can be combined with transient spectroscopic techniques to disclose the order of events and experimentally observe biomolecular interactions in real time. This strategy will further our understanding of how a protein can accommodate, adapt and readjust its structure to answer an incoming signal, revealing more of the dynamical character of proteins.

Comprehensive analysis of SARS-COV-2 drug targets and pharmacological aspects in treating the COVID-19

Corona viruses are enveloped, single-stranded RNA (Ribonucleic acid) viruses and they cause pandemic diseases having a devastating effect on both human healthcare and the global economy. To date, six corona viruses have been identified as pathogenic organisms which are significantly responsible for the infection and also cause severe respiratory diseases. Among them, the novel SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2) caused a major outbreak of corona virus diseases 2019 (COVID-19). Coronaviridae family members can affects both humans and animals. In human, corona viruses cause severe acute respiratory syndrome with mild to severe outcomes. Several structural and genomics have been investigated, and the genome encodes about 28 proteins most of them with unknown function though it shares remarkable sequence identity with other proteins. There is no potent and licensed vaccine against SARS-CoV-2 and several trials are underway to investigate the possible therapeutic agents against viral infection. However, some of the antiviral drugs that have been investigated against SARS-CoV-2 are under clinical trials. In the current review we comparatively emphasize the emergence and pathogenicity of the SARS-CoV-2 and their infection and discuss the various putative drug targets of both viral and host receptors for developing effective vaccines and therapeutic combinations to overcome the viral outbreak.

The Role of Rigid Residues in Modulating TEM-1 β-Lactamase Function and Thermostability

The relationship between protein motions (i.e., dynamics) and enzymatic function has begun to be explored in β-lactamases as a way to advance our understanding of these proteins. In a recent study, we analyzed the dynamic profiles of TEM-1 (a ubiquitous class A β-lactamase) and several ancestrally reconstructed homologues. A chief finding of this work was that rigid residues that were allosterically coupled to the active site appeared to have profound effects on enzyme function, even when separated from the active site by many angstroms. In the present work, our aim was to further explore the implications of protein dynamics on β-lactamase function by altering the dynamic profile of TEM-1 using computational protein design methods. The Rosetta software suite was used to mutate amino acids surrounding either rigid residues that are highly coupled to the active site or to flexible residues with no apparent communication with the active site. Experimental characterization of ten designed proteins indicated that alteration of residues surrounding rigid, highly coupled residues, substantially affected both enzymatic activity and stability; in contrast, native-like activities and stabilities were maintained when flexible, uncoupled residues, were targeted. Our results provide additional insight into the structure-function relationship present in the TEM family of β-lactamases. Furthermore, the integration of computational protein design methods with analyses of protein dynamics represents a general approach that could be used to extend our understanding of the relationship between dynamics and function in other enzyme classes.

Viral PDZ Binding Motifs Influence Cell Behavior Through the Interaction with Cellular Proteins Containing PDZ Domains

Viruses have evolved to interact with their hosts. Some viruses such as human papilloma virus, dengue virus, SARS-CoV, or influenza virus encode proteins including a PBM that interact with cellular proteins containing PDZ domains. There are more than 400 cellular protein isoforms with these domains in the human genome, indicating that viral PBMs have a high potential to influence the behavior of the cell. In this review we analyze the most relevant cellular processes known to be affected by viral PBM-cellular PDZ interactions including the establishment of cell-cell interactions and cell polarity, the regulation of cell survival and apoptosis and the activation of the immune system. Special attention has been provided to coronavirus PBM conservation throughout evolution and to the role of the PBMs of human coronaviruses SARS-CoV and MERS-CoV in pathogenesis.

Emerging strategies on in silico drug development against COVID-19: challenges and opportunities

The importance of coronaviruses as human pathogen has been highlighted by the recent outbreak of SARS-CoV-2 leading to the search of suitable drugs to overcome respiratory infections caused by the virus. Due to the lack of specific drugs against coronavirus, the existing antiviral and antimalarial drugs are currently being administered to the patients infected with SARS-CoV-2. The scientists are also considering repurposing of some of the existing drugs as a suitable option in search of effective drugs against coronavirus till the establishment of a potent drug and/or vaccine. Computer-aided drug discovery provides a promising attempt to enable scientists to develop new and target specific drugs to combat any disease. The discovery of novel targets for COVID-19 using computer-aided drug discovery tools requires knowledge of the structure of coronavirus and various target proteins present in the virus. Targeting viral proteins will make the drug specific against the virus, thereby, increasing the chances of viral mortality. Hence, this review provides the structure of SARS-CoV-2 virus along with the important viral components involved in causing infection. It also focuses on the role of various target proteins in disease, the mechanism by which currently administered drugs act against the virus and the repurposing of few drugs. The gap arising from the absence of specific drugs is addressed by proposing potential antiviral drug targets which might provide insights into structure-based drug development against SARS-CoV-2.

Is There a Link Between the Pathogenic Human Coronavirus Envelope Protein and Immunopathology? A Review of the Literature

Since the severe acute respiratory syndrome (SARS) outbreak in 2003, human coronaviruses (hCoVs) have been identified as causative agents of severe acute respiratory tract infections. Two more hCoV outbreaks have since occurred, the most recent being SARS-CoV-2, the causative agent of coronavirus disease 2019 (COVID-19). The clinical presentation of SARS and MERS is remarkably similar to COVID-19, with hyperinflammation causing a severe form of the disease in some patients. Previous studies show that the expression of the SARS-CoV E protein is associated with the hyperinflammatory response that could culminate in acute respiratory distress syndrome (ARDS), a potentially fatal complication. This immune-mediated damage is largely caused by a cytokine storm, which is induced by significantly elevated levels of inflammatory cytokines interleukin (IL)-1β and IL-6, which are partly mediated by the expression of the SARS-CoV E protein. The interaction between the SARS-CoV E protein and the host protein, syntenin, as well as the viroporin function of SARS-CoV E, are linked to this cytokine dysregulation. This review aims to compare the clinical presentation of virulent hCoVs with a specific focus on the cause of the immunopathology. The review also proposes that inhibition of IL-1β and IL-6 in severe cases can improve patient outcome.

Specificity in PDZ-peptide interaction networks: Computational analysis and review

Globular PDZ domains typically serve as protein-protein interaction modules that regulate a wide variety of cellular functions via recognition of short linear motifs (SLiMs). Often, PDZ mediated-interactions are essential components of macromolecular complexes, and disruption affects the entire scaffold. Due to their roles as linchpins in trafficking and signaling pathways, PDZ domains are attractive targets: both for controlling viral pathogens, which bind PDZ domains and hijack cellular machinery, as well as for developing therapies to combat human disease. However, successful therapeutic interventions that avoid off-target effects are a challenge, because each PDZ domain interacts with a number of cellular targets, and specific binding preferences can be difficult to decipher. Over twenty-five years of research has produced a wealth of data on the stereochemical preferences of individual PDZ proteins and their binding partners. Currently the field lacks a central repository for this information. Here, we provide this important resource and provide a manually curated, comprehensive list of the 271 human PDZ domains. We use individual domain, as well as recent genomic and proteomic, data in order to gain a holistic view of PDZ domains and interaction networks, arguing this knowledge is critical to optimize targeting selectivity and to benefit human health.

Computational Analysis of Energy Landscapes Reveals Dynamic Features That Contribute to Binding of Inhibitors to CFTR-Associated Ligand

The CFTR-associated ligand PDZ domain (CALP) binds to the cystic fibrosis transmembrane conductance regulator (CFTR) and mediates lysosomal degradation of mature CFTR. Inhibition of this interaction has been explored as a therapeutic avenue for cystic fibrosis. Previously, we reported the ensemble-based computational design of a novel peptide inhibitor of CALP, which resulted in the most binding-efficient inhibitor to date. This inhibitor, kCAL01, was designed using osprey and evinced significant biological activity in in vitro cell-based assays. Here, we report a crystal structure of kCAL01 bound to CALP and compare structural features against iCAL36, a previously developed inhibitor of CALP. We compute side-chain energy landscapes for each structure to not only enable approximation of binding thermodynamics but also reveal ensemble features that contribute to the comparatively efficient binding of kCAL01. Finally, we compare the previously reported design ensemble for kCAL01 vs the new crystal structure and show that, despite small differences between the design model and crystal structure, significant biophysical features that enhance inhibitor binding are captured in the design ensemble. This suggests not only that ensemble-based design captured thermodynamically significant features observed in vitro, but also that a design eschewing ensembles would miss the kCAL01 sequence entirely.

Coronavirus envelope protein: current knowledge

BACKGROUND

Coronaviruses (CoVs) primarily cause enzootic infections in birds and mammals but, in the last few decades, have shown to be capable of infecting humans as well. The outbreak of severe acute respiratory syndrome (SARS) in 2003 and, more recently, Middle-East respiratory syndrome (MERS) has demonstrated the lethality of CoVs when they cross the species barrier and infect humans. A renewed interest in coronaviral research has led to the discovery of several novel human CoVs and since then much progress has been made in understanding the CoV life cycle. The CoV envelope (E) protein is a small, integral membrane protein involved in several aspects of the virus' life cycle, such as assembly, budding, envelope formation, and pathogenesis. Recent studies have expanded on its structural motifs and topology, its functions as an ion-channelling viroporin, and its interactions with both other CoV proteins and host cell proteins.

MAIN BODY

This review aims to establish the current knowledge on CoV E by highlighting the recent progress that has been made and comparing it to previous knowledge. It also compares E to other viral proteins of a similar nature to speculate the relevance of these new findings. Good progress has been made but much still remains unknown and this review has identified some gaps in the current knowledge and made suggestions for consideration in future research.

CONCLUSIONS

The most progress has been made on SARS-CoV E, highlighting specific structural requirements for its functions in the CoV life cycle as well as mechanisms behind its pathogenesis. Data shows that E is involved in critical aspects of the viral life cycle and that CoVs lacking E make promising vaccine candidates. The high mortality rate of certain CoVs, along with their ease of transmission, underpins the need for more research into CoV molecular biology which can aid in the production of effective anti-coronaviral agents for both human CoVs and enzootic CoVs.

Emerging Themes in PDZ Domain Signaling: Structure, Function, and Inhibition

Post-synaptic density-95, disks-large and zonula occludens-1 (PDZ) domains are small globular protein-protein interaction domains widely conserved from yeast to humans. They are composed of ∼90 amino acids and form a classical two α-helical/six β-strand structure. The prototypical ligand is the C-terminus of partner proteins; however, they also bind internal peptide sequences. Recent findings indicate that PDZ domains also bind phosphatidylinositides and cholesterol. Through their ligand interactions, PDZ domain proteins are critical for cellular trafficking and the surface retention of various ion channels. In addition, PDZ proteins are essential for neuronal signaling, memory, and learning. PDZ proteins also contribute to cytoskeletal dynamics by mediating interactions critical for maintaining cell-cell junctions, cell polarity, and cell migration. Given their important biological roles, it is not surprising that their dysfunction can lead to multiple disease states. As such, PDZ domain-containing proteins have emerged as potential targets for the development of small molecular inhibitors as therapeutic agents. Recent data suggest that the critical binding function of PDZ domains in cell signaling is more than just glue, and their binding function can be regulated by phosphorylation or allosterically by other binding partners. These studies also provide a wealth of structural and biophysical data that are beginning to reveal the physical features that endow this small modular domain with a central role in cell signaling.

Mutations Utilize Dynamic Allostery to Confer Resistance in TEM-1 β-lactamase

β-lactamases are enzymes produced by bacteria to hydrolyze β-lactam antibiotics as a common mechanism of resistance. Evolution in such enzymes has been rendering a wide variety of antibiotics impotent, therefore posing a major threat. Clinical and in vitro studies of evolution in TEM-1 β-lactamase have revealed a large number of single point mutations that are responsible for driving resistance to antibiotics and/or inhibitors. The distal locations of these mutations from the active sites suggest that these allosterically modulate the antibiotic resistance. We investigated the effects of resistance driver mutations on the conformational dynamics of the enzyme to provide insights about the mechanism of their long-distance interactions. Through all-atom molecular dynamics (MD) simulations, we obtained the dynamic flexibility profiles of the variants and compared those with that of the wild type TEM-1. While the mutational sites in the variants did not have any direct van der Waals interactions with the active site position S70 and E166, we observed a change in the flexibility of these sites, which play a very critical role in hydrolysis. Such long distance dynamic interactions were further confirmed by dynamic coupling index (DCI) analysis as the sites involved in resistance driving mutations exhibited high dynamic coupling with the active sites. A more exhaustive dynamic analysis, using a selection pressure for ampicillin and cefotaxime resistance on all possible types of substitutions in the amino acid sequence of TEM-1, further demonstrated the observed mechanism. Mutational positions that play a crucial role for the emergence of resistance to new antibiotics exhibited high dynamic coupling with the active site irrespective of their locations. These dynamically coupled positions were neither particularly rigid nor particularly flexible, making them more evolvable positions. Nature utilizes these sites to modulate the dynamics of the catalytic sites instead of mutating the highly rigid positions around the catalytic site.

Ancient thioredoxins evolved to modern-day stability-function requirement by altering native state ensemble

Thioredoxins (THRXs)-small globular proteins that reduce other proteins-are ubiquitous in all forms of life, from Archaea to mammals. Although ancestral thioredoxins share sequential and structural similarity with the modern-day (extant) homologues, they exhibit significantly different functional activity and stability. We investigate this puzzle by comparative studies of their (ancient and modern-day THRXs') native state ensemble, as quantified by the dynamic flexibility index (DFI), a metric for the relative resilience of an amino acid to perturbations in the rest of the protein. Clustering proteins using DFI profiles strongly resemble an alternative classification scheme based on their activity and stability. The DFI profiles of the extant proteins are substantially different around the α3, α4 helices and catalytic regions. Likewise, allosteric coupling of the active site with the rest of the protein is different between ancient and extant THRXs, possibly explaining the decreased catalytic activity at low pH with evolution. At a global level, we note that the population of low-flexibility (called hinges) and high-flexibility sites increases with evolution. The heterogeneity (quantified by the variance) in DFI distribution increases with the decrease in the melting temperature typically associated with the evolution of ancient proteins to their modern-day counterparts.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

The maximum penalty criterion for ridge regression: application to the calibration of the force constant in elastic network models

Tikhonov regularization, or ridge regression, is a popular technique to deal with collinearity in multivariate regression. We unveil a formal analogy between ridge regression and statistical mechanics, where the objective function is comparable to a free energy, and the ridge parameter plays the role of temperature. This analogy suggests two novel criteria for selecting a suitable ridge parameter: specific-heat (C_v) and maximum penalty (MP). We apply these fits to evaluate the relative contributions of rigid-body and internal fluctuations, which are typically highly collinear, to crystallographic B-factors. This issue is particularly important for computational models of protein dynamics, such as the elastic network model (ENM), since the amplitude of the predicted internal motion is commonly calibrated using B-factor data. After validation on simulated datasets, our results indicate that rigid-body motions account on average for more than 80% of the amplitude of B-factors. Furthermore, we evaluate the ability of different fits to reproduce the amplitudes of internal fluctuations in X-ray ensembles from the B-factors in the corresponding single X-ray structures. The new ridge criteria are shown to be markedly superior to the commonly used two-parameter fit that neglects rigid-body rotations and to the full fits regularized under generalized cross-validation. In conclusion, the proposed fits ensure a more robust calibration of the ENM force constant and should prove valuable in other applications.

Mechanisms for allosteric activation of protease DegS by ligand binding and oligomerization as revealed from molecular dynamics simulations

A local perturbation of a protein may lead to functional changes at some distal site, a phenomenon denoted as allostery. Here, we study the allosteric control of a protease using molecular dynamics simulations. The system considered is the bacterial protein DegS which includes a protease domain activated on ligand binding to an adjacent PDZ domain. Starting from crystallographic structures of DegS homo-trimers, we perform simulations of the ligand-free and -bound state of DegS at equilibrium. Considering a single protomer only, the trimeric state was mimicked by applying restraints on the residues in contact with other protomers in the DegS trimer. In addition, the bound state was also simulated without any restraints to mimic the monomer. Our results suggest that not only ligand release but also disassembly of a DegS trimer inhibits proteolytic activity. Considering various observables for structural changes, we infer allosteric pathways from the interface with other protomers to the active site. Moreover, we study how ligand release leads to (i) catalytically relevant changes involving residues 199-201 and (ii) a transition from a stretched to a bent conformation for residues 217-219 (which prohibits proper substrate binding). Finally, based on ligand-induced C_α shifts we identify residues in contact with other protomers in the DegS trimer that likely transduce the perturbation from ligand release from a given protomer to adjacent protomers. These residues likely play a key role in the experimentally known effect of ligand release from a protomer on the proteolytic activity of the other protomers. Proteins 2016; 84:1690-1705. © 2016 Wiley Periodicals, Inc.

Characterization of Promiscuous Binding of Phosphor Ligands to Breast-Cancer-Gene 1 (BRCA1) C-Terminal (BRCT): Molecular Dynamics, Free Energy, Entropy and Inhibitor Design

Inhibition of the protein-protein interaction (PPI) mediated by breast-cancer-gene 1 C-terminal (BRCT) is an attractive strategy to sensitize breast and ovarian cancers to chemotherapeutic agents that induce DNA damage. Such inhibitors could also be used for studies to understand the role of this PPI in DNA damage response. However, design of BRCT inhibitors is challenging because of the inherent flexibility associated with this domain. Several studies identified short phosphopeptides as tight BRCT binders. Here we investigated the thermodynamic properties of 18 phosphopeptides or peptide with phosphate mimic and three compounds with phosphate groups binding to BRCT to understand promiscuous molecular recognition and guide inhibitor design. We performed molecular dynamics (MD) simulations to investigate the interactions between inhibitors and BRCT and their dynamic behavior in the free and bound states. MD simulations revealed the key role of loops in altering the shape and size of the binding site to fit various ligands. The mining minima (M2) method was used for calculating binding free energy to explore the driving forces and the fine balance between configuration entropy loss and enthalpy gain. We designed a rigidified ligand, which showed unfavorable experimental binding affinity due to weakened enthalpy. This was because it lacked the ability to rearrange itself upon binding. Investigation of another phosphate group containing compound, C1, suggested that the entropy loss can be reduced by preventing significant narrowing of the energy well and introducing multiple new compound conformations in the bound states. From our computations, we designed an analog of C1 that introduced new intermolecular interactions to strengthen attractions while maintaining small entropic penalty. This study shows that flexible compounds do not always encounter larger entropy penalty, compared with other more rigid binders, and highlights a new strategy for inhibitor design.

Promiscuous proteins are commonly observed in biological systems, such as modular domains that recognize phosphopeptides during signal transduction. The use of phosphopeptides and compounds with phosphate groups as inhibitors to protein–protein interactions have attracted increasing interest for years. By using atomistic molecular dynamics simulations, we are able to perform detailed analyses of the dihedral space to explore protein fluctuation upon ligand binding to better understand promiscuous molecular recognition. Free energy calculation can further provide insights into the mechanism of binding, including both enthalpic and entropic contributions for molecular recognition, which assist in inhibitor design. Our calculation results show that pre-rigidifying a ligand is not always advantageous, suggesting the challenge in retaining optimized intermolecular interactions in pre-rigidified ligand. Instead, certain flexible ligands with multiple binding conformations can reduce entropic penalty, and therefore improves binding affinity. According to our computations, we can introduce new intermolecular interactions to flexible ligand to strengthen attractions while maintaining small entropic penalty by retaining its plasticity in the bound conformation. The study might cast light on a new general strategy for designing inhibitors targeting promiscuous modular domains and protein–protein interactions.

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

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 in the backbone, especially in the β2-β3 loop but also in the β5-α2 and the α2-β6 loop, and are accompanied by changes of inter-residue contacts mainly within the β2-β3 loop as well as between the α2 helix and other segments. The residues showing substantial changes of inter-residue contacts, backbone conformations, or C(α) positions are considered "key residues" for the long-range conformational response of PDZ2. By comparing these residues with various sets of residues highlighted by previous studies of PDZ2, we investigate the statistical correlation of the various approaches. Interestingly, we find a considerable correlation of our findings with several works considering structural changes but no significant correlations with approaches considering energy flow or networks based on inter-residue energies.

Identification of the Mechanisms Causing Reversion to Virulence in an Attenuated SARS-CoV for the Design of a Genetically Stable Vaccine

A SARS-CoV lacking the full-length E gene (SARS-CoV-∆E) was attenuated and an effective vaccine. Here, we show that this mutant virus regained fitness after serial passages in cell culture or in vivo, resulting in the partial duplication of the membrane gene or in the insertion of a new sequence in gene 8a, respectively. The chimeric proteins generated in cell culture increased virus fitness in vitro but remained attenuated in mice. In contrast, during SARS-CoV-∆E passage in mice, the virus incorporated a mutated variant of 8a protein, resulting in reversion to a virulent phenotype. When the full-length E protein was deleted or its PDZ-binding motif (PBM) was mutated, the revertant viruses either incorporated a novel chimeric protein with a PBM or restored the sequence of the PBM on the E protein, respectively. Similarly, after passage in mice, SARS-CoV-∆E protein 8a mutated, to now encode a PBM, and also regained virulence. These data indicated that the virus requires a PBM on a transmembrane protein to compensate for removal of this motif from the E protein. To increase the genetic stability of the vaccine candidate, we introduced small attenuating deletions in E gene that did not affect the endogenous PBM, preventing the incorporation of novel chimeric proteins in the virus genome. In addition, to increase vaccine biosafety, we introduced additional attenuating mutations into the nsp1 protein. Deletions in the carboxy-terminal region of nsp1 protein led to higher host interferon responses and virus attenuation. Recombinant viruses including attenuating mutations in E and nsp1 genes maintained their attenuation after passage in vitro and in vivo. Further, these viruses fully protected mice against challenge with the lethal parental virus, and are therefore safe and stable vaccine candidates for protection against SARS-CoV.

Zoonotic coronaviruses, including SARS-CoV, Middle East respiratory syndrome (MERS-CoV), porcine epidemic diarrhea virus (PEDV) and swine delta coronavirus (SDCoV) have recently emerged causing high morbidity and mortality in human or piglets. No fully protective therapy is still available for these CoVs. Therefore, the development of efficient vaccines is a high priority. Live attenuated vaccines are considered most effective compared to other types of vaccines, as they induce a long-lived, balanced immune response. However, safety is the main concern of this type of vaccines because attenuated viruses can eventually revert to a virulent phenotype. Therefore, an essential feature of any live attenuated vaccine candidate is its stability. In addition, introduction of several safety guards is advisable to increase vaccine safety. In this manuscript, we analyzed the mechanisms by which an attenuated SARS-CoV reverted to a virulent phenotype and describe the introduction of attenuating deletions that maintained virus stability. The virus, engineered with two safety guards, provided full protection against challenge with a lethal SARS-CoV. Understanding the molecular mechanisms leading to pathogenicity and the in vivo evaluation of vaccine genetic stability contributed to a rational design of a promising SARS-CoV vaccine.

Principles of allosteric interactions in cell signaling

Linking cell signaling events to the fundamental physicochemical basis of the conformational behavior of single molecules and ultimately to cellular function is a key challenge facing the life sciences. Here we outline the emerging principles of allosteric interactions in cell signaling, with emphasis on the following points. (1) Allosteric efficacy is not a function of the chemical composition of the allosteric pocket but reflects the extent of the population shift between the inactive and active states. That is, the allosteric effect is determined by the extent of preferred binding, not by the overall binding affinity. (2) Coupling between the allosteric and active sites does not decide the allosteric effect; however, it does define the propagation pathways, the allosteric binding sites, and key on-path residues. (3) Atoms of allosteric effectors can act as "driver" or "anchor" and create attractive "pulling" or repulsive "pushing" interactions. Deciphering, quantifying, and integrating the multiple co-occurring events present daunting challenges to our scientific community.

Identifying key residues for protein allostery through rigid residue scan

Allostery is a ubiquitous process for protein regulatory activity in which a binding event can change a protein's function carried out at a distal site. Despite intensive theoretical and experimental investigation of protein allostery in the past five decades, effective methods have yet to be developed that can systematically identify key residues involved in allosteric mechanisms. In this study, we propose the rigid residue scan as a systematic approach to identify important allosteric residues. The third PDZ domain (PDZ3) in the postsynaptic density 95 protein (PSD-95) is used as a model system, and each amino acid residue is treated as a single rigid body during independent molecular dynamics simulations. Various indices based on cross-correlation matrices are used, which allow for two groups of residues with different functions to be identified. The first group is proposed as "switches" that are needed to "turn on" the binding effect of protein allostery. The second group is proposed as "wire residues" that are needed to propagate energy or information from the binding site to distal locations within the same protein. Among the nine residues suggested as important for PDZ3 intramolecular communication in this study, eight have been reported as critical for allostery in PDZ3. Therefore, the rigid residue scan approach is demonstrated to be an effective method for systemically identifying key residues in protein intramolecular communication and allosteric mechanisms.

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

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 the free-bound transition. A detailed description of the conformational rearrangement following the cis-trans photoisomerization of PDZ2S reveals a series of photoinduced structural changes that propagate from the anchor residues of the photoswitch via intermediate secondary structure segments to the C-terminus of PDZ2S. The changes of the conformational distribution of the C-terminal region is considered as the distal response of the isolated allosteric protein.

Evolution of conformational dynamics determines the conversion of a promiscuous generalist into a specialist enzyme

β-Lactamases are produced by many modern bacteria as a mechanism of resistance toward β-lactam antibiotics, the most common antibiotics in use. β-Lactamases, however, are ancient enzymes that originated billions of years ago. Recently, proteins corresponding to 2- to 3-Gy-old Precambrian nodes in the evolution of Class A β-lactamases have been prepared and shown to be moderately efficient promiscuous catalysts, able to degrade a variety of antibiotics with catalytic efficiency levels similar to those of an average modern enzyme. Remarkably, there are few structural differences (in particular at the active-site regions) between the resurrected enzymes and a penicillin-specialist modern β-lactamase. Here, we propose that the ancestral promiscuity originates from conformational dynamics. We investigate the differences in conformational dynamics of the ancient and extant β-lactamases through MD simulations and quantify the contribution of each position to functionally related dynamics through Dynamic Flexibility Index. The modern TEM-1 lactamase shows a comparatively rigid active-site region, likely reflecting adaptation for efficient degradation of a specific substrate (penicillin), whereas enhanced deformability at the active-site neighborhood in the ancestral resurrected proteins likely accounts for the binding and subsequent degradation of antibiotic molecules of different size and shape. Clustering of the conformational dynamics on the basis of Principal Component Analysis is in agreement with the functional divergence, as the ancient β-lactamases cluster together, separated from their modern descendant. Finally, our analysis leads to testable predictions, as sites of potential relevance for the evolution of dynamics are identified and mutations at those sites are expected to alter substrate-specificity.

The PDZ-binding motif of severe acute respiratory syndrome coronavirus envelope protein is a determinant of viral pathogenesis

A recombinant severe acute respiratory syndrome coronavirus (SARS-CoV) lacking the envelope (E) protein is attenuated in vivo. Here we report that E protein PDZ-binding motif (PBM), a domain involved in protein-protein interactions, is a major determinant of virulence. Elimination of SARS-CoV E protein PBM by using reverse genetics caused a reduction in the deleterious exacerbation of the immune response triggered during infection with the parental virus and virus attenuation. Cellular protein syntenin was identified to bind the E protein PBM during SARS-CoV infection by using three complementary strategies, yeast two-hybrid, reciprocal coimmunoprecipitation and confocal microscopy assays. Syntenin redistributed from the nucleus to the cell cytoplasm during infection with viruses containing the E protein PBM, activating p38 MAPK and leading to the overexpression of inflammatory cytokines. Silencing of syntenin using siRNAs led to a decrease in p38 MAPK activation in SARS-CoV infected cells, further reinforcing their functional relationship. Active p38 MAPK was reduced in lungs of mice infected with SARS-CoVs lacking E protein PBM as compared with the parental virus, leading to a decreased expression of inflammatory cytokines and to virus attenuation. Interestingly, administration of a p38 MAPK inhibitor led to an increase in mice survival after infection with SARS-CoV, confirming the relevance of this pathway in SARS-CoV virulence. Therefore, the E protein PBM is a virulence domain that activates immunopathology most likely by using syntenin as a mediator of p38 MAPK induced inflammation.

SARS-CoV caused a worldwide epidemic infecting 8000 people with a mortality of about 10%. A recombinant SARS-CoV lacking the E protein was attenuated in vivo. The E protein contains a PDZ-binding motif (PBM), a domain potentially involved in the interaction with more than 400 cellular proteins, which highlights its relevance in modulating host-cell behavior. To analyze the contributions of this motif to virulence, recombinant viruses with or without E protein PBM were generated. Recombinant SARS-CoVs lacking E protein PBM caused minimal lung damage and were attenuated, in contrast to viruses containing this motif, indicating that E protein PBM is a virulence determinant. E protein PBM induces the deleterious exacerbated immune response triggered during SARS-CoV infection, and interacts with the cellular protein syntenin, as demonstrated using proteomic analyses. Interestingly, syntenin redistributed from nucleus to cytoplasm during SARS-CoV infection, activating p38 MAPK and triggering the overexpression of inflammatory cytokines. Furthermore, silencing of syntenin using siRNAs led to a decrease in p38 MAPK activation. In addition, administration of a p38 MAPK inhibitor led to an increase in mice survival after SARS-CoV infection. These results indicate that syntenin and p38 MAPK are potential therapeutic targets to reduce the exacerbated immune response during SARS-CoV infection.

Achieving peptide binding specificity and promiscuity by loops: case of the forkhead-associated domain

The regulation of a series of cellular events requires specific protein–protein interactions, which are usually mediated by modular domains to precisely select a particular sequence from diverse partners. However, most signaling domains can bind to more than one peptide sequence. How do proteins create promiscuity from precision? Moreover, these complex interactions typically occur at the interface of a well-defined secondary structure, α helix and β sheet. However, the molecular recognition primarily controlled by loop architecture is not fully understood. To gain a deep understanding of binding selectivity and promiscuity by the conformation of loops, we chose the forkhead-associated (FHA) domain as our model system. The domain can bind to diverse peptides via various loops but only interact with sequences containing phosphothreonine (pThr). We applied molecular dynamics (MD) simulations for multiple free and bound FHA domains to study the changes in conformations and dynamics. Generally, FHA domains share a similar folding structure whereby the backbone holds the overall geometry and the variety of sidechain atoms of multiple loops creates a binding surface to target a specific partner. FHA domains determine the specificity of pThr by well-organized binding loops, which are rigid to define a phospho recognition site. The broad range of peptide recognition can be attributed to different arrangements of the loop interaction network. The moderate flexibility of the loop conformation can help access or exclude binding partners. Our work provides insights into molecular recognition in terms of binding specificity and promiscuity and helpful clues for further peptide design.

BP-Dock: a flexible docking scheme for exploring protein-ligand interactions based on unbound structures

Molecular docking serves as an important tool in modeling protein-ligand interactions. However, it is still challenging to incorporate overall receptor flexibility, especially backbone flexibility, in docking due to the large conformational space that needs to be sampled. To overcome this problem, we developed a novel flexible docking approach, BP-Dock (Backbone Perturbation-Dock) that can integrate both backbone and side chain conformational changes induced by ligand binding through a multi-scale approach. In the BP-Dock method, we mimic the nature of binding-induced events as a first-order approximation by perturbing the residues along the protein chain with a small Brownian kick one at a time. The response fluctuation profile of the chain upon these perturbations is computed using the perturbation response scanning method. These response fluctuation profiles are then used to generate binding-induced multiple receptor conformations for ensemble docking. To evaluate the performance of BP-Dock, we applied our approach on a large and diverse data set using unbound structures as receptors. We also compared the BP-Dock results with bound and unbound docking, where overall receptor flexibility was not taken into account. Our results highlight the importance of modeling backbone flexibility in docking for recapitulating the experimental binding affinities, especially when an unbound structure is used. With BP-Dock, we can generate a wide range of binding site conformations realized in nature even in the absence of a ligand that can help us to improve the accuracy of unbound docking. We expect that our fast and efficient flexible docking approach may further aid in our understanding of protein-ligand interactions as well as virtual screening of novel targets for rational drug design.

New conformational state of NHERF1-CXCR2 signaling complex captured by crystal lattice trapping

NHERF1 is a PDZ adaptor protein that scaffolds the assembly of diverse signaling complexes and has been implicated in many cancers. However, little is known about the mechanism responsible for its scaffolding promiscuity or its ability to bind to multiple targets. Computational studies have indicated that PDZ promiscuity may be attributed to its conformational dynamics, but experimental evidence for this relationship remains very limited. Here we examine the conformational flexibility of the NHERF1 PDZ1 domain using crystal lattice trapping via solving PDZ1 structure of a new crystal form. The structure, together with prior PDZ1 structures of a different space group, reveals that 4 of 11 ligand-interacting residues undergo significant crystal packing-induced structural changes. Most of these residues correspond to the residues involved in allosteric transition when a peptide ligand binds. In addition, a subtle difference in ligand conformations causes the same peptide to bind in slightly different modes in different crystal forms. These findings indicate that substantial structural flexibility is present in the PDZ1 peptide-binding pocket, and the structural substate trapped in the present crystal form can be utilized to represent the conformational space accessible to the protein. Such knowledge will be critical for drug design against the NHERF1 PDZ1 domain, highlighting the continued need for experimentally determined PDZ1-ligand complexes.

Simultaneous prediction of binding free energy and specificity for PDZ domain-peptide interactions

Interactions between protein domains and linear peptides underlie many biological processes. Among these interactions, the recognition of C-terminal peptides by PDZ domains is one of the most ubiquitous. In this work, we present a mathematical model for PDZ domain-peptide interactions capable of predicting both affinity and specificity of binding based on X-ray crystal structures and comparative modeling with ROSETTA. We developed our mathematical model using a large phage display dataset describing binding specificity for a wild type PDZ domain and 91 single mutants, as well as binding affinity data for a wild type PDZ domain binding to 28 different peptides. Structural refinement was carried out through several ROSETTA protocols, the most accurate of which included flexible peptide docking and several iterations of side chain repacking and backbone minimization. Our findings emphasize the importance of backbone flexibility and the energetic contributions of side chain-side chain hydrogen bonds in accurately predicting interactions. We also determined that predicting PDZ domain-peptide interactions became increasingly challenging as the length of the peptide increased in the N-terminal direction. In the training dataset, predicted binding energies correlated with those derived through calorimetry and specificity switches introduced through single mutations at interface positions were recapitulated. In independent tests, our best performing protocol was capable of predicting dissociation constants well within one order of magnitude of the experimental values and specificity profiles at the level of accuracy of previous studies. To our knowledge, this approach represents the first integrated protocol for predicting both affinity and specificity for PDZ domain-peptide interactions.

Effective harmonic potentials: insights into the internal cooperativity and sequence-specificity of protein dynamics

The proper biological functioning of proteins often relies on the occurrence of coordinated fluctuations around their native structure, or on their ability to perform wider and sometimes highly elaborated motions. Hence, there is considerable interest in the definition of accurate coarse-grained descriptions of protein dynamics, as an alternative to more computationally expensive approaches. In particular, the elastic network model, in which residue motions are subjected to pairwise harmonic potentials, is known to capture essential aspects of conformational dynamics in proteins, but has so far remained mostly phenomenological, and unable to account for the chemical specificities of amino acids. We propose, for the first time, a method to derive residue- and distance-specific effective harmonic potentials from the statistical analysis of an extensive dataset of NMR conformational ensembles. These potentials constitute dynamical counterparts to the mean-force statistical potentials commonly used for static analyses of protein structures. In the context of the elastic network model, they yield a strongly improved description of the cooperative aspects of residue motions, and give the opportunity to systematically explore the influence of sequence details on protein dynamics.

The role of flexibility and conformational selection in the binding promiscuity of PDZ domains

In molecular recognition, it is often the case that ligand binding is coupled to conformational change in one or both of the binding partners. Two hypotheses describe the limiting cases involved; the first is the induced fit and the second is the conformational selection model. The conformational selection model requires that the protein adopts conformations that are similar to the ligand-bound conformation in the absence of ligand, whilst the induced-fit model predicts that the ligand-bound conformation of the protein is only accessible when the ligand is actually bound. The flexibility of the apo protein clearly plays a major role in these interpretations. For many proteins involved in signaling pathways there is the added complication that they are often promiscuous in that they are capable of binding to different ligand partners. The relationship between protein flexibility and promiscuity is an area of active research and is perhaps best exemplified by the PDZ domain family of proteins. In this study we use molecular dynamics simulations to examine the relationship between flexibility and promiscuity in five PDZ domains: the human Dvl2 (Dishevelled-2) PDZ domain, the human Erbin PDZ domain, the PDZ1 domain of InaD (inactivation no after-potential D protein) from fruit fly, the PDZ7 domain of GRIP1 (glutamate receptor interacting protein 1) from rat and the PDZ2 domain of PTP-BL (protein tyrosine phosphatase) from mouse. We show that despite their high structural similarity, the PDZ binding sites have significantly different dynamics. Importantly, the degree of binding pocket flexibility was found to be closely related to the various characteristics of peptide binding specificity and promiscuity of the five PDZ domains. Our findings suggest that the intrinsic motions of the apo structures play a key role in distinguishing functional properties of different PDZ domains and allow us to make predictions that can be experimentally tested.

Proteins that are capable of binding to many different ligands are said to have broad specificity. This is sometimes also referred to as promiscuity. Whether a protein is promiscuous or not can sometimes be readily explained by the structure of the protein and the ligand in terms of electrostatic and steric effects. Sometimes however, this simple interpretation can struggle to explain the experimentally observed data. A prominent case in point is the PDZ domains. These small protein domains bind to unstructured regions of other proteins and are involved in many signaling pathways. Some PDZ domains appear to be more promiscuous than others, but this has been difficult to explain purely on the basis of the composition of residues in the binding groove. In this work we examine the dynamics and conformational flexibility of five key PDZ domains and demonstrate that despite similar folds, these proteins can exhibit quite different dynamics. Furthermore the difference in the dynamic behavior appears to correlate with the observed promiscuity. Our findings suggest that knowledge of the dynamic behavior of the PDZs can be used to rationalize the extent of expected promiscuity. Such knowledge will be critical for drug design against PDZ domains.

Dynamics of endoglucanase catalytic domains: implications towards thermostability

Thermostable endoglucanases play a crucial role in the production of biofuels to breakdown plant cellulose. Analyzing their structure-dynamics relationship can inform about the origins of their thermostability. Although tertiary structures of many endoglucanase proteins are available, the relationship between thermostability, structure, and dynamics is not explored fully. We have generated elastic network models for thermostable and mesostable endoglucanases with the (αβ)₈ fold in substrate bound and unbound states. The comparative analyses shed light on the relation between protein dynamics, thermostability, and substrate binding. We observed specific differences in the dynamic behavior of catalytic residues in slow modes: while both the nucleophile and the acid/base donor residues show positively correlated motions in the thermophile, their dynamics is uncoupled in the mesophile. Our proof-of-concept comparison study suggests that global dynamics can be harnessed to further our understanding of thermostability.

Hybrid approaches to molecular simulation

Molecular dynamics (MD) simulation is an established method for studying the conformational changes that are important for protein function. Recent advances in hardware and software have allowed MD simulations over the same timescales as experiment, improving the agreement between theory and experiment to a large extent. However, running such simulations are costly, in terms of resources, storage, and trajectory analysis. There is still a place for techniques that involve short MD simulations. In order to overcome the sampling paucity of short time-scales, hybrid methods that include some form of MD simulation can exploit certain features of the system of interest, often combining experimental information in surprising ways. Here, we review some recent hybrid approaches to the simulation of proteins.

Dissipative electro-elastic network model of protein electrostatics

We propose a dissipative electro-elastic network model to describe the dynamics and statistics of electrostatic fluctuations at active sites of proteins. The model combines the harmonic network of residue beads with overdamped dynamics of the normal modes of the network characterized by two friction coefficients. The electrostatic component is introduced to the model through atomic charges of the protein force field. The overall effect of the electrostatic fluctuations of the network is recorded through the frequency-dependent response functions of the electrostatic potential and electric field at the protein active site. We also consider the dynamics of displacements of individual residues in the network and the dynamics of distances between pairs of residues. The model is tested against loss spectra of residue displacements and the electrostatic potential and electric field at the heme's iron from all-atom molecular dynamics simulations of three hydrated globular proteins.

Plasticity of PDZ domains in ligand recognition and signaling

The PDZ domain is a protein-protein interacting module that plays an important role in the organization of signaling complexes. The recognition of short intrinsically disordered C-terminal peptide motifs is the archetypical PDZ function, but the functional repertoire of this versatile module also includes recognition of internal peptide sequences, dimerization and phospholipid binding. The PDZ function can be tuned by various means such as allosteric effects, changes of physiological buffer conditions and phosphorylation of PDZ domains and/or ligands, which poses PDZ domains as dynamic regulators of cell signaling. This review is focused on the plasticity of the PDZ interactions.

The binding affinities of proteins interacting with the PDZ domain of PICK1

Protein interacting with C kinase (PICK1) is well conserved throughout evolution and plays a critical role in synaptic plasticity by regulating the trafficking and posttranslational modification of its interacting proteins. PICK1 contains a single PSD95/DlgA/Zo-1 (PDZ) protein-protein interaction domain, which is promiscuous and shown to interact with over 60 proteins, most of which play roles in neuronal function. Several reports have suggested the role of PICK1 in disorders such as epilepsy, pain, brain trauma and stroke, drug abuse and dependence, schizophrenia and psychosis. Importantly, lead compounds that block PICK1 interactions are also now becoming available. Here, a new modeling approach was developed to investigate binding affinities of PDZ interactions. Using these methods, the binding affinities of all major PICK1 interacting proteins are reported and the effects of PICK1 mutations on these interactions are described. These modeling methods have important implications in defining the binding properties of proteins interacting with PICK1 as well as the general structural requirements of PDZ interactions. The study also provides modeling methods to support in the drug design of ligands for PDZ domains, which may further aid in development of the family of PDZ domains as a drug target.

Putting into practice domain-linear motif interaction predictions for exploration of protein networks

PDZ domains recognise short sequence motifs at the extreme C-termini of proteins. A model based on microarray data has been recently published for predicting the binding preferences of PDZ domains to five residue long C-terminal sequences. Here we investigated the potential of this predictor for discovering novel protein interactions that involve PDZ domains. When tested on real negative data assembled from published literature, the predictor displayed a high false positive rate (FPR). We predicted and experimentally validated interactions between four PDZ domains derived from the human proteins MAGI1 and SCRIB and 19 peptides derived from human and viral C-termini of proteins. Measured binding intensities did not correlate with prediction scores, and the high FPR of the predictor was confirmed. Results indicate that limitations of the predictor may arise from an incomplete model definition and improper training of the model. Taking into account these limitations, we identified several novel putative interactions between PDZ domains of MAGI1 and SCRIB and the C-termini of the proteins FZD4, ARHGAP6, NET1, TANC1, GLUT7, MARCH3, MAS, ABC1, DLL1, TMEM215 and CYSLTR2. These proteins are localised to the membrane or suggested to act close to it and are often involved in G protein signalling. Furthermore, we showed that, while extension of minimal interacting domains or peptides toward tandem constructs or longer peptides never suppressed their ability to interact, the measured affinities and inferred specificity patterns often changed significantly. This suggests that if protein fragments interact, the full length proteins are also likely to interact, albeit possibly with altered affinities and specificities. Therefore, predictors dealing with protein fragments are promising tools for discovering protein interaction networks but their application to predict binding preferences within networks may be limited.

Change in allosteric network affects binding affinities of PDZ domains: analysis through perturbation response scanning

The allosteric mechanism plays a key role in cellular functions of several PDZ domain proteins (PDZs) and is directly linked to pharmaceutical applications; however, it is a challenge to elaborate the nature and extent of these allosteric interactions. One solution to this problem is to explore the dynamics of PDZs, which may provide insights about how intramolecular communication occurs within a single domain. Here, we develop an advancement of perturbation response scanning (PRS) that couples elastic network models with linear response theory (LRT) to predict key residues in allosteric transitions of the two most studied PDZs (PSD-95 PDZ3 domain and hPTP1E PDZ2 domain). With PRS, we first identify the residues that give the highest mean square fluctuation response upon perturbing the binding sites. Strikingly, we observe that the residues with the highest mean square fluctuation response agree with experimentally determined residues involved in allosteric transitions. Second, we construct the allosteric pathways by linking the residues giving the same directional response upon perturbation of the binding sites. The predicted intramolecular communication pathways reveal that PSD-95 and hPTP1E have different pathways through the dynamic coupling of different residue pairs. Moreover, our analysis provides a molecular understanding of experimentally observed hidden allostery of PSD-95. We show that removing the distal third alpha helix from the binding site alters the allosteric pathway and decreases the binding affinity. Overall, these results indicate that (i) dynamics plays a key role in allosteric regulations of PDZs, (ii) the local changes in the residue interactions can lead to significant changes in the dynamics of allosteric regulations, and (iii) this might be the mechanism that each PDZ uses to tailor their binding specificities regulation.

PDZ domain proteins (PDZs) act as adapters in organizing functional protein complexes. Through dynamic interactions, PDZs play a key role in mediating key cellular functions in the cell, and they are linked to currently challenging diseases including Alzheimer's, Parkinson's and cancer. Moreover, they are associated with allosteric regulations in mediating signaling. Therefore, it is critical to have knowledge of how the allosteric transition occurs in PDZs. We investigate the allosteric response of the two most studied PDZs, PSD-95 and hPTP1E, using the perturbation response scanning (PRS) approach. The method treats the protein as an elastic network and uses linear response theory (LRT) to obtain residue fluctuations upon exerting directed random forces on selected residues. With this efficient and fast approach, we identify the key residues that mediate long-range communication and find the allosteric pathways. Although the structures of PSD-95 and hPTP1E are very similar, our analysis predicts that their allosteric pathways are different. We also observe a significant change in allosteric pathways and a decrease in binding affinity upon removal of the distal α3 helix of PSD-95. This approach enables us to understand how dynamic interactions play an important role in allosteric regulations.

Binding free energy landscape of domain-peptide interactions

Peptide recognition domains (PRDs) are ubiquitous protein domains which mediate large numbers of protein interactions in the cell. How these PRDs are able to recognize peptide sequences in a rapid and specific manner is incompletely understood. We explore the peptide binding process of PDZ domains, a large PRD family, from an equilibrium perspective using an all-atom Monte Carlo (MC) approach. Our focus is two different PDZ domains representing two major PDZ classes, I and II. For both domains, a binding free energy surface with a strong bias toward the native bound state is found. Moreover, both domains exhibit a binding process in which the peptides are mostly either bound at the PDZ binding pocket or else interact little with the domain surface. Consistent with this, various binding observables show a temperature dependence well described by a simple two-state model. We also find important differences in the details between the two domains. While both domains exhibit well-defined binding free energy barriers, the class I barrier is significantly weaker than the one for class II. To probe this issue further, we apply our method to a PDZ domain with dual specificity for class I and II peptides, and find an analogous difference in their binding free energy barriers. Lastly, we perform a large number of fixed-temperature MC kinetics trajectories under binding conditions. These trajectories reveal significantly slower binding dynamics for the class II domain relative to class I. Our combined results are consistent with a binding mechanism in which the peptide C terminal residue binds in an initial, rate-limiting step.

The complex biological processes occurring in living organisms are enabled by numerous networks of interacting proteins. It is therefore of great interest to understand the physical interplay between proteins and, in particular, how this process gives rise to highly specific network connectivities. For a long time, the dominant molecular view of protein-protein interactions was the docking of more or less static folded structures, with specificity obtained from a complementarity in shape and charge distributions. Lately it has been realized that many of the links in protein networks are mediated by interactions between folded domains, on the one hand, and disordered polypeptide segments, on the other. We use an all-atom Monte Carlo based approach which attempts to capture this domain-peptide binding process in full and apply it to representative members of a common domain family. This allows us to examine and compare detailed aspects of the binding free energy landscapes which underlie specificity and affinity. Being able to model domain-peptide binding in a physically sound, yet computationally tractable way is essential for identifying molecular binding mechanisms and opens up possibilities for modifying interaction networks in a controlled way.

Detailed regulatory mechanism of the interaction between ZO-1 PDZ2 and connexin43 revealed by MD simulations

The gap junction protein connexin43 (Cx43) binds to the second PDZ domain of Zonula occludens-1 (ZO-1) through its C-terminal tail, mediating the regulation of gap junction plaque size and dynamics. Biochemical study demonstrated that the very C-terminal 12 residues of Cx43 are necessary and sufficient for ZO-1 PDZ2 binding and phosphorylation at residues Ser (-9) and Ser (-10) of the peptide can disrupt the association. However, only a crystal structure of ZO-1 PDZ2 in complex with a shorter 9 aa peptide of connexin43 was solved experimentally. Here, the interactions between ZO-1 PDZ2 and the short, long and phosphorylated Cx43 peptides were studied using molecular dynamics (MD) simulations and free energy calculation. The short peptide bound to PDZ2 exhibits large structural variations, while the extension of three upstream residues stabilizes the peptide conformation and enhanced the interaction. Phosphorylation at Ser(-9) significantly weakens the binding and results in conformational flexibility of the peptide. Glu210 of ZO-1 PDZ2 was found to be a key regulatory point in Cx43 binding and phosphorylation induced dissociation.

Genome-wide analysis of PDZ domain binding reveals inherent functional overlap within the PDZ interaction network

Binding selectivity and cross-reactivity within one of the largest and most abundant interaction domain families, the PDZ family, has long been enigmatic. The complete human PDZ domain complement (the PDZome) consists of 267 domains and we applied here a Bayesian selectivity model to predict hundreds of human PDZ domain interactions, using target sequences of 22,997 non-redundant proteins. Subsequent analysis of these binding scores shows that PDZs can be divided into two genome-wide clusters that coincide well with the division between canonical class 1 and 2 PDZs. Within the class 1 PDZs we observed binding overlap at unprecedented levels, mediated by two residues at positions 1 and 5 of the second α-helix of the binding pocket. Eight PDZ domains were subsequently selected for experimental binding studies and to verify the basics of our predictions. Overall, the PDZ domain class 1 cross-reactivity identified here implies that auxiliary mechanisms must be in place to overcome this inherent functional overlap and to minimize cross-selectivity within the living cell. Indeed, when we superimpose PDZ domain binding affinities with gene ontologies, network topology data and the domain position within a PDZ superfamily protein, functional overlap is minimized and PDZ domains position optimally in the binding space. We therefore propose that PDZ domain selectivity is achieved through cellular context rather than inherent binding specificity.

The structural and dynamic response of MAGI-1 PDZ1 with noncanonical domain boundaries to the binding of human papillomavirus E6

PDZ domains are protein interaction domains that are found in cytoplasmic proteins involved in signaling pathways and subcellular transport. Their roles in the control of cell growth, cell polarity, and cell adhesion in response to cell contact render this family of proteins targets during the development of cancer. Targeting of these network hubs by the oncoprotein E6 of "high-risk" human papillomaviruses (HPVs) serves to effect the efficient disruption of cellular processes. Using NMR, we have solved the three-dimensional solution structure of an extended construct of the second PDZ domain of MAGI-1 (MAGI-1 PDZ1) alone and bound to a peptide derived from the C-terminus of HPV16 E6, and we have characterized the changes in backbone dynamics and hydrogen bonding that occur upon binding. The binding event induces quenching of high-frequency motions in the C-terminal tail of the PDZ domain, which contacts the peptide upstream of the canonical X-[T/S]-X-[L/V] binding motif. Mutations designed in the C-terminal flanking region of the PDZ domain resulted in a significant decrease in binding affinity for E6 peptides. This detailed analysis supports the notion of a global response of the PDZ domain to the binding event, with effects propagated to distal sites, and reveals unexpected roles for the sequences flanking the canonical PDZ domain boundaries.

Crystallographic and nuclear magnetic resonance evaluation of the impact of peptide binding to the second PDZ domain of protein tyrosine phosphatase 1E

PDZ (PSD95/Discs large/ZO-1) domains are ubiquitous protein interaction motifs found in scaffolding proteins involved in signal transduction. Despite the fact that many PDZ domains show a limited tendency to undergo structural change, the PDZ family has been associated with long-range communication and allostery. One of the PDZ domains studied most in terms of structure and biophysical properties is the second PDZ ("PDZ2") domain from protein tyrosine phosphatase 1E (PTP1E, also known as PTPL1). Previously, we showed through NMR relaxation studies that binding of the RA-GEF2 C-terminal peptide substrate results in long-range propagation of side-chain dynamic changes in human PDZ2 [Fuentes, E. J., et al. (2004) J. Mol. Biol. 335, 1105-1115]. Here, we present the first X-ray crystal structures of PDZ2 in the absence and presence of RA-GEF2 ligand, determined to resolutions of 1.65 and 1.3 Å, respectively. These structures deviate somewhat from previously determined NMR structures and indicate that very minor structural changes in PDZ2 accompany peptide binding. NMR residual dipolar couplings confirm the crystal structures to be accurate models of the time-averaged atomic coordinates of PDZ2. The impact on side-chain dynamics was further tested with a C-terminal peptide from APC, which showed results nearly identical to those of RA-GEF2. Thus, allosteric transmission in PDZ2 induced by peptide binding is conveyed purely and robustly by dynamics. (15)N relaxation dispersion measurements did not detect appreciable populations of a kinetic structural intermediate. Collectively, for ligand binding to PDZ2, these data support a lock-and-key binding model from a structural perspective and an allosteric model from a dynamical perspective, which together suggest a complex energy landscape for functional transitions within the ensemble.

Manipulation of conformational change in proteins by single-residue perturbations

Using the perturbation-response scanning (PRS) technique, we study a set of 25 proteins that display a variety of conformational motions upon ligand binding (e.g., shear, hinge, allosteric). In most cases, PRS determines single residues that may be manipulated to achieve the resulting conformational change. PRS reveals that for some proteins, binding-induced conformational change may be achieved through the perturbation of residues scattered throughout the protein, whereas in others, perturbation of specific residues confined to a highly specific region is necessary. Overlaps between the experimental and PRS-calculated atomic displacement vectors are usually more descriptive of the conformational change than those obtained from a modal analysis of elastic network models. Furthermore, the largest overlaps obtained by the latter approach do not always appear in the most collective modes; there are cases where more than one mode yields comparable overlap sizes. We show that success of the modal analysis depends on an absence of redundant paths in the protein. PRS thus demonstrates that several relevant modes can be induced simultaneously by perturbing a single select residue on the protein. We also illustrate the biological relevance of applying PRS to the GroEL, adenylate kinase, myosin, and kinesin structures in detail by showing that the residues whose perturbation leads to precise conformational changes usually correspond to those experimentally determined to be functionally important.