Search

Physics and physical chemistry are an important thread in computational protein design, complementary to knowledge-based tools. They provide molecular mechanics scoring functions that need little or no ad hoc parameter readjustment, methods to thoroughly sample equilibrium ensembles, and different levels of approximation for conformational flexibility. They led recently to the successful redesign of a small protein using a physics-based folded state energy. Adaptive Monte Carlo or molecular dynamics schemes were discovered where protein variants are populated as per their ligand-binding free energy or catalytic efficiency. Molecular dynamics have been used for backbone flexibility. Implicit solvent models have been refined, polarizable force fields applied, and many physical insights obtained.

Human serum albumin (HSA) is the most prevalent protein in the blood plasma which binds an array of exogenous compounds. Drug binding to HSA is an important consideration when developing new therapeutic molecules, and it also aids in understanding the underlying mechanisms that govern their pharmacological effects. This study aims to investigate the molecular binding of coronavirus disease 2019 (COVID-19) therapeutic candidate molecules to HSA and to identify their putative binding sites. Binding energies and interacting residues were used to evaluate the molecular interaction. Four drug candidate molecules (β-D-N4-hydroxycytidine, Chloroquine, Disulfiram, and Carmofur) demonstrate weak binding to HSA, with binding energies ranging from -5 to -6.7 kcal/mol. Ivermectin, Hydroxychloroquine, Remdesivir, Arbidol, and other twenty drug molecules with binding energies ranging from -6.9 to -9.5 kcal/mol demonstrated moderate binding to HSA. The strong HSA binding drug candidates consist of fourteen molecules (Saquinavir, Ritonavir, Dihydroergotamine, Daclatasvir, Paritaprevir etc.) with binding energies ranging from -9.7 to -12.1 kcal/mol. All these molecules bind to different HSA subdomains (IA, IB, IIA, IIB, IIIA, and IIIB) through molecular forces such as hydrogen bonds and hydrophobic interactions. Various pharmacokinetic properties (gastrointestinal absorption, blood-brain barrier permeation, P-glycoprotein substrate, and cytochrome P450 inhibitor) of each molecule were determined using SwissADME program. Further, the stability of the HSA-ligand complexes was analyzed through 100 ns molecular dynamics simulations considering various geometric properties. The binding free energy between free HSA and compounds were calculated using Molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) and molecular mechanics generalized Born surface area (MM/GBSA) approach. The findings of this study might be useful in understanding the mechanism of COVID-19 drug candidates binding to serum albumin protein, as well as their pharmacodynamics and pharmacokinetics.

Background Dipeptidyl peptidase-IV (DPP-IV) inhibitors, the enhancers of incretin are used for the treatment of diabetes. The non-glycaemic actions of these drugs (under developmental stage) also proved that repurposing of these molecules may be advantageous for other few complicated disorders like cardiovascular diseases, Parkinson's disease, Alzheimer's disease, etc. Objective The present study was aimed to investigate the DPP-IV inhibitory potential of Calebin-A, one of the constituents of Curcuma longa. Material and methods The phytoconstituent was subjected for various in silico studies (using Schrödinger Suite) like, Docking analysis, molecular mechanics combined with generalized Born model and solvent accessibility method (MMGBSA) and Induced fit docking (IFD) after validating the protein using Ramachandran plot. Further, the protein-ligand complex was subjected to molecular dynamic simulation studies for 50 nanoseconds. And finally, the results were confirmed through enzyme inhibition study. Results Insilico results revealed possible inhibitory binding interactions in the catalytic pocket (importantly Glu205, Glu206 and Tyr 662 etc.) and binding affinity in terms of glide g-score and MMGBSA dG bind values were found to be −6.2 kcal/mol and −98.721 kcal/mol. Further, the inhibitory action towards the enzyme was confirmed by an enzyme inhibition assay, in which it showed dose-dependent inhibition, with maximum % inhibition of 55.9 at 26.3 μM. From molecular dynamic studies (50 nanoseconds), it was understood that Calebin A was found to be stable for about 30 nanoseconds in maintaining inhibitory interactions. Conclusion From the in silico and in vitro analysis, the current research emphasizes the consideration of Calebin A to be as a promising or lead compound for the treatment of several ailments where DPP-IV action is culprit., Graphical abstract Image 1

Chromatin structure modeling is a rapidly developing field. Parallel to the enormous growth of available experimental data, there is a growing need of building and visualizing 3D structures of nuclei, chromosomes, chromatin domains, and single loops associated with particular gene loci. Here, we present a tool for chromatin domain modeling; it is available as a webservice and standalone python script. Our tool is based on molecular mechanics and utilizes the OpenMM engine for model generation. In this method the user provides contacts between chromatin regions and obtains a 3D structure that satisfies them. Additional parameters allow for the control of fibre stiffness, initial structure adjustments and simulation resolution, there are also options for structure refinement and modeling in a spherical container. The user may provide contacts in the form of bead indices, or insert interactions in genome coordinates sourced from BEDPE files. After the simulation is complete, the user is able to download the structure in the Protein Data Bank (PDB) format for further analysis. We dedicate this tool to all who are interested in chromatin structures. It is suitable for quick visualization of datasets, studying the impact of structural variants (SVs), inspecting the effects of adding and removing particular contacts, and measuring features such as maximum distances between sites (e.g.promoter-enhancer), or local chromatin density.

