In the present study, we have done a comparative study on the efficacy of some currently used repurposed drugs: Oseltamivir (O), Favipiravir (F) and Hydroxychloroquine (H) in individual and in their combinational mode against CoV-2 infections. The ADME analysis has helped us to identify the inhibitory possibility of the tested drugs towards receptor 3CLpro protein of SARS-CoV-2. Various thermodynamical parameters obtained from Molecular Docking, Molecular dynamics (MD) and MMPBSA simulations like binding affinity, potential energy (Epot), RMSD, RMSF, SASA energy, interaction energies, Gibbs free energy (ΔGbind) etc. also helped us to verify the effectivity of mentioned drugs against CoV-2 protease.

The present study aims at numerically describing to what extent substrate - enzyme complexes in solution may change over time as a natural process of conformational changes for a liganded enzyme in comparison to those movements which occur independently from substrate interaction, i.e. without a ligand. To this end, we selected structurally known pairs of liganded / unliganded CYP450 3A4 enzymes with different geometries hinting at induced fit events. We carried out molecular dynamics simulations (MD) comparing the trajectories in a \"cross-over\" protocol: (i) we added the ligand to the unliganded crystal form which should adopt geometries similar to the known geometry of the liganded crystal structure during MD, and - conversely - (ii) we removed the bound ligand form the known liganded complex to test if a geometry similar to the known unliganded (apo-) form can be adopted during MD. To compare continues changes we measured root means square deviations and frequencies. Results for case (i) hint at larger conformational changes required for accepting the substrate during its approach to final position - in contrast to case (ii) when mobility is fairly reduced by ligand binding (strain energy). In conclusion, a larger conformational sampling prior to ligand binding and the freezing-in (rigidity) of conformations for bound ligands can be interpreted as two conditions linked to induced-fit.

COVID-19 is a viral pandemic caused by SARS-CoV-2. Due to its highly contagious nature, millions of people are getting affected worldwide knocking down the delicate global socio-economic equilibrium. According to the World Health Organization, COVID-19 has affected over 186 million people with a mortality of around 4 million as of July 09, 2021. Currently, there are few therapeutic options available for COVID-19 control. The rapid mutations in SARS-CoV-2 genome and development of new virulent strains with increased infection and mortality among COVID-19 patients, there is a great need to discover more potential drugs for SARS-CoV-2 on a priority basis. One of the key viral enzymes responsible for the replication and maturation of SARS-CoV-2 is Mpro protein. In the current study, structure-based virtual screening was used to identify four potential ligands against SARS-CoV-2 Mpro from a set of 8,722 ASINEX library compounds. These four compounds were evaluated using ADME filter to check their ADME profile and druggability, and all the four compounds were found to be within the current pharmacological acceptable range. They were individually docked to SARS-CoV-2 Mpro protein to assess their molecular interactions. Further, molecular dynamics (MD) simulations was carried out on protein-ligand complex using Desmond at 100 ns to explore their binding conformational stability. Based on RMSD, RMSF and hydrogen bond interactions, it was found that the stability of protein-ligand complex was maintained throughout the entire 100 ns simulations for all the four compounds. Some of the key ligand amino acid residues participated in stabilizing the protein-ligand interactions includes GLN 189, SER 10, GLU 166, ASN 142 with PHE 66 and TRP 132 of SARS-CoV-2 Mpro. Further optimization of these compounds could lead to promising drug candidates for SARS-CoV-2 Mpro target.

