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

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.

Dishevelled2 activates WGEF via its interaction with a unique internal peptide motif of the GEF

The Wnt-planar cell polarity (Wnt-PCP) pathway is crucial in establishing cell polarity during development and tissue homoeostasis. This pathway is found to be dysregulated in many pathological conditions, including cancer and autoimmune disorders. The central event in Wnt-PCP pathway is the activation of Weak-similarity guanine nucleotide exchange factor (WGEF) by the adapter protein Dishevelled (Dvl). The PDZ domain of Dishevelled2 (Dvl2PDZ) binds and activates WGEF by releasing it from its autoinhibitory state. However, the actual Dvl2PDZ binding site of WGEF and the consequent activation mechanism of the GEF have remained elusive. Using biochemical and molecular dynamics studies, we show that a unique "internal-PDZ binding motif" (IPM) of WGEF mediates the WGEF-Dvl2PDZ interaction to activate the GEF. The residues at P2, P0, P-2 and P-3 positions of IPM play an important role in stabilizing the WGEFpep-Dvl2PDZ interaction. Furthermore, MD simulations of modelled Dvl2PDZ-WGEFIPM peptide complexes suggest that WGEF-Dvl2PDZ interaction may differ from the reported Dvl2PDZ-IPM interactions. Additionally, the apo structure of human Dvl2PDZ shows conformational dynamics different from its IPM peptide bound state, suggesting an induced fit mechanism for the Dvl2PDZ-peptide interaction. The current study provides a model for Dvl2 induced activation of WGEF.

Pharmacophore mapping approach to find anti-cancer phytochemicals with metformin-like activities against transforming growth factor (TGF)-beta receptor I kinase: An in silico study

The most frequently prescribed first-line treatment for type II diabetes mellitus is metformin. Recent reports asserted that this diabetes medication can also shield users from cancer. Metformin induces cell cycle arrest in cancer cells. However, the exact mechanism by which this occurs in the cancer system is yet to be elucidated. Here, we investigated the impact of metformin on cell cycle arrest in cancer cells utilizing transforming growth factor (TGF)-beta pathway. TGF-ß pathway has significant effect on cell progression and growth. In order to gain an insight on the underlying molecular mechanism of metformin's effect on TGF beta receptor 1 kinase, molecular docking was performed. Metformin was predicted to interact with transforming growth factor (TGF)-beta receptor I kinase based on molecular docking and molecular dynamics simulations. Furthermore, pharmacophore was generated for metformin-TGF-ßR1 complex to hunt for novel compounds having similar pharmacophore as metformin with enhanced anti-cancer potentials. Virtual screening with 29,000 natural compounds from NPASS database was conducted separately for the generated pharmacophores in Ligandscout® software. Pharmacophore mapping showed 60 lead compounds for metformin-TGF-ßR1 complex. Molecular docking, molecular dynamics simulation for 100 ns and ADMET analysis were performed on these compounds. Compounds with CID 72473, 10316977 and 45140078 showed promising binding affinities and formed stable complexes during dynamics simulation with aforementioned protein and thus have potentiality to be developed into anti-cancer medicaments.

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

General Theory of Specific Binding: Insights from a Genetic-Mechano-Chemical Protein Model

Proteins need to selectively interact with specific targets among a multitude of similar molecules in the cell. However, despite a firm physical understanding of binding interactions, we lack a general theory of how proteins evolve high specificity. Here, we present such a model that combines chemistry, mechanics, and genetics and explains how their interplay governs the evolution of specific protein-ligand interactions. The model shows that there are many routes to achieving molecular discrimination-by varying degrees of flexibility and shape/chemistry complementarity-but the key ingredient is precision. Harder discrimination tasks require more collective and precise coaction of structure, forces, and movements. Proteins can achieve this through correlated mutations extending far from a binding site, which fine-tune the localized interaction with the ligand. Thus, the solution of more complicated tasks is enabled by increasing the protein size, and proteins become more evolvable and robust when they are larger than the bare minimum required for discrimination. The model makes testable, specific predictions about the role of flexibility and shape mismatch in discrimination, and how evolution can independently tune affinity and specificity. Thus, the proposed theory of specific binding addresses the natural question of "why are proteins so big?". A possible answer is that molecular discrimination is often a hard task best performed by adding more layers to the protein.

Refinement of SARS-CoV-2 envelope protein structure in a native-like environment by molecular dynamics simulations

COVID-19 has become an unprecedented threat to human health. The SARS-CoV-2 envelope (E) protein plays a critical role in the viral maturation process and pathogenesis. Despite intensive investigation, its structure in physiological conditions remains mysterious: no high-resolution full-length structure is available and only an NMR structure of the transmembrane (TM) region has been determined. Here, we present a refined E protein structure, using molecular dynamics (MD) simulations to investigate its structure and dynamics in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer system. Our initial homology model based upon the SARS-CoV E protein structure is shown to be unstable in the lipid bilayer, and the H3 helices tend to move away from the membrane center to the membrane-water interface. A more stable model was developed by replacing all H3 helices with the fully equilibrated H3 structure sampled in the MD simulations. This refined model exhibited more favorable contacts with lipids and water than the original homology model and induced local membrane curvature, decreasing local lipid order. Interestingly, the pore radius profiles showed that the channel in both homology and refined models remained in a closed state throughout the simulations. We also demonstrated the utility of this structure to develop anti-SARS-CoV-2 drugs by docking a library of FDA-approved, investigational, and experimental drugs to the refined E protein structure, identifying 20 potential channel blockers. This highlights the power of MD simulations to refine low-resolution structures of membrane proteins in a native-like membrane environment, shedding light on the structural features of the E protein and providing a platform for the development of novel antiviral treatments.