Saffron plant (Crocus sativus L.) is being used as a source of saffron spice and medicine to cure or prevent different types of diseases including cancers. We report the isolation, characterization of bioactive small molecule ([crocetin (β-d-glucosyl) ester] from the leaf biowastes of saffron plant of Kashmir, India. MTTC assay and Bio-autography aided approach were used to assess anti-oxidant activity and anti-cancer properties of crocin (s) against DPPH free radical and breast cancer cell line respectively. Crocetin beta-d-glucosyl ester restrained proliferation of human breast adeno-carcinoma cell model (MCF-7) without significantly affecting normal cell line (L-6). Further studies involving molecular mechanics generalized born surface area and molecular docking showed that crocetin beta-d-glucosyl ester exhibits strong affinity for estrogen receptor alpha and histone deacetylase 2 (crucial receptors involved in breast cancer signalling) as evidenced by the negative docking score and binding free energy (BFE) values. Therefore, crocetin beta-d-glucosyl ester from Crocus sativus biowastes showed antiproliferative effect possibly by inhibiting estrogen receptor alpha and HDAC2 mediated signalling cascade. Keywords: Saffron, Floral biowastes, Antioxidant, Crocetin beta-d-glucosyl ester, Breast cancer, Molecular docking

Solvents play a crucial role in electrocatalytic reactions, yet existing simulation methods for the catalyst-solvent interface are either too slow or too inaccurate. Goddard et. al. report a quantum mechanics/molecular mechanics method for electrocatalysis simulations which they show to be accurate and efficient.

Dehalogenases continue to garner interest of the scientific community due to their potential applications in bioremediation of halogen-contaminated environment and in synthesis of various industrially relevant products. Example of such enzymes is DehL, an L-2-haloacid dehalogenase (EC 3.8.1.2) from Rhizobium sp. RC1 that catalyses the specific cleavage of halide ion from L-2-halocarboxylic acids to produce the corresponding D-2-hydroxycarboxylic acids. Recently, the catalytic residues of DehL have been identified and its catalytic mechanism has been fully elucidated. However, the enantiospecificity determinants of the enzyme remain unclear. This information alongside a well-defined catalytic mechanism are required for rational engineering of DehL for substrate enantiospecificity. Therefore, using quantum mechanics/molecular mechanics and molecular mechanics Poisson-Boltzmann surface area calculations, the current study theoretically investigated the molecular basis of DehL enantiospecificity. The study found that R51L mutation cancelled out the dehalogenation activity of DehL towards it natural substrate, L-2-chloropropionate. The M48R mutation, however introduced a new activity towards D-2-chloropropionate, conveying the possibility of inverting the enantiospecificity of DehL from L-to d-enantiomer with a minimum of two simultaneous mutations. The findings presented here will play important role in the rational design of DehL dehalogenase for improving substrate utility.

Lack of perfect insulin signaling can lead to the insulin resistance, which is the hallmark of diabetes mellitus. Activation of insulin and its binding to the receptor for signaling process initiates via B-chain C-terminal hinge conformational change through an open structure to “wide-open” conformation. Observational studies and basic scientific evidence suggest that vitamin D and E directly and/or indirectly prevent diabetes through improving glucose secretion and tolerance, activating calcium dependent endopeptidases and thus improving insulin exocytosis, antioxidant effect and reducing insulin resistance. On the contrary, clinical trials have yielded inconsistent results about the efficacy of vitamin D supplementations for the control of glucose hemostasis. In this work, best binding modes of vitamin D3 and E on insulin obtained from AutoDock Vina were selected for Molecular Dynamic, MD, study. The binding energy obtained from Molecular Mechanics- Poisson Boltzman Surface Area, MM-PBSA method, revealed that Vitamins D3 and E have good affinity to bind to the insulin and vitamin E has higher binding energy (−46 kj/mol) by engaging more residues in binding site. Distance and angle calculation results illustrated that vitamin E changes the B-chain conformation and it causes the formation of wide-open/active form of insulin. Vitamin E increases the ValB12-TyrB26 distance to ∼15 A and changes the hinge angle to ∼65°. Consequently, essential hydrophobic residues for binding to insulin receptor exposed to surface in the presence of vitamin E. However, our data illustrated that vitamin D3 cannot change B-chain conformation. Thus our MD simulations propose a model for insulin activation through vitamin E interaction for therapeutic approaches.

The fully optimized geometrical structures of the CaMn4Ox (x = 5, 6) clusters in the Si (i = 0–3) states of the Kok cycle of photosynthetic water oxidation by the large-scale quantum mechanics/molecular mechanics (QM/MM) calculations were compared to recent experimental results based on serial femtosecond crystallography (SFX). The Mn Mn and Ca Mn distances obtained by the QM/MM calculations were found to be totally comparable to the SFX experiments, elucidating the entire Kok cycle involving the S4 transition state during the O O bond formation for water oxidation in the oxygen evolving complex (OEC) of photosystem II (PSII).

Amyotrophic lateral sclerosis (ALS) is a calamitous neurodegenerative disorder characterized by denervation of upper and lower motor neurons. Numerous hypotheses suggest that toxic protein misfolding and aggregation cause ALS, similar to that of other neurodegenerative diseases, such as Alzheimers and Parkinsons. Protruding causes of familial ALS are mutations in the gene encoding Cu/Zn superoxide dismutase-1 (SOD1), which decrease protein stability and endorse protein aggregation. Thus, the interference concerning aggregate formation and destabilization in SOD1 is considered to be an impending therapeutic strategy. In this work, we utilized computational chemistry methods to initially study the effect of substitution mutation, His46Arg on SOD1 protein. Further, we described the interaction of two naturally occurring polyphenol compounds, naringin and naringenin with mutant SOD1 that is regarded to hinder the protein aggregation. Subsequently, the use of quantum chemical and molecular mechanics calculations speculated that naringin had a strong binding affinity with mutant SOD1 and impeded the formation of toxic aggregates than that of naringenin. Ultimately, we could conjecture that ingesting of polyphenol-rich foods in ALS patients may be regarded to improvise their living. Moreover, the findings from our study could pave a way in the field of structure-based drug design in developing potential anti-aggregation inhibitors against incurable ALS, affecting the human population.

Multifunctional light-emitting molecular materials have attracted great attention recently, while the amounts are still quite limited. The Cz-AQ molecule with thermally activated delayed fluorescence, aggregation-induced emission and mechanochromic luminescence properties was reported and the bicolor electroluminescence was found in film recently by Jiang's group (J. Mater. Chem. C 5 (2017) 12031–12034). In this paper, the bicolor switching mechanism of the Cz-AQ molecule in solid phase is studied based on first-principles and combined quantum mechanics and molecular mechanics calculations. Our results indicate that the emission of Cz-AQ in solid phase does not come from monomers but dimers. The emission colors have close relationship with intermolecular positions of two molecules in dimers. There is no need of conversion of lattice types for color switching, and slight variation of relative position of molecules in dimers will induce the emission color change between red and yellow. Our theoretical study not only gives reasonable explanation of experimental results but also provides some insights on the design of new type multifunctional light-emitting molecules.

The novel prophylactic agent 7-methoxytacrine-4-pyridinealdoxime is a hybrid compound formerly designed to keep acetylcholinesterase resistant to organophosphates by reactivating it in case of intoxication by such inhibitors. In rational design, a 5-carbon length-spacer hybrid compound was synthesized to evaluate its inhibitory and reactivation capabilities. In this work, theoretical results were achieved through molecular modelling techniques, taking for granted the enzymatic reactivation reaction through nucleophilic substitution. Based on the near attack conformation approach, docking studies were performed to assess the spacer-length from 1 to 10 carbons long of a series of analogues of 7-methoxytacrine-4-pyridinealdoxime. Consequently, the hybrids with length-spacer of 4 and 5 carbons long were the best assessed and subsequently subjected to further molecular dynamics simulations, complemented by Poisson-Boltzmann surface area calculations. As a result, intermolecular interactions with the different binding sites inside human acetylcholinesterase were elucidated. Besides, thermodynamics and kinetics concepts pointed to the 4-carbon linker as optimum for enzymatic reactivation. Further studies, based on quantum mechanics in conjunction with molecular mechanics, were recommended to the presented near attack conformations to achieve more thermodynamics results between the hybrids with 4- and 5-carbon linkers, like values of activation energy for the reactivation reaction. All of those in silico evaluations could be considered as a set of tools for theoretically investigate novel enzymatic reactivators with different shape of spacers.

Transcriptional enhancer activation domain (TEAD) proteins are the downstream transcriptional factor of the Hippo pathway. The transcription co-activators Yes-associated protein (YAP) and its paralog transcription co-activators with PDZ-binding motif (TAZ), binding to TEAD to promote transcription of genes in cell proliferation and anti-apoptosis, are key effectors of the Hippo pathway. TEAD4, one member of TEAD proteins, is specifically required in embryo implantation. The recently reported crystal structure of TEAD4-TAZ complex (PDB Code 5GN0) in mouse reveals that the interactions between the two helices of YAP/TAZ and TEAD4 are highly conserved. Point mutation of the residue Tyr422 of TEAD4 protein would disrupt the relevant hydrogen bond and even abolish the interaction. However, detailed information affected by the mutation at the atom level are still unrevealed. Molecular dynamics (MD) simulations and the molecular mechanics/Generalized-Born surface area (MM/GBSA) free energy calculations were used to explore the effects of mutation Tyr422Ala on the structural flexibility and conformational dynamics. The non-polar interactions play an indispensable role in the binding process of TEAD4 and YAP/TAZ. The helices α1 and α2 of YAP/TAZ provide a primary function to anchor YAP/TAZ well bound to TEAD4. The mutation Tyr422Ala disrupts the hydrogen-bonding network but do not obviously influence the secondary structure stability of TEAD4. The binding conformation of YAP/TAZ distorted by decreased non-polar interaction and the lost hydrogen bonds would lead to reduced interaction activity. The present study would provide important insights into the structure-function relationships of TEAD protein and give a new explanation for the affinity of YAP/TAZ with TEAD.

In order to study the effects of peptide exposure to oxidative attack, we chose a model reaction in which the hydroxyl radical discretely abstracts a hydrogen atom from the α-carbon of each residue of a highly amyloidogenic region in the human calcitonin hormone, hCT15-19. Based on a combined Molecular Mechanics / Quantum Mechanics approach, the extended and folded L- and D-configuration and radical intermediate hCT15-19 peptides were optimized to obtain their compactness, secondary structure and relative thermodynamic data. The results suggest that the epimerization of residues is generally an exergonic process that can explain the cumulative nature of molecular aging. Moreover, the configurational inversion induced conformational changes can cause protein dysfunction. The epimerization of the central residue to the D-configuration induced a hairpin structure in hCT15-19, concomitant with a possible oligomerization of human calcitonin into Aβ(1–42)-like amyloid fibrils present in patients suffering from Alzheimer’s disease.

High production cost, instability, low tumor penetration are some of the shortcomings that have characterized and undermined the use of antibodies as a target for Cytotoxic T-lymphocytes associated protein 4 (CTLA-4). Design and discovery of small molecule inhibitors have therefore become a sine qua non in targeting immune proteins implicated in immune disorders. In this study, we utilized a drug repositioning approach to explore the characteristic feature of unrelated proteins to have similar binding sites and the promiscuity of drugs to repurpose an existing drug to target CTLA-4. CTLA-4 and Kallikrein-7 were found to have similar binding sites, we therefore used 1, 3, 6-trisubstituted 1, 4-diazepane-7-ones (TDSO) which is an inhibitor of Kallikrein-7 as our lead compound. High throughput screening using TDSO as a lead compound resulted in 9 hits with ZINC04515726 and ZINC08985213 having the highest binding score. We went ahead to investigate the interaction of these compounds with CTLA-4 by conducting a molecular dynamic simulation. Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) estimations revealed that TDSO had the highest binding energy value of -28.51Kcal/mol, with ZINC04515726 and ZINC08985213 having -23.76Kcal/mol and -21.03Kcal/mol respectively. The per-residue decomposition highlighted Tyr24, Ala25, Gly28, Ala30, Tyr53 and Asn72 as having significantly high electrostatic energy contributions and the main contributing residues to the binding of TDSO, ZINC04515726 and ZINC08985213 to Cytotoxic T lymphocytes CTLA-4. Summarily, from the results gathered, we proposed that TDSO can be an effective immune check point small molecule inhibitor against the suppression of T-cell activation, proliferation, and tumor cell eradication.