A new series of quinoline derivatives has been designed and synthesized as probable protease inhibitors (PIs) against severe acute respiratory syndrome coronavirus 2. In silico studies using DS v20.1.0.19295 software have shown that these compounds behaved as PIs while interacting at the allosteric site of target Mpro enzyme (6LU7). The designed compounds have shown promising docking results, which revealed that all compounds formed hydrogen bonds with His41, His164, Glu166, Tyr54, Asp187, and showed π-interaction with His41, the highly conserved amino acids in the target protein. Toxicity Prediction by Komputer Assisted Technology results confirmed that the compounds were found to be less toxic than the reference drug. Further, molecular dynamics simulations were performed on compound 5 and remdesivir with protease enzyme. Analysis of conformational stability, residue flexibility, compactness, hydrogen bonding, solvent accessible surface area (SASA), and binding free energy revealed comparable stability of protease:5 complex to the protease: remdesivir complex. The result of hydrogen bonding showed a large number of intermolecular hydrogen bonds formed between protein residues (Glu166 and Gln189) and ligand 5, indicating strong interaction, which validated the docking result. Further, compactness analysis, SASA and interactions like hydrogen-bonding demonstrated inhibitory properties of compound 5 similar to the existing reference drug. Thus, the designed compound 5 might act as a potential inhibitor against the protease enzyme.Communicated by Ramaswamy H. SarmaHighlightsQuinoline derivatives have been designed as protease inhibitors against SARS-CoV-2.The compounds were docked at the allosteric site of SARS-CoV-2-Mpro enzyme (PDB ID: 6LU7) to study the stability of protein-ligand complex.Docking studies indicated the stable ligand-protein complexes for all designed compounds.The Toxicity Prediction by Komputer Assisted Technology protocol in DS v20.1.0.19295 software was used to evaluate the toxicity of the designed quinoline derivatives.Molecular dynamics studies indicated the formation of stable ligand-Mpro complexes.

Antibodies are well known for their high specificity that has enabled them to be of significant use in both therapeutic and diagnostic applications. Antibodies can recognize different antigens, including proteins, carbohydrates, peptides, nucleic acids, lipids, and small molecular weight haptens that are abundantly available as hormones, pharmaceuticals, and pesticides. Here we focus on a structural analysis of hapten-antibody couples and identify potential structural movements originating from the hapten binding by comparison with unbound antibody, utilizing 40 crystal structures from the Protein Data Bank. Our analysis reveals three binding surface trends; S1 where a pocket forms to accommodate the hapten, S2 where a pocket is removed when the hapten binds, and S3 where no pockets changes are found. S1 and S2 are expected for induced-fit binding, whereas S3 indicates that a pre-existing population of optimal binding antibody conformation exists. The structural analysis reveals four classifications of structural reorganization, some of which correlate to S2 but not to the other binding surface changes. These observations demonstrate the complexity of the antibody-antigen interaction, where structural changes can be restricted to the binding sites, or extend through the constant domains to propagate structural changes. This highlights the importance of structural analysis to ensure successful and compatible transformation of small antibody fragments at the early discovery stage into full antibodies during the subsequent development stages, where long-range structural changes are required for an Fc effector response.

Rickettsia rickettsii, the etiological agent of Rocky Mountain spotted fever (RMSF), a life-threatening tick-borne disease that affects humans and various animal species, has been recognized in medicine and science for more than 100 years. Isolate-dependent differences in virulence of R. rickettsii have been documented for many decades; nonetheless, the specific genetic and phenotypic factors responsible for these differences have not been characterized. Using in vivo and in vitro methods, we identified multiple phenotypic differences among six geographically distinct isolates of R. rickettsii, representing isolates from the United States, Costa Rica, and Brazil. Aggregate phenotypic data, derived from growth in Vero E6 cells and from clinical and pathological characteristics following infection of male guinea pigs (Cavia porcellus), allowed separation of these isolates into three categories: nonvirulent (Iowa), mildly virulent (Sawtooth and Gila), and highly virulent (Sheila SmithT, Costa Rica, and Taiaçu). Transcriptional profiles of 11 recognized or putative virulence factors confirmed the isolate-dependent differences between mildly and highly virulent isolates. These data corroborate previous qualitative assessments of strain virulence and suggest further that a critical and previously underappreciated balance between bacterial growth and host immune response could leverage strain pathogenicity. Also, this work provides insight into isolate-specific microbiological factors that contribute to the outcome of RMSF and confirms the hypothesis that distinct rickettsial isolates also differ phenotypically, which could influence the severity of disease in vertebrate hosts.

