1
|
Jones SJ, Perez A. Molecular Modeling of Self-Assembling Peptides. ACS APPLIED BIO MATERIALS 2024; 7:543-552. [PMID: 36795608 DOI: 10.1021/acsabm.2c00921] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
Peptide epitopes mediate as many as 40% of protein-protein interactions and fulfill signaling, inhibition, and activation roles within the cell. Beyond protein recognition, some peptides can self- or coassemble into stable hydrogels, making them a readily available source of biomaterials. While these 3D assemblies are routinely characterized at the fiber level, there are missing atomistic details about the assembly scaffold. Such atomistic detail can be useful in the rational design of more stable scaffold structures and with improved accessibility to functional motifs. Computational approaches can in principle reduce the experimental cost of such an endeavor by predicting the assembly scaffold and identifying novel sequences that adopt said structure. Yet, inaccuracies in physical models and inefficient sampling have limited atomistic studies to short (two or three amino acid) peptides. Given recent developments in machine learning and advances in sampling strategies, we revisit the suitability of physical models for this task. We use the MELD (Modeling Employing Limited Data) approach to drive self-assembly in combination with generic data in cases where conventional MD is unsuccessful. Finally, despite recent developments in machine learning algorithms for protein structure and sequence predictions, we find the algorithms are not yet suited for studying the assembly of short peptides.
Collapse
Affiliation(s)
- Stephen J Jones
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| | - Alberto Perez
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| |
Collapse
|
2
|
Caparotta M, Perez A. When MELD Meets GaMD: Accelerating Biomolecular Landscape Exploration. J Chem Theory Comput 2023; 19:8743-8750. [PMID: 38039424 DOI: 10.1021/acs.jctc.3c01019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2023]
Abstract
We introduce Gaussian accelerated MELD (GaMELD) as a new method for exploring the energy landscape of biomolecules. GaMELD combines the strengths of Gaussian accelerated molecular dynamics (GaMD) and modeling employing limited data (MELD) to navigate complex energy landscapes. MELD uses replica-exchange molecular simulations to integrate limited and uncertain data into simulations via Bayesian inference. MELD has been successfully applied to problems of structure prediction like protein folding and complex structure prediction. However, the computational cost for MELD simulations has limited its broader applicability. The synergy of GaMD and MELD surmounts this limitation efficiently sampling the energy landscape at a lower computational cost (reducing the computational cost by a factor of 2 to six). Effectively, GaMD is used to shift energy distributions along replicas to increase the overlap in energy distributions across replicas, facilitating a random walk in replica space. We tested GaMELD on a benchmark set of 12 small proteins that have been previously studied through MELD and conventional MD. GaMELD consistently achieves accurate predictions with fewer replicas. By increasing the efficacy of replica exchange, GaMELD effectively accelerates convergence in the conformational space, enabling improved sampling across a diverse set of systems.
Collapse
Affiliation(s)
- Marcelo Caparotta
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| | - Alberto Perez
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| |
Collapse
|
3
|
Parui S, Brini E, Dill KA. Computing Free Energies of Fold-Switching Proteins Using MELD x MD. J Chem Theory Comput 2023; 19:6839-6847. [PMID: 37725050 DOI: 10.1021/acs.jctc.3c00679] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/21/2023]
Abstract
Some proteins are conformational switches, able to transition between relatively different conformations. To understand what drives them requires computing the free-energy difference ΔGAB between their stable states, A and B. Molecular dynamics (MD) simulations alone are often slow because they require a reaction coordinate and must sample many transitions in between. Here, we show that modeling employing limited data (MELD) x MD on known endstates A and B is accurate and efficient because it does not require passing over barriers or knowing reaction coordinates. We validate this method on two problems: (1) it gives correct relative populations of α and β conformers for small designed chameleon sequences of protein G; and (2) it correctly predicts the conformations of the C-terminal domain (CTD) of RfaH. Free-energy methods like MELD x MD can often resolve structures that confuse machine-learning (ML) methods.
Collapse
Affiliation(s)
- Sridip Parui
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - Emiliano Brini
- School of Chemistry and Materials Science, 85 Lomb Memorial Drive, Rochester, New York 14623, United States
| | - Ken A Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, United States
| |
Collapse
|
4
|
Parui S, Robertson JC, Somani S, Tresadern G, Liu C, Dill KA. MELD-Bracket Ranks Binding Affinities of Diverse Sets of Ligands. J Chem Inf Model 2023; 63:2857-2865. [PMID: 37093848 DOI: 10.1021/acs.jcim.3c00243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
Affinity ranking of structurally diverse small-molecule ligands is a challenging problem with important applications in structure-based drug discovery. Absolute binding free energy methods can model diverse ligands, but the high computational cost of the current methods limits application to data sets with few ligands. We recently developed MELD-Bracket, a Molecular Dynamics method for efficient affinity ranking of ligands [ JCTC 2022, 18 (1), 374-379]. It utilizes a Bayesian framework to guide sampling to relevant regions of phase space, and it couples this with a bracket-like competition on a pool of ligands. Here we find that 6-competitor MELD-Bracket can rank dozens of diverse ligands that have low structural similarity and different net charges. We benchmark it on four protein systems─PTB1B, Tyk2, BACE, and JAK3─having varied modes of interactions. We also validated 8-competitor and 12-competitor protocols. The MELD-Bracket protocols presented here may have the appropriate balance of accuracy and computational efficiency to be suitable for ranking diverse ligands from typical drug discovery campaigns.
Collapse
Affiliation(s)
- Sridip Parui
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - James C Robertson
- Janssen Research and Development, Spring House, Pennsylvania 19477, United States
| | - Sandeep Somani
- Janssen Research and Development, Spring House, Pennsylvania 19477, United States
| | - Gary Tresadern
- Janssen Research and Development, Turnhoutseweg 30, Beerse B-2340, Belgium
| | - Cong Liu
- Center for the Development of Therapeutics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts 02142, United States
| | - Ken A Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, United States
| |
Collapse
|
5
|
Liu Q, Perez A. Assessing a computational pipeline to identify binding motifs to the α2 β1 integrin. Front Chem 2023; 11:1107400. [PMID: 36860646 PMCID: PMC9968975 DOI: 10.3389/fchem.2023.1107400] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2022] [Accepted: 01/27/2023] [Indexed: 02/16/2023] Open
Abstract
Integrins in the cell surface interact with functional motifs found in the extracellular matrix (ECM) that queue the cell for biological actions such as migration, adhesion, or growth. Multiple fibrous proteins such as collagen or fibronectin compose the ECM. The field of biomechanical engineering often deals with the design of biomaterials compatible with the ECM that will trigger cellular response (e.g., in tissue regeneration). However, there are a relative few number of known integrin binding motifs compared to all the possible peptide epitope sequences available. Computational tools could help identify novel motifs, but have been limited by the challenges in modeling the binding to integrin domains. We revisit a series of traditional and novel computational tools to assess their performance in identifying novel binding motifs for the I-domain of the α2β1 integrin.
Collapse
Affiliation(s)
| | - Alberto Perez
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, FL, United States
| |
Collapse
|
6
|
Kurgan N, Baranowski B, Stoikos J, MacNeil AJ, Fajardo VA, MacPherson REK, Klentrou P. Characterization of sclerostin's response within white adipose tissue to an obesogenic diet at rest and in response to acute exercise in male mice. Front Physiol 2023; 13:1061715. [PMID: 36685192 PMCID: PMC9846496 DOI: 10.3389/fphys.2022.1061715] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2022] [Accepted: 12/14/2022] [Indexed: 01/05/2023] Open
Abstract
Introduction: It is well established that sclerostin antagonizes the anabolic Wnt signalling pathway in bone, however, its physiological role in other tissues remains less clear. This study examined the effect of a high-fat diet (HFD) on sclerostin content and downstream markers of the Wnt signaling pathway (GSK3β and β-catenin) within subcutaneous inguinal white adipose tissue (iWAT), and visceral epididymal WAT (eWAT) depots at rest and in response to acute aerobic exercise. Methods: Male C57BL/6 mice (n = 40, 18 weeks of age) underwent 10 weeks of either a low-fat diet (LFD) or HFD. Within each diet group, mice were assigned to either remain sedentary (SED) or perform 2 h of endurance treadmill exercise at 15 m min-1 with 5° incline (EX), creating four groups: LFD + SED (N = 10), LFD + EX (N = 10), HFD + SED (N = 10), and HFD + EX (N = 10). Serum and WAT depots were collected 2 h post-exercise. Results: Serum sclerostin showed a diet-by-exercise interaction, reflecting HFD + EX mice having higher concentration than HFD + SED (+31%, p = 0.03), and LFD mice being unresponsive to exercise. iWAT sclerostin content decreased post-exercise in both 28 kDa (-31%, p = 0.04) and 30 kDa bands (-36%, main effect for exercise, p = 0.02). iWAT β-catenin (+44%, p = 0.03) and GSK3β content were higher in HFD mice compared to LFD (+128%, main effect for diet, p = 0.005). Monomeric sclerostin content was abolished in eWAT of HFD mice (-96%, main effect for diet, p < 0.0001), was only detectable as a 30 kDa band in LFD mice and was unresponsive to exercise. β-catenin and GSK3β were both unresponsive to diet and exercise within eWAT. Conclusion: These results characterized sclerostin's content to WAT depots in response to acute exercise, which appears to be specific to a reduction in iWAT and identified a differential regulation of sclerostin's form/post-translational modifications depending on diet and WAT depot.
Collapse
Affiliation(s)
- Nigel Kurgan
- Department of Kinesiology, Brock University, St. Catharines, ON, Canada,Centre for Bone and Muscle Health, Brock University, St. Catharines, ON, Canada
| | - Bradley Baranowski
- Department of Health Sciences, Brock University, St. Catharines, ON, Canada
| | - Joshua Stoikos
- Department of Kinesiology, Brock University, St. Catharines, ON, Canada,Centre for Bone and Muscle Health, Brock University, St. Catharines, ON, Canada
| | - Adam J. MacNeil
- Department of Health Sciences, Brock University, St. Catharines, ON, Canada
| | - Val A. Fajardo
- Department of Kinesiology, Brock University, St. Catharines, ON, Canada,Centre for Bone and Muscle Health, Brock University, St. Catharines, ON, Canada
| | | | - Panagiota Klentrou
- Department of Kinesiology, Brock University, St. Catharines, ON, Canada,Centre for Bone and Muscle Health, Brock University, St. Catharines, ON, Canada,*Correspondence: Panagiota Klentrou,
| |
Collapse
|
7
|
Lam K, Kasavajhala K, Gunasekera S, Simmerling C. Accelerating the Ensemble Convergence of RNA Hairpin Simulations with a Replica Exchange Structure Reservoir. J Chem Theory Comput 2022; 18:3930-3947. [PMID: 35502992 PMCID: PMC10658646 DOI: 10.1021/acs.jctc.2c00065] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
RNA is a key participant in many biological processes, but studies of RNA using computer simulations lag behind those of proteins, largely due to less-developed force fields and the slow dynamics of RNA. Generating converged RNA ensembles for force field development and other studies remains a challenge. In this study, we explore the ability of replica exchange molecular dynamics to obtain well-converged conformational ensembles for two RNA hairpin systems in an implicit solvent. Even for these small model systems, standard REMD remains computationally costly, but coupling to a pre-generated structure library using the reservoir REMD approach provides a dramatic acceleration of ensemble convergence for both model systems. Such precise ensembles could facilitate RNA force field development and validation and applications of simulation to more complex RNA systems. The advantages and remaining challenges of applying R-REMD to RNA are investigated in detail.
Collapse
Affiliation(s)
- Kenneth Lam
- Molecular and Cellular Biology, Stony Brook University, Stony Brook, New York 11794, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - Koushik Kasavajhala
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - Sarah Gunasekera
- Program in Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - Carlos Simmerling
- Molecular and Cellular Biology, Stony Brook University, Stony Brook, New York 11794, United States
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| |
Collapse
|
8
|
Esmaeeli R, Andal B, Perez A. Searching for Low Probability Opening Events in a DNA Sliding Clamp. Life (Basel) 2022; 12:life12020261. [PMID: 35207548 PMCID: PMC8876151 DOI: 10.3390/life12020261] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Revised: 02/01/2022] [Accepted: 02/03/2022] [Indexed: 11/27/2022] Open
Abstract
The β subunit of E. coli DNA polymererase III is a DNA sliding clamp associated with increasing the processivity of DNA synthesis. In its free form, it is a circular homodimer structure that can accomodate double-stranded DNA in a nonspecific manner. An open state of the clamp must be accessible before loading the DNA. The opening mechanism is still a matter of debate, as is the effect of bound DNA on opening/closing kinetics. We use a combination of atomistic, coarse-grained, and enhanced sampling strategies in both explicit and implicit solvents to identify opening events in the sliding clamp. Such simulations of large nucleic acid and their complexes are becoming available and are being driven by improvements in force fields and the creation of faster computers. Different models support alternative opening mechanisms, either through an in-plane or out-of-plane opening event. We further note some of the current limitations, despite advances, in modeling these highly charged systems with implicit solvent.
Collapse
|
9
|
Chang L, Perez A, Miranda-Quintana RA. Improving the analysis of biological ensembles through extended similarity measures. Phys Chem Chem Phys 2021; 24:444-451. [PMID: 34897334 DOI: 10.1039/d1cp04019g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
We present new algorithms to classify structural ensembles of macromolecules based on the recently proposed extended similarity measures. Molecular dynamics provides a wealth of structural information on systems of biological interest. As computer power increases, we capture larger ensembles and larger conformational transitions between states. Typically, structural clustering provides the statistical mechanics treatment of the system to identify relevant biological states. The key advantage of our approach is that the newly introduced extended similarity indices reduce the computational complexity of assessing the similarity of a set of structures from O(N2) to O(N). Here we take advantage of this favorable cost to develop several highly efficient techniques, including a linear-scaling algorithm to determine the medoid of a set (which we effectively use to select the most representative structure of a cluster). Moreover, we use our extended similarity indices as a linkage criterion in a novel hierarchical agglomerative clustering algorithm. We apply these new metrics to analyze the ensembles of several systems of biological interest such as folding and binding of macromolecules (peptide, protein, DNA-protein). In particular, we design a new workflow that is capable of identifying the most important conformations contributing to the protein folding process. We show excellent performance in the resulting clusters (surpassing traditional linkage criteria), along with faster performance and an efficient cost-function to identify when to merge clusters.
Collapse
Affiliation(s)
- Liwei Chang
- Department of Chemistry, University of Florida, Gainesville, FL, 32611, USA.
| | - Alberto Perez
- Department of Chemistry, University of Florida, Gainesville, FL, 32611, USA. .,Quantum Theory Project, University of Florida, Gainesville, FL, 32611, USA
| | - Ramón Alain Miranda-Quintana
- Department of Chemistry, University of Florida, Gainesville, FL, 32611, USA. .,Quantum Theory Project, University of Florida, Gainesville, FL, 32611, USA
| |
Collapse
|
10
|
Sharma B, Dill KA. MELD-accelerated molecular dynamics help determine amyloid fibril structures. Commun Biol 2021; 4:942. [PMID: 34354239 PMCID: PMC8342454 DOI: 10.1038/s42003-021-02461-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Accepted: 07/15/2021] [Indexed: 02/07/2023] Open
Abstract
It is challenging to determine the structures of protein fibrils such as amyloids. In principle, Molecular Dynamics (MD) modeling can aid experiments, but normal MD has been impractical for these large multi-molecules. Here, we show that MELD accelerated MD (MELD x MD) can give amyloid structures from limited data. Five long-chain fibril structures are accurately predicted from NMR and Solid State NMR (SSNMR) data. Ten short-chain fibril structures are accurately predicted from more limited restraints information derived from the knowledge of strand directions. Although the present study only tests against structure predictions - which are the most detailed form of validation currently available - the main promise of this physical approach is ultimately in going beyond structures to also give mechanical properties, conformational ensembles, and relative stabilities.
Collapse
Affiliation(s)
- Bhanita Sharma
- grid.36425.360000 0001 2216 9681Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY USA
| | - Ken A. Dill
- grid.36425.360000 0001 2216 9681Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY USA ,grid.36425.360000 0001 2216 9681Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY USA ,grid.36425.360000 0001 2216 9681Departments of Chemistry and Physics, Stony Brook University, Stony Brook, NY USA
| |
Collapse
|
11
|
Identifying hydrophobic protein patches to inform protein interaction interfaces. Proc Natl Acad Sci U S A 2021; 118:2018234118. [PMID: 33526682 DOI: 10.1073/pnas.2018234118] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Interactions between proteins lie at the heart of numerous biological processes and are essential for the proper functioning of the cell. Although the importance of hydrophobic residues in driving protein interactions is universally accepted, a characterization of protein hydrophobicity, which informs its interactions, has remained elusive. The challenge lies in capturing the collective response of the protein hydration waters to the nanoscale chemical and topographical protein patterns, which determine protein hydrophobicity. To address this challenge, here, we employ specialized molecular simulations wherein water molecules are systematically displaced from the protein hydration shell; by identifying protein regions that relinquish their waters more readily than others, we are then able to uncover the most hydrophobic protein patches. Surprisingly, such patches contain a large fraction of polar/charged atoms and have chemical compositions that are similar to the more hydrophilic protein patches. Importantly, we also find a striking correspondence between the most hydrophobic protein patches and regions that mediate protein interactions. Our work thus establishes a computational framework for characterizing the emergent hydrophobicity of amphiphilic solutes, such as proteins, which display nanoscale heterogeneity, and for uncovering their interaction interfaces.
Collapse
|
12
|
Liu C, Brini E, Perez A, Dill KA. Computing Ligands Bound to Proteins Using MELD-Accelerated MD. J Chem Theory Comput 2020; 16:6377-6382. [PMID: 32910647 DOI: 10.1021/acs.jctc.0c00543] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Predicting the poses of small-molecule ligands in protein binding sites is often done by virtual screening algorithms such as DOCK. In principle, molecular dynamics (MD) using atomistic force fields could give better free-energy-based pose selection, but MD is computationally expensive. Here, we ask if modeling employing limited data (MELD)-accelerated MD (MELD × MD) can pick out the best DOCK poses taken as input. We study 30 different ligand-protein pairs. MELD × MD finds native poses, based on best free energies, in 23 out of the 30 cases, 20 of which were previously known DOCK failures. We conclude that MELD × MD can add value for predicting accurate poses of small molecules bound to proteins.
Collapse
Affiliation(s)
- Cong Liu
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794-5252, United States.,Department of Chemistry, Stony Brook University, Stony Brook, New York 11790-3400, United States
| | - Emiliano Brini
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794-5252, United States
| | - Alberto Perez
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Ken A Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794-5252, United States.,Department of Chemistry, Stony Brook University, Stony Brook, New York 11790-3400, United States.,Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794-3800, United States
| |
Collapse
|
13
|
Khramushin A, Marcu O, Alam N, Shimony O, Padhorny D, Brini E, Dill KA, Vajda S, Kozakov D, Schueler-Furman O. Modeling beta-sheet peptide-protein interactions: Rosetta FlexPepDock in CAPRI rounds 38-45. Proteins 2020; 88:1037-1049. [PMID: 31891416 PMCID: PMC7539656 DOI: 10.1002/prot.25871] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2019] [Revised: 12/17/2019] [Accepted: 12/26/2019] [Indexed: 01/09/2023]
Abstract
Peptide-protein docking is challenging due to the considerable conformational freedom of the peptide. CAPRI rounds 38-45 included two peptide-protein interactions, both characterized by a peptide forming an additional beta strand of a beta sheet in the receptor. Using the Rosetta FlexPepDock peptide docking protocol we generated top-performing, high-accuracy models for targets 134 and 135, involving an interaction between a peptide derived from L-MAG with DLC8. In addition, we were able to generate the only medium-accuracy models for a particularly challenging target, T121. In contrast to the classical peptide-mediated interaction, in which receptor side chains contact both peptide backbone and side chains, beta-sheet complementation involves a major contribution to binding by hydrogen bonds between main chain atoms. To establish how binding affinity and specificity are established in this special class of peptide-protein interactions, we extracted PeptiDBeta, a benchmark of solved structures of different protein domains that are bound by peptides via beta-sheet complementation, and tested our protocol for global peptide-docking PIPER-FlexPepDock on this dataset. We find that the beta-strand part of the peptide is sufficient to generate approximate and even high resolution models of many interactions, but inclusion of adjacent motif residues often provides additional information necessary to achieve high resolution model quality.
Collapse
Affiliation(s)
- Alisa Khramushin
- Department of Microbiologyand Molecular Genetics, Institute
for Medical Research Israel-Canada, Faculty of Medicine, The Hebrew University,
Jerusalem, Israel
| | - Orly Marcu
- Department of Microbiologyand Molecular Genetics, Institute
for Medical Research Israel-Canada, Faculty of Medicine, The Hebrew University,
Jerusalem, Israel
| | - Nawsad Alam
- Department of Microbiologyand Molecular Genetics, Institute
for Medical Research Israel-Canada, Faculty of Medicine, The Hebrew University,
Jerusalem, Israel
| | - Orly Shimony
- Department of Microbiologyand Molecular Genetics, Institute
for Medical Research Israel-Canada, Faculty of Medicine, The Hebrew University,
Jerusalem, Israel
| | - Dzmitry Padhorny
- Department of Applied Mathematics and Statistics, Stony
Brook University, New York, New York
- Laufer Center for Physical and Quantitative Biology, Stony
Brook University, New York, New York
| | - Emiliano Brini
- Laufer Center for Physical and Quantitative Biology, Stony
Brook University, New York, New York
| | - Ken A. Dill
- Laufer Center for Physical and Quantitative Biology, Stony
Brook University, New York, New York
- Department of Physics and Astronomy, Stony Brook
University, New York, New York
- Department of Chemistry, Stony Brook University, New York,
New York
| | - Sandor Vajda
- Department of Biomedical Engineering, Boston University,
Boston, Massachusetts
- Department of Chemistry, Boston University, Boston,
Massachusetts
| | - Dima Kozakov
- Department of Applied Mathematics and Statistics, Stony
Brook University, New York, New York
- Laufer Center for Physical and Quantitative Biology, Stony
Brook University, New York, New York
| | - Ora Schueler-Furman
- Department of Microbiologyand Molecular Genetics, Institute
for Medical Research Israel-Canada, Faculty of Medicine, The Hebrew University,
Jerusalem, Israel
| |
Collapse
|