Molecular dynamics simulations on the critical states of the farnesyltransferase enzyme

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

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

Mechanistic insight into impact of phosphorylation on the enzymatic steps of farnesyltransferase

Farnesyltransferase (FTase) is a heterodimeric enzyme, which catalyzes covalent attachment of the farnesyl group to target proteins, thus coordinating their trafficking in the cell. FTase has been demonstrated to be highly expressed in cancer and neurological diseases; hence considered as a hot target for therapeutic purposes. However, due to the nonspecific inhibition, there has been only one inhibitor that could be translated into the clinic. Importantly, it has been shown that phosphorylation of the α-subunit of FTase increases the activity of the enzyme in certain diseases. As such, understanding the impact of phosphorylation on dynamics of FTase provides a basis for targeting a specific state of the enzyme that emerges under pathological conditions. To this end, we performed 18 μs molecular dynamics (MD) simulations using complexes of (non)-phosphorylated FTase that are representatives of the farnesylation reaction. We demonstrated that phosphorylation modulated the catalytic site by rearranging interactions between farnesyl pyrophosphate (FPP)/peptide substrate, catalytic Zn2+ ion/coordinating residues and hot-spot residues at the interface of the subunits, all of which led to the stabilization of the substrate and facilitation of the release of the product, thus collectively expediting the reaction rate. Importantly, we also identified a likely allosteric pocket on the phosphorylated FTase, which might be used for specific targeting of the enzyme. To the best of our knowledge, this is the first study that systematically examines the impact of phosphorylation on the enzymatic reaction steps, hence opens up new avenues for drug discovery studies that focus on targeting phosphorylated FTase.

Mutational scanning of spike RBD protein for enhanced ACE2 affinity emerging Southeast Asia in the late transmission phase

The COVID-19 pandemic has changed the quality of life and economic systems all over the world, as the virus can be transmitted from human to human via air-droplets. Since the SARS-CoV-2 virus was first identified in 2019, the virus has naturally mutated over time. Southeast Asia is one of the areas in the world that has implemented various procedures and measures to slow down the disease outbreaks. The first cluster of COVID-19 was identified from the tourist-travel history, and then the diversity of coronavirus victims has posed a serious issue of human security on a massive scale. To evaluate whether or not naturally occurring mutations have strengthened the infectivity of SARS-CoV-2, we computed in silico the structural dynamics of the RBD-spike protein mutation enhancing ACE2-binding. When considering emerging variations in Southeast Asia, 14 dominant mutations were analyzed by applying the structural and energetic characterization using MD simulations. The ones in the RBD region displayed higher affinity to ACE2 due to the improved interfacial stability of the RBD β-strand surrounding the ACE2 across salt bridge hotspots. The binding hotspots and structurally conserved conformational-epitopes have been identified, which are deleterious for RBD mutation and ACE2 binding. We present an interactive visualization to facilitate the development of effective neutralizing agents for vaccination, prevention and treatment.

Drug Repurposing Targeting Pseudomonas aeruginosa MvfR Using Docking, Virtual Screening, Molecular Dynamics, and Free-Energy Calculations

Pseudomonas aeruginosa is an opportunistic Gram-negative bacterium responsible for acute and chronic infections in planktonic state or in biofilms. The sessile structures are known to confer physical stability, increase virulence, and work as a protective armor against antimicrobial compounds. P. aeruginosa can control the expression of genes, population density, and biofilm formation through a process called quorum sensing (QS), a rather complex and hierarchical system of communication. A recent strategy to try and overcome bacterial resistance is to target QS proteins. In this study, a combined multi-level computational approach was applied to find possible inhibitors against P. aeruginosa QS regulator protein MvfR, also known as PqsR, using a database of approved FDA drugs, as a repurposing strategy. Fifteen compounds were identified as highly promising putative MvfR inhibitors. On those 15 MvfR ligand complexes, molecular dynamic simulations and MM/GBSA free-energy calculations were performed to confirm the docking predictions and elucidate on the mode of interaction. Ultimately, the five compounds that presented better binding free energies of association than the reference molecules (a known antagonist, M64 and a natural inducer, 2-nonyl-4-hydroxyquinoline) were highlighted as very promising MvfR inhibitors.

Synthesis, computational and nanoencapsulation studies on eugenol-derived insecticides

A new set of alkoxy alcohols were synthesised by reaction of eugenol oxirane with aliphatic and aromatic alcohols. These eugenol derivatives were evaluated against their effect upon the viability of...

