Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the ongoing global pandemic that has infected more than 31 million people in more than 180 countries worldwide. Like ...other coronaviruses, SARS-CoV-2 is thought to have been transmitted to humans from wild animals. Given the scale and widespread geographical distribution of the current pandemic and confirmed cases of cross-species transmission, the question of the extent to which this transmission is possible emerges, as well as what molecular features distinguish susceptible from non-susceptible animal species. Here, we investigated the structural properties of several ACE2 orthologs bound to the SARS-CoV-2 spike protein. We found that species known not to be susceptible to SARS-CoV-2 infection have non-conservative mutations in several ACE2 amino acid residues that disrupt key polar and charged contacts with the viral spike protein. Our models also allow us to predict affinity-enhancing mutations that could be used to design ACE2 variants for therapeutic purposes. Finally, our study provides a blueprint for modeling viral-host protein interactions and highlights several important considerations when designing these computational studies and analyzing their results.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The pdb-tools are a collection of Python scripts for working with molecular structure data in the Protein Data Bank (PDB) format. They allow users to edit, convert, and validate PDB files, from the ...command-line, in a simple but efficient manner. The pdb-tools are implemented in Python, without any external dependencies, and are freely available under the open-source Apache License at https://github.com/haddocking/pdb-tools/ and on
PyPI.
Protein-protein interactions are fundamental to many biological processes. Experimental screens have identified tens of thousands of interactions, and structural biology has provided detailed ...functional insight for select 3D protein complexes. An alternative rich source of information about protein interactions is the evolutionary sequence record. Building on earlier work, we show that analysis of correlated evolutionary sequence changes across proteins identifies residues that are close in space with sufficient accuracy to determine the three-dimensional structure of the protein complexes. We evaluate prediction performance in blinded tests on 76 complexes of known 3D structure, predict protein-protein contacts in 32 complexes of unknown structure, and demonstrate how evolutionary couplings can be used to distinguish between interacting and non-interacting protein pairs in a large complex. With the current growth of sequences, we expect that the method can be generalized to genome-wide elucidation of protein-protein interaction networks and used for interaction predictions at residue resolution.
Salmeterol is a partial agonist for the β2 adrenergic receptor (β2AR) and the first long-acting β2AR agonist to be widely used clinically for the treatment of asthma and chronic obstructive pulmonary ...disease. Salmeterol’s safety and mechanism of action have both been controversial. To understand its unusual pharmacological action and partial agonism, we obtained the crystal structure of salmeterol-bound β2AR in complex with an active-state-stabilizing nanobody. The structure reveals the location of the salmeterol exosite, where sequence differences between β1AR and β2AR explain the high receptor-subtype selectivity. A structural comparison with the β2AR bound to the full agonist epinephrine reveals differences in the hydrogen-bond network involving residues Ser2045.43 and Asn2936.55. Mutagenesis and biophysical studies suggested that these interactions lead to a distinct active-state conformation that is responsible for the partial efficacy of G-protein activation and the limited β-arrestin recruitment for salmeterol.
The Protein Data Bank (PDB) file format remains a popular format used and supported by many software to represent coordinates of macromolecular structures. It however suffers from drawbacks such as ...error‐prone manual editing. Because of that, various software toolkits have been developed to facilitate its editing and manipulation, but, to date, there is no online tool available for this purpose. Here we present PDB‐Tools Web, a flexible online service for manipulating PDB files. It offers a rich and user‐friendly graphical user interface that allows users to mix‐and‐match more than 40 individual tools from the pdb‐tools suite. Those can be combined in a few clicks to perform complex pipelines, which can be saved and uploaded. The resulting processed PDB files can be visualized online and downloaded. The web server is freely available at https://wenmr.science.uu.nl/pdbtools.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
The inflammatory chemokine CCL5, which binds the chemokine receptor CCR5 in a two‐step mechanism so as to activate signaling pathways in hematopoetic cells, plays an important role in immune ...surveillance, inflammation, and development as well as in several immune system pathologies. The recently published crystal structure of CCR5 bound to a high‐affinity variant of CCL5 lacks the N‐terminal segment of the receptor that is post‐translationally sulfated and is known to be important for high‐affinity binding. Here, we report the NMR solution structure of monomeric CCL5 bound to a synthetic doubly sulfated peptide corresponding to the missing first 27 residues of CCR5. Our structures show that two sulfated tyrosine residues, sY10 and sY14, as well as the unsulfated Y15 form a network of strong interactions with a groove on a surface of CCL5 that is formed from evolutionarily conserved basic and hydrophobic amino acids. We then use our NMR structures, in combination with available crystal data, to create an atomic model of full‐length wild‐type CCR5:CCL5. Our findings reveal the structural determinants involved in the recognition of CCL5 by the CCR5 N terminus. These findings, together with existing structural data, provide a complete structural framework with which to understand the specificity of receptor:chemokine interactions.
Database
Structural data are available in the PDB under the accession number 6FGP
The inflammatory chemokine CCL5 binds the chemokine receptor CCR5 in hematopoietic cells, activating signalling pathways with important roles in immune surveillance, inflammation, development and several immune system pathologies. In this work, Meital Abayev and colleagues report the NMR solution structure of a doubly sulfated N‐terminal peptide of CCR5 bound to CCL5. Their findings reveal that two sulfated tyrosines on CCR5 interact with evolutionarily conserved residues on the chemokine. Bioinformatics analyses suggest that such contacts may be conserved across the chemokine family, explaining both the specificity and affinity of receptor:chemokine interactions.
Full text
Available for:
BFBNIB, DOBA, FZAB, GIS, IJS, IZUM, KILJ, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBMB, SIK, UILJ, UKNU, UL, UM, UPUK
Predicting the 3D structure of protein interactions remains a challenge in the field of computational structural biology. This is in part due to difficulties in sampling the complex energy landscape ...of multiple interacting flexible polypeptide chains. Coarse-graining approaches, which reduce the number of degrees of freedom of the system, help address this limitation by smoothing the energy landscape, allowing an easier identification of the global energy minimum. They also accelerate the calculations, allowing for modeling larger assemblies. Here, we present the implementation of the MARTINI coarse-grained force field for proteins into HADDOCK, our integrative modeling platform. Docking and refinement are performed at the coarse-grained level, and the resulting models are then converted back to atomistic resolution through a distance restraints-guided morphing procedure. Our protocol, tested on the largest complexes of the protein docking benchmark 5, shows an overall ∼7-fold speed increase compared to standard all-atom calculations, while maintaining a similar accuracy and yielding substantially more near-native solutions. To showcase the potential of our method, we performed simultaneous 7 body docking to model the 1:6 KaiC-KaiB complex, integrating mutagenesis and hydrogen/deuterium exchange data from mass spectrometry with symmetry restraints, and validated the resulting models against a recently published cryo-EM structure.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
Free energy perturbation (FEP) remains an indispensable method for computationally assaying prospective compounds in advance of synthesis. However, before FEP can be deployed prospectively, it must ...demonstrate retrospective recapitulation of known experimental data where the subtle details of the atomic ligand–receptor model are consequential. An open question is whether AlphaFold models can serve as useful initial models for FEP in the regime where there exists a congeneric series of known chemical matter but where no experimental structures are available either of the target or of close homologues. As AlphaFold structures are provided without a bound ligand, we employ induced fit docking to refine the AlphaFold models in the presence of one or more congeneric ligands. In this work, we first validate the performance of our latest induced fit docking technology, IFD-MD, on a retrospective set of public experimental GPCR structures with 95% of cross-docks producing a pose with a ligand RMSD ≤ 2.5 Å in the top two predictions. We then apply IFD-MD and FEP on AlphaFold models of the somatostatin receptor family of GPCRs. We use AlphaFold models produced prior to the availability of any experimental structure from this family. We arrive at FEP-validated models for SSTR2, SSTR4, and SSTR5, with RMSE around 1 kcal/mol, and explore the challenges of model validation under scenarios of limited ligand data, ample ligand data, and categorical data.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
The KoBaMIN web server provides an online interface to a simple, consistent and computationally efficient protein structure refinement protocol based on minimization of a knowledge-based potential of ...mean force. The server can be used to refine either a single protein structure or an ensemble of proteins starting from their unrefined coordinates in PDB format. The refinement method is particularly fast and accurate due to the underlying knowledge-based potential derived from structures deposited in the PDB; as such, the energy function implicitly includes the effects of solvent and the crystal environment. Our server allows for an optional but recommended step that optimizes stereochemistry using the MESHI software. The KoBaMIN server also allows comparison of the refined structures with a provided reference structure to assess the changes brought about by the refinement protocol. The performance of KoBaMIN has been benchmarked widely on a large set of decoys, all models generated at the seventh worldwide experiments on critical assessment of techniques for protein structure prediction (CASP7) and it was also shown to produce top-ranking predictions in the refinement category at both CASP8 and CASP9, yielding consistently good results across a broad range of model quality values. The web server is fully functional and freely available at http://csb.stanford.edu/kobamin.