1
|
Zeng Q, Meng G, Zhao B, Lin H, Guan Y, Qin X, Yuan Y, Li Y, Wang Q. Peptide Drug Design Using Alchemical Free Energy Calculation: An Application and Validation on Agonists of Ghrelin Receptor. J Chem Inf Model 2024; 64:4863-4876. [PMID: 38836743 DOI: 10.1021/acs.jcim.4c00414] [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: 06/06/2024]
Abstract
With recent large-scale applications and validations, the relative binding free energy (RBFE) calculated using alchemical free energy methods has been proven to be an accurate measure to probe the binding of small-molecule drug candidates. On the other hand, given the flexibility of peptides, it is of great interest to find out whether sufficient sampling could be achieved within the typical time scale of such calculation, and a similar level of accuracy could be reached for peptide drugs. However, the systematic evaluation of such calculations on protein-peptide systems has been less reported. Most reported studies of peptides were restricted to a limited number of data points or lacking experimental support. To demonstrate the applicability of the alchemical free energy method for protein-peptide systems in a typical real-world drug discovery project, we report an application of the thermodynamic integration (TI) method to the RBFE calculation of ghrelin receptor and its peptide agonists. Along with the calculation, the synthesis and in vitro EC50 activity of relamorelin and 17 new peptide derivatives were also reported. A cost-effective criterion to determine the data collection time was proposed for peptides in the TI simulation. The average of three TI repeats yielded a mean absolute error of 0.98 kcal/mol and Pearson's correlation coefficient (R) of 0.77 against the experimental free energy derived from the in vitro EC50 activity, showing good repeatability of the proposed method and a slightly better agreement than the results obtained from the arbitrary time frames up to 20 ns. Although it is limited by having one target and a deduced binding pose, we hope that this study can add some insights into alchemical free energy calculation of protein-peptide systems, providing theoretical assistance to the development of peptide drugs.
Collapse
Affiliation(s)
- Qin Zeng
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu 610041, China
| | - Guangpeng Meng
- Chengdu Sintanovo Biotechnology Co., Ltd., Chengdu 610000, China
| | - Bingyu Zhao
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu 610041, China
| | - Haodian Lin
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu 610041, China
| | - Yuqing Guan
- Chengdu Sintanovo Biotechnology Co., Ltd., Chengdu 610000, China
| | - Xiaobin Qin
- Chengdu Sintanovo Biotechnology Co., Ltd., Chengdu 610000, China
| | - Yu Yuan
- Chengdu Sintanovo Biotechnology Co., Ltd., Chengdu 610000, China
| | - Yuanbo Li
- Chengdu Sintanovo Biotechnology Co., Ltd., Chengdu 610000, China
| | - Qiantao Wang
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu 610041, China
| |
Collapse
|
2
|
Zhang L, Li J. Molecular Dynamics Simulations on Spike Protein Mutants Binding with Human β Defensin Type 2. J Phys Chem B 2024; 128:415-428. [PMID: 38189674 DOI: 10.1021/acs.jpcb.3c05460] [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: 01/09/2024]
Abstract
Human β defensin type 2 (hBD-2), a cationic cysteine-rich peptide secreted from the human innate immune system, can bind Spike-RBD at the same site as receptor ACE2, thus blocking viral entry into ACE2-expressing cells. In order to find out the impact of CoV-2 mutations on hBD-2's antiviral activity, it is important to investigate the binding and interaction of hBD-2 with RBD mutants. All-atom molecular dynamics simulations were conducted on typical RBD mutants, including N501Y, E484K, P479S, T478I, S477N, N439K, K417N, and N501Y-E484K-K417N, binding with hBD-2. Starting from the stable binding structure of hBD-2 and wt-RBD and ClusPro and HADDOCK docking-predicted initial structures, the RBD variants bound with hBD-2 simulations were set up, and NAMD simulations were conducted. Based on the structure and dynamics analysis, it was found that most RBD variants can still form a similar number of hydrogen bonds with hBD-2, in addition to having a similar-sized buried surface area (BSA) and a similar binding interface to the RBD wildtype. However, the RBD triple mutant formed a less stable binding structure with hBD-2 than other variants. Additionally, the free energy perturbation (FEP) method was applied to calculate the contribution of key mutant residues to the binding and the free energy change caused by the mutations. The result shows that N439K, K417N, and the trimutation increase the binding free energy of RBD with hBD-2; thus, RBD should bind less stably with hBD-2. E484K decreases the binding free energy, thus it should bind more stably with hBD-2, while N501Y, S477N, T478I, and P479S almost do not change the binding free energy with hBD-2. The MM-GBSA method predicted the binding interaction energy which shows that the trimutant should be able to escape the binding with hBD-2 but N501Y should not. The result can provide insight into understanding the functional mechanism of hBD-2 combating SARS-CoV-2 mutants.
Collapse
Affiliation(s)
- Liqun Zhang
- Chemical Engineering Department, University of Rhode Island, Kingston, Rhode Island 02881, United States
| | - Jadeson Li
- Newton North High School, 457 Walnut Street, Newton, Massachusetts 02460, United States
| |
Collapse
|
3
|
Papadourakis M, Sinenka H, Matricon P, Hénin J, Brannigan G, Pérez-Benito L, Pande V, van Vlijmen H, de Graaf C, Deflorian F, Tresadern G, Cecchini M, Cournia Z. Alchemical Free Energy Calculations on Membrane-Associated Proteins. J Chem Theory Comput 2023; 19:7437-7458. [PMID: 37902715 PMCID: PMC11017255 DOI: 10.1021/acs.jctc.3c00365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Indexed: 10/31/2023]
Abstract
Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.
Collapse
Affiliation(s)
- Michail Papadourakis
- Biomedical
Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| | - Hryhory Sinenka
- Institut
de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg Cedex, France
| | - Pierre Matricon
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Jérôme Hénin
- Laboratoire
de Biochimie Théorique UPR 9080, CNRS and Université Paris Cité, 75005 Paris, France
| | - Grace Brannigan
- Center
for Computational and Integrative Biology, Rutgers University−Camden, Camden, New Jersey 08103, United States of America
- Department
of Physics, Rutgers University−Camden, Camden, New Jersey 08102, United States
of America
| | - Laura Pérez-Benito
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Vineet Pande
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Herman van Vlijmen
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Chris de Graaf
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Francesca Deflorian
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Gary Tresadern
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Marco Cecchini
- Institut
de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg Cedex, France
| | - Zoe Cournia
- Biomedical
Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| |
Collapse
|
4
|
Rucker G, Qin H, Zhang L. Structure, dynamics and free energy studies on the effect of point mutations on SARS-CoV-2 spike protein binding with ACE2 receptor. PLoS One 2023; 18:e0289432. [PMID: 37796794 PMCID: PMC10553274 DOI: 10.1371/journal.pone.0289432] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2023] [Accepted: 09/11/2023] [Indexed: 10/07/2023] Open
Abstract
The ongoing COVID-19 pandemic continues to infect people worldwide, and the virus continues to evolve in significant ways which can pose challenges to the efficiency of available vaccines and therapeutic drugs and cause future pandemic. Therefore, it is important to investigate the binding and interaction of ACE2 with different RBD variants. A comparative study using all-atom MD simulations was conducted on ACE2 binding with 8 different RBD variants, including N501Y, E484K, P479S, T478I, S477N, N439K, K417N and N501Y-E484K-K417N on RBD. Based on the RMSD, RMSF, and DSSP results, overall the binding of RBD variants with ACE2 is stable, and the secondary structure of RBD and ACE2 are consistent after the point mutation. Besides that, a similar buried surface area, a consistent binding interface and a similar amount of hydrogen bonds formed between RBD and ACE2 although the exact residue pairs on the binding interface were modified. The change of binding free energy from point mutation was predicted using the free energy perturbation (FEP) method. It is found that N501Y, N439K, and K417N can strengthen the binding of RBD with ACE2, while E484K and P479S weaken the binding, and S477N and T478I have negligible effect on the binding. Point mutations modified the dynamic correlation of residues in RBD based on the dihedral angle covariance matrix calculation. Doing dynamic network analysis, a common intrinsic network community extending from the tail of RBD to central, then to the binding interface region was found, which could communicate the dynamics in the binding interface region to the tail thus to the other sections of S protein. The result can supply unique methodology and molecular insight on studying the molecular structure and dynamics of possible future pandemics and design novel drugs.
Collapse
Affiliation(s)
- George Rucker
- Chemical Engineering Department, Tennessee Technological University, Cookeville, TN, United States of America
| | - Hong Qin
- Computer Science Department, University of Tennessee Chattanooga, Chattanooga, TN, United States of America
| | - Liqun Zhang
- Chemical Engineering Department, University of Rhode Island, Kingston, RI, United States of America
| |
Collapse
|
5
|
Zheng J, Nurit H, Zhang L, Wei T, Shao Q. Editorial: Computational data-driven design and modeling of biomolecules and biomimetics. Biophys Chem 2022; 289:106877. [PMID: 36007359 DOI: 10.1016/j.bpc.2022.106877] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Affiliation(s)
- Jie Zheng
- Dept. of Chemical, Biomolecular, and Corrosion Engineering, University of Akron, United States of America.
| | - Haspel Nurit
- Dept. of Computer Science, University of Massachusetts Boston, United States of America
| | - Liqun Zhang
- Dept. of Chemical Engineering, Tennessee Technological University, United States of America
| | - Tao Wei
- Dept. of Chemical Engineering, Howard University, United States of America
| | - Qing Shao
- Dept. of Chemical and Materials Engineering, University of Kentucky, United States of America
| |
Collapse
|
6
|
Zhang L. Interaction of Human β Defensin Type 3 (hBD-3) with Different PIP2-Containing Membranes, a Molecular Dynamics Simulation Study. J Chem Inf Model 2021; 61:4670-4686. [PMID: 34473496 DOI: 10.1021/acs.jcim.1c00805] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Human β defensin type 3 (hBD-3) is a cysteine-rich small antibacterial peptide. It belongs to the human innate immune system. hBD-3 has six cysteine residues, which form three pairs of disulfide bonds, and those bonds break in the reducing condition. It is known that hBD-3 can interact with bacterial membrane, and even eukaryotic cell membrane, which has a low concentration of phosphatidylinositol 4,5-bisphosphate (PIP2) lipids. PIP2 is a vital component in cell membranes and has been found to play important roles during antimicrobial peptide (AMP) interaction with membranes. To understand the functional mechanism of hBD-3 interacting with PIP2-containing membranes, the binding structures of hBD-3 on 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers mixed with 10% of PIP2 were predicted using two kinds of methods. The first one is by placing the hBD-3 monomer in different orientations above the POPC + 10%PIP2 membrane to set up five different initial simulation systems and performing long-term simulations on each to predict the most stable binding structure. It was found that hBD-3 analogue binds on the mixed lipid membrane on the two loop regions. The second method is by running long-term simulations on one or nine hBD-3 dimers binding on POPC mixed with 10%PIP2 lipid bilayer starting from the solid-state NMR (ssNMR)-suggested orientation. The dimer dissociated, and the most stable binding of hBD-3 in wild-type on the mixed membrane is also through the two loop regions, which agrees with the prediction from both the first method and the lipid self-assembly result. The PIP2 lipids can form long-lasting hydrogen bonds with positively charged residues such as Arg and Lys on hBD-3, thus forming clusters with hBD-3. As a comparison, hBD-3 dimers binding with a combined bilayer having 1,2-palmitoyl-oleoyl-sn-glycero-3-phosphoserine (POPS) on the upper and POPC on the lower leaflets and the combined POPS + POPC bilayer mixing with 10%PIP2 were also studied. The long-term simulation result shows that hBD-3 can bind with the heads of negatively charged POPS and PIP2 lipids and form hydrogen bonds. The stable binding sites of hBD-3 on PIP2 or POPS mixed bilayers are still on the two loop regions. On the combined POPS + POPC mixed with 10%PIP2 bilayer, the binding of hBD-3 with PIP2 lipids became less stable and fewer because of the competition of binding with the POPS lipids. Besides that, binding with hBD-3 can decrease the membrane thickness of the POPC + PIP2, POPS + POPC, and POPS + POPC + PIP2 bilayers and make POPS and PIP2 lipids more flexible based on the order parameter calculations. Our results supply molecular insight on AMP binding with different membranes and can help understand the functional mechanism of hBD-3 disrupting PIP2-containing membranes.
Collapse
Affiliation(s)
- Liqun Zhang
- Department of Chemical Engineering, Tennessee Technological University, Cookeville, Tennessee 38505, United States
| |
Collapse
|