Amino Alcohols from Eugenol as Potential Semisynthetic Insecticides: Chemical, Biological, and Computational Insights

A series of β-amino alcohols were prepared by the reaction of eugenol epoxide with aliphatic and aromatic amine nucleophiles. The synthesized compounds were fully characterized and evaluated as potential insecticides through the assessment of their biological activity against Sf9 insect cells, compared with a commercial synthetic pesticide (chlorpyrifos, CHPY). Three derivatives bearing a terminal benzene ring, either substituted or unsubstituted, were identified as the most potent molecules, two of them displaying higher toxicity to insect cells than CHPY. In addition, the most promising molecules were able to increase the activity of serine proteases (caspases) pivotal to apoptosis and were more toxic to insect cells than human cells. Structure-based inverted virtual screening and molecular dynamics simulations demonstrate that these molecules likely target acetylcholinesterase and/or the insect odorant-binding proteins and are able to form stable complexes with these proteins. Encapsulation assays in liposomes of DMPG and DPPC/DMPG (1:1) were performed for the most active compound, and high encapsulation efficiencies were obtained. A thermosensitive formulation was achieved with the compound release being more efficient at higher temperatures.

New carvacrol and thymol derivatives as potential insecticides: synthesis, biological activity, computational studies and nanoencapsulation

New compounds with potential insecticide activity were synthesized by structural modifications performed in the monoterpenoid phenolic moieties of carvacrol and thymol, resulting in a set of derivatives with the ether function containing the propyl, chloropropyl or hydroxypropyl chains, as well as a bicyclic ether with an unsaturated chain containing a carboxylic acid terminal. In addition, an analogue of carvacrol and thymol isomers bearing methoxyl, 1-hydroxyethyl and (3-chlorobenzoyl)oxy, instead of the three original methyl groups, was also synthesized. Several structural changes that resulted in diminished insecticide activity have been identified, but two significantly active molecules have been synthesized, one of them being less toxic to human cells than the naturally-derived starting materials. Structure-based inverted virtual screening and molecular dynamics simulations demonstrate that these active molecules likely target the insect odorant binding proteins and/or acetylcholinesterase and are able to form stable complexes. For the most promising compounds, nanoencapsulation assays were carried out in liposomes of egg phosphatidylcholine/cholesterol (7 : 3) prepared by both thin film hydration and ethanolic injection methods. The compound-loaded liposomes were generally monodisperse and with sizes smaller than or around 200 nm. The thin film hydration method allowed high encapsulation efficiencies (above 85%) for both compounds and a delayed release, while for the systems prepared by ethanolic injection the encapsulation efficiency is lower than 50%, but the release is almost complete in two days.

Dynamics and binding interactions of peptide inhibitors of dengue virus entry

In this study, we investigate the binding interactions of two synthetic antiviral peptides (DET2 and DET4) on type II dengue virus (DENV2) envelope protein domain III. These two antiviral peptides are designed based on the domain III of the DENV2 envelope protein, which has shown significant inhibition activity in previous studies and can be potentially modified further to be active against all dengue strains. Molecular docking was performed using AutoDock Vina and the best-ranked peptide-domain III complex was further explored using molecular dynamics simulations. Molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) was used to calculate the relative binding free energies and to locate the key residues of peptide-protein interactions. The predicted binding affinity correlated well with the previous experimental studies. DET4 outperformed DET2 and is oriented within the binding site through favorable vdW and electrostatic interactions. Pairwise residue decomposition analysis has revealed several key residues that contribute to the binding of these peptides. Residues in DET2 interact relatively lesser with the domain III compared to DET4. Dynamic cross-correlation analysis showed that both the DET2 and DET4 trigger different dynamic patterns on the domain III. Correlated motions were seen between the residue pairs of DET4 and the binding site while binding of DET2 results in anti-correlated motion on the binding site. This work showcases the use of computational study in elucidating and explaining the experiment observation on an atomic level.

A spectroscopic and molecular dynamics simulation approach towards the stabilizing effect of ammonium-based ionic liquids on bovine serum albumin

Experimental and theoretical evidence in support of the stabilizing effect of ammonium-based ionic liquids on thermal unfolding/refolding of bovine serum albumin is provided in this article.

Molecular dynamic simulations and structure-based pharmacophore development for farnesyltransferase inhibitors discovery