Proteostasis Regulators in Cystic Fibrosis: Current Development and Future Perspectives

In cystic fibrosis (CF), the deletion of phenylalanine 508 (F508del) in the CF transmembrane conductance regulator (CFTR) leads to misfolding and premature degradation of the mutant protein. These defects can be targeted with pharmacological agents named potentiators and correctors. During the past years, several efforts have been devoted to develop and approve new effective molecules. However, their clinical use remains limited, as they fail to fully restore F508del-CFTR biological function. Indeed, the search for CFTR correctors with different and additive mechanisms has recently increased. Among them, drugs that modulate the CFTR proteostasis environment are particularly attractive to enhance therapy effectiveness further. This Perspective focuses on reviewing the recent progress in discovering CFTR proteostasis regulators, mainly describing the design, chemical structure, and structure-activity relationships. The opportunities, challenges, and future directions in this emerging and promising field of research are discussed, as well.

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.

Structural communication between the GTPase Sec4p and its activator Sec2p: Determinants of GEF activity and early deformations to nucleotide release

Ras GTPases are molecular switches that cycle between OFF and ON states depending on the bound nucleotide (i.e. GDP-bound and GTP-bound, respectively).

The Rab GTPase, Sec4p, plays regulatory roles in multiple steps of intracellular vesicle trafficking. Nucleotide release is catalyzed by the Guanine Nucleotide Exchange Factor (GEF) Sec2p.

Here, the integration of structural information with molecular dynamics (MD) simulations addressed a number of questions concerning the intrinsic and stimulated dynamics of Sec2p and Sec4p as well as the chain of structural deformations leading to GEF-assisted activation of the Rab GTPase.

Sec2p holds an intrinsic ability to adopt the conformation found in the crystallographic complexes with Sec4p, thus suggesting that the latter selects and shifts the conformational equilibrium towards a pre-existing bound-like conformation of Sec2p.

The anchoring of Sec4p to a suitable conformation of Sec2p favors the Sec2p-assisted pulling on itself of the α1/switch 1 (SWI) loop and of SWI, which loose any contact with GDP. Those deformations of Sec4p would occur earlier. Formation of the final Sec2p-Sec4p hydrophobic interface, accomplishes later. Disruption of the nucleotide cage would cause firstly loss of interactions with the guanine ring and secondly loss of interactions with the phosphates.

The ease in sampling the energy landscape and adopting a bound-like conformation likely favors the catalyzing ability of GEFs for Ras GTPases.

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.

Targeting PDZ domains as potential treatment for viral infections, neurodegeneration and cancer

The interaction between proteins is a fundamental event for cellular life that is generally mediated by specialized protein domains or modules. PDZ domains are the largest class of protein-protein interaction modules, involved in several cellular pathways such as signal transduction, cell-cell junctions, cell polarity and adhesion, and protein trafficking. Because of that, dysregulation of PDZ domain function often causes the onset of pathologies, thus making this family of domains an interesting pharmaceutical target. In this review article we provide an overview of the structural and functional features of PDZ domains and their involvement in the cellular and molecular pathways at the basis of different human pathologies. We also discuss some of the strategies that have been developed with the final goal to hijack or inhibit the interaction of PDZ domains with their ligands. Because of the generally low binding selectivity of PDZ domain and the scarce efficiency of small molecules in inhibiting PDZ binding, this task resulted particularly difficult to pursue and still demands increasing experimental efforts in order to become completely feasible and successful in vivo.

Discovery of potent and selective Cdc25 phosphatase inhibitors via rapid assembly and in situ screening of Quinonoid-focused libraries

Cell division cycle 25 (Cdc25) phosphatase is an attractive target for drug discovery. The rapid assembly and in situ screening of focused combinatorial fragment libraries using efficient modular reactions is a highly robust strategy for discovering bioactive molecules. In this study, we have utilized miniaturized synthesis to generate several quinonoid-focused libraries, by standard CuAAC reaction and HBTU-based amide coupling chemistry. Then the enzyme inhibition screening afforded some potent and selective Cdc25s inhibitors. Compound M5N36 (Cdc25A: IC₅₀ = 0.15 ± 0.05 μM; Cdc25B: IC₅₀ = 0.19 ± 0.06 μM; Cdc25C: IC₅₀ = 0.06 ± 0.04 μM) exhibited higher inhibitory activity than the initial lead NSC663284 (Cdc25A: IC₅₀ = 0.27 ± 0.02 μM; Cdc25B: IC₅₀ = 0.42 ± 0.01 μM; Cdc25C: IC₅₀ = 0.23 ± 0.01 μM). Moreover, M5N36 displayed about three-fold more potent against Cdc25C than Cdc25A and B, indicating that M5N36 could act as a relatively selective Cdc25C inhibitor. Cell viability evaluation, western blotting and molecular simulations provided a mechanistic understanding of the activity of M5N36. It showed promising anti-growth activity against the MDA-MB-231 cell line and desirable predicted physicochemical properties. Overall, M5N36 was proven to be a promising novel Cdc25C inhibitor.

Drug targets, mechanisms of drug action, and therapeutics against SARS-CoV-2

