Binding interaction of laccases from Bacillus Subtilis after industrial dyes exposure: Molecular docking and molecular dynamics simulation studies

Dyes are becoming more widely used around the world wide, but there is no effective bioremediation approach for removing them completely from the environment. Several dyes are mentioned to be degraded through bacteria; however, it's still unknown how the particular enzymes act throughout the dye degradation. The behavior and function of these enzymes in the biodegradation of azo dyes (Textile dyes) had been investigated experimentally by the numbers of the researchers, however, the molecular mechanisms remain unclear. Therefore, the interaction mechanisms of textile dye (methyl orange) with laccase from B. subtilis were explored through molecular docking and molecular dynamics simulations, the three selected dyes (methyl orange, malachite green, and acid blue 62) that interact positively with laccase on the basis of their maximum binding energy, molecular docking results indicate that one of the three dyes is more stable as a target for degradation through Bacillus subtilis laccase. Therefore, subsequent research focused solely on one substrate: methyl orange. Molecular Dynamics simulation study was applied after the molecular docking to determine the interaction between laccases and methyl orange dyes. The trajectory was proved with root mean square deviation and root mean square fluctuation analysis. According to the molecular dynamics simulation results, laccase-methyl orange complexes remain stable during the catalytic reaction. So, this study demonstrates how laccase is involved in methyl orange bioremediation.

Cryptosporidiosis, a public health challenge: A combined 3D shape-based virtual screening, docking study, and molecular dynamics simulation approach to identify inhibitors with novel scaffolds for the treatment of cryptosporidiosis

Cryptosporidiosis is a neglected tropical disease caused by the protozoan parasite Cryptosporidium parvum. Limited therapeutic options, limitation in in vitro parasite culture, and lack of a reliable animal model of parasite for replication of in vivo life cycle and drug testing demand alternative methods for drug development. The in silico methods of drug discovery prove a crucial process in such conditions.Recent research reported a limited number of small molecules for drug development. Purine nucleotide biosynthesis in Cryptosporidium species is dependent on the IMPDH (CpIMPDH) enzyme, so distortion of parasite IMPDH has been pursued as a compelling strategy for curbing Cryptosporidium infection due to its different kinetics from the host enzyme. Our study's primary aim was to discover novel ligand molecules with noticeable activity against Cryptosporidium parvum IMPDH. For this purpose, we selected 18 previously discovered ligands to understand the interaction feature between ligand and receptor, and their shape and electronic features are employed as a template for shape-based virtual screening of the ZINC database (drug-like subset) search approach via Schrodinger-2019 (Maestro 11.9). The obtained hits were subsequently subjected to structure-based screening, quantum polarized ligand docking (QPLD), and molecular dynamics simulations to fetch potential small molecules with the highest binding affinity for CpIMPDH protein. Further ligand binding energy and pharmacokinetic analysis were also taken into consideration as filtering criteria for selecting the most promising drug-like compounds. On this experimentation analysis, three top-ranked (ZINC24855054, ZINC58171263, and ZINC08000072) molecules were found to have appropriate pharmacokinetic properties along with surpassing in silico inhibitory potential towards the CpIMPDH compared to known inhibitors. The molecular docking and molecular dynamics simulation analysis results satisfactorily confirmed the inhibitory action. Therefore, these new scaffolds deduced by the presented computational methodology could recommend lead molecules for designing promising anti-cryptosporidial drugs targeting CpIMPDH protein.

Computational approaches towards the discovery and optimisation of cruzain inhibitors

The need to develop safer and more efficacious drugs to treat Chagas disease has motivated the search for cruzain inhibitors. Cruzain is the recombinant, truncated version of cruzipain, a cysteine protease from Trypanosoma cruzi with important roles during the parasite life cycle. Several computational techniques have been applied to discover and optimise cruzain inhibitors, providing a molecular basis to guide this process. Here, we review some of the most recent computational studies that provided important information for the design of cruzain inhibitors. Moreover, we highlight the diversity of applications of in silico techniques and their impact.