Farnesyltransferase is one of the enzyme targets for the development of drugs for diseases, including cancer, malaria, progeria, etc. In the present study, the structure-based pharmacophore models have been developed from five complex structures (1LD7, 1NI1, 2IEJ, 2ZIR and 2ZIS) obtained from the protein data bank. Initially, molecular dynamic (MD) simulations were performed for the complexes for 10 ns using AMBER 12 software. The conformers of the complexes (75) generated from the equilibrated protein were undergone protein-ligand interaction fingerprint (PLIF) analysis. The results showed that some important residues, such as LeuB96, TrpB102, TrpB106, ArgB202, TyrB300, AspB359 and TyrB361, are predominantly present in most of the complexes for interactions. These residues form side chain acceptor and surface (hydrophobic or π-π) kind of interactions with the ligands present in the complexes. The structure-based pharmacophore models were generated from the fingerprint bits obtained from PLIF analysis. The pharmacophore models have 3-4 pharmacophore contours consist of acceptor and metal ligation (Acc & ML), hydrophobic (HydA) and extended acceptor (Acc2) features with the radius ranging between 1-3 Å for Acc & ML and 1-2 Å for HydA. The excluded volumes of the pharmacophore contours radius are between 1-2 Å. Further, the distance between the interacting groups, root mean square deviation (RMSD), root mean square fluctuation (RMSF) and radial distribution function (RDF) analysis were performed for the MD-simulated proteins using PTRAJ module. The generated pharmacophore models were used to screen a set of natural compounds and database compounds to select significant HITs. We conclude that the developed pharmacophore model can be a significant model for the identification of HITs as FTase inhibitors.

Insights into the mechanistic dichotomy of the protein farnesyltransferase peptide substrates CVIM and CVLS

Protein farnesyltransferase (FTase) catalyzes farnesylation of a variety of peptide substrates. (3)H α-secondary kinetic isotope effect (α-SKIE) measurements of two peptide substrates, CVIM and CVLS, are significantly different and have been proposed to reflect a rate-limiting S(N)2-like transition state with dissociative characteristics for CVIM, while, due to the absence of an isotope effect, CVLS was proposed to have a rate-limiting peptide conformational change. Potential of mean force quantum mechanical/molecular mechanical studies coupled with umbrella sampling techniques were performed to further probe this mechanistic dichotomy. We observe the experimentally proposed transition state (TS) for CVIM but find that CVLS has a symmetric S(N)2 TS, which is also consistent with the absence of a (3)H α-SKIE. These calculations demonstrate facile substrate-dependent alterations in the transition state structure catalyzed by FTase.

In Silico–Based Structural Analysis of Arylthiophene Derivatives for FTase Inhibitory Activity, hERG, and Other Toxic Effects

In the present investigation, the authors have performed an in silico–based analysis on a series of arylthiophene derivatives for the determination of their structural features responsible for farnesyltransferase (FTase) inhibitory activity, hERG blocking activity, and toxicity by quantitative structure–activity relationship and pharmacophore analysis techniques. The statistically significant models derived through multiple linear regression analysis were validated by different validation methods. The applicability of the descriptors contributed in the selected models show that the polar and polarizable properties on the van der Waals (vdW) surface area of the molecules are important for the FTase inhibitory and hERG blocking activities, while being detrimental for the toxicity of the molecules. It is interesting to note that the topological properties, molecular flexibility, and connectivity of the molecules are positively correlated to all the activities (FTase inhibition, hERG blocking, and toxicity). This implies that the flexibility of the molecules is the common feature for interaction in all targets, whereas the presence of polar groups on the molecular surface (vdW) is a determinant for the favorable (FTase inhibition) or unwanted effect (hERG blocking and toxicity) of the molecules. The pharmacophore analysis of the molecules demonstrated that the aromatic/hydrophobicity and polarizability features are important pharmacophore contours favorable for these activities.

Finding a needle in the haystack: computational modeling of Mg2+ binding in the active site of protein farnesyltransferase

