Molecular Mechanism and Energy Basis of Conformational Diversity of Antibody SPE7 Revealed by Molecular Dynamics Simulation and Principal Component Analysis

Unveiling the ESR1 Conformational Stability and Screening Potent Inhibitors for Breast Cancer Treatment

BACKGROUND

The current study recognizes the significance of estrogen receptor alpha (ERα) as a member of the nuclear receptor protein family, which holds a central role in the pathophysiology of breast cancer. ERα serves as a valuable prognostic marker, with its established relevance in predicting disease outcomes and treatment responses.

METHODS

In this study, computational methods are utilized to search for suitable drug-like compounds that demonstrate analogous ligand binding kinetics to ERα.

RESULTS

Docking-based simulation screened out the top 5 compounds - ZINC13377936, NCI35753, ZINC35465238, ZINC14726791, and NCI663569 against the targeted protein. Further, their dynamics studies reveal that the compounds ZINC13377936 and NCI35753 exhibit the highest binding stability and affinity.

CONCLUSION

Anticipating the competitive inhibition of ERα protein expression in breast cancer, we envision that both ZINC13377936 and NCI35753 compounds hold substantial promise as potential therapeutic agents. These candidates warrant thorough consideration for rigorous In vitro and In vivo evaluations within the context of clinical trials. The findings from this current investigation carry significant implications for the advancement of future diagnostic and therapeutic approaches for breast cancer.

Computational assessment of herbal medicine-derived compounds as potential inhibitors of SARS-CoV-2 main protease

Since the main protease (Mpro) is crucial for the COVID-19 virus replication and transcription, searching for Mpro inhibitors is one possible treatment option. In our study, 258 small molecules were collected from lung-related herbal medicines, and their structures were optimized with the B3LYP-D3/6-31G* method. After the molecular docking with Mpro, we selected the top 20 compounds for the further geometry optimization with the larger basis sets. After the further molecular docking, the top eight compounds were screened out. Then we performed molecular dynamics simulations and binding free energy calculations to determine stability of the complexes. Our results show that mulberrofuran G, Xambioona, and kuwanon D can bind Mpro well. In quantum chemistry studies, such as ESP and CDFT analyses, the compounds properties are predicted. Additionally, the drug-likeness analyses and ADME studies on these three candidate compounds verified that all of them conform to Libinski's rule and may be drug-like compounds.

Computational biophysics approach towards the discovery of multi-kinase blockers for the management of MAPK pathway dysregulation

The MAPK pathway is important in human lung cancer and is improperly activated in a substantial proportion through number of ways. Strategies on dual-targeting RAF and MEK are an alternative option to diminish the limitations in this pathway inhibition. Hence, we implemented parallel pharmacophore screening of 11,808 DrugBank compounds against RAF and MEK. ADHRR and DHHRR were modeled as a pharmacophore hypothesis for RAF and MEK respectively. Importantly, these hypotheses resulted an AUC value of > 0.90 with the external data set. As a result of phase screening, glide docking, and prime-MM/GBSA scoring, it is determined that DB08424 and DB08907 have the best chances of acting as multi-kinase inhibitors. The pi-cation interaction with key amino acid residues of both target receptors may responsible for the stronger binding with these kinases. Cumulative 600 ns MD simulation studies validate the binding ability of these compounds. Significantly, the hit compounds resulted higher number of stable conformational state with less atomic movements than the reference compound against both targets. The anti-cancer efficacy of the lead compounds was validated through machine learning-based approaches. These findings suggest that DB08424 and DB08907 might be novel molecules to be explored further experimentally to block the MAPK signaling in lung cancer patients.

Affinity enhancement of CR3022 binding to RBD; in silico site directed mutagenesis using molecular dynamics simulation approaches

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a disease which caused by a novel beta coronavirus. Structural and non-structural proteins are expressed by the virus gene fragments. The RBD of the S1 protein of the virus has the ability to interact with potent antibodies including CR3022, which was characterized to target the S protein of the virus which can efficiently neutralize the SARS-CoV in vitro and in vivo. In current study, we aimed to design CR3022 based antibody with high affinity compared with wild-type CR3022 using MD simulation method. Two variants were designed based on the amino acid binding conformation and the free binding energy of the critical amino acids which involved in CR3022-RBD interactions were evaluated. In this study three complexes were evaluated; CR3022-RBD, V1-RBD and V2-RBD using molecular dynamics simulations carried out for 100 ns in each case. Then, all the complexes were simulated for 100 ns. In the next step, to calculate the free binding affinity of the wild CR3022 and mutant antibody (V1 and V2) with RBD, the PMF method was performed. The RMSD profile demonstrated that all three complexes were equilibrated after 85 ns. Furthermore, the free binding energy results indicated that the V2-RBD complex has the higher binding affinity than V1-RBD and CR3022-RBD complexes. It should be noted that in above variants, the electrostatic energy and the number of H-bonds between the antibody and RBD increased. Thus, it is suggested that both designed antibodies could be considered as appropriate candidates for covid-19 disease treatment.Communicated by Ramaswamy H. Sarma.