Computation-aided design of the flexible region of zearalenone hydrolase improves its thermal stability / 生物工程学报

The zearalenone hydrolase (ZHD101) derived from Clonostachys rosea can effectively degrade the mycotoxin zearalenone (ZEN) present in grain by-products and feed. However, the low thermal stability of ZHD101 hampers its applications. High throughput screening of variants using spectrophotometer is challenging because the reaction of hydrolyzing ZEN does not change absorbance. In this study, we used ZHD101 as a model enzyme to perform computation-aided design followed by experimental verification. By comparing the molecular dynamics simulation trajectories of ZHD101 at different temperatures, 32 flexible sites were selected. 608 saturated mutations were introduced into the 32 flexible sites virtually, from which 12 virtual mutants were screened according to the position specific score and enzyme conformation free energy calculation. Three of the mutants N156F, S194T and T259F showed an increase in thermal melting temperature (ΔTm>4 °C), and their enzyme activities were similar to or even higher than that of the wild type (relative enzyme activity 95.8%, 131.6% and 169.0%, respectively). Molecular dynamics simulation analysis showed that the possible mechanisms leading to the improved thermal stability were NH-π force, salt bridge rearrangement, and hole filling on the molecular surface. The three mutants were combined iteratively, and the combination of N156F/S194T showed the highest thermal stability (ΔTm=6.7 °C). This work demonstrated the feasibility of engineering the flexible region to improve enzyme performance by combining virtual computational mutations with experimental verification.

Analyzing The Role Of Phytochemicals In Targeting Drug Transporter Protein Abcc6 Using Molecular Docking And Molecular Dynamics Simulations

Triple negative breast cancer (TNBC) is an aggressive breast cancer subtype that lacks hormonal receptors. This reduces the therapeutic options for TNBC patients creating more focus on chemotherapy. Drug resistance has posed as a major hurdle in treating TNBC patients. Deregulation of drug transporter proteins is one of major factors that cause resistance to chemotherapeutic drugs. In this study, ABCC6 a drug transporter protein that is found dysregulated in several resistant cancer cells has been docked with natural compounds or phytochemicals with known anti-cancer activities. Subtrifloralactone G, a withanolide extracted from Deprea subtriflora is found to show highest binding energy with ABCC6 protein. Molecular dynamics simulations further prove the stability of the ABCC6 protein- Subtrifloralactone G ligand complex. ADMET analysis shows that phytochemical Subtrifloralactone G can be used as an anti-cancer therapeutic drug in treating resistant cancer cells. The study mainly focuses on the role of phytochemicals in treating resistant TNBC cells.

An in silico investigation of phytochemicals as potential inhibitors against non-structural protein 1 from dengue virus 4

Dengue fever has emerged as a big threat to human health since the last decade owing to high morbidity with considerable mortalities. The proposed study aims at the in silico investigation of the inhibitory action against DENV4-NS1 of phytochemicals from two local medicinal plants of Pakistan. Non-Structural Protein 1 of Dengue Virus 4 (DENV4-NS1) is known to be involved in the replication and maturation of viron in the host cells. A total of 129 phytochemicals (50 from Tanacetum parthenium and 79 from Silybum marianum) were selected for this study. The tertiary structure of DENV4-NS1 was predicted based on homology modelling using Modeller 9.18 and the structural stability was evaluated using molecular dynamics simulations. Absorption, distribution, metabolism, excretion and toxicity (ADMET) along with the drug-likeness was also predicted for these phytochemicals using SwissADME and PreADMET servers. The results of ADMET and drug-likeness predictions exhibited that 54 phytochemicals i.e. 25 from Tanacetum parthenium and 29 from Silybum marianum showed effective druglikeness. These phytochemicals were docked against DENV4-NS1 using AutoDock Vina and 18 most suitable phytochemicals with binding affinities ≤ -6.0 kcal/mol were selected as potential inhibitors for DENV4-NS1. Proposed study also exploits the novel inhibitory action of Jaceidin, Centaureidin, Artecanin, Secotanaparthenolide, Artematin, Schizolaenone B, Isopomiferin, 6, 8-Diprenyleriodictyol, and Anthraxin against dengue virus. It is concluded that the screened 18 phytochemicals have strong inhibition potential against Dengue Virus 4.