Catechol-containing polyphenols present in coffee and tea, while serving as excellent substrates for catechol-O-methyltransferase (COMT)-catalyzed O-methylation, can also operate as COMT inhibitors. However, little is known about the relationship between COMT and the characteristic phenolics present in extra virgin olive oil (EVOO). We here selected the EVOO dihydroxy-phenol oleacein for a computational study of COMT-driven methylation using classic molecular docking/molecular dynamics simulations and hybrid quantum mechanical/molecular mechanics, which were supported by in vitro activity studies using human COMT. Oleacein could be superimposed onto the catechol-binding site of COMT, maintaining the interactions with the atomic positions involved in methyl transfer from the S-adenosyl-L-methionine cofactor. The transition state structure for the meta-methylation in the O5 position of the oleacein benzenediol moiety was predicted to occur preferentially. Enzyme analysis of the conversion ratio of catechol to O-alkylated guaiacol confirmed the inhibitory effect of oleacein on human COMT, which remained unaltered when tested against the protein version encoded by the functional Val158Met polymorphism of the COMT gene. Our study provides a theoretical determination of how EVOO dihydroxy-phenols can be metabolized via COMT. The ability of oleacein to inhibit COMT adds a new dimension to the physiological and therapeutic utility of EVOO secoiridoids.

A molecular mechanics model for covalent and ionic double-network polymer composites was developed in this study to investigate mechanisms of mechanochemically induced unzipping and scission of chemical bonds. Morse potential function was firstly applied to investigate mechanical unzipping of the covalent bonds, and then stress-dependent mechanical energy for the interatomic covalent bonds was discussed. A new mechanochemical model was formulated for describing the mechanically induced ionic bond scissions based on the Morse potential model and equations for electrostatic forces. Based on this newly proposed model, mechanochemical behaviors of several common interatomic interaction types (e.g., A+B−, A2+B2−/A2+2B−/2A+B2− and A3+B3−/A3+3B−/3A+B3−) of the ionic bonds have been quantitatively described and analyzed. Finally, mechanochemical unzipping of the covalent bonds and dissociation of the ionic bonds have been characterized and analyzed based on the molecular mechanics model, which has well predicted the chemical and mechanochemical activations in the covalent and ionic double-network polymer composites.

The field of molecular mechanics studies of proteins has developed enormously since its origin in the 1970’s, and many applications and methodologies have branched from the original idea of the force field. The applications of such methodologies are far spread and commonplace in neuroscience research today. In this mini-review, we outline the main methodologies applied when studying events ranging from ligands binding within small binding sites, through overall large-scale conformational changes, to the even larger-scale oligomerization events of neurological membrane proteins. The limitations and caveats of the methods are discussed, while examples of recent applications are described and their implications discussed. We have chosen to focus on the monoamine transporters throughout, with a few examples from neurological membrane proteins such as ionotropic and metabotropic neurotransmitter receptors.

Plastic materials are directly or indirectly emitted into the environment arousing increasing environmental concerns. This work focused on investigating the binding and the degradation mechanism of biodegradable plastics polyvinyl alcohol by oxidized polyvinyl alcohol hydrolase (OPH) purified from Pseudomonas. Molecular dynamics (MD) and quantum mechanics/molecular mechanics calculations (QM/MM) were performed for the investigation. The formation of hydrogen bonds between Gly65, Ser66, Val67 with ligand exert extremely important promotion effect on the reaction. Besides, significant electrostatic influence of Glu198 also contributes to the catalysis. These results provide valuable guidance to future mutation studies and enzyme engineering for efficient degradation of plastics.