Effect of phosphorylation of protamine-like cationic peptide on the binding affinity to DNA

Protamines are more arginine-rich and more basic than histones and are responsible for providing a highly compacted shape to the sperm heads in the testis. Phosphorylation and dephosphorylation are two events that occur in the late phase of spermatogenesis before the maturation of sperms. In this work, we have studied the effect of phosphorylation of protamine-like cationic peptides using all-atom molecular dynamics simulations. Through thermodynamic analyses, we found that phosphorylation reduces the binding efficiency of such cationic peptides on DNA duplexes. Peptide phosphorylation leads to a less efficient DNA condensation, due to a competition between DNA-peptide and peptide-peptide interactions. We hypothesize that the decrease of peptide bonds between DNA together with peptide self-assembly might allow an optimal re-organization of chromatin and an efficient condensation through subsequent peptide dephosphorylation. Based on the globular and compact conformations of phosphorylated peptides mediated by arginine-phosphoserine H-bonding, we furthermore postulate that phosphorylated protamines could more easily intrude into chromatin and participate to histone release through disruption of histone-histone and histone-DNA binding during spermatogenesis.

Molecular dynamics study on the strengthening behavior of Delta and Omicron SARS-CoV-2 spike RBD improved receptor-binding affinity

The COVID-19 pandemic caused by a virus that can be transmitted from human to human via air droplets has changed the quality of life and economic systems all over the world. The viral DNA has mutated naturally over time leading to the diversity of coronavirus victims which has posed a serious threat to human security on a massive scale. The current variants have developed in a dominant way and are considered "Variants of Concern" by the World Health Organization (WHO). In this work, Kappa (B.1.617.1), Delta (B.1.617.2), and Omicron (B.1.1.529) variants were obtained to evaluate whether naturally occurring mutations have strengthened viral infectivity. We apply reliable in silico structural dynamics and energetic frameworks of the mutated S-RBD protein for ACE2-binding to analyze and compare the structural information related to the wild-type. In particular, the hotspot residues at Q493, Q498, and N501 on the S-RBD protein were determined as contributing factors to the employment stability of the relevant binding interface. The L452R mutation induces an increment of the hydrogen bonds formed by changing the Q493 environment for ACE2 binding. Moreover, the Q493K exchange in Omicron enables the formation of two additional salt bridges, leading to a strong binding affinity by increased electrostatic interaction energy. These results could be used in proposing concrete informative data for a structure-based design engaged in finding better therapeutics against novel variants.

Remdesivir MD Simulations Suggest a More Favourable Binding to SARS-CoV-2 RNA Dependent RNA Polymerase Mutant P323L Than Wild-Type

SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) protein is the target for the antiviral drug Remdesivir (RDV). With RDV clinical trials on COVID-19 patients showing a reduced hospitalisation time. During the spread of the virus, the RdRp has developed several mutations, with the most frequent being A97V and P323L. The current study sought to investigate whether A97V and P323L mutations influence the binding of RDV to the RdRp of SARS-CoV-2 compared to wild-type (WT). The interaction of RDV with WT-, A97V-, and P323L-RdRp were measured using molecular dynamic (MD) simulations, and the free binding energies were extracted. Results showed that RDV that bound to WT- and A97V-RdRp had a similar dynamic motion and internal residue fluctuations, whereas RDV interaction with P323L-RdRp exhibited a tighter molecular conformation, with a high internal motion near the active site. This was further corroborated with RDV showing a higher binding affinity to P323L-RdRp (-24.1 kcal/mol) in comparison to WT-RdRp (-17.3 kcal/mol). This study provides insight into the potential significance of administering RDV to patients carrying the SARS-CoV-2 P323L-RdRp mutation, which may have a more favourable chance of alleviating the SARS-CoV-2 illness in comparison to WT-RdRp carriers, thereby suggesting further scientific consensus for the usage of Remdesivir as clinical candidate against COVID-19.

Covalent Stabilization of Antibody Recruitment Enhances Immune Recognition of Cancer Targets