The recent outbreak of COVID-19 has affected human lives severely. The human-to-human transmission of this viral disease has become deadly due to the unavailability of COVID-19 specific drugs. Here, an overview of various attempts made to design different therapeutic agents against various structural and non-structural proteins of SARS-CoV-2 has been summarized. Emphasis has been made to highlight the mechanisms of drug action and ways to design better inhibitors of these proteins. The roles of anti-oxidants and vitamins in suppressing COVID-19 are also discussed.

Phosphorylation of SNX27 by MAPK11/14 links cellular stress-signaling pathways with endocytic recycling

Endocytosed proteins can be delivered to lysosomes for degradation or recycled to either the trans-Golgi network or the plasma membrane. It remains poorly understood how the recycling versus degradation of cargoes is determined. Here, we show that multiple extracellular stimuli, including starvation, LPS, IL-6, and EGF treatment, can strongly inhibit endocytic recycling of multiple cargoes through the activation of MAPK11/14. The stress-induced kinases in turn directly phosphorylate SNX27, a key regulator of endocytic recycling, at serine 51 (Ser51). Phosphorylation of SNX27 at Ser51 alters the conformation of its cargo-binding pocket and decreases the interaction between SNX27 and cargo proteins, thereby inhibiting endocytic recycling. Our study indicates that endocytic recycling is highly dynamic and can crosstalk with cellular stress–signaling pathways. Suppression of endocytic recycling and enhancement of receptor lysosomal degradation serve as new mechanisms for cells to cope with stress and save energy.

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.

Development of Peptide-Based PDZ Domain Inhibitors

Over the past decades, peptide-based drugs have gained increasing interest in a wide range of treatment applications, primarily because of high potency and selectivity, as well as good efficacy, tolerability, and safety often achieved with peptides. Attempts to target postsynaptic density protein of 95 (PSD-95) PSD-95/Discs large/Zonula occludens-1 (PDZ) domains, which mediate the formation of a ternary complex with the N-methyl-D-aspartate (NMDA) receptor and neuronal nitric oxide synthase (nNOS) responsible for excitotoxicity in ischemic stroke, by high-affinity small molecules have failed in the past. In this chapter, we focus on the discovery of peptide-based drugs targeting PSD-95, using AVLX-144 as an example, from the synthesis, over binding assays to its target, to further in vitro experiments based on the development of AVLX-144, a potential stroke treatment, which is planned to enter clinical trials in 2020.

Protonation-Induced Dynamic Allostery in PDZ Domain: Evidence of Perturbation-Independent Universal Response Network

Dynamic allostery is a relatively new paradigm where certain external perturbations may lead to modulation of conformational dynamics at a distant part of a protein without significant changes in the overall structure. While most well-characterized examples of dynamic allostery involve binding with other entities like small molecules, peptides, or nucleic acids, in this work we demonstrate that chemical modifications like protonation may lead to significant dynamical allosteric response in a PDZ domain protein. Tuning the protonation states of two histidine residues (H317 and H372), we identify the allosteric pathways responsible for the dynamic response. Interestingly, the same set of residues that constitute the allosteric response network upon ligand binding seem to be responsible for protonation-induced dynamic allostery. Thus, we propose the existence of an inherent universal response network in signaling proteins, where the same set of residues can respond to varying types of external perturbations in terms of rearrangement of hydrogen-bonded network and redistribution of electrostatic interaction energies.

Managing the COVID-19 Pandemic: Research Strategies Based on the Evolutionary and Molecular Characteristics of Coronaviruses

Coronavirus disease 2019 (COVID-19), an ongoing global health emergency, is a highly transmittable and pathogenic viral infection caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Emerging in Wuhan, China, in December 2019, it spread widely across the world causing panic-worst ever economic depression is visibly predictable. Coronaviruses (CoVs) have emerged as a major public health concern having caused three zoonotic outbreaks; severe acute respiratory syndrome-CoV (SARS-CoV) in 2002-2003, Middle East respiratory syndrome-CoV (MERS-CoV) in 2012, and currently this devastating COVID-19. Research strategies focused on understanding the evolutionary origin, transmission, and molecular basis of SARS-CoV-2 and its pathogenesis need to be urgently formulated to manage the current and possible future coronaviral outbreaks. Current response to the COVID-19 outbreak has been largely limited to monitoring/containment. Although frantic global efforts for developing safe and effective prophylactic and therapeutic agents are on, no licensed antiviral treatment or vaccine exists till date. In this review, research strategies for coping with COVID-19 based on evolutionary and molecular aspects of coronaviruses have been proposed.

Dynamics and structural communication in the ternary complex of fully phosphorylated V2 vasopressin receptor, vasopressin, and β-arrestin 1

