The molecular interaction of six single-stranded DNA aptamers to cardiac troponin I revealed by docking and molecular dynamics simulation.

Cardiac troponin I (cTnI) is a cardiac biomarker for diagnosing ischemic heart disease and acute myocardial infarction. Current biochemical assays use antibodies (Abs) due to their high specificity and sensitivity. However, there are some limitations, such as the high-cost production of Abs due to complex instruments, reagents, and steps; the variability of Abs quality from batch to batch; the low stability at high temperatures; and the difficulty of chemical modification. Aptamer overcomes the limitations of antibodies, such as relatively lower cost, high reproducibility, high stability, and ease of being chemically modified. Aptamers are three-dimensional architectures of single-stranded RNA or DNA that bind to targets such as proteins. Six aptamers (Tro1-Tro6) with higher binding affinity than an antibody have been identified, but the molecular interaction has not been studied. In this study, six DNA aptamers were modeled and docked to cTnI protein. Molecular docking revealed that the interaction between all aptamer and cTnI happened in the similar cTnI region. The interaction between aptamer and cTnI involved hydrophobic interaction, hydrogen bonds, π-cation interactions, π-stack interactions, and salt-bridge formation. The calculated binding energy of all complexes was negative, which means that the complex formation was thermodynamically favorable. The electrostatic energy term was the main driving force of the interaction between all aptamer and cTnI. This study could be used to predict the behavior of further modified aptamer to improve aptamer performance.

Designing a Novel Multiepitope Vaccine from the Human Papilloma Virus E1 and E2 Proteins for Indonesia with Immunoinformatics and Molecular Dynamics Approaches.

One of the deadliest malignant cancer in women globally is cervical cancer. Specifically, cervical cancer is the second most common type of cancer in Indonesia. The main infectious agent of cervical cancer is the human papilloma virus (HPV). Although licensed prophylactic vaccines are available, cervical cancer cases are on the rise. Therapy using multiepitope-based vaccines is a very promising therapy for cervical cancer. This study aimed to develop a multiepitope vaccine based on the E1 and E2 proteins of HPV 16, 18, 45, and 52 using in silico. In this study, we develop a novel multiepitope vaccine candidate using an immunoinformatic approach. We predicted the epitopes of the cytotoxic T lymphocyte (CTL) and helper T lymphocyte (HTL) and evaluated their immunogenic properties. Population coverage analysis of qualified epitopes was conducted to determine the successful use of the vaccine worldwide. The epitopes were constructed into a multiepitope vaccine by using AAY linkers between the CTL epitopes and GPGPG linkers between the HTL epitopes. The tertiary structure of the multiepitope vaccine was modeled with AlphaFold and was evaluated by Prosa-web. The results of vaccine construction were analyzed for B-cell epitope prediction, molecular docking with Toll like receptor-4 (TLR4), and molecular dynamics simulation. The results of epitope prediction obtained 4 CTL epitopes and 7 HTL epitopes that are eligible for construction of multiepitope vaccines. Prediction of the physicochemical properties of multiepitope vaccines obtained good results for recombinant protein production. The interaction showed that the interaction of the multiepitope vaccine-TLR4 complex is stable based on the binding free energy value -106.5 kcal/mol. The results of the immune response simulation show that multiepitope vaccine candidates could activate the adaptive and humoral immune systems and generate long-term B-cell memory. According to these results, the development of a multiepitope vaccine with a reverse vaccinology approach is a breakthrough to develop potential cervical cancer therapeutic vaccines.

The SARS-CoV-2 Mpro Dimer-Based Screening System: A Synthetic Biology Tool for Identifying Compounds with Dimerization Inhibitory Potential.

The COVID-19 endemic remains a global concern. The search for effective antiviral candidates is still needed to reduce disease risk. However, the availability of high biosafety level laboratory facilities for drug screening is limited in number. To address this issue, a screening system that could be utilized at lower biosafety levels remains essential. This study aimed to develop a novel SARS-CoV-2 main protease (Mpro) dimer-based screening system (DBSS) utilizing synthetic biology in Escherichia coli BL21(DE3). We linked the SARS-CoV-2 Mpro with the DNA-binding domain of AraC regulatory protein, which regulates the reporter gene expression. Protein modeling and molecular docking showed that saquinavir could bind to AraC-Mpro both in its monomer and dimer forms. The constructed DBSS assay indicated the screening system could detect saquinavir inhibitory activity at a concentration range of 4-10 µg/mL compared to the untreated control (P ≤ 0.05). The Vero E6 cell assay validated the DBSS result that saquinavir at 4-10 µg/mL exhibited antiviral activity against SARS-CoV-2. Our DBSS could be used for preliminary screening of numerous drug candidates that possess a dimerization inhibitor activity of SARS-CoV-2 Mpro and also minimize the use of a high biosafety level laboratory.