Aplicação de modelagem molecular e de formalismo do CAMD (Computer-Aided Molecular Design) na elucidação do mecanismo de ação de inibidores de metalopropteinases de matriz / Molecular modeling methods and computer-aided molecular design (CAMD) formalisms for elucidating the mechanism of action of matrix metalloproteinases inhibitors

As metaloproteinases de matriz (MMP) são enzimas superexpressas em quase todos os tumores humanos, sendo que os subtipos MMP-2 e MMP-9 têm sido associados ao potencial metastático e prognóstico desfavorável em neoplasias malignas como, por exemplo, melanoma metastático e glioma. Compostos capazes de inibir a atividade destas enzimas podem representar potenciais agentes terapêuticos. O composto 4-nerolidilcatecol (4-NC), isolado de plantas do gênero Pothomorphe, apresentou resultados promissores para o tratamento do melanoma e glioma e foi capaz de atuar em várias etapas bioquímicas importantes envolvidas na progressão dessas patologias, inclusive inibindo MMP-2 e MMP-9. No entanto, o mecanismo de ação do 4-NC não está completamente elucidado. O presente estudo envolveu a aplicação de métodos de modelagem molecular e de formalismos do planejamento de novas moléculas auxiliado por computador, CAMD (Computer-Aided Molecular Design) a fim de explorar a interação entre esta molécula e as enzimas MMP-2 e MMP-9, além de planejar novos inibidores para estes alvos. Análise exploratória de dados, que compreende a análise de agrupamentos hierárquicos e de componentes principais. foi desenvolvida para um conjunto de hidroxamatos (N=64) descritos como inibidores de MMP-2 e MMP-9, a fim de identificar as propriedades moleculares que mais influenciavam o processo de discriminação dos compostos. As propriedades termodinâmicas, eletrônicas e estéricas foram importantes para descrever os compostos mais ativos no conjunto de dados da MMP-2. Para a MMP-9, o coeficiente de distribuição (ClogD) em pH 1,5 foi relevante no processo de discriminação do conjunto. A presença de substituintes volumosos na porção R3 parece ser crucial para o conjunto de inibidores investigados. Esta região está envolvida em interações moleculares com a cavidade S1 de ambas as enzimas, mas há um limite de volume a ser considerado para estes substituintes. O formalismo QSAR-4D independente do receptor (IR) foi aplicado ao mesmo conjunto de dados e permitiu estabelecer o mapeamento do farmacóforo, além de explorar diferentes alinhamentos para a obtenção da hipótese de conformação bioativa prevista pelo melhor modelo de QSAR. OS modelos QSAR apresentaram boa capacidade de previsão, auxiliaram na proposição de novos inibidores e estimaram a atividade do 4-NC. Com o melhor modelo QSAR para MMP-9 (N=64), a atividade prevista para o 4-NC foi classificada na faixa dos inibidores com atividade moderada. Entretanto, o melhor modelo QSAR obtido para MMP-2 (N=38) não foi capaz de prever, de forma adequada, a atividade de compostos com arcabouço químico diferente daqueles utilizados na construção dos modelos. Estudos de ancoramento molecular foram desenvolvidos para investigar a orientação do 4-NC no sitio catalítico das duas enzimas e as interações que poderiam ser estabelecidas nestes complexos. Duas conformações favoráveis foram encontradas. Simulações computacionais de dinâmica molecular foram desenvolvidas com os complexos mais promissores selecionados nos estudos de ancoramento, a fim de obter informações mais detalhadas e de maior confiabilidade. sobre suas interações intermoleculares. O 4-NC tende a se orientar no sítio de forma a acomodar sua cadeia lateral no bolso S1 adjacente ao sítio catalítico em ambas as enzimas. Ensaios de zimografia também foram realizados com o objetivo de elucidar possíveis contribuições da cadeia lateral e do núcleo catecólico do 4-NC na atividade inibitória frente às enzimas em estudo. O núcleo catecólico parece ser o responsável por sua atividade, pois o composto 1,2dimetoxibenzeno, que possui as hidroxilas bloqueadas por grupos metil, não foi capaz de exercer atividade inibitória significante frente à MMP-2 e MMP-9. Estudos de voltametria reforçaram a hipótese de que o 4-NC tem a capacidade de quelar os íons zinco presentes no tampão de incubação