G protein-coupled receptors (GPCRs) are critically regulated by arrestins, which not only desensitize G-protein signaling but also initiate a G protein-independent wave of signaling. The information from structure determination was herein exploited to build a structural model of the ternary complex, comprising fully phosphorylated V2 vasopressin receptor (V2R), the agonist arginine vasopressin (AVP), and β-arrestin 1 (β-arr1). Molecular simulations served to explore dynamics and structural communication in the ternary complex. Flexibility and mechanical profiles reflect fold of V2R and β-arr1. Highly conserved amino acids tend to behave as hubs in the structure network and contribute the most to the mechanical rigidity of V2R seven-helix bundle and of β-arr1. Two structurally and dynamically distinct receptor-arrestin interfaces assist the twist of the N- and C-terminal domains (ND and CD, respectively) of β-arr1 with respect to each other, which is linked to arrestin activation. While motion of the ND is essentially assisted by the fully phosphorylated C-tail of V2R (V2RCt), that of CD is assisted by the second and third intracellular loops and the cytosolic extensions of helices 5 and 6. In the presence of the receptor, the β-arr1 inter-domain twist angle correlates with the modes describing the essential subspace of the ternary complex. β-arr1 motions are also influenced by the anchoring to the membrane of the C-edge-loops in the β-arr1-CD. Overall fluctuations reveal a coupling between motions of the agonist binding site and of β-arr1-ND, which are in allosteric communication between each other. Mechanical rigidity points, often acting as hubs in the structure network and distributed along the main axis of the receptor helix bundle, contribute to establish a preferential communication pathway between agonist ligand and the ND of arrestin. Such communication, mediated by highly conserved amino acids, involves also the first amino acid in the arrestin C-tail, which is highly dynamic and is involved in clathrin-mediated GPCR internalization.

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.

Interconnecting Flexibility, Structural Communication, and Function in RhoGEF Oncoproteins

Dbl family Rho guanine nucleotide exchange factors (RhoGEFs) play a central role in cell biology by catalyzing the exchange of guanosine 5'-triphosphate for guanosine 5'-diphosphate (GDP) on RhoA. Insights into the oncogenic constitutive activity of the Lbc RhoGEF were gained by analyzing the structure and dynamics of the protein in different functional states and in comparison with a close homologue, leukemia-associated RhoGEF. Higher intrinsic flexibility, less dense and extended structure network, and less stable allosteric communication pathways in Lbc, compared to the nonconstitutively active homologue, emerged as major determinants of the constitutive activity. Independent of the state, the essential dynamics of the two RhoGEFs is contributed by the last 10 amino acids of Dbl homology (DH) and the whole pleckstrin homology (PH) domains and tends to be equalized by the presence of RhoA. The catalytic activity of the RhoGEF relies on the scaffolding action of the DH domain that primarily turns the switch I (SWI) of RhoA on itself through highly conserved amino acids participating in the stability core and essential for function. Changes in the conformation of SWI and disorganization of the RhoA regions deputed to nucleotide binding are among the major RhoGEF effects leading to GDP release. Binding of RhoA reorganizes the allosteric communication on RhoGEF, strengthening the communication among the canonical RhoA binding site on DH, a secondary RhoA binding site on PH, and the binding site for heterotrimeric G proteins, suggesting dual roles for RhoA as a catalysis substrate and as a regulatory protein. The structure network-based analysis tool employed in this study proved to be useful for predicting potentially druggable regulatory sites in protein structures.

Conformational Dynamics and Cooperativity Drive the Specificity of a Protein-Ligand Interaction

Molecular recognition is critical for the fidelity of signal transduction in biology. Conversely, the disruption of protein-protein interactions can lead to disease. Thus, comprehension of the molecular determinants of specificity is essential for understanding normal biological signaling processes and for the development of precise therapeutics. Although high-resolution structures have provided atomic details of molecular interactions, much less is known about the influence of cooperativity and conformational dynamics. Here, we used the Tiam2 PSD-95/Dlg/ZO-1 (PDZ) domain and a quadruple mutant (QM), engineered by swapping the identity of four residues important for specificity in the Tiam1 PDZ into the Tiam2 PDZ domain, as a model system to investigate the role of cooperativity and dynamics in PDZ ligand specificity. Surprisingly, equilibrium binding experiments found that the ligand specificity of the Tiam2 QM was switched to that of the Tiam1 PDZ. NMR-based studies indicated that Tiam2 QM PDZ, but not other mutants, had extensive microsecond to millisecond motions distributed throughout the entire domain suggesting structural cooperativity between the mutated residues. Thermodynamic analyses revealed energetic cooperativity between residues in distinct specificity subpockets that was dependent upon the identity of the ligand, indicating a context-dependent binding mechanism. Finally, isothermal titration calorimetry experiments showed distinct entropic signatures along the mutational trajectory from the Tiam2 wild-type to the QM PDZ domain. Collectively, our studies provide unique insights into how structure, conformational dynamics, and thermodynamics combine to modulate ligand-binding specificity and have implications for the evolution, regulation, and design of protein-ligand interactions.

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.

Allosteric Modulation of Binding Specificity by Alternative Packing of Protein Cores

Hydrophobic cores are often viewed as tightly packed and rigid, but they do show some plasticity and could thus be attractive targets for protein design. Here we explored the role of different functional pressures on the core packing and ligand recognition of the SH3 domain from human Fyn tyrosine kinase. We randomized the hydrophobic core and used phage display to select variants that bound to each of three distinct ligands. The three evolved groups showed remarkable differences in core composition, illustrating the effect of different selective pressures on the core. Changes in the core did not significantly alter protein stability, but were linked closely to changes in binding affinity and specificity. Structural analysis and molecular dynamics simulations revealed the structural basis for altered specificity. The evolved domains had significantly reduced core volumes, which in turn induced increased backbone flexibility. These motions were propagated from the core to the binding surface and induced significant conformational changes. These results show that alternative core packing and consequent allosteric modulation of binding interfaces could be used to engineer proteins with novel functions.

Role of Severe Acute Respiratory Syndrome Coronavirus Viroporins E, 3a, and 8a in Replication and Pathogenesis

