Enabling structure-based drug discovery utilizing predicted models.

High-quality predicted structures enable structure-based approaches to an expanding number of drug discovery programs. We propose that by utilizing free energy perturbation (FEP), predicted structures can be confidently employed to achieve drug design goals. We use structure-based modeling of hERG inhibition to illustrate this value of FEP.

Using AlphaFold and Experimental Structures for the Prediction of the Structure and Binding Affinities of GPCR Complexes via Induced Fit Docking and Free Energy Perturbation.

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.

An Electrostatically-steered Conformational Selection Mechanism Promotes SARS-CoV-2 Spike Protein Variation.

After two years since the outbreak, the COVID-19 pandemic remains a global public health emergency. SARS-CoV-2 variants with substitutions on the spike (S) protein emerge increasing the risk of immune evasion and cross-species transmission. Here, we analyzed the evolution of the S protein as recorded in 276,712 samples collected before the start of vaccination efforts. Our analysis shows that most variants destabilize the S protein trimer, increase its conformational heterogeneity and improve the odds of the recognition by the host cell receptor. Most frequent substitutions promote overall hydrophobicity by replacing charged amino acids, reducing stabilizing local interactions in the unbound S protein trimer. Moreover, our results identify "forbidden" regions that rarely show any sequence variation, and which are related to conformational changes occurring upon fusion. These results are significant for understanding the structure and function of SARS-CoV-2 related proteins which is a critical step in vaccine development and epidemiological surveillance.

Predicting the zoonotic capacity of mammals to transmit SARS-CoV-2.

Back and forth transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) between humans and animals will establish wild reservoirs of virus that endanger long-term efforts to control COVID-19 in people and to protect vulnerable animal populations. Better targeting surveillance and laboratory experiments to validate zoonotic potential requires predicting high-risk host species. A major bottleneck to this effort is the few species with available sequences for angiotensin-converting enzyme 2 receptor, a key receptor required for viral cell entry. We overcome this bottleneck by combining species' ecological and biological traits with three-dimensional modelling of host-virus protein-protein interactions using machine learning. This approach enables predictions about the zoonotic capacity of SARS-CoV-2 for greater than 5000 mammals-an order of magnitude more species than previously possible. Our predictions are strongly corroborated by in vivo studies. The predicted zoonotic capacity and proximity to humans suggest enhanced transmission risk from several common mammals, and priority areas of geographic overlap between these species and global COVID-19 hotspots. With molecular data available for only a small fraction of potential animal hosts, linking data across biological scales offers a conceptual advance that may expand our predictive modelling capacity for zoonotic viruses with similarly unknown host ranges.

Predicting the Specificity- Determining Positions of Receptor Tyrosine Kinase Axl.

Owing to its clinical significance, modulation of functionally relevant amino acids in protein-protein complexes has attracted a great deal of attention. To this end, many approaches have been proposed to predict the partner-selecting amino acid positions in evolutionarily close complexes. These approaches can be grouped into sequence-based machine learning and structure-based energy-driven methods. In this work, we assessed these methods' ability to map the specificity-determining positions of Axl, a receptor tyrosine kinase involved in cancer progression and immune system diseases. For sequence-based predictions, we used SDPpred, Multi-RELIEF, and Sequence Harmony. For structure-based predictions, we utilized HADDOCK refinement and molecular dynamics simulations. As a result, we observed that (i) sequence-based methods overpredict partner-selecting residues of Axl and that (ii) combining Multi-RELIEF with HADDOCK-based predictions provides the key Axl residues, covered by the extensive molecular dynamics simulations. Expanding on these results, we propose that a sequence-structure-based approach is necessary to determine specificity-determining positions of Axl, which can guide the development of therapeutic molecules to combat Axl misregulation.

Predicting the zoonotic capacity of mammals to transmit SARS-CoV-2.

Back and forth transmission of SARS-CoV-2 between humans and animals may lead to wild reservoirs of virus that can endanger efforts toward long-term control of COVID-19 in people, and protecting vulnerable animal populations that are particularly susceptible to lethal disease. Predicting high risk host species is key to targeting field surveillance and lab experiments that validate host zoonotic potential. A major bottleneck to predicting animal hosts is the small number of species with available molecular information about the structure of ACE2, a key cellular receptor required for viral cell entry. We overcome this bottleneck by combining species' ecological and biological traits with 3D modeling of virus and host cell protein interactions using machine learning methods. This approach enables predictions about the zoonotic capacity of SARS-CoV-2 for over 5,000 mammals - an order of magnitude more species than previously possible. The high accuracy predictions achieved by this approach are strongly corroborated by in vivo empirical studies. We identify numerous common mammal species whose predicted zoonotic capacity and close proximity to humans may further enhance the risk of spillover and spillback transmission of SARS-CoV-2. Our results reveal high priority areas of geographic overlap between global COVID-19 hotspots and potential new mammal hosts of SARS-CoV-2. With molecular sequence data available for only a small fraction of potential host species, predictive modeling integrating data across multiple biological scales offers a conceptual advance that may expand our predictive capacity for zoonotic viruses with similarly unknown and potentially broad host ranges.

PDB-tools web: A user-friendly interface for the manipulation of PDB files.

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.

Insights on cross-species transmission of SARS-CoV-2 from structural modeling.

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.