Epigenetic changes, histone acetylation and deacetylation in chromatin have been intensively studied due to their significance in regulating the gene expression. According to the type of tumor, the levels of histone deacetylases (HDAC) are varied. HDAC inhibitors are a new promising class of compounds that inhibit the proliferation of tumor cells. In this study, the inhibitory efficacy of some HDAC inhibitors such as vorinostat, panobinostat, abexinostat, belinostat, resminostat, dacinostat and pracinostat was studied using molecular dynamics simulation. The inhibitory efficacy was estimated by computing the enzyme\'s stability, positional stability of the individual amino acids and interaction energies of HDLP-inhibitor complexes. It is hoped that this investigation may improve our understanding of the atomic-level description of the inhibitor binding site and how the HDAC inhibitors change the environment of the enzyme\'s active site. The results obtained from the root-mean-square deviation, the radius of gyration, solvent-accessible surface area, root-mean-square fluctuation, stride server and Ramachandran plot have revealed that the stability of HDLP enzyme with vorinostat, panobinostat and abexinostat is higher than the other studied complexes. According to the calculated values for MM-PBSA, LIE, semi-LIE binding free energies and interaction energies, the stability of the HDLP enzyme varies as panobinostat > abexinostat > vorinostat where resminostat complex showed relatively low stability. The ligandability and drugability values also give the same trend as above. The findings revealed that the panobinostat and abexinostat are potential lead compounds as reference inhibitor vorinostat. Therefore, it is possible to use these drugs as HDAC inhibitors in clinical practices. Also, the outcomes of this study could be utilized to identify new inhibitors for clinical research.Communicated by Ramaswamy H. Sarma.

Prenatal COVID infection is one of the worst affected and least attended aspects of the COVID-19 disease. Like other coronaviruses, CoV2 infection is anticipated to affect fetal development by maternal inflammatory response on the fetus and placenta. Studies showed that higher prenatal choline level in mother\'s body can safeguard the developing brain of the fetus from the adverse effects of CoV2 infection. Choline is commonly used as food supplement. By virtual screening, molecular docking and molecular dynamics techniques, we have established a strong inhibitory possibility of choline for SARS 3CLpro protease which may provide a lead for prenatal COVID-19 treatment.

Protein 3D structures support their biological functions. As the number of protein structures is negligible in regards to the number of available protein sequences, prediction methodologies relying only on protein sequences are essential tools. In this field, protein secondary structure prediction (PSSPs) is a mature area, and is considered to have reached a plateau. Nonetheless, proteins are highly dynamical macromolecules, a property that could impact the PSSP methods. Indeed, in a previous study, the stability of local protein conformations was evaluated demonstrating that some regions easily changed to another type of secondary structure. The protein sequences of this dataset were used by PSSPs and their results compared to molecular dynamics to investigate their potential impact on the quality of the secondary structure prediction. Interestingly, a direct link is observed between the quality of the prediction and the stability of the assignment to the secondary structure state. The more stable a local protein conformation is, the better the prediction will be. The secondary structure assignment not taken from the crystallized structures but from the conformations observed during the dynamics slightly increase the quality of the secondary structure prediction. These results show that evaluation of PSSPs can be done differently, but also that the notion of dynamics can be included in development of PSSPs and other approaches such as de novo approaches.

Alzheimer\'s disease (AD) is the most common type of dementia and usually manifests as diminished episodic memory and cognitive functions. Caspases are crucial mediators of neuronal death in a number of neurodegenerative diseases, and caspase 8 is considered a major therapeutic target in the context of AD. In the present study, we performed a virtual screening of 200 natural compounds by molecular docking with respect to their abilities to bind with caspase 8. Among them, rutaecarpine was found to have the highest (negative) binding energy (-6.5 kcal/mol) and was further subjected to molecular dynamics (MD) simulation analysis. Caspase 8 was determined to interact with rutaecarpine through five amino acid residues, specifically Thr337, Lys353, Val354, Phe355, and Phe356, and two hydrogen bonds (ligand: H35-A: LYS353:O and A:PHE355: N-ligand: N5). Furthermore, a 50 ns MD simulation was conducted to optimize the interaction, to predict complex flexibility, and to investigate the stability of the caspase 8-rutaecarpine complex, which appeared to be quite stable. The obtained results propose that rutaecarpine could be a lead compound that bears remarkable anti-Alzheimer\'s potential against caspase 8.