Viroporins are viral proteins with ion channel (IC) activity that play an important role in several processes, including virus replication and pathogenesis. While many coronaviruses (CoVs) encode two viroporins, severe acute respiratory syndrome CoV (SARS-CoV) encodes three: proteins 3a, E, and 8a. Additionally, proteins 3a and E have a PDZ-binding motif (PBM), which can potentially bind over 400 cellular proteins which contain a PDZ domain, making them potentially important for the control of cell function. In the present work, a comparative study of the functional motifs included within the SARS-CoV viroporins was performed, mostly focusing on the roles of the IC and PBM of E and 3a proteins. Our results showed that the full-length E and 3a proteins were required for maximal SARS-CoV replication and virulence, whereas viroporin 8a had only a minor impact on these activities. A virus missing both the E and 3a proteins was not viable, whereas the presence of either protein with a functional PBM restored virus viability. E protein IC activity and the presence of its PBM were necessary for virulence in mice. In contrast, the presence or absence of the homologous motifs in protein 3a did not influence virus pathogenicity. Therefore, dominance of the IC and PBM of protein E over those of protein 3a was demonstrated in the induction of pathogenesis in mice.IMPORTANCE Collectively, these results demonstrate key roles for the ion channel and PBM domains in optimal virus replication and pathogenesis and suggest that the viral viroporins and PBMs are suitable targets for antiviral therapy and for mutation in attenuated SARS-CoV vaccines.

Destabilizing mutations encode nongenetic variation that drives evolutionary innovation

Evolutionary innovations are often achieved by repurposing existing genes to perform new functions; however, the mechanisms enabling the transition from old to new remain controversial. We identified mutations in bacteriophage λ's host-recognition gene J that confer enhanced adsorption to λ's native receptor, LamB, and the ability to access a new receptor, OmpF. The mutations destabilize λ particles and cause conformational bistability of J, which yields progeny of multiple phenotypic forms, each proficient at different receptors. This work provides an example of how nongenetic protein variation can catalyze an evolutionary innovation. We propose that cases where a single genotype can manifest as multiple phenotypes may be more common than previously expected and offer a general mechanism for evolutionary innovation.

Structure function relations in PDZ-domain-containing proteins: Implications for protein networks in cellular signalling

Protein scaffolds as essential backbones for organization of supramolecular signalling complexes are a recurrent theme in several model systems. Scaffold proteins preferentially employ linear peptide binding motifs for recruiting their interaction partners. PDZ domains are one of the more commonly encountered peptide binding domains in several proteins including those involved in scaffolding functions. This domain is known for its promiscuity both in terms of ligand selection, mode of interaction with its ligands as well as its association with other protein interaction domains. PDZ domains are subject to several means of regulations by virtue of their functional diversity. Additionally, the PDZ domains are refractive to the effect of mutations and maintain their three-dimensional architecture under extreme mutational load. The biochemical and biophysical basis for this selectivity as well as promiscuity has been investigated and reviewed extensively. The present review focuses on the plasticity inherent in PDZ domains and its implications for modular organization as well as evolution of cellular signalling pathways in higher eukaryotes.

Activation processes in ligand-activated G protein-coupled receptors: A case study of the adenosine A2A receptor

Here we review concepts related to an ensemble description of G-protein-coupled receptors (GPCRs). The ensemble is characterized by both inactive and active states, whose equilibrium populations and exchange rates depend sensitively on ligand, environment, and allosteric factors. This review focuses on the adenosine A₂ receptor (A_2A R), a prototypical class A GPCR. ¹⁹ F Nuclear Magnetic Resonance (NMR) studies show that apo A_2A R is characterized by a broad ensemble of conformers, spanning inactive to active states, and resembling states defined earlier for rhodopsin. In keeping with ideas associated with a conformational selection mechanism, addition of agonist serves to allosterically restrict the overall degrees of freedom at the G protein binding interface and bias both states and functional dynamics to facilitate G protein binding and subsequent activation. While the ligand does not necessarily "induce" activation, it does bias sampling of states, increase the cooperativity of the activation process and thus, the lifetimes of functional activation intermediates, while restricting conformational dynamics to that needed for activation.

Nucleotide Dependent Switching in Rho GTPase: Conformational Heterogeneity and Competing Molecular Interactions

Ras superfamily of GTPases regulate myriad cellular processes through a conserved nucleotide (GTP/GDP) dependent switching mechanism. Unlike Ras family of GTPases, for the Rho GTPases, there is no clear evidence for the existence of “sub-states” such as state 1 & state 2 in the GTP bound form. To explore the nucleotide dependent conformational space of the Switch I loop and also to look for existence of state 1 like conformations in Rho GTPases, atomistic molecular dynamics and metadynamics simulations on RhoA were performed. These studies demonstrate that both the nucleotide-free state and the GDP bound “OFF” state have very similar conformations, whereas the GTP bound “ON” state has unique conformations with signatures of two intermediate states. The conformational free energy landscape for these systems suggests the presence of multiple intermediate states. Interestingly, the energetic penalty of exposing the non-polar residues in the GTP bound form is counter balanced by the favourable hydrogen bonded interactions between the γ-phosphate group of GTP with the highly conserved Tyr34 and Thr37 residues. These competing molecular interactions lead to a tuneable energy landscape of the Switch I conformation, which can undergo significant changes based on the local environment including changes upon binding to effectors.

Structural Determinants of Constitutive Activation of Gα Proteins: Transducin as a Paradigm