The initial Mn2+ binding site in photoassembly of the Mn4CaO5 cluster, the catalytic center of water oxidation in photosystem II, was studied using quantum mechanics/molecular mechanics calculations. Among the five metal sites (Mn1–Mn4 and Ca) of the Mn4CaO5 cluster, the Mn1 site involving D1-H332 was most stable for Mn2+ binding, while the Mn2 site showed a slightly higher energy comparable to the thermal fluctuation energy at room temperature. It is thus suggested that Mn2+ binding at the Mn1 site, likely in equilibrium with that at the Mn2 site, is the initial step of photoassembly of the Mn4CaO5 cluster.

Survival of the Acinetobacter baumannii inside host requires different micronutrients such as iron, but their bioavailability is limited because of nutritional immunity created by host. A. baumannii has to develop mechanisms to acquire nutrient iron during infection. The present study is an attempt to identify membrane proteins involved in iron sequestration mechanism of A. baumannii using two-dimensional electrophoresis and LC-MS/MS analysis. The identified iron-regulated membrane protein (IRMP) of A. baumannii was used for its interaction studies with different siderophores, and designing of the inhibitor against A. baumannii targeting this IRMP. Membrane proteomic results identified over-expression of four membrane proteins (Fhu-E receptor, ferric-acinetobactin receptor, ferrienterochelin receptor, and ferric siderophore receptor) under iron-limited condition. A. baumannii produces siderophores that have good interaction with the FhuE receptor. Result also showed that FhuE receptor has interaction with siderophores produced by other bacteria. Interaction of FhuE receptor and siderophores helps in iron sequestration and survival of Acinetobacter under nutritional immunity imposed by the host. Hence it becomes essential to find a potential inhibitor for the FhuE receptor that can inhibit the survival of A. baumannii in the host. In-silico screening, and molecular mechanics studies identified ZINC03794794 and ZINC01530652 as a likely lead to design inhibitor against the FhuE receptor of A. baumannii. The designed inhibitor is experimentally validated for its antibacterial activity on the A. baumannii. Therefore, designed inhibitor interferes with the iron acquisition mechanism of Acinetobacter hence may prove useful for preventing infection caused by A. baumannii by limiting nutrient availability.

Scaffold known as designed ankyrin repeat proteins (DARPins) have been utilized by scientists extensively in protein binding, as they acquire remarkable binding properties. Trimodular ankyrin repeat protein (Ank1D4) is a potential macromolecular and antiviral intracellular inhibitor that aims the N-terminal domain capsid protein (NTDCA) of viral Human Immunodeficiency Virus type-1 (HIV-1) at the virus assembly and budding machinery process. The protein-protein interaction between Ank1D4 and the viral capsid CA alpha helixes H1 and H7 have been studied in depth and point mutation strategy at the S451D4 key residue position was conducted to further improve the binding interaction and efficacy on the wildtype Ank1D4 structure with NTDCA. A computational approach via molecular dynamics (MD) simulation has been fully utilized in this study to simulate both the wildtype and potential novel complexes. Free energy contributions on the modified Ank1D4-NTDCA (H1, H7) complexes were acquired by using vigorous scoring functions of implicit solvent models; Molecular Mechanics/Poisson-Boltzmann surface area (MM-PBSA) and Molecular Mechanics/Generalized Born surface area (MM-GBSA), respectively. The wildtype structure of Ank1D4-NTDCA H7 has been proven to possess better binding stability than the H1 complex. From there, two mutants of H7 complex were generated and the computed ∆Gbind values based on both models have revealed that S45Y-H7 has the lowest binding free energy compared to both S45W-H7 and S45-H7 (Wildtype). The binding strength predictions given by both MM-PBSA and MM-GBSA in the following ranking order for the three systems in terms of the value of ∆Gbind: S45Y-H7 > S45 W-H7 > S45-H7. These findings provide a good basis in directing us a step further to discover better novel treatment modalities to treat millions of HIV-1 patients.

Tuberculosis (TB) is one of the world's deadliest diseases resulting from infection by the bacterium, Mycobacterium tuberculosis (M.tb). The L,D-transpeptidase enzymes catalyze the synthesis of 3 → 3 transpeptide linkages which are predominant in the peptidoglycan of the M.tb cell wall. Carbapenems is class of β-lactams that inactivate L,D-transpeptidases by acylation, although differences in antibiotic side chains modulate drug binding and acylation rates. Herein, we used a two-layered our Own N-layer integrated Molecular Mechanics ONIOM method to investigate the catalytic mechanism of L,D-transpeptidase 5 (LdtMt5) by β-lactam derivatives. LdtMt5 complexes with six β-lactams, ZINC03788344 (1), ZINC02462884 (2), ZINC03791246 (3), ZINC03808351 (4), ZINC03784242 (5) and ZINC02475683 (6) were simulated. The QM region (high-level) comprises the β-lactam, one water molecule and the Cys360 catalytic residue, while the rest of the LdtMt5 residues were treated with AMBER force field. The activation energies (ΔG#) were calculated with B3LYP, M06-2X and ωB97X density functionals with 6–311++G(2d, 2p) basis set. The ΔG# for the acylation of LdtMt5 by the selected β-lactams were obtained as 13.67, 20.90, 22.88, 24.29, 27.86 and 28.26 kcal mol−1respectively. Several of the compounds showed an improved ΔG# when compared to the previously calculated energies for imipenem and meropenem for the acylation step for LdtMt5. This model provides further validation of the catalytic inhibition mechanism of LDTs with atomistic detail.

A series of new indole-3-acetic acid (IAA)-tacrine hybrids as dual acetylcholinesterase (AChE)/butyrylcholinesterase (BChE) inhibitors were designed and prepared based on the molecular docking mode of AChE with an IAA derivative (1a), a moderate AChE inhibitor identified by screening our compound library for anti-Alzheimer's disease (AD) drug leads. The enzyme assay results revealed that some hybrids, e.g. 5d and 5e, displayed potent dual in vitro inhibitory activities against AChE/BChE with IC50 values in low nanomolar range. Molecular modeling studies in tandem with kinetic analysis suggest that these hybrids target both catalytic active site and peripheral anionic site of cholinesterase (ChE). Molecular dynamic simulations and Molecular Mechanics/Poisson-Boltzmann Surface Area (MM-PBSA) calculations indicate that 5e has more potent binding affinity than hit 1a, which may explain the stronger inhibitory effect of 5e on AChE. Furthermore, their predicted pharmacokinetic properties and in vitro influences on mouse brain neural network electrical activity were discussed. Taken together, compound 5e can be highlighted as a lead compound worthy of further optimization for designing new anti-AD drugs.

The ecdysone receptor (EcR) is an insect nuclear receptor that is activated by the molting hormone, 20-hydroxyecdysone. Because synthetic EcR ligands disrupt the normal growth of insects, they are attractive candidates for new insecticides. In this study, the Molecular Mechanics/Poisson–Boltzmann Surface Area (MM/PBSA) method was used to predict the binding activity of EcR ligands. Validity analyses using 40 known EcR ligands showed that the binding activity was satisfactorily predicted when the ligand conformational free energy term was introduced. Subsequently, this MM/PBSA method was applied to structure-based hierarchical virtual screening, and 12 candidate compounds were selected from a database of 3.8 million compounds. Five of these compounds were active in a cell-based competitive binding assay. The most potent compound is a simple proline derivative with low micromolar binding activity, representing a valuable lead compound for further structural optimization.