Matrix metalloproteinases (MMP) enzymes are overexpressed in almost all human tumors, and MMP-2 and MMP-9 subtypes have been associated with metastatic potential and poor prognosis in malignant tumors, such as metastatic melanoma and glioma. Compounds capable of inhibiting the activity of theses enzymes would be considered as potential therapeutic agents. The 4-nerolidylcatechol compound (4-NC), isolated from plants of genus Pothomorphe, has showed promising results in the treatment of melanoma and glioma, and was able to act in several important biochemical steps involved in the progression of these diseases, as well as inhibiting MMP-2 and MMP-9. However, the 4-NC mechanism of action is not completely understood. This study has involved the application of molecular modeling methods and formalisms of computer-aided molecular design (CAMD) in order to explore the interaction between 4-NC and MMP-2/MMP-9, and to design new inhibitors for these targets. Exploratory data analysis, which comprises hierarchical cluster analysis and principal components analysis, was performed to a set of hydroxamates (N=64). previously reported as MMP-2 and MMP-9 inhibitors, in order lo identify the molecular properties that is most critical for the discrimination process regarding the investigated compounds. The thermodynamic, electronic, and steric properties were: quite important to describe the highly active compounds in the data set of MMP-2, whereas the apparent partition coefficient (ClogD) at pH 1.5 was the property more relevant for MMP-9 data set. The presence of bulky substituents on the R3 moiety seems to be crucial for this set of inhibitors due to the molecular interaction with the S1 subsite of both enzymes. However, there is a limit regarding the substituents volume in this region. Receptor independent (RI) 4D-QSAR analysis was applied lo the same data set and it was possible to establish the pharmacophore mapping, besides to explore different alignments in order to generate the hypothesized bioactive conformation through the best QSAR model. The QSAR models have presented good predictability, assisted in proposing new inhibitors, and estimated the activity of 4-NC. Regarding the best QSAR model for MMP-9 (N=64), the 4-NC predicted activity was classified in the range of the moderate active inhibitors. The best QSAR model obtained for MMP-2 (N=38), however was not able to properly predict the activity for compounds with different chemical scaffold from those used to build up the QSAR model. Molecular docking studies have been developed to investigate the 4-NC binding mode into the catalytic site of the two enzymes and the interactions that could be established in those complexes. The results have shown two favorable conformers regarding the MMP inhibition. Molecular dynamics computational simulation were combined to molecular docking studies in order to obtain more detailed and reliable information regarding the intermolecular interactions of each complex. The 4-NC molecule tends to accommodate the side chain in the S1 pocket adjacent to the catalytic site in both enzymes. Experimental zymography assays were also performed to elucidate the possible contribution of the side chain and the catechol core in the 4-NC inhibitory activity against the MMP-2 and MMP-9 enzymes. The catechol core seems to be responsible for its activity, since the 1,2 dimethoxybenzene compound, which has the hydroxyl blocked by a methyl group, was not able to exert any significant inhibition on enzymes. Voltametric assays confirmed the hypothesis that 4-NC chelates zinc ions present in the incubation buffer

Conformational behavior of phenylglycines and hydroxyphenylglycines and non-planarity of phenyl rings.