Studies aimed at elucidating the unknown Mg2+ binding site in protein farnesyltransferase (FTase) are reported. FTase catalyzes the transfer of a farnesyl group to a conserved cysteine residue (Cys1p) on a target protein, an important step for proteins in the signal transduction pathways (e.g., Ras). Mg2+ ions accelerate the protein farnesylation reaction by up to 700-fold. The exact function of Mg2+ in catalysis and the structural characteristics of its binding remain unresolved to date. Molecular dynamics (MD) simulations addressing the role of magnesium ions in FTase are presented, and relevant octahedral binding motifs for Mg2+ in wild-type (WT) FTase and the Dβ352A mutant are explored. Our simulations suggest that the addition of Mg2+ ions causes a conformational change to occur in the FTase active site, breaking interactions known to keep FPP in its inactive conformation. Two relevant Mg2+ ion binding motifs were determined in WT FTase. In the first binding motif, WT1, the Mg2+ ion is coordinated to D352β, zinc-bound D297β, two water molecules, and one oxygen atom from the α- and β-phosphates of farnesyl diphosphate (FPP). The second binding motif, WT2, is identical with the exception of the zinc-bound D297β being replaced by a water molecule in the Mg2+ coordination complex. In the Dβ352A mutant Mg2+ binding motif, D297β, three water molecules, and one oxygen atom from the α- and β-phosphates of FPP complete the octahedral coordination sphere of Mg2+. Simulations of WT FTase, in which Mg2+ was replaced by water in the active site, recreated the salt bridges and hydrogen-bonding patterns around FPP, validating these simulations. In all Mg2+ binding motifs, a key hydrogen bond was identified between a magnesium-bound water and Cys1p, bridging the two metallic binding sites and, thereby, reducing the equilibrium distance between the reacting atoms of FPP Cys1p. The free energy profiles calculated for these systems provide a qualitative understanding of experimental results. They demonstrate that the two reactive atoms approach each other more readily in the presence of Mg2+ in WT FTase and mutant. The flexible WT2 model was found to possess the lowest barrier toward the conformational change, suggesting it is the preferred Mg2+ binding motif in WT FTase. In the mutant, the absence of D352β makes the transition toward a conformational change harder. Our calculations find support for the proposal that D352β performs a critical role in Mg2+ binding and Mg2+ plays an important role in the conformational transition step.

Virtual screening and pharmacophore studies for ftase inhibitors using Indian plant anticancer compounds database

Farnesyl transferase (FTase) is an enzyme responsible for post-translational modification in proteins having a carboxy-terminal CaaX motif in human. It catalyzes the attachment of a lipid group in proteins of RAS superfamily, which is essential in signal transduction. FTase has been recognized as an important target for anti cancer therapeutics. In this work, we performed virtual screening against FTase with entire 125 compounds from Indian Plant Anticancer Database using AutoDock 3.0.5 software. All compounds were docked within binding pocket containing Lys164, Tyr300, His248 and Tyr361 residues in crystal structure of FTase. These complexes were ranked according to their docking score, using methodology that was shown to achieve maximum accuracy. Finally we got three potent compounds with the best Autodock docking Score (Vinorelbine: -21.28 Kcal/mol, Vincristine: -21.74 Kcal/mol and Vinblastine: -22.14 Kcal/mol) and their energy scores were better than the FTase bound co-crystallized ligand (L- 739: ―7.9 kcal/mol). These three compounds belong to Vinca alkaloids were analyzed through Python Molecular Viewer for their interaction studies. It predicted similar orientation and binding modes for these compounds with L-739 in FTase.Thus from the complex scoring and binding ability it is concluded that these Vinca alkaloids could be promising inhibitors for FTase. A 2-D pharmacophore was generated for these alkaloids using LigandScout to confirm it. A shared feature pharmacophore was also constructed that shows four common features (one hydogen bond Donar, Two hydrogen bond Acceptor and one ionizable area) help compounds to interact with this enzyme.

Determination of ensemble-average pairwise root mean-square deviation from experimental B-factors

Root mean-square deviation (RMSD) after roto-translational least-squares fitting is a measure of global structural similarity of macromolecules used commonly. On the other hand, experimental x-ray B-factors are used frequently to study local structural heterogeneity and dynamics in macromolecules by providing direct information about root mean-square fluctuations (RMSF) that can also be calculated from molecular dynamics simulations. We provide a mathematical derivation showing that, given a set of conservative assumptions, a root mean-square ensemble-average of an all-against-all distribution of pairwise RMSD for a single molecular species, <RMSD(2)>(1/2), is directly related to average B-factors (<B>) and <RMSF(2)>(1/2). We show this relationship and explore its limits of validity on a heterogeneous ensemble of structures taken from molecular dynamics simulations of villin headpiece generated using distributed-computing techniques and the Folding@Home cluster. Our results provide a basis for quantifying global structural diversity of macromolecules in crystals directly from x-ray experiments, and we show this on a large set of structures taken from the Protein Data Bank. In particular, we show that the ensemble-average pairwise backbone RMSD for a microscopic ensemble underlying a typical protein x-ray structure is approximately 1.1 A, under the assumption that the principal contribution to experimental B-factors is conformational variability.