Coordination between plant development and nutrient availability ensures a suitable supply of macromolecules for growth and developmental programs. Nitrate is an important source of nitrogen (N) that acts as a signal molecule to modulate gene expression, physiological, growth and developmental responses throughout the life of the plant. New key players in the nitrate signaling pathway have been described and knowledge of the molecular mechanics of how it impacts growth and developmental processes is increasing fast. Importantly, mechanisms for nitrate-control of growth and developmental processes have been proposed for both local as well as systemic responses. This article provides a synthesis of recent insights into molecular mechanisms by which nitrate impacts growth and development over Arabidopsis life-cycle.

We here precisely analyzed the electronic structures during the oxygen-oxygen (O-O) bond formation of the oxygen-evolving complex (OEC) in photosystem II (PSII) based on the quantum mechanics/molecular mechanics (QM/MM) calculations. Changes in the spin density and the chemical indices derived from the natural orbitals suggested that the O(5)-O(6) bond formation proceeds with a concerted bond breaking of the Mn4-O(5) coupled with the through-bond one-electron transfer from the O(6) site to the Mn4 site (Mn4(IV) → Mn4(III)). This mechanism can be best described as a Concerted Bond Switching (CBS) mechanism, which is a non-radical reaction of the O(5)-O(6) bond formation.

Coronaviruses are known to infect respiratory tract and intestine. These viruses possess highly conserved viral macro domain A1pp having adenosine diphosphate (ADP)-ribose binding and phosphatase activity sites. A1pp inhibits adenosine diphosphate (ADP)-ribosylation in the host and promotes viral infection and pathogenesis. We performed in silico screening of FDA approved drugs and nucleoside analogue library against the recently reported crystal structure of SARS-CoV-2 A1pp domain. Docking scores and interaction profile analyses exhibited strong binding affinity of eleven FDA approved drugs and five nucleoside analogues NA1 (−13.84), nadide (−13.65), citicholine (−13.54), NA2 (−12.42), and NA3 (−12.27). The lead compound NA1 exhibited significant hydrogen bonding and hydrophobic interaction at the natural substrate binding site. The root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent accessible surface (SASA), hydrogen bond formation, principle component analysis, and free energy landscape calculations for NA1 bound protein displayed stable complex formation in 100 ns molecular dynamics simulation, compared to unbound macro domain and natural substrate adenosine-5-diphosphoribose bound macro domain that served as a positive control. The molecular mechanics Poisson–Boltzmann surface area analysis of NA1 demonstrated binding free energy of −175.978 ± 0.401 kJ/mol in comparison to natural substrate which had binding free energy of −133.403 ± 14.103 kJ/mol. In silico analysis by modelling tool ADMET and prediction of biological activity of these compounds further validated them as putative therapeutic molecules against SARS-CoV-2. Taken together, this study offers NA1 as a lead SARS-CoV-2 A1pp domain inhibitor for future testing and development as therapeutics against human coronavirus., Highlights • SARS-CoV-2 ADP-Ribose phosphatase is a drugable target. • FDA approved drugs act as potential inhibitors of ADP-Ribose phosphatase. • Nucleoside inhibitor binds and stabilizes the ADP-Ribose phosphatase.

To assess a tumor biomarker like the cyclooxygenase-2 enzyme (COX-2), non-invasive imaging techniques are powerful tools. The (non-) linear optical properties of activatable fluorescent probes which are selectively bound to the biomarker can therefore be exploited. The here presented molecular modelling results based on multi-scale modelling techniques highlight the importance of the conformational versatility and of changes in the electronic interactions of such probes when they are embedded in water or in the COX-2 homodimer enzyme. The ANQ-IMC-6 probe, which combines the binding domain/scaffold of indomethacin (IMC) on COX-2 with the optical properties of acenaphtho[1,2–b]quinoxaline (ANQ), is found to be folded in the solvent and unfolded in the enzyme. A concerted movement of the probe and the protein is seen, while the rotational autocorrelation function exhibits also the intrinsic properties of the probe. Hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) calculations are used to simulate the one-photon and two-photon absorption spectra along with the first hyperpolarizability. The transition has a local character in vacuum, but changes to a charge transfer one in the presence of the microenvironment of the enzyme. This is also visible through a change of the shape of the absorption spectrum, while at the same time the simulated signals of second harmonic generation experiments are strongly enhanced. The results of this work prove that an environment sensitive probe with an anchoring group and an optical active part can be constructed for use in absorption spectroscopy, without the need to revert to fluorescence experiments.

Angiotensin-converting enzyme 2 (ACE2) is a membrane-bound zinc metallopeptidase that generates the vasodilatory peptide angiotensin 1–7 and thus performs a protective role in heart disease. It is considered an important therapeutic target in controlling the COVID-19 outbreak, since SARS-CoV-2 enters permissive cells via an ACE2-mediated mechanism. The present in silico study attempted to repurpose existing drugs for use as prospective viral-entry inhibitors targeting human ACE2. Initially, a clinically approved drug library of 7,173 ligands was screened against the receptor using molecular docking, followed by energy minimization and rescoring of docked ligands. Finally, potential binders were inspected to ensure molecules with different scaffolds were engaged in favorable contacts with both the metal cofactor and the critical residues lining the receptor’s active site. The results of the calculations suggest that lividomycin, burixafor, quisinostat, fluprofylline, pemetrexed, spirofylline, edotecarin, and diniprofylline emerge as promising repositionable drug candidates for stabilizing the closed (substrate/inhibitor-bound) conformation of ACE2, thereby shifting the relative positions of the receptor’s critical exterior residues recognized by SARS-CoV-2. This study is among the rare ones in the relevant scientific literature to search for potential ACE2 inhibitors. In practical terms, the drugs, unmodified as they are, may be introduced into the therapeutic armamentarium of the ongoing fight against COVID-19 now, or their scaffolds may serve as rich skeletons for designing novel ACE2 inhibitors in the near future., Graphical abstract Image 1, Highlights • A drug library of 7173 ligands was screened against ACE2 using molecular docking. • The predicted ACE2–ligand complexes were optimized via molecular mechanics. • On closer examination, eight promising repositionable drug candidates were proposed. • The drugs, unmodified as they are, may serve as entry inhibitors of SARS-CoV-2. • Their scaffolds may serve as rich skeletons for designing novel ACE2 inhibitors.