Heterotrimeric guanine nucleotide-binding proteins (Gα proteins) are intracellular nanomachines deputed to signal transduction. The switch-on process requires the release of bound GDP from a site at the interface between GTPase and helical domains. Nucleotide release is catalyzed by G protein Coupled Receptors (GPCRs). Here we investigate the functional dynamics of wild type (WT) and six constitutively active mutants (CAMs) of the Gα protein transducin (Gt) by combining atomistic molecular dynamics (MD) simulations with Maxwell-Demod discrete MD (MDdMD) simulations of the receptor-catalyzed transition between GDP-bound and nucleotide-free states. Compared to the WT, Gt CAMs increase the overall fluctuations of nucleotide and its binding site. This is accompanied by weakening of native links involving GDP, α1, the G boxes, β1-β3, and α5. Collectively, constitutive activation by the considered mutants seems to associate with weakening of the interfaces between α5 and the surrounding portions and the interface between GTPase (G) and helical (H) domains. These mutational effects associate with increases in the overall fluctuations of the G and H domains, which reflect on the collective motions of the protein. Gt CAMs, with prominence to G56P, T325A, and F332A, prioritize collective motions of the H domain overlapping with the collective motions associated with receptor-catalyzed nucleotide release. In spite of different local perturbations, the mechanisms of nucleotide exchange catalyzed by activating mutations and by receptor are expected to employ similar molecular switches in the nucleotide binding site and to share the detachment of the H domain from the G domain.

Determinants of Macromolecular Specificity from Proteomics-Derived Peptide Substrate Data

BACKGROUND

Recent advances in proteomics methodologies allow for high throughput profiling of proteolytic cleavage events. The resulting substrate peptide distributions provide deep insights in the underlying macromolecular recognition events, as determinants of biomolecular specificity identified by proteomics approaches may be compared to structure-based analysis of corresponding protein-protein interfaces.

METHOD

Here, we present an overview of experimental and computational methodologies and tools applied in the area and provide an outlook beyond the protein class of proteases.

RESULTS AND CONCLUSION

We discuss here future potential, synergies and needs of the emerging overlap disciplines of proteomics and structure-based modelling.

Flexible Backbone Methods for Predicting and Designing Peptide Specificity

Protein-protein interactions play critical roles in essentially every cellular process. These interactions are often mediated by protein interaction domains that enable proteins to recognize their interaction partners, often by binding to short peptide motifs. For example, PDZ domains, which are among the most common protein interaction domains in the human proteome, recognize specific linear peptide sequences that are often at the C-terminus of other proteins. Determining the set of peptide sequences that a protein interaction domain binds, or it's "peptide specificity," is crucial for understanding its cellular function, and predicting how mutations impact peptide specificity is important for elucidating the mechanisms underlying human diseases. Moreover, engineering novel cellular functions for synthetic biology applications, such as the biosynthesis of biofuels or drugs, requires the design of protein interaction specificity to avoid crosstalk with native metabolic and signaling pathways. The ability to accurately predict and design protein-peptide interaction specificity is therefore critical for understanding and engineering biological function. One approach that has recently been employed toward accomplishing this goal is computational protein design. This chapter provides an overview of recent methodological advances in computational protein design and highlights examples of how these advances can enable increased accuracy in predicting and designing peptide specificity.

New open conformation of SMYD3 implicates conformational selection and allostery

SMYD3 plays a key role in cancer cell viability, adhesion, migration and invasion. SMYD3 promotes formation of inducible regulatory T cells and is involved in reducing autoimmunity. However, the nearly “closed” substrate-binding site and poor in vitro H3K4 methyltransferase activity have obscured further understanding of this oncogenically related protein. Here we reveal that SMYD3 can adopt an “open” conformation using molecular dynamics simulation and small-angle X-ray scattering. This ligand-binding-capable open state is related to the crystal structure-like closed state by a striking clamshell-like inter-lobe dynamics. The two states are characterized by many distinct structural and dynamical differences and the conformational transition pathway is mediated by a reversible twisting motion of the C-terminal domain (CTD). The spontaneous transition from the closed to open states suggests two possible, mutually non-exclusive models for SMYD3 functional regulation and the conformational selection mechanism and allostery may regulate the catalytic or ligand binding competence of SMYD3. This study provides an immediate clue to the puzzling role of SMYD3 in epigenetic gene regulation.

Structural and evolutionary analyses reveal determinants of DNA binding specificities of nucleoid-associated proteins HU and IHF

Nucleoid-associated proteins (NAPs) are chromosome-organizing factors, which affect the transcriptional landscape of a bacterial cell. HU is an NAP, which binds to DNA with a broad specificity while homologous IHF (Integration Host Factor), binds DNA with moderately higher specificity. Specificity and differential binding affinity of HU/IHF proteins towards their target binding sites play a crucial role in their regulatory dynamics. Decades of biochemical and genomic studies have been carried out for HU and IHF like proteins. Yet, questions related to their DNA binding specificity, and differential ability to bend DNA thus affecting the binding site length remained unanswered. In addition, the problem has not been investigated from an evolutionary perspective. Our phylogenetic analysis revealed three major clades belonging to HU, IHFα and IHFβ like proteins with reference to E. coli. We carried out a comparative analysis of three-dimensional structures of HU/IHF proteins to gain insight into the structural basis of clade division. The present study revealed three major features which contribute to differential DNA binding specificity of HU/IHF proteins, (I) conformational restriction of DNA binding residues due to salt-bridge formation, (II) the enrichment of alanine in the DNA binding site increasing conformational space of flexible side chains in its vicinity and (III) nature of DNA binding residue (Arg to Lys bias in different clades) which interacts differentially to DNA bases. We observed an extended electropositive surface at the DNA draping site for IHF clade proteins compared to HU, which stabilizes the DNA bend. Differences in the dimer stabilization strategies between HU and IHF were also observed. Our analysis reveals a comprehensive evolutionary picture, which rationalizes the origin of multi-specificity of HU/IHF proteins using sequence and structure-based determinants, which could also be applied to understand differences in binding specificities of other nucleic acid binding proteins.