Antibody recruiting molecules (ARMs) represent an important class of "proximity-inducing" chemical tools with therapeutic potential. ARMs function by simultaneously binding to a hapten-specific serum antibody (Ab) (e.g., anti-dinitrophenyl (DNP)) and a cancer cell surface protein, enforcing their proximity. ARM anticancer efficacy depends on the formation of ARM:Ab complexes on the cancer cell surface, which activate immune cell recognition and elimination of the cancer cell. Problematically, ARM function in human patients may be limited by conditions that drive the dissociation of ARM:Ab complexes, namely, intrinsically low binding affinity and/or low concentrations of anti-hapten antibodies in human serum. To address this potential limitation, we previously developed a covalent ARM (cARM) chemical tool that eliminates the ARM:antibody equilibrium through a covalent linkage. In the current study, we set out to determine to what extent maximizing the stability of ARM:antibody complexes via cARMs enhances target immune recognition. We observe cARMs significantly increase target immune recognition relative to ARMs across a range of therapeutically relevant antibody concentrations. These results demonstrate that ARM therapeutic function can be dramatically enhanced by increasing the kinetic stability of ARM:antibody complexes localized on cancer cells. Our findings suggest that a) high titres/concentrations of target antibody in human serum are not neccessary and b) saturative antibody recruitment to cancer cells not sufficient, to achieve maximal ARM therapeutic function.

Enfuvirtide, an HIV-1 fusion inhibitor peptide, can act as a potent SARS-CoV-2 fusion inhibitor: an in silico drug repurposing study

Regarding the urgency of therapeutic measures for coronavirus disease 2019 (COVID-19) pandemic, the use of available drugs with FDA approval is preferred because of the less time and cost required for their development. In silico drug repurposing is an accurate way to speed up the screening of the existing FDA-approved drugs to find a therapeutic option for COVID-19. The similarity in SARS-CoV-2 and HIV-1 fusion mechanism to host cells can be a key point for Inhibit SARS-CoV-2 entry into host cells by HIV fusion inhibitors. Accordingly, in this study, an HIV-1 fusion inhibitor called Enfuvirtide (Enf) was selected. The affinity and essential residues involving in the Enf binding to the S2 protein of SARS-CoV-2, HIV-1 gp41 protein and angiotensin-converting enzyme 2 (ACE-2) as a negative control, was evaluated using molecular docking. Eventually, Enf-S2 and Enf-gp41 protein complexes were simulated by molecular dynamics (MD) in terms of binding affinity and stability. Based on the most important criteria such as docking score, cluster size, energy and dissociation constant, the strongest interaction was observed between Enf with the S2 protein. In addition, MD results confirmed that Enf-S2 protein interaction was remarkably stable and caused the S2 protein residues to undergo the fewest fluctuations. In conclusion, it can be stated that Enf can act as a strong SARS-CoV-2 fusion inhibitor and demonstrates the potential to enter the clinical trial phase of COVID-19.

Communicated by Ramaswamy H. Sarma

Interplay among Structural Stability, Plasticity, and Energetics Determined by Conformational Attuning of Flexible Loops in PD-1

The dynamics and plasticity of the PD-1/PD-L1 axis are the bottlenecks for the discovery of small-molecule antagonists to perturb this interaction interface significantly. Understanding the process of this protein-protein interaction (PPI) is of fundamental biological interest in structure-based drug designing. Food and Drug Administration (FDA)-approved anti-PD-1 monoclonal antibodies (mAbs) are the first-in-class with distinct binding modes to access this axis clinically; however, their mechanistic aspects remain elusive. Here, we have unveiled the interactive interfaces with PD-L1 and mAbs to investigate the native plasticity of PD-1 at global (structural and dynamical) and local (residue side-chain orientations) levels. We found that the structural stability and coordinated C_α movements are increased in the presence of PD-1's binding partners. The rigorous analysis of these PPIs using computational biophysical approaches revealed PD-1's intrinsic plasticity, its concerted loops' movement (BC, FG, and CC'), distal side-chain motions, and the thermodynamic landscape, which are perturbed remarkably from its unbound to bound states. Based on intra-/inter-residues' contact networks and energetics, the hot-spots have been identified that were found to be essential to arrest the dynamical motions of PD-1 significantly for the rational design of therapeutic agents by mimicking the mAbs mechanism.

Revealing the Inhibition Mechanism of RNA-Dependent RNA Polymerase (RdRp) of SARS-CoV-2 by Remdesivir and Nucleotide Analogues: A Molecular Dynamics Simulation Study