The recombinant human keratinocyte growth factor (rhKGF) is clinically applied to decrease the incidence and duration of cancer therapeutic agents. Particularly, it is extensively used for oral mucositis after chemotherapy-induced damage of different human cancers. However, the usage of rhKGF in treatment is limited owing to its short half-life, poor stability, immunogenicity, tendency to aggregate, and side effects. Therefore, there is a need to enhance the stability and to reduce immunogenicity of rhKGF for therapeutic applications. In this study, the stability, activity, and immunogenicity of rhKGF were improved using computational methods. The several mutations were generated based on sequence alignment, amino acids physic-chemical properties, and the structure simulation. The 3D structure of rhKGF and proposed mutants were predicted by Modeller v9.15 program, and then were evaluated using PROSESS, PROCHECK, and ProSA web tools. Afterwards, the effect of these mutants on rhKGF structure, stability, activity, and its interaction with fibroblast growth factor receptor2-IIb (FGFR2-IIb) was analyzed through utilizing GROMACS molecular dynamics simulations and docking tools, respectively. Also, binding free energies were calculated by the Molecular Mechanics/Poisson–Boltzmann Surface Area (MM/PBSA) method. We found that F63Y, R121K, and combine1 (K38R, F63Y, K72E, N105S) mutants lead to reduction of the number of T-cell epitopes. However, all of the selected mutants, except for R121K, could considerably increase stability and affinity of the rhKGF to FGFR2-IIb, in silico. In conclusion, this study, for the first time, offered that the combine1 and F63Y mutants could highly improve the stability and activity of rhKGF and even reduce immunogenicity without having any significant effect on the biological functions of rhKGF.

Although a multitude of theoretical studies exist on light-harvesting complex 2 (LH2), less is known about the light-harvesting complex 3 (LH3) of similar ring-like structure. In a comparative study of three system, i.e., the LH2 protein-pigment aggregate of the purple bacterium Rhodospirillum molischianum as well as for the LH2 and LH3 complexes of Rhodoblastus acidophilus the similarities and the differences in the excitonic system were analyzed. To this end, the systems have been studied in a multi-scale approach that combines molecular dynamics simulations with quantum chemistry calculations and quantum dynamics. Along the ground-state molecular dynamics trajectories, the excited energy gaps were determined in a quantum mechanics/molecular mechanics hybrid fashion. Based on the simulations, spectral densities, absorption spectra and exciton dynamics have been determined. After correcting for some shortcoming in the absorption spectra, the exciton dynamics within and between the ring systems have been determined and discussed.

As one of the first-in-class inhibitor, TH588 was found to be efficient in the suppression of MutT homolog1 (MTH1). A recent work shows that the inhibitory potency of TH588 against human MTH1 (hsMTH1) is approximately 20-fold over that of mouse MTH1 (mmMTH1) and identifies residue in position 116 in MTH1 has an important contribution to TH588 affinity. But the effect of residue Leu or Met in position 116 on the binding affinity remains unclear. In this study, molecular dynamics (MD) simulations and free energy calculations were used to elucidate the mechanism about the effect of residue 116 to the different inhibitory potency of TH588 against MTH1. The binding free energy of TH588 in M116 complexes predicated by the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) is much lower than that in L116 complexes, which is consistent with the experiment results. The analysis of the individual energy terms suggests that the non-polar interactions are important for distinguishing the binding of TH588. The MD results show that the Leu116 disrupts the interactions between Asn33 and TH588, thus induces the conformational changes of Asn33 as well as TH588. The altered interactions between TH588 and mmMTH1 change the flexibility of TH588, which could induce the remarkable conformational fluctuation of mmMTH1. The conformations of the two loops covering the binding pocket have obvious influence on the opening or closure of the active site. The more open binding site may explain the lower inhibitor potency of TH588 against mmMTH1 than hsMTH1. Our results provide mechanistic insight into the effect of different residue Leu or Met in position 116 on the binding affinity of TH588 for MTH1, which is expected to contribute to the further rational design of more potent inhibitors.