Biosimilar structural comparability assessment by NMR: from small proteins to monoclonal antibodies

Biosimilar drug products must have a demonstrated similarity with respect to the reference product's molecules in order to ensure both the effectiveness of the drug and the patients' safety. In this paper the fusion framework of a highly sensitive NMR fingerprinting approach for conformational changes and mathematically-based biosimilarity metrics is introduced. The final goal is to translate the complex spectral information into biosimilarity scores, which are then used to estimate the degree of similarity between the biosimilar and the reference product. The proposed method was successfully applied to a small protein, i.e., filgrastim (neutropenia treatment), which is the first biosimilar approved in the United States, and a relatively large protein, i.e., monoclonal antibody rituximab (lymphoma treatment). This innovative approach introduces a new level of sensitivity to structural changes that are induced by, e.g., a small pH shift or other changes in the protein formulation.

Role of sequence evolution and conformational dynamics in the substrate specificity and oligomerization mode of thymidylate kinases

Thymidylate kinase (TMK) is a key enzyme for the synthesis of DNA, making it an important target for the development of anticancer, antibacterial, and antiparasitic drugs. TMK homologs exhibit significant variations in sequence, residue conformation, substrate specificity, and oligomerization mode. However, the influence of sequence evolution and conformational dynamics on its quaternary structure and function has not been studied before. Based on extensive sequence and structure analyses, our study detected several non-conserved residues which are linked by co-evolution and are implicated in the observed variations in flexibility, oligomeric assembly, and substrate specificity among the homologs. These lead to differences in the pattern of interactions at the active site in TMKs of different specificity. The method was further tested on TMK from Sulfolobus tokodaii (StTMK) which has substantial differences in sequence and structure compared to other TMKs. Our analyses pointed to a more flexible dTMP-binding site in StTMK compared to the other homologs. Binding assays proved that the protein can accommodate both purine and pyrimidine nucleotides at the dTMP binding site with comparable affinity. Additionally, the residues responsible for the narrow specificity of Brugia malayi TMK, whose three-dimensional structure is unavailable, were detected. Our study provides a residue-level understanding of the differences observed among TMK homologs in previous experiments. It also illustrates the correlation among sequence evolution, conformational dynamics, oligomerization mode, and substrate recognition in TMKs and detects co-evolving residues that affect binding, which should be taken into account while designing novel inhibitors.

Frizzled 7 and PIP2 binding by syntenin PDZ2 domain supports Frizzled 7 trafficking and signalling

PDZ domain-containing proteins work as intracellular scaffolds to control spatio-temporal aspects of cell signalling. This function is supported by the ability of their PDZ domains to bind other proteins such as receptors, but also phosphoinositide lipids important for membrane trafficking. Here we report a crystal structure of the syntenin PDZ tandem in complex with the carboxy-terminal fragment of Frizzled 7 and phosphatidylinositol 4,5-bisphosphate (PIP2). The crystal structure reveals a tripartite interaction formed via the second PDZ domain of syntenin. Biophysical and biochemical experiments establish co-operative binding of the tripartite complex and identify residues crucial for membrane PIP2-specific recognition. Experiments with cells support the importance of the syntenin-PIP2 interaction for plasma membrane targeting of Frizzled 7 and c-jun phosphorylation. This study contributes to our understanding of the biology of PDZ proteins as key players in membrane compartmentalization and dynamics.

Interdomain salt-bridges in the Ebola virus protein VP40 and their role in domain association and plasma membrane localization

The Ebola virus protein VP40 is a transformer protein that possesses an extraordinary ability to accomplish multiple functions by transforming into various oligomeric conformations. The disengagement of the C-terminal domain (CTD) from the N-terminal domain (NTD) is a crucial step in the conformational transformations of VP40 from the dimeric form to the hexameric form or octameric ring structure. Here, we use various molecular dynamics (MD) simulations to investigate the dynamics of the VP40 protein and the roles of interdomain interactions that are important for the domain-domain association and dissociation, and report on experimental results of the behavior of mutant variants of VP40. The MD studies find that various salt-bridge interactions modulate the VP40 domain dynamics by providing conformational specificity through interdomain interactions. The MD simulations reveal a novel salt-bridge between D45-K326 when the CTD participates in a latch-like interaction with the NTD. The D45-K326 salt-bridge interaction is proposed to help domain-domain association, whereas the E76-K291 interaction is important for stabilizing the closed-form structure. The effects of the removal of important VP40 salt-bridges on plasma membrane (PM) localization, VP40 oligomerization, and virus like particle (VLP) budding assays were investigated experimentally by live cell imaging using an EGFP-tagged VP40 system. It is found that the mutations K291E and D45K show enhanced PM localization but D45K significantly reduced VLP formation.

Promiscuity in the Enzymatic Catalysis of Phosphate and Sulfate Transfer