The assessment of molecular dynamics results of three-dimensional RNA aptamer structure prediction.

Aptamers are single-stranded DNA or RNA that bind to specific targets such as proteins, thus having similar characteristics to antibodies. It can be synthesized at a lower cost, with no batch-to-batch variations, and is easier to modify chemically than antibodies, thus potentially being used as therapeutic and biosensing agents. The current method for RNA aptamer identification in vitro uses the SELEX method, which is considered inefficient due to its complex process. Computational models of aptamers have been used to predict and study the molecular interaction of modified aptamers to improve affinity. In this study, we generated three-dimensional models of five RNA aptamers from their sequence using mFold, RNAComposer web server, and molecular dynamics simulation. The model structures were then evaluated and compared with the experimentally determined structures. This study showed that the combination of mFold, RNAComposer, and molecular dynamics simulation could generate 14-16, 28, or 29 nucleotides length of 3D RNA aptamer with similar geometry and topology to the experimentally determined structures. The non-canonical basepair structure of the aptamer loop was formed through the MD simulation, which also improved the three-dimensional RNA aptamers model. Clustering analysis was recommended to choose the more representative model.

Evaluating the potency of Sulawesi propolis compounds as ACE-2 inhibitors through molecular docking for COVID-19 drug discovery preliminary study.

Coronavirus disease (COVID-19) is a global pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Up to date, there has been no specific cure to treat the disease. Indonesia is one of the countries that is still fighting to control virus transmission. Yet, at the same time, Indonesia has a rich biodiversity of natural medicinal products that potentially become an alternative cure. Thus, this study examined the potency of a natural medicinal product, Sulawesi propolis compounds produced by Tetragonula sapiens, inhibiting angiotensin-converting activity enzyme-2 (ACE-2), a receptor of SARS-CoV-2 in the human body. In this study, molecular docking was done to analyze the docking scores as the representation of binding affinity and the interaction profiles of propolis compounds toward ACE-2. The results illustrated that by considering the docking score and the presence of interaction with targeted sites, five compounds, namely glyasperin A, broussoflavonol F, sulabiroins A, (2S)-5,7-dihydroxy-4'-methoxy-8-prenylflavanone and isorhamnetin are potential to inhibit the binding of ACE-2 and SARS-CoV-2, with the docking score of -10.8, -9.9, -9.5, -9.3 and -9.2 kcal/mol respectively. The docking scores are considered to be more favorable compared to MLN-4760 as a potent inhibitor.

Molecular interaction analysis of Sulawesi propolis compounds with SARS-CoV-2 main protease as preliminary study for COVID-19 drug discovery.

Coronavirus disease 2019 (COVID-19), a respiratory disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a global health concern, as the World Health Organization declared this outbreak to be a global pandemic in March 2020. The need for an effective treatment is urgent because the development of an effective vaccine may take years given the complexity of the virus and its rapid mutation. One promising treatment target for COVID-19 is SARS-CoV-2 main protease. Thus, this study was aimed to examine whether Sulawesi propolis compounds produced by Tetragonula sapiens inhibit the enzymatic activity of SARS-CoV-2 main protease. In this study, molecular docking was performed to analyze the interaction profiles of propolis compounds with SARS-CoV-2 main protease. The results illustrated that two compounds, namely glyasperin A and broussoflavonol F, are potential drug candidates for COVID-19 based on their binding affinity of -7.8 kcal/mol and their ability to interact with His41 and Cys145 as catalytic sites. Both compounds also displayed favorable interaction profiles with SARS-CoV-2 main protease with binding similarities compared to inhibitor 13b as positive control 63% and 75% respectively.

Designing a less immunogenic nattokinase from Bacillus subtilis subsp. natto: a computational mutagenesis.

Nattokinase is an enzyme produced by Bacillus subtilis subsp. natto that contains strong fibrinolytic activity. It has potential to treat cardiovascular diseases. In silico analysis revealed that nattokinase is considered as an antigen, thus hindering its application for injectable therapeutic protein. Various web servers were used to predict B-cell epitopes of nattokinase both continuously and discontinuously to determine which amino acid residues had been responsible for the immunogenicity. With the exclusion of the predicted conserved amino acids, four amino acids such as S18, Q19, T242, and Q245 were allowed for mutation. Substitution mutation was done to lower the immunogenicity of native nattokinase. Through the stability of the mutated protein with the help of Gibbs free energy difference, the proposed mutein was S18D, Q19I, T242Y, and Q245W. The 3D model of the mutated nattokinase was modeled and validated with various tools. Physicochemical properties and stability analysis of the protein indicated that the mutation brought higher stability without causing any changes in the catalytic site of nattokinase. Molecular dynamics simulation implied that the mutation indicated similar stability, conformation, and behavior compared to the native nattokinase. These results are highly likely to contribute to the wet lab experiment to develop safer nattokinase.