Allergen induced IgE dependent type I hypersensitivity is the main cause of the allergy, which would be a burden on medical setup in coming years. Allergens of Glycine max have been isolated, and their disease relationships are documented. Therefore, it becomes important to investigate the interaction of different allergens of Glycine max with IgE and also screen suitable therapeutics to prevent this interaction. The amino acid sequences of all allergens of Glycine max and their isoallergens have been taken, and 3D structure of allergens (Gly m 3, Gly m 4, Gly m 5, Gly m 6 and Gly m 8) and their isoallergens were generated using Modeller v9.17. The modeled structures were further validated using PSVS, ProSA, RAMPAGE, and PDBsum. HL domain of Fab region of human IgE (PDBID: 2R56 ) was generated using UCSFchimera. The HL domain was minimized by Schrodinger software using the OPLS_2005 force field. SiteMap identified epitope binding site of the minimized domain. All the predicted epitopes of different allergens were docked to the binding site of HL domain using the Patchdock server. We have also designed a peptidomimetics based inhibitor targeted at interaction interface of Gly m8 and IgE, using in-silico virtual screening, molecular mechanics, and molecular dynamics simulation studies. These studies identified BDE32166344 ((N-(1-{[1-(1-aminocyclopentanecarbonyl)-3-hydroxypyrrolidin-3-yl]methyl}piperidin-4-yl)acetamide) as a peptidomimetics based lead with binding energy of −72.77 kcal/mol. Therefore, the present study investigates the interaction between different Gly m allergens and IgE antibody and identifies peptidomimetics based lead that might be developed as a suitable therapeutics against allergy caused by allergen of Glycine max.

Survivin is the smallest member of IAP (inhibitor of apoptosis protein) family, which plays important roles in both mitosis and apoptosis. It has become an attractive drug target due to its overexpression in many human cancers. Survivin has been proven to bind to Smac/DIABLO protein that indirectly inhibits apoptosis. Meanwhile, it is the key subunit of chromosome passenger complex (CPC) which bind to the N-terminal tail of phosphorylated histone H3T3ph during mitosis. Up to now, Survivin directly targeting inhibitor has yet to merge since the difficulty of disrupting the protein-protein interactions (PPIs) between Survivin and its substrate proteins. Nevertheless, currently known binding partners of Survivin provide crucial information about conserved recognition mechanism, which can assist in the detection of some uncharted substrates and also the Survivin inhibitors. Herein, we adopted a method that using four substrates to analyze the common binding mode of Survivin. To accomplish this, conventional molecular dynamics (MD) simulations, molecular mechanics/generalized born surface area (MM-GBSA) binding free energy calculations and energy decomposition were carried out to assess the binding affinity and per-residue contributions. We found that there are two anchor sites of Survivin responsible for maintaining the binding conformation and one sub-pocket for intermolecular recognition. The results of this study synthetically describe the binding mechanism and provide valuable guidance for rational drug design of PPI inhibitor.

Aim Lactoperoxidase (LPO) is an essential protein with broad spectrum antimicrobial activity present in mammalian milk. It imparts immunity to infants against wide range of pathogenic infections. Several in vitro studies have shown inhibition of LPO activity by pharmaceutical compounds including commonly used antibiotics such as ampicillin and gentamicin, and molecules like prednisolone, norepinephrine, etc. Prescription of such drugs to lactating mothers might have adverse health effects on infants. The aim of our study was the elucidation of the structural aspects of the inhibitory mechanism of ampicillin, gentamicin, amoxicillin, prednisolone and norepinephrine on LPO. Material and methods Three dimensional structure of camel LPO (cLPO) was developed using homology modeling and used for in silico experimental studies. The Schrodinger induced fit docking along with binding affinity estimation experiments were performed. The cLPO and Ligands were prepared using Protein Preparation Wizard and Ligprep modules available in Schrodinger suite. For estimating Binding affinity Prime Molecular Mechanics with Generalized Born and Surface Area (MMGB-SA) module was used. Key results The five drug ligands formed three to five hydrogen bonding interactions with cLPO. Amino acids Arg-231, Asp-232, Ser-370, Arg-371 and Glu-374 of cLPO were crucial for these interactions. The binding affinity values for gentamicin were highest and for norepinephrine were the lowest. Significance This study concludes that the five drug molecules show potential ability to inhibit the LPO activity. Further, a very high sequence similarity of cLPO with human LPO imparts high significance to these conclusions in relation to human health especially in new born infants.

The synthesis and characterization of two hybrid N-methylated carbazole derivatives containing a thiazolyl or a thienyl ring is reported. The thiazolyl derivative has been also characterised by X-ray diffraction analysis. The study of its reactivity in front of [MCl2(dmso)2] (M = Pd or Pt) or Na2[PdCl4] in methanol has allowed us to isolate and characterize its complexes. However, for the thienyl analogue, the formation of any Pd(II) or Pt(II) complex was not detected, indicating that it is less prone to bind to the M(II) ions than its thiazolyl analogue. Density Functional Theory (DFT) and Time-Dependent Density Functional Theory (TD-DFT) calculations have also been carried out in order to rationalize the influence of the nature of the thiazolyl or thienyl group on the electronic delocalization. Molecular mechanics calculations show that the free rotation of the thiazolyl in relation to the carbazole requires a greater energy income than for its thienyl analogue. Studies of the cytotoxic activity of the new compounds on colon (HCT116) and breast (MDA-MB231 and MCF7) cancer cell lines show that the thiazolyl carbazole ligand and its Pt(II) complex are the most active agents of the series and in the MCF7 line their potency is higher than that of cisplatin. In the non-tumoral human skin fibroblast BJ cell line, all the compounds were less toxic than cisplatin. Their potential ability to modify the electrophoretic mobility of pBluescript SK+ plasmid DNA and to act as inhibitors of Topoisomerases I and IIα or cathepsin B has also been investigated.

There are over 150 currently known, highly diverse chemically modified RNAs, which are dynamic, reversible, and can modulate RNA-protein interactions. Yet, little is known about the wealth of such interactions. This can be attributed to the lack of tools that allow the rapid study of all the potential RNA modifications that might mediate RNA-protein interactions. As a promising step toward this direction, here we present a computational protocol for the characterization of interactions between proteins and RNA containing post-transcriptional modifications. Given an RNA-protein complex structure, potential RNA modified ribonucleoside positions, and molecular mechanics parameters for capturing energetics of RNA modifications, our protocol operates in two stages. In the first stage, a decision-making tool, comprising short simulations and interaction energy calculations, performs a fast and efficient search in a high-throughput fashion, through a list of different types of RNA modifications categorized into trees according to their structural and physicochemical properties, and selects a subset of RNA modifications prone to interact with the target protein. In the second stage, RNA modifications that are selected as recognized by the protein are examined in-detail using all-atom simulations and free energy calculations. We implement and experimentally validate this protocol in a test case involving the study of RNA modifications in complex with Escherichia coli (E. coli) protein Polynucleotide Phosphorylase (PNPase), depicting the favorable interaction between 8-oxo-7,8-dihydroguanosine (8-oxoG) RNA modification and PNPase. Further advancement of the protocol can broaden our understanding of protein interactions with all known RNA modifications in several systems.

The ubiquitin-proteasome system plays an important role in protein quality control. Currently, inhibition of the proteasome has been validated as a promising approach in anticancer therapy. The 20S core particle of the proteasome harbors β5 subunit which is a crucial active site in proteolysis. Targeting proteasome β5 subunit which is responsible for the chymotrypsin-like activity of small molecules has been regarded as an important way for achieving therapeutics target. In the present study, a series of porphyrin derivatives bearing either pyridine or pyrazole rings as meso-substituents were designed and evaluated as an inhibitor for the β5 subunit of the proteasome by employing molecular docking and dynamics simulations. The molecular docking was performed with the help of AutoDock 4.2, while molecular dynamics simulation was done using AMBER 14. All compounds bound to the proteasome with similar binding modes, and each porphyrin-proteasome complex was stable during 30 ns MD simulation as indicated by root-mean-square-deviation (RMSD) value. An analysis on protein residue fluctuation of porphyrin binding demonstrates that in all complexes, porphyrin binding produces minor fluctuation on amino acid residues. The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free energy calculation shows that the binding affinities of mono-H2PyP, bis-H2PzP, and tetra-H2PyP were comparable with that of the potential inhibitor, HU10. It is noted that the electrostatic interaction increases with the number of meso-substituents, which was favourable for porphyrin binding. The present study shows that both electrostatic and van der Waals interaction are the main force which controls the interaction of porphyrin compounds with the proteasome.

Computational investigation of orotate phosphoribosyltransferase (OPRT) action, an enzymatic reaction between phosphoribosyl pyrophosphate (PRPP) and orotic acid (OA) to yield orotidine 5′-monophosphate (OMP), was carried out. Insights into the pathways of the substrate attack step of the reaction were developed under the quantum mechanics/molecular mechanics framework with S. cerevisiae strain as the representative enzyme bearer. Four pathways were proposed for PRPP and OA binding differing in the sequence of PRPP, OA and Mg2+ ion complexation with OPRT. The formation of Mg2+-OPRT complex was accompanied by a small energy change while the largest stabilization was observed for the formation of Mg2+-PRPP complex supporting the experimental observation of Mg2+-PRPP complex as the true substrate for the reaction. Formation of PRPP-OPRT complex was found to be energetically not probable rendering the pathway requiring Mg2+-OA complex not probable. Further, PRPP migration towards the active site was found to be energetically not favoured rendering the pathway involving Mg2+-OA complexation improbable. Migration of OA and Mg2+-PRPP complex towards the active site was found to be energetically probable with a large stabilization of the system when Mg2+-PRPP complex bound to the OA-OPRT complex. This conclusively proved the sequential binding of OA and Mg2+-PRPP complexes during OPRT action.

The linear and nonlinear properties of isomeric forms of pyridinium-N-phenoxide betaine dye were investigated in protic and aprotic solvents using atomistic simulations. We employed the sequential Quantum Mechanics/Molecular Mechanics (S-QM/MM) and the free energy gradient (FEG) methods to optimize the geometry of each isomer in chloroform, acetonitrile, methanol and water. The results show a complex dependence of the first hyperpolarizability with respect to the solvent nature and isomeric form, with a marked effect of conformational changes for para-betaine. Large contrasts of the first hyperpolarizability show a clear distinction between isomeric forms in solution that could be experimentally detected.

Here we propose a 3D-molecular structural model for dextran sulphate sodium (DSS) in a neutral aqueous environment based on the results of a molecular modelling study. The DSS structure is dominated by the stereochemistry of the 1,6-linked α-glucose units and the presence of two sulphate groups on each α-glucose unit. The structure of DSS can be best described as a helix with various patterns of di-sulphate substitution on the glucose rings. The presence of a side chain does not alter the 3D-structure of the linear main chain much, but affects the overall spatial dimension of the polymer. The simulated polymers have a diameter similar to or in some cases even larger than model α-hemolysin nano-pores for macromolecule transport in many biological processes, indicating a size-limited translocation through such pores. All results of the molecular modelling study are in line with previously reported experimental data. This study establishes the three-dimensional structure of DSS and summarizes the spatial dimension of the polymer, serving as the basis for a better understanding on the molecular level of DSS-involved electrostatic interaction processes with biological components like proteins and cell pores.

The small GTPase Ras transmits signals in a variety of cellular signaling pathways, most prominently in cell proliferation. GTP hydrolysis in the active center of Ras acts as a prototype for many GTPases and is the key to the understanding of several diseases, including cancer. Therefore, Ras has been the focus of intense research over the last decades. A recent neutron diffraction crystal structure of Ras indicated a protonated γ-guanylyl imidodiphosphate (γ-GppNHp) group, which has put the protonation state of GTP in question. A possible protonation of GTP was not considered in previously published mechanistic studies. To determine the detailed prehydrolysis state of Ras, we calculated infrared and NMR spectra from quantum mechanics/molecular mechanics (QM/MM) simulations and compared them with those from previous studies. Furthermore, we measured infrared spectra of GTP and several GTP analogs bound to lipidated Ras on a membrane system under near-native conditions. Our findings unify results from previous studies and indicate a structural model confirming the hypothesis that γ-GTP is fully deprotonated in the prehydrolysis state of Ras.

The GA codon box incorporates the two-fold degeneracy of aspartic acid and of glutamic acid. Using the molecular mechanics approach of the AMBER suite, the four codons of the GA box are paired via H-bonding with two aspartic acid anticodons and two glutamic acid anticodons to yield 8 cognate and 11 non-cognate codon-anticodon duplexes. In addition four select non-cognate duplexes between the GA box codons and three alanine anticodons are also studied. These 23 duplexes display a variety of base-pairing possibilities at the wobble position. Cognate duplexes are differentiated from non-cognate duplexes on the grounds of structure and stability (chiefly the former). The results are in line with Crick's wobble hypothesis, and corroborate the observed reading properties of the aspartic acid anticodons GUC and QUC and of the glutamic acid anticodons CUC and SmnUC.

The relative performance of MM-PBSA and MM-3D-RISM methods to estimate the binding free energy of protein–ligand complexes is investigated by applying these to three proteins (Dihydrofolate Reductase, Catechol-O-methyltransferase, and Stromelysin-1) differing in the number of metal ions they contain. None of the computational methods could distinguish all the ligands based on their calculated binding free energies (as compared to experimental values). The difference between the two comes from both polar and non-polar part of solvation. For charged ligand case, MM-PBSA and MM-3D-RISM give a qualitatively different result for the polar part of solvation.

The CBP (CREB (cAMP responsive element binding protein) binding protein) bromodomain (BRD) could recognize and bind with acetyl K382 of human tumor suppressor protein p53 which the mutation of encoding gene might cause human cancers. CBP-BRD serves as a promising drug target for several disease pathways and a series of effective drug have been discovered. In this study, molecular dynamics (MD) simulations and molecular mechanics generalized born surface area (MM-GB/SA) approaches were performed to investigate the different binding modes between five inhibitors with CBP-BRD. Based on the energy and conformation analyses, a potent core fragment is chosen to act as the starting point for new inhibitor design by means of LUDI and rational drug design approaches. Then, T.E.S.T and molinspirition were applied to evaluate oral bioavailability and drug promiscuity of the new molecules. These results shed light on the idea for further inhibitor design.

Acinetobacter baumannii, an ESKAPE pathogen, causes various nosocomial infections and has capacity to produce biofilm. Biofilm produced by this bacterium is highly tolerant to environmental factors and different antibiotics. Biofilm-associated protein (Bap) plays a significant role in the biofilm formation by A. baumannii and found in the extra cellular matrix of the biofilm. Therefore, it becomes essential to find a potential drug against Bap that has capacity to inhibit biofilm formation by A. baumannii. In-silico screening, molecular mechanics and molecular dynamics studies identified ZINC00039089 (L-Adrenaline) as an inhibitor for Bap of A. baumannii. Recently, it is reported that Bap can form amyloid like structure; hence we have created dimer of Bap protein. This inhibitor can bind to dimeric Bap with good affinity. It confirms that ZINC00039089 (L-Adrenaline) can bind with Bap monomer as well as oligomeric Bap, responsible for amyloid formation and biofilm formation. Hence, we have tested Adrenaline as an anti-biofilm molecule and determined its IC50 value against biofilm. The result showed Adrenaline has anti-biofilm activity with IC50 value of 75μg/ml. Therefore; our finding suggests that L-Adrenaline can be developed to inhibit biofilm formation by carbapenem resistant strain of Acinetobacter baumannii.

The steered molecular dynamics (SMD), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and free energy perturbation (FEP) methods were used to determine the binding affinity of the L-742,001 inhibitor to the endonuclease domain of the A/H1N1/PA influenza viruses (including wild type (WT) and three mutations I79L, E119D and F105S for both pH1N1 PA and PR8 PA viruses). Calculated results showed that the L-742,001 inhibitor not only binds to the PR8 PAs (1934 A influenza virus) better than to the pH1N1 PAs (2009 A influenza virus) but also more strongly interacts with the WT endonuclease domain than with three mutant variants for both pH1N1 PA and PR8 PA viruses. The binding affinities obtained by the SMD, MM-PBSA and FEP methods attain high correlation with available experimental data. Here the FEP method appears to provide a more accurate determination of the binding affinity than the SMD and MM-PBSA counterparts.

Filamentous temperature-sensitive Z (FtsZ) protein plays essential role in bacteria cell division, and its inhibition prevents Mycobacteria reproduction. Here we adopted curcumin derivatives as candidates of novel inhibitors and investigated their specific interactions with FtsZ, using ab initio molecular simulations based on protein–ligand docking, classical molecular mechanics and ab initio fragment molecular orbital (FMO) calculations. Based on FMO calculations, we specified the most preferable site of curcumin binding to FtsZ and highlighted the key amino acid residues for curcumin binding at an electronic level. The result will be useful for proposing novel inhibitors against FtsZ based on curcumin derivatives.

Quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics simulations were performed to study the hydrogen bonding and polar interactions of an inhibitor namely 2-{[(5- bromothiophen-2-yl) sulfonyl] amino}- 4-chlorobenzoic acid with HCV NS5B polymerase. Using the initial poses, our aim is to identify the minimum energy positions and calculate binding affinities of inhibitor to HCV. Results suggest that the accuracy of binding affinity prediction can be significantly affected by the proper treatment of the hydrogen bonding and polar interactions. Using this protocol, we successfully identified important poses of low RMSD with respect to crystallographic structure from the initially docked poses. This work shows that the most important energetic component contributing to binding for this particular protein-ligand system is the conformational energy. In addition, formation of hydrogen bonds can produce favorable van der Waals interactions with some important residues. This systematic study can provide guidance for the choice of QM/MM MD methods and design of new inhibitors targeting HCV.