The enzymes that facilitate phosphate and sulfate hydrolysis are among the most proficient natural catalysts known to date. Interestingly, a large number of these enzymes are promiscuous catalysts that exhibit both phosphatase and sulfatase activities in the same active site and, on top of that, have also been demonstrated to efficiently catalyze the hydrolysis of other additional substrates with varying degrees of efficiency. Understanding the factors that underlie such multifunctionality is crucial both for understanding functional evolution in enzyme superfamilies and for the development of artificial enzymes. In this Current Topic, we have primarily focused on the structural and mechanistic basis for catalytic promiscuity among enzymes that facilitate both phosphoryl and sulfuryl transfer in the same active site, while comparing this to how catalytic promiscuity manifests in other promiscuous phosphatases. We have also drawn on the large number of experimental and computational studies of selected model systems in the literature to explore the different features driving the catalytic promiscuity of such enzymes. Finally, on the basis of this comparative analysis, we probe the plausible origins and determinants of catalytic promiscuity in enzymes that catalyze phosphoryl and sulfuryl transfer.

Beyond Membrane Protein Structure: Drug Discovery, Dynamics and Difficulties

Most of the previous content of this book has focused on obtaining the structures of membrane proteins. In this chapter we explore how those structures can be further used in two key ways. The first is their use in structure based drug design (SBDD) and the second is how they can be used to extend our understanding of their functional activity via the use of molecular dynamics. Both aspects now heavily rely on computations. This area is vast, and alas, too large to consider in depth in a single book chapter. Thus where appropriate we have referred the reader to recent reviews for deeper assessment of the field. We discuss progress via the use of examples from two main drug target areas; G-protein coupled receptors (GPCRs) and ion channels. We end with a discussion of some of the main challenges in the area.

Mechanisms of Differential Allosteric Modulation in Homologous Proteins: Insights from the Analysis of Internal Dynamics and Energetics of PDZ Domains

Allostery is a general phenomenon in proteins whereby a perturbation at one site reverberates into a functional change at another one, through modulation of its conformational dynamics. Herein, we address the problem of how the molecular signal encoded by a ligand is differentially transmitted through the structures of two homologous PDZ proteins: PDZ2, which responds to binding with structural and dynamical changes in regions distal from the ligand site, and PDZ3, which is characterized by less-intense dynamical variations. We use novel methods of analysis of MD simulations in the unbound and bound states to investigate the determinants of the differential allosteric behavior of the two proteins. The analysis of the correlations between the redistribution of stabilization energy and local fluctuation patterns highlights the nucleus of residues responsible for the stabilization of the 3D fold, the stability core, as the substructure that defines the difference in the allosteric response: in PDZ2, it undergoes a consistent dynamic and energetic reorganization, whereas in PDZ3, it remains largely unperturbed. Specifically, we observe for PDZ2 a significant anticorrelation between the motions of distal loops and residues of the stability core and differences in the correlation patterns between the bound and unbound states. Such variation is not observed in PDZ3, indicating that its energetics and internal dynamics are less affected by the presence/absence of the ligand. Finally, we propose a model with a direct link between the modulation of the structural, energetic and dynamic properties of a protein, and its allosteric response to a perturbation.

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.

Dynamics Govern Specificity of a Protein-Protein Interface: Substrate Recognition by Thrombin

Biomolecular recognition is crucial in cellular signal transduction. Signaling is mediated through molecular interactions at protein-protein interfaces. Still, specificity and promiscuity of protein-protein interfaces cannot be explained using simplistic static binding models. Our study rationalizes specificity of the prototypic protein-protein interface between thrombin and its peptide substrates relying solely on binding site dynamics derived from molecular dynamics simulations. We find conformational selection and thus dynamic contributions to be a key player in biomolecular recognition. Arising entropic contributions complement chemical intuition primarily reflecting enthalpic interaction patterns. The paradigm “dynamics govern specificity” might provide direct guidance for the identification of specific anchor points in biomolecular recognition processes and structure-based drug design.

Comparison of the Internal Dynamics of Metalloproteases Provides New Insights on Their Function and Evolution

Metalloproteases have evolved in a vast number of biological systems, being one of the most diverse types of proteases and presenting a wide range of folds and catalytic metal ions. Given the increasing understanding of protein internal dynamics and its role in enzyme function, we are interested in assessing how the structural heterogeneity of metalloproteases translates into their dynamics. Therefore, the dynamical profile of the clan MA type protein thermolysin, derived from an Elastic Network Model of protein structure, was evaluated against those obtained from a set of experimental structures and molecular dynamics simulation trajectories. A close correspondence was obtained between modes derived from the coarse-grained model and the subspace of functionally-relevant motions observed experimentally, the later being shown to be encoded in the internal dynamics of the protein. This prompted the use of dynamics-based comparison methods that employ such coarse-grained models in a representative set of clan members, allowing for its quantitative description in terms of structural and dynamical variability. Although members show structural similarity, they nonetheless present distinct dynamical profiles, with no apparent correlation between structural and dynamical relatedness. However, previously unnoticed dynamical similarity was found between the relevant members Carboxypeptidase Pfu, Leishmanolysin, and Botulinum Neurotoxin Type A, despite sharing no structural similarity. Inspection of the respective alignments shows that dynamical similarity has a functional basis, namely the need for maintaining proper intermolecular interactions with the respective substrates. These results suggest that distinct selective pressure mechanisms act on metalloproteases at structural and dynamical levels through the course of their evolution. This work shows how new insights on metalloprotease function and evolution can be assessed with comparison schemes that incorporate information on protein dynamics. The integration of these newly developed tools, if applied to other protein families, can lead to more accurate and descriptive protein classification systems.

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.