1
|
Chuntakaruk H, Hengphasatporn K, Shigeta Y, Aonbangkhen C, Lee VS, Khotavivattana T, Rungrotmongkol T, Hannongbua S. FMO-guided design of darunavir analogs as HIV-1 protease inhibitors. Sci Rep 2024; 14:3639. [PMID: 38351065 PMCID: PMC10864397 DOI: 10.1038/s41598-024-53940-1] [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] [Received: 06/03/2023] [Accepted: 02/06/2024] [Indexed: 02/16/2024] Open
Abstract
The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.
Collapse
Affiliation(s)
- Hathaichanok Chuntakaruk
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, 10330, Thailand
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Kowit Hengphasatporn
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8577, Japan.
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8577, Japan
| | - Chanat Aonbangkhen
- Center of Excellence in Natural Products Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Vannajan Sanghiran Lee
- Chemistry Department, Faculty of Science, University Malaya, Kuala Lumpur, 50603, Malaysia
| | - Tanatorn Khotavivattana
- Center of Excellence in Natural Products Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| | - Thanyada Rungrotmongkol
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, 10330, Thailand.
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| | - Supot Hannongbua
- Center of Excellence in Computational Chemistry (CECC), Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| |
Collapse
|
2
|
Yunaini L, Erlina L, Fadilah F, Pujianto DA. In silico docking analysis of beta-defensin 20 against cation channel sperm-associated protein 1-4 to predict its role in the sperm maturation. Asian J Androl 2023; 25:528-532. [PMID: 36571327 PMCID: PMC10411257 DOI: 10.4103/aja2022103] [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] [Received: 07/08/2022] [Accepted: 11/16/2022] [Indexed: 12/27/2022] Open
Abstract
Beta-defensin 20 (DEFB20) is widely expressed in the epididymis with gene features involved in epididymal sperm maturation. However, the action mechanism and function of DEFB20 in sperm maturation are still unclear. One of the important roles of beta-defensin is the ion channel activity. The cation channel sperm-associated protein (CatSper) alpha is an ion channel protein found on the sperm surface. This study aimed to investigate the interaction between DEFB20 and CatSper1-4 protein in relation to the sperm maturation process. Protein sequences were obtained from the National Center for Biotechnology Information (NCBI). Protein modeling and validation were carried out by using the Robetta modeling server and the Ramachandran plot method. Rosetta web server was used for the docking analysis. The results revealed a natural interaction between DEFB20 and CatSper1-4. The interaction occurred at the cation channel (close to the casein kinase II), ion transport protein, and kinase c phosphorylation of the CatSper1-4 active site. The DEFB20 region interacting with CatSper2-4 was the beta-defensin domain, while with CatSper1 was the non-beta-defensin domain. Based on the analysis, DEFB20 may interact with CatSper α subunits, particularly CatsSper1, to affect ion channel activity during sperm maturation.
Collapse
Affiliation(s)
- Luluk Yunaini
- Doctoral Program for Biomedical Sciences, Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
- Department of Medical Biology, Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
| | - Linda Erlina
- Department of Medical Chemistry, Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
- Bioinformatics Core Facilities, Indonesia Medical Education and Research Institute (IMERI), Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
| | - Fadilah Fadilah
- Department of Medical Chemistry, Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
- Bioinformatics Core Facilities, Indonesia Medical Education and Research Institute (IMERI), Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
| | - Dwi Ari Pujianto
- Department of Medical Biology, Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
| |
Collapse
|
3
|
Rana N, Singh AK, Shuaib M, Gupta S, Habiballah MM, Alkhanani MF, Haque S, Reshi MS, Kumar S. Drug Resistance Mechanism of M46I-Mutation-Induced Saquinavir Resistance in HIV-1 Protease Using Molecular Dynamics Simulation and Binding Energy Calculation. Viruses 2022; 14:v14040697. [PMID: 35458427 PMCID: PMC9031992 DOI: 10.3390/v14040697] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Revised: 03/05/2022] [Accepted: 03/07/2022] [Indexed: 02/06/2023] Open
Abstract
Drug-resistance-associated mutation in essential proteins of the viral life cycle is a major concern in anti-retroviral therapy. M46I, a non-active site mutation in HIV-1 protease has been clinically associated with saquinavir resistance in HIV patients. A 100 ns molecular dynamics (MD) simulation and MM-PBSA calculations were performed to study the molecular mechanism of M46I-mutation-based saquinavir resistance. In order to acquire deeper insight into the drug-resistance mechanism, the flap curling, closed/semi-open/open conformations, and active site compactness were studied. The M46I mutation significantly affects the energetics and conformational stability of HIV-1 protease in terms of RMSD, RMSF, Rg, SASA, and hydrogen formation potential. This mutation significantly decreased van der Waals interaction and binding free energy (∆G) in the M46I–saquinavir complex and induced inward flap curling and a wider opening of the flaps for most of the MD simulation period. The predominant open conformation was reduced, but inward flap curling/active site compactness was increased in the presence of saquinavir in M46I HIV-1 protease. In conclusion, the M46I mutation induced structural dynamics changes that weaken the protease grip on saquinavir without distorting the active site of the protein. The produced information may be utilized for the discovery of inhibitor(s) against drug-resistant HIV-1 protease.
Collapse
Affiliation(s)
- Nilottam Rana
- Molecular Signaling & Drug Discovery Laboratory, Department of Biochemistry, Central University of Punjab, Bathinda 151401, Punjab, India; (N.R.); (A.K.S.); (M.S.)
| | - Atul Kumar Singh
- Molecular Signaling & Drug Discovery Laboratory, Department of Biochemistry, Central University of Punjab, Bathinda 151401, Punjab, India; (N.R.); (A.K.S.); (M.S.)
| | - Mohd Shuaib
- Molecular Signaling & Drug Discovery Laboratory, Department of Biochemistry, Central University of Punjab, Bathinda 151401, Punjab, India; (N.R.); (A.K.S.); (M.S.)
| | - Sanjay Gupta
- Department of Urology, Pharmacology and Pathology, Case Western Reserve University, Cleveland, OH 44106, USA;
| | - Mahmoud M. Habiballah
- Medical Laboratory Technology Department, Jazan University, Jazan 45142, Saudi Arabia;
- SMIRES for Consultation in Specialized Medical Laboratories, Jazan University, Jazan 45142, Saudi Arabia
| | - Mustfa F. Alkhanani
- Emergency Service Department, College of Applied Sciences, AlMaarefa University, Riyadh 11597, Saudi Arabia;
| | - Shafiul Haque
- Research and Scientific Studies Unit, College of Nursing and Allied Health Sciences, Jazan University, Jazan 45142, Saudi Arabia;
| | - Mohd Salim Reshi
- Toxicology and Pharmacology Lab., Department of Zoology, School of Biosciences and Biotechnology, Baba Ghulam Shah Badshah University, Rajouri 185234, Jammu & Kashmir, India;
| | - Shashank Kumar
- Molecular Signaling & Drug Discovery Laboratory, Department of Biochemistry, Central University of Punjab, Bathinda 151401, Punjab, India; (N.R.); (A.K.S.); (M.S.)
- Correspondence:
| |
Collapse
|
4
|
Varga JK, Diffley K, Welker Leng KR, Fierke CA, Schueler-Furman O. Structure-based prediction of HDAC6 substrates validated by enzymatic assay reveals determinants of promiscuity and detects new potential substrates. Sci Rep 2022; 12:1788. [PMID: 35110592 PMCID: PMC8810773 DOI: 10.1038/s41598-022-05681-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Accepted: 01/17/2022] [Indexed: 01/25/2023] Open
Abstract
Histone deacetylases play important biological roles well beyond the deacetylation of histone tails. In particular, HDAC6 is involved in multiple cellular processes such as apoptosis, cytoskeleton reorganization, and protein folding, affecting substrates such as ɑ-tubulin, Hsp90 and cortactin proteins. We have applied a biochemical enzymatic assay to measure the activity of HDAC6 on a set of candidate unlabeled peptides. These served for the calibration of a structure-based substrate prediction protocol, Rosetta FlexPepBind, previously used for the successful substrate prediction of HDAC8 and other enzymes. A proteome-wide screen of reported acetylation sites using our calibrated protocol together with the enzymatic assay provide new peptide substrates and avenues to novel potential functional regulatory roles of this promiscuous, multi-faceted enzyme. In particular, we propose novel regulatory roles of HDAC6 in tumorigenesis and cancer cell survival via the regulation of EGFR/Akt pathway activation. The calibration process and comparison of the results between HDAC6 and HDAC8 highlight structural differences that explain the established promiscuity of HDAC6.
Collapse
Affiliation(s)
- Julia K Varga
- Department of Microbiology and Molecular Genetics, Institute for Medical Research Israel-Canada (IMRIC), The Hebrew University of Jerusalem, Faculty of Medicine, POB 12272, 9112102, Jerusalem, Israel
| | - Kelsey Diffley
- Department of Chemistry, University of Michigan, 930 North University Avenue, Ann Arbor, MI, 48109, USA
| | - Katherine R Welker Leng
- Department of Chemistry, University of Michigan, 930 North University Avenue, Ann Arbor, MI, 48109, USA
| | - Carol A Fierke
- Department of Chemistry, University of Michigan, 930 North University Avenue, Ann Arbor, MI, 48109, USA
- Department of Biochemistry, Brandeis University, 415 South Street, Waltham, MA, 02453, USA
| | - Ora Schueler-Furman
- Department of Microbiology and Molecular Genetics, Institute for Medical Research Israel-Canada (IMRIC), The Hebrew University of Jerusalem, Faculty of Medicine, POB 12272, 9112102, Jerusalem, Israel.
| |
Collapse
|
5
|
Mahajan SP, Srinivasan Y, Labonte JW, DeLisa MP, Gray JJ. Structural basis for peptide substrate specificities of glycosyltransferase GalNAc-T2. ACS Catal 2021; 11:2977-2991. [PMID: 34322281 DOI: 10.1021/acscatal.0c04609] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The polypeptide N-acetylgalactosaminyl transferase (GalNAc-T) enzyme family initiates O-linked mucin-type glycosylation. The family constitutes 20 isoenzymes in humans. GalNAc-Ts exhibit both redundancy and finely tuned specificity for a wide range of peptide substrates. In this work, we deciphered the sequence and structural motifs that determine the peptide substrate preferences for the GalNAc-T2 isoform. Our approach involved sampling and characterization of peptide-enzyme conformations obtained from Rosetta Monte Carlo-minimization-based flexible docking. We computationally scanned 19 amino acid residues at positions -1 and +1 of an eight-residue peptide substrate, which comprised a dataset of 361 (19x19) peptides with previously characterized experimental GalNAc-T2 glycosylation efficiencies. The calculations recapitulated experimental specificity data, successfully discriminating between glycosylatable and non-glycosylatable peptides with a probability of 96.5% (ROC-AUC score), a balanced accuracy of 85.5% and a false positive rate of 7.3%. The glycosylatable peptide substrates viz. peptides with proline, serine, threonine, and alanine at the -1 position of the peptide preferentially exhibited cognate sequon-like conformations. The preference for specific residues at the -1 position of the peptide was regulated by enzyme residues R362, K363, Q364, H365 and W331, which modulate the pocket size and specific enzyme-peptide interactions. For the +1 position of the peptide, enzyme residues K281 and K363 formed gating interactions with aromatics and glutamines at the +1 position of the peptide, leading to modes of peptide-binding sub-optimal for catalysis. Overall, our work revealed enzyme features that lead to the finely tuned specificity observed for a broad range of peptide substrates for the GalNAc-T2 enzyme. We anticipate that the key sequence and structural motifs can be extended to analyze specificities of other isoforms of the GalNAc-T family and can be used to guide design of variants with tailored specificity.
Collapse
Affiliation(s)
- Sai Pooja Mahajan
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Yashes Srinivasan
- Department of Bioengineering, University of California—Los Angeles, Los Angeles, California 90095, United States
| | - Jason W. Labonte
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Department of Chemistry, Franklin & Marshall College, Lancaster, Pennsylvania 17604, United States
| | - Matthew P. DeLisa
- Robert Frederick Smith School of Chemical and Biomolecular Engineering, Department of Microbiology, and Nancy E. and Peter C. Meinig School of Biomedical Engineering, Biochemistry, Molecular and Cell Biology, Cornell University, Ithaca, New York 14853, United States
| | - Jeffrey J. Gray
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins School of Medicine, Baltimore, Maryland 21224, United States
| |
Collapse
|
6
|
Fischer A, Sellner M, Mitusińska K, Bzówka M, Lill MA, Góra A, Smieško M. Computational Selectivity Assessment of Protease Inhibitors against SARS-CoV-2. Int J Mol Sci 2021; 22:2065. [PMID: 33669738 PMCID: PMC7922391 DOI: 10.3390/ijms22042065] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Revised: 02/08/2021] [Accepted: 02/11/2021] [Indexed: 12/27/2022] Open
Abstract
The pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses a serious global health threat. Since no specific therapeutics are available, researchers around the world screened compounds to inhibit various molecular targets of SARS-CoV-2 including its main protease (Mpro) essential for viral replication. Due to the high urgency of these discovery efforts, off-target binding, which is one of the major reasons for drug-induced toxicity and safety-related drug attrition, was neglected. Here, we used molecular docking, toxicity profiling, and multiple molecular dynamics (MD) protocols to assess the selectivity of 33 reported non-covalent inhibitors of SARS-CoV-2 Mpro against eight proteases and 16 anti-targets. The panel of proteases included SARS-CoV Mpro, cathepsin G, caspase-3, ubiquitin carboxy-terminal hydrolase L1 (UCHL1), thrombin, factor Xa, chymase, and prostasin. Several of the assessed compounds presented considerable off-target binding towards the panel of proteases, as well as the selected anti-targets. Our results further suggest a high risk of off-target binding to chymase and cathepsin G. Thus, in future discovery projects, experimental selectivity assessment should be directed toward these proteases. A systematic selectivity assessment of SARS-CoV-2 Mpro inhibitors, as we report it, was not previously conducted.
Collapse
Affiliation(s)
- André Fischer
- Computational Pharmacy, Departement of Pharmaceutical Sciences, University of Basel, 4056 Basel, Switzerland; (A.F.); (M.S.)
| | - Manuel Sellner
- Computational Pharmacy, Departement of Pharmaceutical Sciences, University of Basel, 4056 Basel, Switzerland; (A.F.); (M.S.)
| | - Karolina Mitusińska
- Tunneling Group, Biotechnology Centre, ul. Krzywoustego 8, Silesian University of Technology, 44-100 Gliwice, Poland; (K.M.); (M.B.)
| | - Maria Bzówka
- Tunneling Group, Biotechnology Centre, ul. Krzywoustego 8, Silesian University of Technology, 44-100 Gliwice, Poland; (K.M.); (M.B.)
| | - Markus A. Lill
- Computational Pharmacy, Departement of Pharmaceutical Sciences, University of Basel, 4056 Basel, Switzerland; (A.F.); (M.S.)
| | - Artur Góra
- Tunneling Group, Biotechnology Centre, ul. Krzywoustego 8, Silesian University of Technology, 44-100 Gliwice, Poland; (K.M.); (M.B.)
| | - Martin Smieško
- Computational Pharmacy, Departement of Pharmaceutical Sciences, University of Basel, 4056 Basel, Switzerland; (A.F.); (M.S.)
| |
Collapse
|
7
|
Yu XX, Liang WY, Yin JY, Zhou Q, Chen DM, Zhang YH. Combining experimental techniques with molecular dynamics to investigate the impact of different enzymatic hydrolysis of β-lactoglobulin on the antigenicity reduction. Food Chem 2021; 350:129139. [PMID: 33588281 DOI: 10.1016/j.foodchem.2021.129139] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2020] [Revised: 12/07/2020] [Accepted: 01/16/2021] [Indexed: 11/26/2022]
Abstract
β-Lactoglobulin (β-LG) is one of the major food allergens. Enzymatic hydrolysis is a promising strategy to reduce the antigenicity of β-LG in industrial production. The relationship between the cleavage sites of β-LG by protease and its antigenic active sites were explored in this study. Molecular docking and molecular dynamics (MD) were used to analyze the active sites and interaction force of β-LG and IgG antibody. Whey protein was hydrolyzed by four specific enzymes and the antigenicity of the hydrolysates were determined by ELISA. The results of MD showed that the amino acid residue Gln155 (-4.48 kcal mol-1) played the most important roles in the process of binding. Hydrolysates produced by AY-10, which was the only one with specificity towards cleavage sites next to a Gln, had the lowest antigenicity at the same hydrolysis degree. Antigenicity decrease was related to the energy contribution of the cleavage site in the active sites.
Collapse
Affiliation(s)
- Xin-Xin Yu
- Key Laboratory of Dairy Science, Ministry of Education, Northeast Agricultural University, Harbin 150030, PR China; Department of Food Science, Northeast Agricultural University, Harbin 150030, PR China
| | - Wei-Yue Liang
- Key Laboratory of Dairy Science, Ministry of Education, Northeast Agricultural University, Harbin 150030, PR China; Department of Food Science, Northeast Agricultural University, Harbin 150030, PR China
| | - Jia-Yi Yin
- Key Laboratory of Dairy Science, Ministry of Education, Northeast Agricultural University, Harbin 150030, PR China; Department of Food Science, Northeast Agricultural University, Harbin 150030, PR China
| | - Qian Zhou
- Key Laboratory of Dairy Science, Ministry of Education, Northeast Agricultural University, Harbin 150030, PR China; Department of Food Science, Northeast Agricultural University, Harbin 150030, PR China
| | - Dong-Mei Chen
- Key Laboratory of Dairy Science, Ministry of Education, Northeast Agricultural University, Harbin 150030, PR China; Department of Food Science, Northeast Agricultural University, Harbin 150030, PR China
| | - Ying-Hua Zhang
- Key Laboratory of Dairy Science, Ministry of Education, Northeast Agricultural University, Harbin 150030, PR China; Department of Food Science, Northeast Agricultural University, Harbin 150030, PR China.
| |
Collapse
|
8
|
On the cutting edge: protease-based methods for sensing and controlling cell biology. Nat Methods 2020; 17:885-896. [PMID: 32661424 DOI: 10.1038/s41592-020-0891-z] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2018] [Accepted: 06/09/2020] [Indexed: 02/06/2023]
Abstract
Sequence-specific proteases have proven to be versatile building blocks for tools that report or control cellular function. Reporting methods link protease activity to biochemical signals, whereas control methods rely on engineering proteases to respond to exogenous inputs such as light or chemicals. In turn, proteases have inherent control abilities, as their native functions are to release, activate or destroy proteins by cleavage, with the irreversibility of proteolysis allowing sustained downstream effects. As a result, protease-based synthetic circuits have been created for diverse uses such as reporting cellular signaling, tuning protein expression, controlling viral replication and detecting cancer states. Here, we comprehensively review the development and application of protease-based methods for reporting and controlling cellular function in eukaryotes.
Collapse
|
9
|
Abstract
This chapter describes the current status of development of the fragment molecular orbital (FMO) method for analyzing the electronic state and intermolecular interactions of biomolecular systems in solvent. The orbital energies and the inter-fragment interaction energies (IFIEs) for a specific molecular structure can be obtained directly by performing FMO calculations by exposing water molecules and counterions around biomolecular systems. Then, it is necessary to pay attention to the thickness of the water shell surrounding the biomolecules. The single-point calculation for snapshots from MD trajectory does not incorporate the effects of temperature and configurational fluctuation, but the SCIFIE (statistically corrected IFIE) method is proposed as a many-body correlated method that partially compensates for this deficiency. Furthermore, implicit continuous dielectric models have been developed as effective approaches to incorporating the screening effect of the solvent in thermal equilibrium, and we illustrate their usefulness for theoretical evaluation of IFIEs and ligand-binding free energy on the basis of the FMO-PBSA (Poisson-Boltzmann surface area) method and other computational methods.
Collapse
|
10
|
Voshavar C. Protease Inhibitors for the Treatment of HIV/AIDS: Recent Advances and Future Challenges. Curr Top Med Chem 2019; 19:1571-1598. [PMID: 31237209 DOI: 10.2174/1568026619666190619115243] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2019] [Revised: 05/22/2019] [Accepted: 05/27/2019] [Indexed: 02/07/2023]
Abstract
Acquired Immunodeficiency Syndrome (AIDS) is a chronic disease characterized by multiple life-threatening illnesses caused by a retro-virus, Human Immunodeficiency Virus (HIV). HIV infection slowly destroys the immune system and increases the risk of various other infections and diseases. Although, there is no immediate cure for HIV infection/AIDS, several drugs targeting various cruxes of HIV infection are used to slow down the progress of the disease and to boost the immune system. One of the key therapeutic strategies is Highly Active Antiretroviral Therapy (HAART) or ' AIDS cocktail' in a general sense, which is a customized combination of anti-retroviral drugs designed to combat the HIV infection. Since HAART's inception in 1995, this treatment was found to be effective in improving the life expectancy of HIV patients over two decades. Among various classes of HAART treatment regimen, Protease Inhibitors (PIs) are known to be widely used as a major component and found to be effective in treating HIV infection/AIDS. For the past several years, a variety of protease inhibitors have been reported. This review outlines the drug design strategies of PIs, chemical and pharmacological characteristics of some mechanism-based inhibitors, summarizes the recent developments in small molecule based drug discovery with HIV protease as a drug target. Further discussed are the pharmacology, PI drug resistance on HIV PR, adverse effects of HIV PIs and challenges/impediments in the successful application of HIV PIs as an important class of drugs in HAART regimen for the effective treatment of AIDS.
Collapse
Affiliation(s)
- Chandrashekhar Voshavar
- College of Pharmacy and Pharmaceutical Sciences, Florida A&M University, Tallahassee, FL 32307, United States
| |
Collapse
|
11
|
Paulsen JL, Leidner F, Ragland DA, Kurt Yilmaz N, Schiffer CA. Interdependence of Inhibitor Recognition in HIV-1 Protease. J Chem Theory Comput 2017; 13:2300-2309. [PMID: 28358514 PMCID: PMC5425943 DOI: 10.1021/acs.jctc.6b01262] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
![]()
Molecular recognition
is a highly interdependent process. Subsite
couplings within the active site of proteases are most often revealed
through conditional amino acid preferences in substrate recognition.
However, the potential effect of these couplings on inhibition and
thus inhibitor design is largely unexplored. The present study examines
the interdependency of subsites in HIV-1 protease using a focused
library of protease inhibitors, to aid in future inhibitor design.
Previously a series of darunavir (DRV) analogs was designed to systematically
probe the S1′ and S2′ subsites. Co-crystal structures
of these analogs with HIV-1 protease provide the ideal opportunity
to probe subsite interdependency. All-atom molecular dynamics simulations
starting from these structures were performed and systematically analyzed
in terms of atomic fluctuations, intermolecular interactions, and
water structure. These analyses reveal that the S1′ subsite
highly influences other subsites: the extension of the hydrophobic
P1′ moiety results in 1) reduced van der Waals contacts in
the P2′ subsite, 2) more variability in the hydrogen bond frequencies
with catalytic residues and the flap water, and 3) changes in the
occupancy of conserved water sites both proximal and distal to the
active site. In addition, one of the monomers in this homodimeric
enzyme has atomic fluctuations more highly correlated with DRV than
the other monomer. These relationships intricately link the HIV-1
protease subsites and are critical to understanding molecular recognition
and inhibitor binding. More broadly, the interdependency of subsite
recognition within an active site requires consideration in the selection
of chemical moieties in drug design; this strategy is in contrast
to what is traditionally done with independent optimization of chemical
moieties of an inhibitor.
Collapse
Affiliation(s)
- Janet L Paulsen
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School , Worcester, Massachusetts 01605, United States
| | - Florian Leidner
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School , Worcester, Massachusetts 01605, United States
| | - Debra A Ragland
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School , Worcester, Massachusetts 01605, United States
| | - Nese Kurt Yilmaz
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School , Worcester, Massachusetts 01605, United States
| | - Celia A Schiffer
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School , Worcester, Massachusetts 01605, United States
| |
Collapse
|
12
|
Modeling Peptide-Protein Structure and Binding Using Monte Carlo Sampling Approaches: Rosetta FlexPepDock and FlexPepBind. Methods Mol Biol 2017; 1561:139-169. [PMID: 28236237 DOI: 10.1007/978-1-4939-6798-8_9] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Many signaling and regulatory processes involve peptide-mediated protein interactions, i.e., the binding of a short stretch in one protein to a domain in its partner. Computational tools that generate accurate models of peptide-receptor structures and binding improve characterization and manipulation of known interactions, help to discover yet unknown peptide-protein interactions and networks, and bring into reach the design of peptide-based drugs for targeting specific systems of medical interest.Here, we present a concise overview of the Rosetta FlexPepDock protocol and its derivatives that we have developed for the structure-based characterization of peptide-protein binding. Rosetta FlexPepDock was built to generate precise models of protein-peptide complex structures, by effectively addressing the challenge of the considerable conformational flexibility of the peptide. Rosetta FlexPepBind is an extension of this protocol that allows characterizing peptide-binding affinities and specificities of various biological systems, based on the structural models generated by Rosetta FlexPepDock. We provide detailed descriptions and guidelines for the usage of these protocols, and on a specific example, we highlight the variety of different challenges that can be met and the questions that can be answered with Rosetta FlexPepDock.
Collapse
|
13
|
Pethe MA, Rubenstein AB, Khare SD. Large-Scale Structure-Based Prediction and Identification of Novel Protease Substrates Using Computational Protein Design. J Mol Biol 2016; 429:220-236. [PMID: 27932294 DOI: 10.1016/j.jmb.2016.11.031] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2016] [Revised: 11/23/2016] [Accepted: 11/30/2016] [Indexed: 12/16/2022]
Abstract
Characterizing the substrate specificity of protease enzymes is critical for illuminating the molecular basis of their diverse and complex roles in a wide array of biological processes. Rapid and accurate prediction of their extended substrate specificity would also aid in the design of custom proteases capable of selectively and controllably cleaving biotechnologically or therapeutically relevant targets. However, current in silico approaches for protease specificity prediction, rely on, and are therefore limited by, machine learning of sequence patterns in known experimental data. Here, we describe a general approach for predicting peptidase substrates de novo using protein structure modeling and biophysical evaluation of enzyme-substrate complexes. We construct atomic resolution models of thousands of candidate substrate-enzyme complexes for each of five model proteases belonging to the four major protease mechanistic classes-serine, cysteine, aspartyl, and metallo-proteases-and develop a discriminatory scoring function using enzyme design modules from Rosetta and AMBER's MMPBSA. We rank putative substrates based on calculated interaction energy with a modeled near-attack conformation of the enzyme active site. We show that the energetic patterns obtained from these simulations can be used to robustly rank and classify known cleaved and uncleaved peptides and that these structural-energetic patterns have greater discriminatory power compared to purely sequence-based statistical inference. Combining sequence and energetic patterns using machine-learning algorithms further improves classification performance, and analysis of structural models provides physical insight into the structural basis for the observed specificities. We further tested the predictive capability of the model by designing and experimentally characterizing the cleavage of four novel substrate motifs for the hepatitis C virus NS3/4 protease using an in vivo assay. The presented structure-based approach is generalizable to other protease enzymes with known or modeled structures, and complements existing experimental methods for specificity determination.
Collapse
Affiliation(s)
- Manasi A Pethe
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Aliza B Rubenstein
- Computational Biology & Molecular Biophysics Program, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Sagar D Khare
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Computational Biology & Molecular Biophysics Program, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA.
| |
Collapse
|
14
|
Governa P, Giachetti D, Biagi M, Manetti F, De Vico L. Hypothesis on Serenoa repens (Bartram) small extract inhibition of prostatic 5 α-reductase through an in silico approach on 5 β-reductase x-ray structure. PeerJ 2016; 4:e2698. [PMID: 27904805 PMCID: PMC5126621 DOI: 10.7717/peerj.2698] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2016] [Accepted: 10/18/2016] [Indexed: 11/20/2022] Open
Abstract
Benign prostatic hyperplasia is a common disease in men aged over 50 years old, with an incidence increasing to more than 80% over the age of 70, that is increasingly going to attract pharmaceutical interest. Within conventional therapies, such as α-adrenoreceptor antagonists and 5α-reductase inhibitor, there is a large requirement for treatments with less adverse events on, e.g., blood pressure and sexual function: phytotherapy may be the right way to fill this need. Serenoa repens standardized extract has been widely studied and its ability to reduce lower urinary tract symptoms related to benign prostatic hyperplasia is comprehensively described in literature. An innovative investigation on the mechanism of inhibition of 5α-reductase by Serenoa repens extract active principles is proposed in this work through computational methods, performing molecular docking simulations on the crystal structure of human liver 5β-reductase. The results confirm that both sterols and fatty acids can play a role in the inhibition of the enzyme, thus, suggesting a competitive mechanism of inhibition. This work proposes a further confirmation for the rational use of herbal products in the management of benign prostatic hyperplasia, and suggests computational methods as an innovative, low cost, and non-invasive process for the study of phytocomplex activity toward proteic targets.
Collapse
Affiliation(s)
- Paolo Governa
- Department of Physical Sciences, Earth and Environment, University of Siena, Siena, Italy; Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Daniela Giachetti
- Department of Physical Sciences, Earth and Environment, University of Siena , Siena , Italy
| | - Marco Biagi
- Department of Physical Sciences, Earth and Environment, University of Siena , Siena , Italy
| | - Fabrizio Manetti
- Department of Biotechnology, Chemistry and Pharmacy, University of Siena , Siena , Italy
| | - Luca De Vico
- Department of Chemistry, University of Copenhagen , Copenhagen , Denmark
| |
Collapse
|
15
|
Abstract
The HIV genome encodes a small number of viral proteins (i.e., 16), invariably establishing cooperative associations among HIV proteins and between HIV and host proteins, to invade host cells and hijack their internal machineries. As a known example, the HIV envelope glycoprotein GP120 is closely associated with GP41 for viral entry. From a genome-wide perspective, a hypothesis can be worked out to determine whether 16 HIV proteins could develop 120 possible pairwise associations either by physical interactions or by functional associations mediated via HIV or host molecules. Here, we present the first systematic review of experimental evidence on HIV genome-wide protein associations using a large body of publications accumulated over the past 3 decades. Of 120 possible pairwise associations between 16 HIV proteins, at least 34 physical interactions and 17 functional associations have been identified. To achieve efficient viral replication and infection, HIV protein associations play essential roles (e.g., cleavage, inhibition, and activation) during the HIV life cycle. In either a dispensable or an indispensable manner, each HIV protein collaborates with another viral protein to accomplish specific activities that precisely take place at the proper stages of the HIV life cycle. In addition, HIV genome-wide protein associations have an impact on anti-HIV inhibitors due to the extensive cross talk between drug-inhibited proteins and other HIV proteins. Overall, this study presents for the first time a comprehensive overview of HIV genome-wide protein associations, highlighting meticulous collaborations between all viral proteins during the HIV life cycle.
Collapse
Affiliation(s)
- Guangdi Li
- Department of Metabolism and Endocrinology, Metabolic Syndrome Research Center, Key Laboratory of Diabetes Immunology, Ministry of Education, National Clinical Research Center for Metabolic Diseases, The Second Xiangya Hospital, Central South University, Changsha, Hunan, China KU Leuven-University of Leuven, Rega Institute for Medical Research, Department of Microbiology and Immunology, Leuven, Belgium
| | - Erik De Clercq
- KU Leuven-University of Leuven, Rega Institute for Medical Research, Department of Microbiology and Immunology, Leuven, Belgium
| |
Collapse
|
16
|
HIV Genome-Wide Protein Associations: a Review of 30 Years of Research. Microbiol Mol Biol Rev 2016; 80:679-731. [PMID: 27357278 DOI: 10.1128/mmbr.00065-15] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The HIV genome encodes a small number of viral proteins (i.e., 16), invariably establishing cooperative associations among HIV proteins and between HIV and host proteins, to invade host cells and hijack their internal machineries. As a known example, the HIV envelope glycoprotein GP120 is closely associated with GP41 for viral entry. From a genome-wide perspective, a hypothesis can be worked out to determine whether 16 HIV proteins could develop 120 possible pairwise associations either by physical interactions or by functional associations mediated via HIV or host molecules. Here, we present the first systematic review of experimental evidence on HIV genome-wide protein associations using a large body of publications accumulated over the past 3 decades. Of 120 possible pairwise associations between 16 HIV proteins, at least 34 physical interactions and 17 functional associations have been identified. To achieve efficient viral replication and infection, HIV protein associations play essential roles (e.g., cleavage, inhibition, and activation) during the HIV life cycle. In either a dispensable or an indispensable manner, each HIV protein collaborates with another viral protein to accomplish specific activities that precisely take place at the proper stages of the HIV life cycle. In addition, HIV genome-wide protein associations have an impact on anti-HIV inhibitors due to the extensive cross talk between drug-inhibited proteins and other HIV proteins. Overall, this study presents for the first time a comprehensive overview of HIV genome-wide protein associations, highlighting meticulous collaborations between all viral proteins during the HIV life cycle.
Collapse
|
17
|
Bayden AS, Gomez EF, Audie J, Chakravorty DK, Diller DJ. A combined cheminformatic and bioinformatic approach to address the proteolytic stability challenge in peptide-based drug discovery. Biopolymers 2015; 104:775-89. [PMID: 26270398 DOI: 10.1002/bip.22711] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2015] [Revised: 07/22/2015] [Accepted: 08/09/2015] [Indexed: 11/10/2022]
Abstract
We have created models to predict cleavage sites for several human proteases including caspase-1, caspase-3, caspase-6, caspase-7, cathepsin B, cathepsin D, cathepsin G, cathepsin K, cathepsin L, elastase-2, granzyme A, granzyme B, matrix metallopeptidase-2 (MMP2), MMP7, MMP9, thrombin, and trypsin-1. Rather than representing the sequence pattern around the potential cleavage site through a series of flags with each flag representing one of the 20 standard amino acids, we first represent each amino acid by its calculated properties. For these calculated properties, we use validated cheminformatic descriptors, such as molecular weight, logP, and polar surface area, of the individual amino acids. Finally, the cleavage site-specific descriptors are calculated through various combinations of the individual amino acid descriptors for the residues surrounding the cleavage site. Some of these combinations do not take into account the location of the residue, as long as it is in a prescribed neighborhood of the potential cleavage site, whereas others are sensitive to the precise order of the residues in the sequence. The key advantage of this approach is that it allows one to perform meaningful calculations with nonstandard amino acids for which little or no data exists. Finally, using both docking and molecular dynamics simulations, we examine the potential for and limitations of protease crystal structures to impact the design of proteolytically stable peptides.
Collapse
Affiliation(s)
| | - Edwin F Gomez
- Department of Chemistry, University of New Orleans, New Orleans, LA
| | - Joseph Audie
- CMDBioscience Inc., 5 Science Park, New Haven, CT
| | | | | |
Collapse
|
18
|
Sellers MS, Hurley MM. XPairIt Docking Protocolfor peptide docking and analysis. MOLECULAR SIMULATION 2015. [DOI: 10.1080/08927022.2015.1025267] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
19
|
Liu S, Liu S, Wang Y, Liao Z. The P2/P2′ sites affect the substrate cleavage of TNF-α converting enzyme (TACE). Mol Immunol 2014; 62:122-8. [DOI: 10.1016/j.molimm.2014.05.017] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2014] [Revised: 05/28/2014] [Accepted: 05/29/2014] [Indexed: 01/08/2023]
|
20
|
Jensen JH, Willemoës M, Winther JR, De Vico L. In silico prediction of mutant HIV-1 proteases cleaving a target sequence. PLoS One 2014; 9:e95833. [PMID: 24796579 PMCID: PMC4010418 DOI: 10.1371/journal.pone.0095833] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2013] [Accepted: 03/31/2014] [Indexed: 11/17/2022] Open
Abstract
HIV-1 protease represents an appealing system for directed enzyme re-design, since it has various different endogenous targets, a relatively simple structure and it is well studied. Recently Chaudhury and Gray (Structure (2009) 17: 1636–1648) published a computational algorithm to discern the specificity determining residues of HIV-1 protease. In this paper we present two computational tools aimed at re-designing HIV-1 protease, derived from the algorithm of Chaudhuri and Gray. First, we present an energy-only based methodology to discriminate cleavable and non cleavable peptides for HIV-1 proteases, both wild type and mutant. Secondly, we show an algorithm we developed to predict mutant HIV-1 proteases capable of cleaving a new target substrate peptide, different from the natural targets of HIV-1 protease. The obtained in silico mutant enzymes were analyzed in terms of cleavability and specificity towards the target peptide using the energy-only methodology. We found two mutant proteases as best candidates for specificity and cleavability towards the target sequence.
Collapse
Affiliation(s)
- Jan H Jensen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Martin Willemoës
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Jakob R Winther
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Luca De Vico
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
21
|
London N, Raveh B, Schueler-Furman O. Peptide docking and structure-based characterization of peptide binding: from knowledge to know-how. Curr Opin Struct Biol 2013; 23:894-902. [PMID: 24138780 DOI: 10.1016/j.sbi.2013.07.006] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2013] [Revised: 07/04/2013] [Accepted: 07/08/2013] [Indexed: 11/25/2022]
Abstract
Peptide-mediated interactions are gaining increased attention due to their predominant roles in the many regulatory processes that involve dynamic interactions between proteins. The structures of such interactions provide an excellent starting point for their characterization and manipulation, and can provide leads for targeted inhibitor design. The relatively few experimentally determined structures of peptide-protein complexes can be complemented with an outburst of modeling approaches that have been introduced in recent years, with increasing accuracy and applicability to ever more systems. We review different methods to address the considerable challenges in modeling the binding of a short yet highly flexible peptide to its partner. These methods apply an array of sampling strategies and draw from a recent amassing of knowledge about the biophysical nature of peptide-protein interactions. We elaborate on applications of these structure-based approaches and in particular on the characterization of peptide binding specificity to different peptide-binding domains and enzymes. Such applications can identify new biological targets and thus complement our current view of protein-protein interactions in living organisms. Accurate peptide-protein docking is of particular importance in the light of increased appreciation of the crucial functional roles of disordered regions and the many linear binding motifs embedded within.
Collapse
Affiliation(s)
- Nir London
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA
| | | | | |
Collapse
|
22
|
Insight into structural and biochemical determinants of substrate specificity of PFI1625c: Correlation analysis of protein-peptide molecular models. J Mol Graph Model 2013; 43:21-30. [DOI: 10.1016/j.jmgm.2013.03.008] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2012] [Revised: 03/18/2013] [Accepted: 03/28/2013] [Indexed: 11/21/2022]
|
23
|
Recent work in the development and application of protein–peptide docking. Future Med Chem 2012; 4:1619-44. [DOI: 10.4155/fmc.12.99] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Interest in the development of novel peptide-based drugs is growing. There is, thus, a pressing need for the development of effective methods to enable novel peptide-based drug discovery. A cogent case can be made for the development and application of computational or in silico methods to assist with peptide discovery. In particular, there is a need for the development of effective protein–peptide docking methods. Here, recent work in the area of protein–peptide docking method development is reviewed and several drug-discovery projects that benefited from protein–peptide docking are discussed. In the present review, special attention is given to the search and scoring problems, the use of peptide docking to enable hit identification, and the use of peptide docking to help rationalize experimental results, and generate and test structure-based hypotheses. Finally, some recommendations are made for improving the future development and application of protein–peptide docking.
Collapse
|
24
|
London N, Gullá S, Keating AE, Schueler-Furman O. In silico and in vitro elucidation of BH3 binding specificity toward Bcl-2. Biochemistry 2012; 51:5841-50. [PMID: 22702834 DOI: 10.1021/bi3003567] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Interactions between Bcl-2-like proteins and BH3 domains play a key role in the regulation of apoptosis. Despite the overall structural similarity of their interaction with helical BH3 domains, Bcl-2-like proteins exhibit an intricate spectrum of binding specificities whose underlying basis is not well understood. Here, we characterize these interactions using Rosetta FlexPepBind, a protocol for the prediction of peptide binding specificity that evaluates the binding potential of different peptides based on structural models of the corresponding peptide-receptor complexes. For two prominent players, Bcl-xL and Mcl-1, we obtain good agreement with a large set of experimental SPOT array measurements and recapitulate the binding specificity of peptides derived by yeast display in a previous study. We extend our approach to a third member of this family, Bcl-2: we test our blind prediction of the binding of 180 BIM-derived peptides with a corresponding experimental SPOT array. Both prediction and experiment reveal a Bcl-2 binding specificity pattern that resembles that of Bcl-xL. Finally, we extend this application to accurately predict the specificity pattern of additional human BH3-only derived peptides. This study characterizes the distinct patterns of binding specificity of BH3-only derived peptides for the Bcl-2 like proteins Bcl-xL, Mcl-1, and Bcl-2 and provides insight into the structural basis of determinants of specificity.
Collapse
Affiliation(s)
- Nir London
- Department of Microbiology and Molecular Genetics, Institute for Medical Research Israel-Canada, Hadassah Medical School, The Hebrew University, POB 12272, Jerusalem 91120, Israel
| | | | | | | |
Collapse
|
25
|
London N, Lamphear CL, Hougland JL, Fierke CA, Schueler-Furman O. Identification of a novel class of farnesylation targets by structure-based modeling of binding specificity. PLoS Comput Biol 2011; 7:e1002170. [PMID: 21998565 PMCID: PMC3188499 DOI: 10.1371/journal.pcbi.1002170] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2011] [Accepted: 07/01/2011] [Indexed: 11/19/2022] Open
Abstract
Farnesylation is an important post-translational modification catalyzed by farnesyltransferase (FTase). Until recently it was believed that a C-terminal CaaX motif is required for farnesylation, but recent experiments have revealed larger substrate diversity. In this study, we propose a general structural modeling scheme to account for peptide binding specificity and recapitulate the experimentally derived selectivity profile of FTase in vitro. In addition to highly accurate recovery of known FTase targets, we also identify a range of novel potential targets in the human genome, including a new substrate class with an acidic C-terminal residue (CxxD/E). In vitro experiments verified farnesylation of 26/29 tested peptides, including both novel human targets, as well as peptides predicted to tightly bind FTase. This study extends the putative range of biological farnesylation substrates. Moreover, it suggests that the ability of a peptide to bind FTase is a main determinant for the farnesylation reaction. Finally, simple adaptation of our approach can contribute to more accurate and complete elucidation of peptide-mediated interactions and modifications in the cell.
Collapse
Affiliation(s)
- Nir London
- Department of Microbiology and Molecular Genetics, Institute for Medical Research Israel-Canada, Hadassah Medical School, The Hebrew University, Jerusalem, Israel
| | - Corissa L. Lamphear
- Department of Biological Chemistry, University of Michigan, Ann Arbor, Michigan, United States of America
| | - James L. Hougland
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Carol A. Fierke
- Department of Biological Chemistry, University of Michigan, Ann Arbor, Michigan, United States of America
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Ora Schueler-Furman
- Department of Microbiology and Molecular Genetics, Institute for Medical Research Israel-Canada, Hadassah Medical School, The Hebrew University, Jerusalem, Israel
- * E-mail:
| |
Collapse
|
26
|
Donsky E, Wolfson HJ. PepCrawler: a fast RRT-based algorithm for high-resolution refinement and binding affinity estimation of peptide inhibitors. Bioinformatics 2011; 27:2836-42. [PMID: 21880702 DOI: 10.1093/bioinformatics/btr498] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
MOTIVATION Design of protein-protein interaction (PPI) inhibitors is a key challenge in structural bioinformatics and computer-aided drug design. Peptides, which partially mimic the interface area of one of the interacting proteins, are natural candidates to form protein-peptide complexes competing with the original PPI. The prediction of such complexes is especially challenging due to the high flexibility of peptide conformations. RESULTS In this article, we present PepCrawler, a new tool for deriving binding peptides from protein-protein complexes and prediction of peptide-protein complexes, by performing high-resolution docking refinement and estimation of binding affinity. By using a fast path planning approach, PepCrawler rapidly generates large amounts of flexible peptide conformations, allowing backbone and side chain flexibility. A newly introduced binding energy funnel 'steepness score' was applied for the evaluation of the protein-peptide complexes binding affinity. PepCrawler simulations predicted high binding affinity for native protein-peptide complexes benchmark and low affinity for low-energy decoy complexes. In three cases, where wet lab data are available, the PepCrawler predictions were consistent with the data. Comparing to other state of the art flexible peptide-protein structure prediction algorithms, our algorithm is very fast, and takes only minutes to run on a single PC. AVAILABILITY http://bioinfo3d.cs.tau.ac.il/PepCrawler/ CONTACT eladdons@tau.ac.il; wolfson@tau.ac.il.
Collapse
Affiliation(s)
- Elad Donsky
- Blavatnik School of Computer Science, Tel Aviv University, Tel Aviv 69978, Israel.
| | | |
Collapse
|
27
|
Chaudhury S, Berrondo M, Weitzner BD, Muthu P, Bergman H, Gray JJ. Benchmarking and analysis of protein docking performance in Rosetta v3.2. PLoS One 2011; 6:e22477. [PMID: 21829626 PMCID: PMC3149062 DOI: 10.1371/journal.pone.0022477] [Citation(s) in RCA: 223] [Impact Index Per Article: 17.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2011] [Accepted: 06/22/2011] [Indexed: 11/30/2022] Open
Abstract
RosettaDock has been increasingly used in protein docking and design strategies in order to predict the structure of protein-protein interfaces. Here we test capabilities of RosettaDock 3.2, part of the newly developed Rosetta v3.2 modeling suite, against Docking Benchmark 3.0, and compare it with RosettaDock v2.3, the latest version of the previous Rosetta software package. The benchmark contains a diverse set of 116 docking targets including 22 antibody-antigen complexes, 33 enzyme-inhibitor complexes, and 60 ‘other’ complexes. These targets were further classified by expected docking difficulty into 84 rigid-body targets, 17 medium targets, and 14 difficult targets. We carried out local docking perturbations for each target, using the unbound structures when available, in both RosettaDock v2.3 and v3.2. Overall the performances of RosettaDock v2.3 and v3.2 were similar. RosettaDock v3.2 achieved 56 docking funnels, compared to 49 in v2.3. A breakdown of docking performance by protein complex type shows that RosettaDock v3.2 achieved docking funnels for 63% of antibody-antigen targets, 62% of enzyme-inhibitor targets, and 35% of ‘other’ targets. In terms of docking difficulty, RosettaDock v3.2 achieved funnels for 58% of rigid-body targets, 30% of medium targets, and 14% of difficult targets. For targets that failed, we carry out additional analyses to identify the cause of failure, which showed that binding-induced backbone conformation changes account for a majority of failures. We also present a bootstrap statistical analysis that quantifies the reliability of the stochastic docking results. Finally, we demonstrate the additional functionality available in RosettaDock v3.2 by incorporating small-molecules and non-protein co-factors in docking of a smaller target set. This study marks the most extensive benchmarking of the RosettaDock module to date and establishes a baseline for future research in protein interface modeling and structure prediction.
Collapse
Affiliation(s)
- Sidhartha Chaudhury
- Program in Molecular Biophysics, Johns Hopkins University, Baltimore, Maryland, United States of America
| | - Monica Berrondo
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland, United States of America
| | - Brian D. Weitzner
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland, United States of America
| | - Pravin Muthu
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland, United States of America
| | - Hannah Bergman
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland, United States of America
| | - Jeffrey J. Gray
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland, United States of America
- * E-mail:
| |
Collapse
|
28
|
Pantazes RJ, Grisewood MJ, Maranas CD. Recent advances in computational protein design. Curr Opin Struct Biol 2011; 21:467-72. [PMID: 21600758 DOI: 10.1016/j.sbi.2011.04.005] [Citation(s) in RCA: 73] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2011] [Accepted: 04/28/2011] [Indexed: 11/30/2022]
Affiliation(s)
- Robert J Pantazes
- The Pennsylvania State University, Department of Chemical Engineering, 112 Fenske Lab, University Park, PA 16802, USA
| | | | | |
Collapse
|
29
|
Abstract
Peptide-protein interactions are prevalent in the living cell and form a key component of the overall protein-protein interaction network. These interactions are drawing increasing interest due to their part in signaling and regulation, and are thus attractive targets for computational structural modeling. Here we report an overview of current techniques for the high resolution modeling of peptide-protein complexes. We dissect this complicated challenge into several smaller subproblems, namely: modeling the receptor protein, predicting the peptide binding site, sampling an initial peptide backbone conformation and the final refinement of the peptide within the receptor binding site. For each of these conceptual stages, we present available tools, approaches, and their reported performance. We summarize with an illustrative example of this process, highlighting the success and current challenges still facing the automated blind modeling of peptide-protein interactions. We believe that the upcoming years will see considerable progress in our ability to create accurate models of peptide-protein interactions, with applications in binding-specificity prediction, rational design of peptide-mediated interactions and the usage of peptides as therapeutic agents.
Collapse
Affiliation(s)
- Nir London
- Department of Microbiology and Molecular Genetics, Institute for Medical Research Israel-Canada, Hadassah Medical School, The Hebrew University, Jerusalem, Israel
| | | | | |
Collapse
|
30
|
Raveh B, London N, Schueler-Furman O. Sub-angstrom modeling of complexes between flexible peptides and globular proteins. Proteins 2010; 78:2029-40. [PMID: 20455260 DOI: 10.1002/prot.22716] [Citation(s) in RCA: 319] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
A wide range of regulatory processes in the cell are mediated by flexible peptides that fold upon binding to globular proteins. Computational efforts to model these interactions are hindered by the large number of rotatable bonds in flexible peptides relative to typical ligand molecules, and the fact that different peptides assume different backbone conformations within the same binding site. In this study, we present Rosetta FlexPepDock, a novel tool for refining coarse peptide-protein models that allows significant changes in both peptide backbone and side chains. We obtain high resolution models, often of sub-angstrom backbone quality, over an extensive and general benchmark that is based on a large nonredundant dataset of 89 peptide-protein interactions. Importantly, side chains of known binding motifs are modeled particularly well, typically with atomic accuracy. In addition, our protocol has improved modeling quality for the important application of cross docking to PDZ domains. We anticipate that the ability to create high resolution models for a wide range of peptide-protein complexes will have significant impact on structure-based functional characterization, controlled manipulation of peptide interactions, and on peptide-based drug design.
Collapse
Affiliation(s)
- Barak Raveh
- Department of Microbiology and Molecular Genetics, Insitute for Medical Research Israel-Canada, Hadassah Medical School, The Hebrew University, Jerusalem, 91120 Israel
| | | | | |
Collapse
|