Antiviral drug therapy against SARS-CoV-2 is not yet established and posing a serious global health issue. Remdesivir is the first antiviral compound approved by the US FDA for the SARS-CoV-2 treatment for emergency use, targeting RNA-dependent RNA polymerase (RdRp) enzyme. In this work, we have examined the action of remdesivir and other two ligands screened from the library of nucleotide analogues using docking and molecular dynamics (MD) simulation studies. The MD simulations have been performed for all the ligand-bound RdRp complexes for the 30 ns time scale. This is one of the earlier reports to perform the MD simulations studies using the SARS-CoV-2 RdRp crystal structure (PDB ID 7BTF). The MD trajectories were analyzed and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) calculations were performed to calculate the binding free energy. The binding energy data reveal that compound-17 (-59.6 kcal/mol) binds more strongly as compared to compound-8 (-46.3 kcal/mol) and remdesivir (-29.7 kcal/mol) with RdRp. The detailed analysis of trajectories shows that the remdesivir binds in the catalytic site and forms a hydrogen bond with the catalytic residues from 0 to 0.46 ns. Compound-8 binds in the catalytic site but does not form direct hydrogen bonds with catalytic residues. Compound-17 showed the formation of hydrogen bonds with catalytic residues throughout the simulation process. The MD simulation results such as hydrogen bonding, the center of mass distance analysis, snapshots at a different time interval, and binding energy suggest that compound-17 binds strongly with RdRp of SARS-CoV-2 and has the potential to develop as a new antiviral against COVID-19. Further, the frontier molecular orbital analysis and molecular electrostatic potential (MESP) iso-surface analysis using DFT calculations shed light on the superior binding of compound-17 with RdRp compared to remdesivir and compound-8. The computed as well as the experimentally reported pharmacokinetics and toxicity parameters of compound-17 is encouraging and therefore can be one of the potential candidates for the treatment of COVID-19.

Immunoinformatics and molecular dynamics studies to predict T-cell-specific epitopes of four Klebsiella pneumoniae fimbriae antigens

Klebsiella pneumoniae (K. pneumoniae) is a causative agent of severe infections in humans. There is no publically available vaccine for K. pneumoniae infections yet. Here, using comprehensive immunoinformatics methods, T-cell-specific epitopes of four type 1 fimbriae antigens of K. pneumoniae were predicted and evaluated as potential vaccine candidates. Both CD8+ (class I) and CD4+ (class II) T-cell-specific epitopes were predicted and the epitopes similar to human proteome were excluded. Subsequently, the windows of class-II epitopes containing class-I epitopes were determined. The immunogenicity, IFN-γ production and population coverage were also estimated. Using the 3D structure of HLA and epitopes, molecular docking was carried out. Two best epitopes were selected for molecular dynamics studies. Our prediction and analyses resulted in the several dominant epitopes for each antigen. The docking results showed that all selected epitopes can bind to their restricted HLA molecules with high affinity. The molecular dynamics results indicated the stability of system with minimum possible deviation, suggesting the selected epitopes can be promising candidates for stably binding to HLA molecules. Altogether, our results suggest that the selected T-cell-specific epitopes of K. pneumoniae fimbriae antigens, particularly the two epitopes confirmed by molecular dynamics, can be applied for vaccine development. However, the in vitro and in vivo studies are required to authenticate the results of the present study.Communicated by Ramaswamy H. Sarma.

Effects of unsaturated fatty acids (Arachidonic/Oleic Acids) on stability and structural properties of Calprotectin using molecular docking and molecular dynamics simulation approach

Calprotectin is a heterodimeric protein complex with two subunits called S100A8/A9. The protein has an essential role in inflammation process and various human diseases. It has the ability to bind to unsaturated fatty acids including Arachidonic acid, Oleic acid and etc., which could be considered as a major carrier for fatty acids. In this study we aimed to appraise the thermodynamics and structural changes of Calprotectin in presence of Arachidonic acid/Oleic acid) using docking and molecular dynami simulation method. To create the best conformation of Calprotectin-Oleic acid/Arachidonic acid complexes, the docking process was performed. The complexes with the best binding energy were selected as the models for molecular dynamics simulation process. Furthermore, the structural and thermodynamics properties of the complexes were evaluated too. The Root Mean Square Deviation and Root Mean Square Fluctuation results showed that the binding of Arachidonic acid/Oleic acid to Calprotectin can cause the protein structural changes which was confirmed by Define Secondary Structure of Proteins results. Accordingly, the binding free energy results verified that binding of Oleic acid to Calprotectin leads to instability of S100A8/A9 subunits in the protein. Moreover, the electrostatic energy contribution of the complexes (Calprotectin-Oleic acid/Arachidonic acid) was remarkably higher than van der Waals energy. Thus, the outcome of this study confirm that Oleic acid has a stronger interaction with Calprotectin in comparison with Arachidonic acid. Our findings indicated that binding of unsaturated fatty acids to Calprotectin leads to structural changes of the S100A8/A9 subunits which could be beneficial to play a biological role in inflammation process.