The non-proteinogenic amino acids — phenylglycine (PG) and hydroxyphenylglycine (HPG) are crucial components of certain peptidic natural products and are important for the preparation of various medicines. In this, study, the conformation of model dipeptides Ac-X-NHMe of PG, p-HPG and 3, 5-di-hydroxyphenylglycine (3, 5-DHPG) was studied both in R and S form by quantum mechanical (QM) and molecular dynamics approaches. On the energy scale, the conformational states of these molecules in both the R and S were found to be degenerate by QM studies, stabilized by non-covalent interactions like carbonyl--carbonyl interactions, carbonyl-lp··π (aromatic ring) interactions etc. These interactions disappeared/weakened due to interaction of water molecules with carbonyl groups of backbone in simulation and water was found to interact with the aromatic ring through Ow-H··π or Owlp··π interactions. The degeneracy of conformational states was lifted in favor of R-form of PG and DHPG and water molecules interactions with aromatic ring led to non-planarity of the aromatic ring. In simulation studies, irrespective of the starting geometry, the , values for the R form correspond to inverse b/inverse collagen region and for the S-form, the , values correspond to b/collagen region i.e., adopt single conformation. The obtained results were in conformity with the CD spectroscopic data on D-PG and D-p-HPG. The conformational behavior of the unusual amino acids might be of great help in designing of bioactive peptides/peptide based drugs to be realized in single conformation – an essential requirement.

Explicit and implicit modeling of nanobubbles in hydrophobic confinement

Water at normal conditions is a fluid thermodynamically close to the liquid-vapor phase coexistence and features a large surface tension. This combination can lead to interesting capillary phenomena on microscopic scales. Explicit water molecular dynamics (MD) computer simulations of hydrophobic solutes, for instance, give evidence of capillary evaporation on nanometer scales, i.e., the formation of nanometer-sized vapor bubbles (nanobubbles) between confining hydrophobic surfaces. This phenomenon has been exemplified for solutes with varying complexity, e.g., paraffin plates, coarse-grained homopolymers, biological and solid-state channels, and atomistically resolved proteins. It has been argued that nanobubbles strongly impact interactions in nanofluidic devices, translocation processes, and even in protein stability, function, and folding. As large-scale MD simulations are computationally expensive, the efficient multiscale modeling of nanobubbles and the prediction of their stability poses a formidable task to the'nanophysical' community. Recently, we have presented a conceptually novel and versatile implicit solvent model, namely, the variational implicit solvent model (VISM), which is based on a geometric energy functional. As reviewed here, first solvation studies of simple hydrophobic solutes using VISM coupled with the numerical level-set scheme show promising results, and, in particular, capture nanobubble formation and its subtle competition to local energetic potentials in hydrophobic confinement.

Água em condições normais consiste de um fluido termodinamicamente próximo à fase líquida-vapor exibindo alta tensão superficial. Esta combinação conduz a fenômenos capilares interessantes na escala microscópica. Simulações computacionais baseadas em técnicas de Dinâmica Molecular em solutos hidrofóbicos por exemplo fornecem evidências do fenômeno de evaporação capilar em escalas nanométricas dando origem à formação de bolhas nanométricas confinadas entre superfícies hidrofóbicas. Este fenômeno tem sido exemplificado em solutos de complexidade variável, ex placas de parafina, polímeros, canais biológicos e de estado sólido bem como proteínas tratadas atomicamente. Tem sido arguido quebolhas nanométricas possuem impacto significativo nos mecanismos fluídicos nanométricos e na estabilidade e dobramento de proteínas. Dado que simulações baseadas em dinâmica molecular são custosas computacionalmente, o desenvolvimento de uma modelagem multiescala eficiente sob o ponto de vista computacional impõe uma tarefa formidável à comunidadenano-física. Recentemente apresentamos um novo e versátil modelo de solvente denominado modelo implícito variacional de solvente (VISM) o qual é baseado em um funcional deenergia geométrica. Tal como apresentado aqui os primeiros estudos de solvatação de solutos hidrofóbicos simples usando VISM acoplados com esquemas numéricos de conjunto de níveis mostraram resultados promissores e em particular capturaram a formação de nano-bolhas e a sua competição com potenciais localmente energéticos em condições de confinamento hidrofóbico.