1
|
Banerjee T, Gosai A, Yousefi N, Garibay OO, Seal S, Balasubramanian G. Examining sialic acid derivatives as potential inhibitors of SARS-CoV-2 spike protein receptor binding domain. J Biomol Struct Dyn 2024; 42:6342-6358. [PMID: 37424217 DOI: 10.1080/07391102.2023.2234044] [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: 05/02/2023] [Accepted: 07/01/2023] [Indexed: 07/11/2023]
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS CoV-2) has been the primary reason behind the COVID-19 global pandemic which has affected millions of lives worldwide. The fundamental cause of the infection is the molecular binding of the viral spike protein receptor binding domain (SP-RBD) with the human cell angiotensin-converting enzyme 2 (ACE2) receptor. The infection can be prevented if the binding of RBD-ACE2 is resisted by utilizing certain inhibitors or drugs that demonstrate strong binding affinity towards the SP RBD. Sialic acid based glycans found widely in human cells and tissues have notable propensity of binding to viral proteins of the coronaviridae family. Recent experimental literature have used N-acetyl neuraminic acid (Sialic acid) to create diagnostic sensors for SARS-CoV-2, but a detailed interrogation of the underlying molecular mechanisms is warranted. Here, we perform all atom molecular dynamics (MD) simulations for the complexes of certain Sialic acid-based molecules with that of SP RBD of SARS CoV-2. Our results indicate that Sialic acid not only reproduces a binding affinity comparable to the RBD-ACE2 interactions, it also assumes the longest time to dissociate completely from the protein binding pocket of SP RBD. Our predictions corroborate that a combination of electrostatic and van der Waals energies as well the polar hydrogen bond interactions between the RBD residues and the inhibitors influence free energy of binding.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Tanumoy Banerjee
- Department of Mechanical Engineering and Mechanics, Lehigh University, Bethlehem, PA, USA
| | | | - Niloofar Yousefi
- Industrial Engineering and Management Systems, University of Central Florida, Orlando, FL, USA
| | - Ozlem Ozmen Garibay
- Industrial Engineering and Management Systems, University of Central Florida, Orlando, FL, USA
| | - Sudipta Seal
- College of Medicine, Bionix Cluster, University of Central Florida, Orlando, FL, USA
- Advanced Materials Processing and Analysis Center, Dept. of Materials Science and Engineering, University of Central Florida, Orlando, FL, USA
| | - Ganesh Balasubramanian
- Department of Mechanical Engineering and Mechanics, Lehigh University, Bethlehem, PA, USA
- Institute of Functional Materials & Devices and College of Health, Lehigh University, Bethlehem, PA, USA
| |
Collapse
|
2
|
Ciura P, Smardz P, Spodzieja M, Sieradzan AK, Krupa P. Multilayered Computational Framework for Designing Peptide Inhibitors of HVEM-LIGHT Interaction. J Phys Chem B 2024; 128:6770-6785. [PMID: 38958133 DOI: 10.1021/acs.jpcb.4c02255] [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: 07/04/2024]
Abstract
The herpesvirus entry mediator (HVEM) and its ligand LIGHT play crucial roles in immune system regulation, including T-cell proliferation, B-cell differentiation, and immunoglobulin secretion. However, excessive T-cell activation can lead to chronic inflammation and autoimmune diseases. Thus, inhibiting the HVEM-LIGHT interaction emerges as a promising therapeutic strategy for these conditions and in preventing adverse reactions in organ transplantation. This study focused on designing peptide inhibitors, targeting the HVEM-LIGHT interaction, using molecular dynamics (MD) simulations of 65 peptides derived from HVEM. These peptides varied in length and disulfide-bond configurations, crucial for their interaction with the LIGHT trimer. By simulating 31 HVEM domain variants, including the full-length protein, we assessed conformational changes upon LIGHT binding to understand the influence of HVEM segments and disulfide bonds on the binding mechanism. Employing multitrajectory microsecond-scale, all-atom MD simulations and molecular mechanics with generalized Born and surface area (MM-GBSA) binding energy estimation, we identified promising CRD2 domain variants with high LIGHT affinity. Notably, point mutations in these variants led to a peptide with a single disulfide bond (C58-C73) and a K54E substitution, exhibiting the highest binding affinity. The importance of the CRD2 domain and Cys58-Cys73 disulfide bond for interrupting HVEM-LIGHT interaction was further supported by analyzing truncated CRD2 variants, demonstrating similar binding strengths and mechanisms. Further investigations into the binding mechanism utilized steered MD simulations at various pulling speeds and umbrella sampling to estimate the energy profile of HVEM-based inhibitors with LIGHT. These comprehensive analyses revealed key interactions and different binding mechanisms, highlighting the increased binding affinity of selected peptide variants. Experimental circular dichroism techniques confirmed the structural properties of these variants. This study not only advances our understanding of the molecular basis of HVEM-LIGHT interactions but also provides a foundation for developing novel therapeutic strategies for immune-related disorders. Furthermore, it sets a gold standard for peptide inhibitor design in drug development due to its systematic approach.
Collapse
Affiliation(s)
- Piotr Ciura
- Faculty of Chemistry, Fahrenheit Union of Universities in Gdańsk, University of Gdańsk, Baż̇yńskiego 8, 80-309 Gdansḱ, Poland
| | - Pamela Smardz
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland
| | - Marta Spodzieja
- Faculty of Chemistry, Fahrenheit Union of Universities in Gdańsk, University of Gdańsk, Baż̇yńskiego 8, 80-309 Gdansḱ, Poland
| | - Adam K Sieradzan
- Faculty of Chemistry, Fahrenheit Union of Universities in Gdańsk, University of Gdańsk, Baż̇yńskiego 8, 80-309 Gdansḱ, Poland
| | - Pawel Krupa
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
3
|
Liu T, Simine L. DeltaGzip: Computing Biopolymer-Ligand Binding Affinity via Kolmogorov Complexity and Lossless Compression. J Chem Inf Model 2024. [PMID: 38980667 DOI: 10.1021/acs.jcim.4c00461] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
The design of biosequences for biosensing and therapeutics is a challenging multistep search and optimization task. In principle, computational modeling may speed up the design process by virtual screening of sequences based on their binding affinities to target molecules. However, in practice, existing machine-learned models trained to predict binding affinities lack the flexibility with respect to reaction conditions, and molecular dynamics simulations that can incorporate reaction conditions suffer from high computational costs. Here, we describe a computational approach called DeltaGzip that evaluates the free energy of binding in biopolymer-ligand complexes from ultrashort equilibrium molecular dynamics simulations. The entropy of binding is evaluated using the Kolmogorov complexity definition of entropy and approximated using a lossless compression algorithm, Gzip. We benchmark the method on a well-studied data set of protein-ligand complexes comparing the predictions of DeltaGzip to the free energies of binding obtained using Jarzynski equality and experimental measurements.
Collapse
Affiliation(s)
- Tao Liu
- Department of Chemistry, McGill University, Montreal, Quebec H3A 0B8, Canada
| | - Lena Simine
- Department of Chemistry, McGill University, Montreal, Quebec H3A 0B8, Canada
| |
Collapse
|
4
|
Truong DT, Ho K, Nhi HTY, Nguyen VH, Dang TT, Nguyen MT. Imidazole[1,5-a]pyridine derivatives as EGFR tyrosine kinase inhibitors unraveled by umbrella sampling and steered molecular dynamics simulations. Sci Rep 2024; 14:12218. [PMID: 38806555 PMCID: PMC11133355 DOI: 10.1038/s41598-024-62743-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2024] [Accepted: 05/21/2024] [Indexed: 05/30/2024] Open
Abstract
Although the use of the tyrosine kinase inhibitors (TKIs) has been proved that it can save live in a cancer treatment, the currently used drugs bring in many undesirable side-effects. Therefore, the search for new drugs and an evaluation of their efficiency are intensively carried out. Recently, a series of eighteen imidazole[1,5-a]pyridine derivatives were synthetized by us, and preliminary analyses pointed out their potential to be an important platform for pharmaceutical development owing to their promising actions as anticancer agents and enzyme (kinase, HIV-protease,…) inhibitors. In the present theoretical study, we further analyzed their efficiency in using a realistic scenario of computational drug design. Our protocol has been developed to not only observe the atomistic interaction between the EGFR protein and our 18 novel compounds using both umbrella sampling and steered molecular dynamics simulations, but also determine their absolute binding free energies. Calculated properties of the 18 novel compounds were in detail compared with those of two known drugs, erlotinib and osimertinib, currently used in cancer treatment. Inspiringly the simulation results promote three imidazole[1,5-a]pyridine derivatives as promising inhibitors into a further step of clinical trials.
Collapse
Affiliation(s)
- Duc Toan Truong
- Laboratory for Chemical Computation and Modeling, Institute for Computational Science and Artificial Intelligence, Van Lang University, Ho Chi Minh City, 70000, Vietnam
- Faculty of Applied Technology, School of Technology, Van Lang University, Ho Chi Minh City, 70000, Vietnam
| | - Kiet Ho
- Institute for Computational Science and Technology (ICST), Quang Trung Software City, Ho Chi Minh City, 70000, Vietnam
| | - Huynh Thi Yen Nhi
- Laboratory for Chemical Computation and Modeling, Institute for Computational Science and Artificial Intelligence, Van Lang University, Ho Chi Minh City, 70000, Vietnam
| | - Van Ha Nguyen
- Faculty of Chemistry, VNU University of Science, Vietnam National University, Hanoi, 19 Le Thanh Tong, Hoan Kiem, Hanoi, 11021, Vietnam
| | - Tuan Thanh Dang
- Faculty of Chemistry, VNU University of Science, Vietnam National University, Hanoi, 19 Le Thanh Tong, Hoan Kiem, Hanoi, 11021, Vietnam
| | - Minh Tho Nguyen
- Laboratory for Chemical Computation and Modeling, Institute for Computational Science and Artificial Intelligence, Van Lang University, Ho Chi Minh City, 70000, Vietnam.
- Faculty of Applied Technology, School of Technology, Van Lang University, Ho Chi Minh City, 70000, Vietnam.
| |
Collapse
|
5
|
Ravichandran A, Das R. The Thermodynamic Properties of Fat10ylated Proteins Are Regulated by the Fat10ylation Site. ACS OMEGA 2024; 9:22265-22276. [PMID: 38799324 PMCID: PMC11112694 DOI: 10.1021/acsomega.4c01396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/13/2024] [Revised: 04/24/2024] [Accepted: 04/29/2024] [Indexed: 05/29/2024]
Abstract
Degradation of proteins by the proteasome is crucial in regulating their levels in the cell. Post-translational modifications, such as ubiquitylation and Fat10ylation, trigger proteasomal degradation of the substrate proteins. While ubiquitylation regulates multiple cellular pathways, Fat10ylation functions explicitly in the inflammatory response pathway. At the proteasome, ubiquitin is recycled after being cleaved from the substrate, while Fat10 is degraded simultaneously with its substrate. Although the thermodynamic properties of the substrate are critical for effective proteasomal degradation, they remain poorly understood for the Fat10-proteasome pathway. We studied the thermodynamic properties of the Fat10∼substrate conjugate to uncover mechanistic details of the pathway. First, the mechanical unfolding of Fat10∼substrate was studied by molecular dynamics simulations, which suggested that the unfolding pathway and unfolding energy of the substrate depend on the site of Fat10 modification. We also investigated different pathways for the entry of the Fat10∼substrate into the proteasome core. Our analysis supports a model where the entry of Fat10, followed by the substrate, is the energetically preferred pathway. Further, we studied Fat10's effect on the thermodynamic properties of distinct substrates, considering their size, flexibility, and surface properties. The results uncovered significant entropic destabilization of substrates due to Fat10ylation, particularly in smaller substrates. For larger substrates, multi-monoFat10ylation is necessary to induce destabilization. Our study further reveals that Fat10 modification at negative patches on substrate surfaces is essential for optimal destabilization and subsequent degradation. These findings provide atomistic insights into the degradation mechanisms in the Fat10 proteasome pathway with potential implications for therapeutic interventions.
Collapse
Affiliation(s)
- Aravind Ravichandran
- National
Center for Biological Sciences, Tata Institute
of Fundamental Research, Bangalore 560065, India
- SASTRA
University, Thirumalaisamudram, Thanjavur 613401, India
| | - Ranabir Das
- National
Center for Biological Sciences, Tata Institute
of Fundamental Research, Bangalore 560065, India
| |
Collapse
|
6
|
Truong DT, Ho K, Pham DQH, Chwastyk M, Nguyen-Minh T, Nguyen MT. Treatment of flexibility of protein backbone in simulations of protein-ligand interactions using steered molecular dynamics. Sci Rep 2024; 14:10475. [PMID: 38714683 PMCID: PMC11076533 DOI: 10.1038/s41598-024-59899-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2024] [Accepted: 04/16/2024] [Indexed: 05/10/2024] Open
Abstract
To ensure that an external force can break the interaction between a protein and a ligand, the steered molecular dynamics simulation requires a harmonic restrained potential applied to the protein backbone. A usual practice is that all or a certain number of protein's heavy atoms or Cα atoms are fixed, being restrained by a small force. This present study reveals that while fixing both either all heavy atoms and or all Cα atoms is not a good approach, while fixing a too small number of few atoms sometimes cannot prevent the protein from rotating under the influence of the bulk water layer, and the pulled molecule may smack into the wall of the active site. We found that restraining the Cα atoms under certain conditions is more relevant. Thus, we would propose an alternative solution in which only the Cα atoms of the protein at a distance larger than 1.2 nm from the ligand are restrained. A more flexible, but not too flexible, protein will be expected to lead to a more natural release of the ligand.
Collapse
Affiliation(s)
- Duc Toan Truong
- Laboratory for Chemical Computation and Modeling, Institute for Computational Science and Artificial Intelligence, Van Lang University, Ho Chi Minh City, 70000, Vietnam
- Faculty of Applied Technology, School of Technology, Van Lang University, Ho Chi Minh City, 70000, Vietnam
| | - Kiet Ho
- Institute for Computational Science and Technology (ICST), Quang Trung Software City, Ho Chi Minh City, 70000, Vietnam
| | | | - Mateusz Chwastyk
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland
| | - Thai Nguyen-Minh
- University of Medicine and Pharmacy at Ho Chi Minh City, Ho Chi Minh City, 70000, Vietnam
| | - Minh Tho Nguyen
- Laboratory for Chemical Computation and Modeling, Institute for Computational Science and Artificial Intelligence, Van Lang University, Ho Chi Minh City, 70000, Vietnam.
- Faculty of Applied Technology, School of Technology, Van Lang University, Ho Chi Minh City, 70000, Vietnam.
| |
Collapse
|
7
|
Iida S, Kameda T. Dissociation Rate Calculation via Constant-Force Steered Molecular Dynamics Simulation. J Chem Inf Model 2023. [PMID: 37188657 DOI: 10.1021/acs.jcim.2c01529] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
Steered molecular dynamics (SMD) simulations are used to study molecular dissociation events by applying a harmonic force to the molecules and pulling them at a constant velocity. Instead of constant-velocity pulling, we use a constant force: the constant-force SMD (CF-SMD) simulation. The CF-SMD simulation employs a constant force to reduce the activation barrier of molecular dissociation, thereby enhancing the dissociation event. Here, we present the capability of the CF-SMD simulation to estimate the dissociation time at equilibrium. We performed all-atom CF-SMD simulations for NaCl and protein-ligand systems, producing dissociation time at various forces. We extrapolated these values to the dissociation rate without a constant force using Bell's model or the Dudko-Hummer-Szabo model. We demonstrate that the CF-SMD simulations with the models predicted the dissociation time in equilibrium. A CF-SMD simulation is a powerful tool for estimating the dissociation rate in a direct and computationally efficient manner.
Collapse
Affiliation(s)
- Shinji Iida
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, Tokyo 135-0064, Japan
| | - Tomoshi Kameda
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, Tokyo 135-0064, Japan
| |
Collapse
|
8
|
Hepatoprotective Effect of Millettia dielsiana: In Vitro and In Silico Study. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27248978. [PMID: 36558112 PMCID: PMC9787921 DOI: 10.3390/molecules27248978] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/26/2022] [Revised: 12/07/2022] [Accepted: 12/08/2022] [Indexed: 12/24/2022]
Abstract
In silico docking studies of 50 selected compounds from Millettia dielsiana Harms ex Diels (family Leguminosae) were docked into the binding pocket of the PI3K/mTOR protein. In there, compounds trans-3-O-p-hydroxycinnamoyl ursolic acid (1) and 5,7,4'-trihydroxyisoflavone 7-O-β-D-apiofuranosyl-(1→6)-β-D-glucopyranoside (2) are predicted to be very promising inhibitors against PI3K/mTOR. They direct their cytotoxic activity against Hepatocellular carcinoma with binding affinity (BA) values, the pulling work spent to the co-crystallized ligand from the binding site of PI3K/mTOR (W and Fmax), and the non-equilibrium binding free energy (∆GneqJar) as BA values = -9.237 and -9.083 kcal/mol, W = 83.5 ± 10.6 kcal/mol with Fmax = 336.2 ± 45.3 pN and 126.6 ± 21.7 kcal/mol with Fmax = 430.3 ± 84.0 pN, and ∆GneqJar = -69.86074 and -101.2317 kcal/mol, respectively. In molecular dynamic simulation, the RMSD value of the PI3K/mTOR complex with compounds (1 and 2) was in the range of 0.3 nm to the end of the simulation. Therefore, the compounds (1 and 2) are predicted to be very promising inhibitors against PI3K/mTOR. The crude extract, ethyl acetate fraction and compounds (1 and 2) from Millettia dielsiana exhibited moderate to potent in vitro cytotoxicity on Hepatocellular carcinoma cell line with IC50 values of 81.2 µg/mL, 60.4 µg/mL, 23.1 μM, and 16.3 μM, respectively, and showed relatively potent to potent in vitro antioxidant activity on mouse hepatocytes with ED50 values of 24.4 µg/mL, 19.3 µg/mL, 30.7 μM, and 20.5 μM, respectively. In conclusion, Millettia dielsiana and compounds (1 and 2) are predicted to have very promising cytotoxic activity against Hepatocellular carcinoma and have a hepatoprotective effect.
Collapse
|
9
|
Iida S, Tomoshi K. Free energy and kinetic rate calculation via non-equilibrium molecular simulation: application to biomolecules. Biophys Rev 2022; 14:1303-1314. [PMID: 36659997 PMCID: PMC9842846 DOI: 10.1007/s12551-022-01036-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Accepted: 11/26/2022] [Indexed: 12/30/2022] Open
Abstract
Non-equilibrium molecular dynamics (NEMD) simulation has been recognized as a powerful tool for examining biomolecules and provides fruitful insights into not only non-equilibrium but also equilibrium processes. We review recent advances in NEMD simulation and relevant, fundamental results of non-equilibrium statistical mechanics. We first introduce Crooks fluctuation theorem and Jarzynski equality that relate free energy difference to work done on a physical system during a non-equilibrium process. The theorems are beneficial for the analysis of NEMD trajectories. We then describe rate theory, a framework to calculate molecular kinetics from a non-equilibrium process; this theoretical framework enables us to calculate a reaction time-mean-first passage time-from NEMD trajectories. We, in turn, present recent NEMD techniques that apply an external force to a system to enhance molecular dissociation and introduce their application to biomolecules. Lastly, we show the current status of an appropriate selection of reaction coordinates for NEMD simulation.
Collapse
Affiliation(s)
- Shinji Iida
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-Ku, Tokyo, 135-0064 Japan
| | - Kameda Tomoshi
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-Ku, Tokyo, 135-0064 Japan
| |
Collapse
|
10
|
Thu TTM, Li MS. Protein aggregation rate depends on mechanical stability of fibrillar structure . J Chem Phys 2022; 157:055101. [DOI: 10.1063/5.0088689] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The formation of the fibrillar structure of amyloid proteins/peptides is believed to be associated with neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, etc. Since the rate of aggregation can influence neurotoxicity, finding the key factors that control this rate is of paramount importance. It was recently found evidence that the rate of protein aggregation is related to the mechanical stability of the fibrillar structure, such that the higher the mechanical stability, the faster the fibril is formed. However, this conclusion was supported by a limited dataset. In this work, we expand the previous study to a larger dataset, including the wild type of Aβ42 peptide and its 20 mutants, the aggregation rate of which was measured experimentally. By using all-atom steered molecular dynamics (SMD) simulations we can access the mechanical stability of the fibril structure, which is characterized by the rupture force, pulling work and unbinding free energy barrier. Our result confirms that mechanical stability is indeed related to the aggregation rate. Since estimation of the aggregation rate using all-atom simulations is almost forbidden by the current computational capabilities, our result is useful for predicting it based on information obtained from fast SMD simulations for fibrils.
Collapse
Affiliation(s)
| | - Mai Suan Li
- Theoretical Physics, Institute of Physics, Polish Academy of Sciences, Poland
| |
Collapse
|
11
|
Nguyen H, Li MS. Antibody-nanobody combination increases their neutralizing activity against SARS-CoV-2 and nanobody H11-H4 is effective against Alpha, Kappa and Delta variants. Sci Rep 2022; 12:9701. [PMID: 35690632 PMCID: PMC9188278 DOI: 10.1038/s41598-022-14263-1] [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: 03/02/2022] [Accepted: 06/03/2022] [Indexed: 11/14/2022] Open
Abstract
The global spread of COVID-19 is devastating health systems and economies worldwide. While the use of vaccines has yielded encouraging results, the emergence of new variants of SARS-CoV-2 shows that combating COVID-19 remains a big challenge. One of the most promising treatments is the use of not only antibodies, but also nanobodies. Recent experimental studies revealed that the combination of antibody and nanobody can significantly improve their neutralizing ability through binding to the SARS-CoV-2 spike protein, but the molecular mechanisms underlying this observation remain largely unknown. In this work, we investigated the binding affinity of the CR3022 antibody and H11-H4 nanobody to the SARS-CoV-2 receptor binding domain (RBD) using molecular modeling. Both all-atom steered molecular dynamics simulations and coarse-grained umbrella sampling showed that, consistent with the experiment, CR3022 associates with RBD more strongly than H11-H4. We predict that the combination of CR3022 and H11-H4 considerably increases their binding affinity to the spike protein. The electrostatic interaction was found to control the association strength of CR3022, but the van der Waals interaction dominates in the case of H11-H4. However, our study for a larger set of nanobodies and antibodies showed that the relative role of these interactions depends on the specific complex. Importantly, we showed Beta, Gamma, Lambda, and Mu variants reduce the H11-H4 activity while Alpha, Kappa and Delta variants increase its neutralizing ability, which is in line with experiment reporting that the nanobody elicited from the llama is very promising for fighting against the Delta variant.
Collapse
Affiliation(s)
- Hung Nguyen
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668, Warsaw, Poland
| | - Mai Suan Li
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668, Warsaw, Poland. .,Life Science Lab, Institute for Computational Science and Technology, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam.
| |
Collapse
|
12
|
Reif MM, Zacharias M. Improving the Potential of Mean Force and Nonequilibrium Pulling Simulations by Simultaneous Alchemical Modifications. J Chem Theory Comput 2022; 18:3873-3893. [PMID: 35653503 DOI: 10.1021/acs.jctc.1c01194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present an approach combining alchemical modifications and physical-pathway methods to calculate absolute binding free energies. The employed physical-pathway method is either a stratified umbrella sampling to calculate a potential of mean force or nonequilibrium pulling. We devised two basic approaches: the simultaneous approach (S-approach), where, along the physical unbinding pathway, an alchemical transformation of ligand-protein interactions is installed and deinstalled, and the prior-plus-simultaneous approach (PPS-approach), where, prior to the physical-pathway simulation, an alchemical transformation of ligand-protein interactions is installed in the binding site and deinstalled during the physical-pathway simulation. Using a mutant of T4 lysozyme with a benzene ligand as an example, we show that installation and deinstallation of soft-core interactions concurrent with physical ligand unbinding (S-approach) allow successful potential of mean force calculations and nonequilibrium pulling simulations despite the problems posed by the occluded nature of the lysozyme binding pocket. Good agreement between the potential of the mean-force-based S-approach and double decoupling simulations as well as a remarkable efficiency and accuracy of the nonequilibrium-pulling-based S-approach is found. The latter turned out to be more compute-efficient than the potential of mean force calculation by approximately 70%. Furthermore, we illustrate the merits of reducing ligand-protein interactions prior to potential of mean force calculations using the murine double minute homologue protein MDM2 with a p53-derived peptide ligand (PPS-approach). Here, the problem of breaking strong interactions in the binding pocket is transferred to a prior alchemical transformation that reduces the free-energy barrier between the bound and unbound state in the potential of mean force. Besides, disentangling physical ligand displacement from the deinstallation of ligand-protein interactions was seen to allow a more uniform sampling of distance histograms in the umbrella sampling. In the future, physical ligand unbinding combined with simultaneous alchemical modifications may prove useful in the calculation of protein-protein binding free energies, where sampling problems posed by multiple, possibly sticky interactions and potential steric clashes can thus be reduced.
Collapse
Affiliation(s)
- Maria M Reif
- Center for Protein Assemblies (CPA), Physics Department, Chair of Theoretical Biophysics (T38), Technical University of Munich, Ernst-Otto-Fischer-Str. 8, Garching 85748, Germany
| | - Martin Zacharias
- Center for Protein Assemblies (CPA), Physics Department, Chair of Theoretical Biophysics (T38), Technical University of Munich, Ernst-Otto-Fischer-Str. 8, Garching 85748, Germany
| |
Collapse
|
13
|
Nguyen HL, Thai NQ, Li MS. Determination of Multidirectional Pathways for Ligand Release from the Receptor: A New Approach Based on Differential Evolution. J Chem Theory Comput 2022; 18:3860-3872. [PMID: 35512104 PMCID: PMC9202309 DOI: 10.1021/acs.jctc.1c01158] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
![]()
Steered molecular
dynamics (SMD) simulation is a powerful method
in computer-aided drug design as it can be used to access the relative
binding affinity with high precision but with low computational cost.
The success of SMD depends on the choice of the direction along which
the ligand is pulled from the receptor-binding site. In most simulations,
the unidirectional pathway was used, but in some cases, this choice
resulted in the ligand colliding with the complex surface of the exit
tunnel. To overcome this difficulty, several variants of SMD with
multidirectional pulling have been proposed, but they are not completely
devoid of disadvantages. Here, we have proposed to determine the direction
of pulling with a simple scoring function that minimizes the receptor–ligand
interaction, and an optimization algorithm called differential evolution
is used for energy minimization. The effectiveness of our protocol
was demonstrated by finding expulsion pathways of Huperzine A and
camphor from the binding site of Torpedo California acetylcholinesterase
and P450cam proteins, respectively, and comparing them with the previous
results obtained using memetic sampling and random acceleration molecular
dynamics. In addition, by applying this protocol to a set of ligands
bound with LSD1 (lysine specific demethylase 1), we obtained a much
higher correlation between the work of pulling force and experimental
data on the inhibition constant IC50 compared to that obtained using
the unidirectional approach based on minimal steric hindrance.
Collapse
Affiliation(s)
- Hoang Linh Nguyen
- Life Science Lab, Institute for Computational Science and Technology, QuangTrung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City 729110, Vietnam.,Ho Chi Minh City University of Technology (HCMUT), Ho Chi Minh City 740500, Vietnam.,Vietnam National University, Ho Chi Minh City 71300, Vietnam
| | - Nguyen Quoc Thai
- Life Science Lab, Institute for Computational Science and Technology, QuangTrung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City 729110, Vietnam.,Dong Thap University, 783 Pham Huu Lau Street, Ward 6, Cao Lanh City, Dong Thap 81100, Vietnam
| | - Mai Suan Li
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, Warsaw 02-668, Poland
| |
Collapse
|
14
|
Nguyen H, Lan PD, Nissley DA, O’Brien EP, Li MS. Cocktail of REGN Antibodies Binds More Strongly to SARS-CoV-2 Than Its Components, but the Omicron Variant Reduces Its Neutralizing Ability. J Phys Chem B 2022; 126:2812-2823. [PMID: 35403431 PMCID: PMC9016775 DOI: 10.1021/acs.jpcb.2c00708] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
A promising approach to combat Covid-19 infections is the development of effective antiviral antibodies that target the SARS-CoV-2 spike protein. Understanding the structures and molecular mechanisms underlying the binding of antibodies to SARS-CoV-2 can contribute to quickly achieving this goal. Recently, a cocktail of REGN10987 and REGN10933 antibodies was shown to be an excellent candidate for the treatment of Covid-19. Here, using all-atom steered molecular dynamics and coarse-grained umbrella sampling, we examine the interactions of the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein with REGN10987 and REGN10933 separately as well as together. Both computational methods show that REGN10933 binds to RBD more strongly than REGN10987. Importantly, the cocktail binds to RBD (simultaneous binding) more strongly than its components. The dissociation constants of REGN10987-RBD and REGN10933-RBD complexes calculated from the coarse-grained simulations are in good agreement with the experimental data. Thus, REGN10933 is probably a better candidate for treating Covid-19 than REGN10987, although the cocktail appears to neutralize the virus more efficiently than REGN10933 or REGN10987 alone. The association of REGN10987 with RBD is driven by van der Waals interactions, while electrostatic interactions dominate in the case of REGN10933 and the cocktail. We also studied the effectiveness of these antibodies on the two most dangerous variants Delta and Omicron. Consistent with recent experimental reports, our results confirmed that the Omicron variant reduces the neutralizing activity of REGN10933, REGN10987, and REGN10933+REGN10987 with the K417N, N440K, L484A, and Q498R mutations playing a decisive role, while the Delta variant slightly changes their activity.
Collapse
Affiliation(s)
- Hung Nguyen
- Institute
of Physics, Polish Academy of Sciences, al. Lotnikow 32/46, 02-668 Warsaw, Poland
| | - Pham Dang Lan
- Life
Science Lab, Institute for Computational
Science and Technology, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, 729110 Ho Chi Minh City, Vietnam,Faculty
of Physics and Engineering Physics, VNUHCM-University
of Science, 227, Nguyen
Van Cu Street, District 5, 749000 Ho Chi Minh City, Vietnam
| | - Daniel A. Nissley
- Department
of Statistics, University of Oxford, Oxford
Protein Bioinformatics Group, OX1 2JD Oxford, United Kingdom
| | - Edward P. O’Brien
- Department
of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States,Bioinformatics
and Genomics Graduate Program, The Huck Institutes of the Life Sciences, Penn State University, University Park, Pennsylvania 16802, United States,Institute
for Computational and Data Sciences, Penn
State University, University Park, Pennsylvania 16802, United States
| | - Mai Suan Li
- Institute
of Physics, Polish Academy of Sciences, al. Lotnikow 32/46, 02-668 Warsaw, Poland,
| |
Collapse
|
15
|
Pham HA, Truong DT, Li MS. Dependence of Work on the Pulling Speed in Mechanical Ligand Unbinding. J Phys Chem B 2021; 125:8325-8330. [PMID: 34292743 PMCID: PMC8389893 DOI: 10.1021/acs.jpcb.1c01818] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In single-molecule force spectroscopy, the rupture force Fmax required for mechanical unfolding of a biomolecule or for pulling a ligand out of a binding site depends on the pulling speed V and, in the linear Bell-Evans regime, Fmax ∼ ln(V). Recently, it has been found that non-equilibrium work W is better than Fmax in describing relative ligand binding affinity, but the dependence of W on V remains unknown. In this paper, we developed an analytical theory showing that in the linear regime, W ∼ c1 ln(V) + c2 ln2(V), where c1 and c2 are constants. This quadratic dependence was also confirmed by all-atom steered molecular dynamics simulations of protein-ligand complexes. Although our theory was developed for ligand unbinding, it is also applicable to other processes, such as mechanical unfolding of proteins and other biomolecules, due to its universality.
Collapse
Affiliation(s)
- Hong An Pham
- Institute for Computational Science and Technology, QuangTrung Software City, Tan Chanh, Hiep Ward, District 12, Ho Chi Minh City 700000, Vietnam
| | - Duc Toan Truong
- Institute for Computational Science and Technology, QuangTrung Software City, Tan Chanh, Hiep Ward, District 12, Ho Chi Minh City 700000, Vietnam
| | - Mai Suan Li
- Institute of Physics, Polish Academy Science, Al. Lotnikow 32/46, Warsaw 02-668, Poland
| |
Collapse
|
16
|
Nguyen H, Lan PD, Nissley DA, O’Brien EP, Li MS. Electrostatic Interactions Explain the Higher Binding Affinity of the CR3022 Antibody for SARS-CoV-2 than the 4A8 Antibody. J Phys Chem B 2021; 125:7368-7379. [PMID: 34228472 PMCID: PMC8276604 DOI: 10.1021/acs.jpcb.1c03639] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2021] [Revised: 05/30/2021] [Indexed: 12/23/2022]
Abstract
A structural understanding of the mechanism by which antibodies bind SARS-CoV-2 at the atomic level is highly desirable as it can tell the development of more effective antibodies to treat Covid-19. Here, we use steered molecular dynamics (SMD) and coarse-grained simulations to estimate the binding affinity of the monoclonal antibodies CR3022 and 4A8 to the SARS-CoV-2 receptor-binding domain (RBD) and SARS-CoV-2 N-terminal domain (NTD). Consistent with experiments, our SMD and coarse-grained simulations both indicate that CR3022 has a higher affinity for SARS-CoV-2 RBD than 4A8 for the NTD, and the coarse-grained simulations indicate the former binds three times stronger to its respective epitope. This finding shows that CR3022 is a candidate for Covid-19 therapy and is likely a better choice than 4A8. Energetic decomposition of the interaction energies between these two complexes reveals that electrostatic interactions explain the difference in the observed binding affinity between the two complexes. This result could lead to a new approach for developing anti-Covid-19 antibodies in which good candidates must contain charged amino acids in the area of contact with the virus.
Collapse
Affiliation(s)
- Hung Nguyen
- Institute
of Physics, Polish Academy of Sciences, al. Lotnikow 32/46, 02-668 Warsaw, Poland
| | - Pham Dang Lan
- Life
Science Lab, Institute for Computational Science and Technology, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam
- Faculty
of Physics and Engineering Physics, VNUHCM-University
of Science, 227, Nguyen
Van Cu Street, District 5, Ho Chi Minh City, Vietnam
| | - Daniel A. Nissley
- Department
of Statistics, University of Oxford, Oxford
Protein Bioinformatics Group, Oxford OX1 2JD, United Kingdom
| | - Edward P. O’Brien
- Department
of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States
- Bioinformatics
and Genomics Graduate Program, The Huck
Institutes of the Life Sciences, Penn State University, University Park, Pennsylvania 16802, United States
- Institute
for Computational and Data Sciences, Penn
State University, University Park, Pennsylvania 16802, United States
| | - Mai Suan Li
- Institute
of Physics, Polish Academy of Sciences, al. Lotnikow 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
17
|
Nguyen H, Thai NQ, Truong DT, Li MS. Remdesivir Strongly Binds to Both RNA-Dependent RNA Polymerase and Main Protease of SARS-CoV-2: Evidence from Molecular Simulations. J Phys Chem B 2020; 124:11337-11348. [PMID: 33264025 PMCID: PMC7724981 DOI: 10.1021/acs.jpcb.0c07312] [Citation(s) in RCA: 56] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2020] [Revised: 11/12/2020] [Indexed: 02/07/2023]
Abstract
The outbreak of a new coronavirus SARS-CoV-2 (severe acute respiratory syndrome-coronavirus 2) has caused a global COVID-19 (coronavirus disease 2019) pandemic, resulting in millions of infections and thousands of deaths around the world. There is currently no drug or vaccine for COVID-19, but it has been revealed that some commercially available drugs are promising, at least for treating symptoms. Among them, remdesivir, which can block the activity of RNA-dependent RNA polymerase (RdRp) in old SARS-CoV and MERS-CoV viruses, has been prescribed to COVID-19 patients in many countries. A recent experiment showed that remdesivir binds to SARS-CoV-2 with an inhibition constant of μM, but the exact target has not been reported. In this work, combining molecular docking, steered molecular dynamics, and umbrella sampling, we examined its binding affinity to two targets including the main protease (Mpro), also known as 3C-like protease, and RdRp. We showed that remdesivir binds to Mpro slightly weaker than to RdRp, and the corresponding inhibition constants, consistent with the experiment, fall to the μM range. The binding mechanisms of remdesivir to two targets differ in that the electrostatic interaction is the main force in stabilizing the RdRp-remdesivir complex, while the van der Waals interaction dominates in the Mpro-remdesivir case. Our result indicates that remdesivir can target not only RdRp but also Mpro, which can be invoked to explain why this drug is effective in treating COVID-19. We have identified residues of the target protein that make the most important contribution to binding affinity, and this information is useful for drug development for this disease.
Collapse
Affiliation(s)
- Hoang
Linh Nguyen
- Life Science Lab, Institute for Computational Science and Technology, Quang Trung Software City, Tan Chanh
Hiep Ward, District 12, Ho Chi Minh City 700000, Vietnam
| | - Nguyen Quoc Thai
- Life Science Lab, Institute for Computational Science and Technology, Quang Trung Software City, Tan Chanh
Hiep Ward, District 12, Ho Chi Minh City 700000, Vietnam
- Dong
Thap University, 783 Pham Huu Lau Street, Ward 6, Cao Lanh
City 870000, Dong Thap, Vietnam
| | - Duc Toan Truong
- Life Science Lab, Institute for Computational Science and Technology, Quang Trung Software City, Tan Chanh
Hiep Ward, District 12, Ho Chi Minh City 700000, Vietnam
| | - Mai Suan Li
- Institute of Physics, Polish Academy of Sciences, al. Lotnikow 32/46, Warsaw 02-668, Poland
| |
Collapse
|
18
|
Nguyen HL, Lan PD, Thai NQ, Nissley DA, O'Brien EP, Li MS. Does SARS-CoV-2 Bind to Human ACE2 More Strongly Than Does SARS-CoV? J Phys Chem B 2020; 124:7336-7347. [PMID: 32790406 PMCID: PMC7433338 DOI: 10.1021/acs.jpcb.0c04511] [Citation(s) in RCA: 90] [Impact Index Per Article: 22.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
![]()
The
2019 novel coronavirus (SARS-CoV-2) epidemic, which was first
reported in December 2019 in Wuhan, China, was declared a pandemic
by the World Health Organization in March 2020. Genetically, SARS-CoV-2
is closely related to SARS-CoV, which caused a global epidemic with
8096 confirmed cases in more than 25 countries from 2002 to 2003.
Given the significant morbidity and mortality rate, the current pandemic
poses a danger to all of humanity, prompting us to understand the
activity of SARS-CoV-2 at the atomic level. Experimental studies have
revealed that spike proteins of both SARS-CoV-2 and SARS-CoV bind
to angiotensin-converting enzyme 2 (ACE2) before entering the cell
for replication. However, the binding affinities reported by different
groups seem to contradict each other. Wrapp et al. (Science2020, 367, 1260–1263) showed
that the spike protein of SARS-CoV-2 binds to the ACE2 peptidase domain
(ACE2-PD) more strongly than does SARS-CoV, and this fact may be associated
with a greater severity of the new virus. However, Walls et al. (Cell2020, 181, 281–292)
reported that SARS-CoV-2 exhibits a higher binding affinity, but the
difference between the two variants is relatively small. To understand
the binding mechnism and experimental results, we investigated how
the receptor binding domain (RBD) of SARS-CoV (SARS-CoV-RBD) and SARS-CoV-2
(SARS-CoV-2-RBD) interacts with a human ACE2-PD using molecular modeling.
We applied a coarse-grained model to calculate the dissociation constant
and found that SARS-CoV-2 displays a 2-fold higher binding affinity.
Using steered all-atom molecular dynamics simulations, we demonstrate
that, like a coarse-grained simulation, SARS-CoV-2-RBD was associated
with ACE2-PD more strongly than was SARS-CoV-RBD, as evidenced by
a higher rupture force and larger pulling work. We show that the binding
affinity of both viruses to ACE2 is driven by electrostatic interactions.
Collapse
Affiliation(s)
- Hoang Linh Nguyen
- Life Science Lab, Institute for Computational Science and Technology, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam
| | - Pham Dang Lan
- Life Science Lab, Institute for Computational Science and Technology, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam.,Faculty of Physics and Engineering Physics, VNUHCM-University of Science, 227, Nguyen Van Cu Street, District 5, Ho Chi Minh City, Vietnam
| | - Nguyen Quoc Thai
- Life Science Lab, Institute for Computational Science and Technology, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam.,Dong Thap University, 783 Pham Huu Lau Street, Ward 6, Cao Lanh City, Dong Thap, Vietnam
| | - Daniel A Nissley
- Department of Statistics, University of Oxford, Oxford OX1 2JD, United Kingdom
| | - Edward P O'Brien
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States.,Bioinformatics and Genomics Graduate Program, The Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, Pennsylvania 16802, United States.,Institute for Computational and Data Sciences, Pennsylvania State University, University Park, Pennsylvania 16802,United States
| | - Mai Suan Li
- Institute of Physics, Polish Academy of Sciences, al. Lotnikow 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
19
|
Abstract
![]()
Accurate determination
of the binding affinity of the ligand to
the receptor remains a difficult problem in computer-aided drug design.
Here, we study and compare the efficiency of Jarzynski’s equality
(JE) combined with steered molecular dynamics and the linear interaction
energy (LIE) method by assessing the binding affinity of 23 small
compounds to six receptors, including β-lactamase, thrombin,
factor Xa, HIV-1 protease (HIV), myeloid cell leukemia-1, and cyclin-dependent
kinase 2 proteins. It was shown that Jarzynski’s nonequilibrium
binding free energy ΔGneqJar correlates with the available
experimental data with the correlation levels R =
0.89, 0.86, 0.83, 0.80, 0.83, and 0.81 for six data sets, while for
the binding free energy ΔGLIE obtained
by the LIE method, we have R = 0.73, 0.80, 0.42,
0.23, 0.85, and 0.01. Therefore, JE is recommended to be used for
ranking binding affinities as it provides accurate and robust results.
In contrast, LIE is not as reliable as JE, and it should be used with
caution, especially when it comes to new systems.
Collapse
Affiliation(s)
- Kiet Ho
- Institute for Computational Sciences and Technology, Quang Trung Software City, SBI Building, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam
| | - Duc Toan Truong
- Institute for Computational Sciences and Technology, Quang Trung Software City, SBI Building, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam.,Department of Theoretical Physics, Faculty of Physics and Engineering Physics, Ho Chi Minh University of Science, Ho Chi Minh City, Vietnam
| | - Mai Suan Li
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
20
|
Hacisuleyman A, Erman B. ModiBodies: A computational method for modifying nanobodies in nanobody-antigen complexes to improve binding affinity and specificity. J Biol Phys 2020; 46:189-208. [PMID: 32418062 DOI: 10.1007/s10867-020-09548-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2020] [Accepted: 04/28/2020] [Indexed: 11/26/2022] Open
Abstract
Nanobodies are special derivatives of antibodies, which consist of single domain fragments. They have become of considerable interest as next-generation biotechnological tools for antigen recognition. They can be easily engineered due to their high stability and compact size. Nanobodies have three complementarity-determining regions, CDRs, which are enlarged to provide a similar binding surface to that of human immunoglobulins. Here, we propose a benchmark testing algorithm that uses 3D structures of already existing protein-nanobody complexes as initial structures followed by successive mutations on the CDR domains. The aim is to find optimum binding amino acids for hypervariable residues of CDRs. We use molecular dynamics simulations to compare the binding energies of the resulting complexes with that of the known complex and accept those that are improved by mutations. We use the MDM4-VH9 complex, (PDB id 2VYR), fructose-bisphosphate aldolase from Trypanosoma congolense (PDB id 5O0W) and human lysozyme (PDB id 4I0C) as benchmark complexes. By using this algorithm, better binding nanobodies can be generated in a short amount of time. We suggest that this method can complement existing immune and synthetic library-based methods, without a need for extensive experimentation or large libraries.
Collapse
Affiliation(s)
- Aysima Hacisuleyman
- Department of Chemical and Biological Engineering, Koc University, Istanbul, Turkey.
| | - Burak Erman
- Department of Chemical and Biological Engineering, Koc University, Istanbul, Turkey
| |
Collapse
|
21
|
Pham T, Nguyen HL, Phan-Toai T, Nguyen H. Investigation of Binding Affinity between Potential Antiviral Agents and PB2 Protein of Influenza A: Non-equilibrium Molecular Dynamics Simulation Approach. Int J Med Sci 2020; 17:2031-2039. [PMID: 32788882 PMCID: PMC7415388 DOI: 10.7150/ijms.46231] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/22/2020] [Accepted: 07/09/2020] [Indexed: 11/30/2022] Open
Abstract
The PB2 protein of the influenza virus RNA polymerase is a major virulence determinant of influenza viruses. It binds to the cap structure at the 5' end of host mRNA to generate short capped RNA fragments that are used as primers for viral transcription named cap-snatching. A large number of the compounds were shown to bind the minimal cap-binding domain of PB2 to inhibit the cap-snatching machinery. However, their binding in the context of an extended form of the PB2 protein has remained elusive. A previous study reported some promising compounds including azaindole and hydroxymethyl azaindole, which were analyzed here to predict binding affinity to PB2 protein using the steered molecular dynamics (SMD) and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) methods. The results show that the rupture force (Fmax) value of three complexes is in agreement with the binding free energy value (ΔGbind) estimated by the MM-PBSA method, whereas for the non-equilibrium pulling work (Wpull) value a small difference between A_PB2-4 and A_PB2-12 was observed. The binding affinity results indicate the A_PB2-12 complex is more favorable than the A_PB2-4 and A_PB2-16 complexes, which means the inhibitor (12) has the potential to be further developed as anti-influenza agents in the treatment of influenza A.
Collapse
Affiliation(s)
- Tri Pham
- Institute for Computational Science and Technology, Ho Chi Minh City, Vietnam.,VNUHCM-University of Technology, Ho Chi Minh City, Vietnam
| | - Hoang Linh Nguyen
- Institute for Computational Science and Technology, Ho Chi Minh City, Vietnam.,VNUHCM-University of Technology, Ho Chi Minh City, Vietnam
| | - Tuyn Phan-Toai
- Institute for Computational Science and Technology, Ho Chi Minh City, Vietnam
| | - Hung Nguyen
- Institute for Computational Science and Technology, Ho Chi Minh City, Vietnam
| |
Collapse
|
22
|
Detailed potential of mean force studies on host-guest systems from the SAMPL6 challenge. J Comput Aided Mol Des 2018; 32:1013-1026. [PMID: 30143917 DOI: 10.1007/s10822-018-0153-7] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2018] [Accepted: 08/11/2018] [Indexed: 12/14/2022]
Abstract
Accurately predicting receptor-ligand binding free energies is one of the holy grails of computational chemistry with many applications in chemistry and biology. Many successes have been reported, but issues relating to sampling and force field accuracy remain significant issues affecting our ability to reliably calculate binding free energies. In order to explore these issues in more detail we have examined a series of small host-guest complexes from the SAMPL6 blind challenge, namely octa-acids (OAs)-guest complexes and Curcurbit[8]uril (CB8)-guest complexes. Specifically, potential of mean force studies using umbrella sampling combined with the weighted histogram method were carried out on both systems with both known and unknown binding affinities. We find that using standard force fields and straightforward simulation protocols we are able to obtain satisfactory results, but that simply scaling our results allows us to significantly improve our predictive ability for the unknown test sets: the overall RMSD of the binding free energy versus experiment is reduced from 5.59 to 2.36 kcal/mol; for the CB8 test system, the RMSD goes from 8.04 to 3.51 kcal/mol, while for the OAs test system, the RSMD goes from 2.89 to 0.95 kcal/mol. The scaling approach was inspired by studies on structurally related known benchmark sets: by simply scaling, the RMSD was reduced from 6.23 to 1.19 kcal/mol and from 2.96 to 0.62 kcal/mol for the CB8 benchmark system and the OA benchmark system, respectively. We find this scaling procedure to correct absolute binding affinities to be highly effective especially when working across a "congeneric" series with similar charge states. It is less successful when applied to mixed ligands with varied charges and chemical characteristics, but improvement is still realized in the present case. This approach suggests that there are large systematic errors in absolute binding free energy calculations that can be straightforwardly accounted for using a scaling procedure. Random errors are still an issue, but near chemical accuracy can be obtained using the present strategy in select cases.
Collapse
|