Enhanced metastable state models of TAM kinase binding to cabozantinib explains the dynamic nature of receptor tyrosine kinases

Receptor tyrosine kinases (RTKs) are essential proteins in the regulation of cell signaling. Tyro3, Axl and Mer are members of TAM RTKs and are overexpressed in several cancer forms. Kinase inhibitors such as cabozantinib, foretinib are reported to inhibit TAM kinases at nanomolar concentrations. The atomistic details of structure and mechanism of functional regulation is required to understand their normal physiological process and when bound to an inhibitor. The docking of cabozantinib into the active state conformations of TAM kinases (crystal structure and computational models) revealed the best binding pose and the complex formation that is mediated through non-bonding interactions involving the hinge region residues. The alterations in the conformations and the regions of flexibility in apo and complexed TAM kinases as a course of time are studied using 250 ns molecular dynamics (MD) simulations. The post-MD trajectory analysis using Python libraries like ProDy, MDTraj and PyEMMA revealed encrypted protein dynamic motions in active kinetic metastable states. Comparison between Principal component analysis and Anisotropic mode analysis deciphered structural residue interactions and salt bridge contacts between apo and inhibitor bound TAM kinases. Various structural changes occurred in αC-helix and activation loop involving hydrogen bonding between residues from Lys-(β3 sheet), Glu-(αC-helix) and Asp-(DFG-motif) resulting in higher RMSD. Mechanical stiffness plots revealed that similar regions in apo and cabozantinib bound Axl fluctuated during MD simulations whereas different regions in Tyro3 and Mer kinases. The residue interaction network plots revealed important salt bridges that lead to constrained domain motions in the TAM kinases.Communicated by Ramaswamy H. Sarma.

In silico assessment of human Calprotectin subunits (S100A8/A9) in presence of sodium and calcium ions using Molecular Dynamics simulation approach

Calprotectin is a heterodimeric protein complex which consists of two subunits including S100A8 and S100A9. This protein has a major role in different inflammatory disease and various types of cancers. In current study we aimed to evaluate the structural and thermodynamic changes of the subunits and the complex in presence of sodium and calcium ions using molecular dynamics (MD) simulation. Therefore, the residue interaction network (RIN) was visualized in Cytoscape program. In next step, to measure the binding free energy, the potential of mean force (PMF) method was performed. Finally, the molecular mechanics Poisson-Boltzmann surface area (MMPBSA) method was applied as an effective tool to calculate the molecular model affinities. The MD simulation results of the subunits represented their structural changes in presence of Ca2+. Moreover, the RIN and Hydrogen bond analysis demonstrated that cluster interactions between Calprotectin subunits in presence of Ca2+ were greater in comparison with Na+. Our findings indicated that the binding free energy of the subunits in presence of Ca2+ was significantly greater than Na+. The results revealed that Ca2+ has the ability to induce structural changes in subunits in comparison with Na+ which lead to create stronger interactions between. Hence, studying the physical characteristics of the human proteins could be considered as a powerful tool in theranostics and drug design purposes.

Zinc ion-induced conformational changes in new Delphi metallo-β-lactamase 1 probed by molecular dynamics simulations and umbrella sampling

The hydrolysis of a β-lactam core ring caused by new Delphi metallo-β-lactamase 1 (NDM-1) with the help of two zinc cofactors induces significant resistance toward β-lactam antibiotics. Molecular dynamics (MD) simulations and the umbrella sampling method are integrated to study the conformational change mechanism of NDM-1 mediated by zinc ion binding. The statistical analyses of interaction contacts of the antibiotic ampicillin (AMP) with residues based on MD trajectories suggest that two Zn ions are essential for maintaining the binding of AMP with NDM-1. Umbrella sampling simulations further reveal that double-Zn coordination exerts strong restriction on the motions of loop L10 relative to loops L3 and L4. Principal component (PC) analysis also demonstrates that zinc ion binding totally inhibits the motion extent of NDM-1 and changes internal motion modes in NDM-1. We expect that the current study can provide significant dynamical information involving conformational changes of NDM-1 for the development of efficient inhibitors to decrease drug resistance of NDM-1 toward antibiotics.

Mutation L1196M-induced conformational changes and the drug resistant mechanism of anaplastic lymphoma kinase studied by free energy perturbation and umbrella sampling

Anaplastic lymphoma kinase (ALK) has been regarded as a promising drug target in the treatment of tumors and the mutation L1196M induces different levels of drug resistance toward the existing inhibitors.