Structural models of human ACE2 variants with SARS-CoV-2 Spike protein for structure-based drug design.

Emergence of coronaviruses poses a threat to global health and economy. The current outbreak of SARS-CoV-2 has infected more than 28,000,000 people and killed more than 915,000. To date, there is no treatment for coronavirus infections, making the development of therapies to prevent future epidemics of paramount importance. To this end, we collected information regarding naturally-occurring variants of the Angiotensin-converting enzyme 2 (ACE2), an epithelial receptor that both SARS-CoV and SARS-CoV-2 use to enter the host cells. We built 242 structural models of variants of human ACE2 bound to the receptor binding domain (RBD) of the SARS-CoV-2 surface spike glycoprotein (S protein) and refined their interfaces with HADDOCK. Our dataset includes 140 variants of human ACE2 representing missense mutations found in genome-wide studies, 39 mutants with reported effects on the recognition of the RBD, and 63 predictions after computational alanine scanning mutagenesis of ACE2-RBD interface residues. This dataset will help accelerate the design of therapeutics against SARS-CoV-2, as well as contribute to prevention of possible future coronaviruses outbreaks.

Bacterial flagellar motor PL-ring disassembly subcomplexes are widespread and ancient.

The bacterial flagellum is an amazing nanomachine. Understanding how such complex structures arose is crucial to our understanding of cellular evolution. We and others recently reported that in several Gammaproteobacterial species, a relic subcomplex comprising the decorated P and L rings persists in the outer membrane after flagellum disassembly. Imaging nine additional species with cryo-electron tomography, here, we show that this subcomplex persists after flagellum disassembly in other phyla as well. Bioinformatic analyses fail to show evidence of any recent horizontal transfers of the P- and L-ring genes, suggesting that this subcomplex and its persistence is an ancient and conserved feature of the flagellar motor. We hypothesize that one function of the P and L rings is to seal the outer membrane after motor disassembly.

Less Is More: Coarse-Grained Integrative Modeling of Large Biomolecular Assemblies with HADDOCK.

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.

Publisher Correction: Structural insights into binding specificity, efficacy and bias of a ß2AR partial agonist.

In the version of this paper originally published, the structure for epinephrine shown in Figure 1a was redrawn with an extra carbon. The structure has been replaced in the HTML and PDF versions of the article. The original and corrected versions of the structure are shown below.

Structural insights into binding specificity, efficacy and bias of a ß2AR partial agonist.

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 solution structure of monomeric CCL5 in complex with a doubly sulfated N-terminal segment of CCR5.

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.

pdb-tools: a swiss army knife for molecular structures.

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.

Augmenting Research, Education, and Outreach with Client-Side Web Programming.

The evolution of computing and web technologies over the past decade has enabled the development of fully fledged scientific applications that run directly on web browsers. Powered by JavaScript, the lingua franca of web programming, these 'web apps' are starting to revolutionize and democratize scientific research, education, and outreach.

M3: an integrative framework for structure determination of molecular machines.

We present a broadly applicable, user-friendly protocol that incorporates sparse and hybrid experimental data to calculate quasi-atomic-resolution structures of molecular machines. The protocol uses the HADDOCK framework, accounts for extensive structural rearrangements both at the domain and atomic levels and accepts input from all structural and biochemical experiments whose data can be translated into interatomic distances and/or molecular shapes.

Supramolecular Organization and Functional Implications of K+ †Channel Clusters in Membranes.

The segregation of cellular surfaces in heterogeneous patches is considered to be a common motif in bacteria and eukaryotes that is underpinned by the observation of clustering and cooperative gating of signaling membrane proteins such as receptors or channels. Such processes could represent an important cellular strategy to shape signaling activity. Hence, structural knowledge of the arrangement of channels or receptors in supramolecular assemblies represents a crucial step towards a better understanding of signaling across membranes. We herein report on the supramolecular organization of clusters of the K+ channel KcsA in bacterial membranes, which was analyzed by a combination of DNP-enhanced solid-state NMR experiments and MD simulations. We used solid-state NMR spectroscopy to determine the channel-channel interface and to demonstrate the strong correlation between channel function and clustering, which suggests a yet unknown mechanism of communication between K+ channels.

Structural and Functional Analysis of a ß2-Adrenergic Receptor Complex with GRK5.

The phosphorylation of agonist-occupied G-protein-coupled receptors (GPCRs) by GPCR kinases (GRKs) functions to turn off G-protein signaling and turn on arrestin-mediated signaling. While a structural understanding of GPCR/G-protein and GPCR/arrestin complexes has emerged in recent years, the molecular architecture of a GPCR/GRK complex remains poorly defined. We used a comprehensive integrated approach of cross-linking, hydrogen-deuterium exchange mass spectrometry (MS), electron microscopy, mutagenesis, molecular dynamics simulations, and computational docking to analyze GRK5 interaction with the ß2-adrenergic receptor (ß2AR). These studies revealed a dynamic mechanism of complex formation that involves large conformational changes in the GRK5 RH/catalytic domain interface upon receptor binding. These changes facilitate contacts between intracellular loops 2 and 3 and the C terminus of the ß2AR with the GRK5 RH bundle subdomain, membrane-binding surface, and kinase catalytic cleft, respectively. These studies significantly contribute to our understanding of the mechanism by which GRKs regulate the function of activated GPCRs. PAPERCLIP.