The psychrophilic enzyme is an interesting subject to study due to its special ability to adapt to extreme temperatures, unlike typical enzymes. Utilizing computer-aided software, the predicted structure and function of the enzyme lipase AMS8 (LipAMS8) (isolated from the psychrophilic Pseudomonas sp., obtained from the Antarctic soil) are studied. The enzyme shows significant sequence similarities with lipases from Pseudomonas sp. MIS38 and Serratia marcescens. These similarities aid in the prediction of the 3D molecular structure of the enzyme. In this study, 12 ns MD simulation is performed at different temperatures for structural flexibility and stability analysis. The results show that the enzyme is most stable at 0°C and 5°C. In terms of stability and flexibility, the catalytic domain (N-terminus) maintained its stability more than the noncatalytic domain (C-terminus), but the non-catalytic domain showed higher flexibility than the catalytic domain. The analysis of the structure and function of LipAMS8 provides new insights into the structural adaptation of this protein at low temperatures. The information obtained could be a useful tool for low temperature industrial applications and molecular engineering purposes, in the near future.
Glycolipids form materials of considerable potential for a wide range of surfactant and thin film applications. Understanding the effect of glycolipid covalent structure on the properties of their thermotropic and lyotropic assemblies is a key step toward rational design of new glycolipid-based materials. Here, we perform molecular dynamics simulations of anhydrous bilayers of dodecyl β-maltoside, dodecyl β-cellobioside, dodecyl β-isomaltoside, and a C(12)C(10) branched β-maltoside. Specifically, we examine the consequences of chain branching and headgroup identity on the structure and dynamics of the lamellar assemblies. Chain branching of the glycolipid leads to measurable differences in the dimensions and interactions of the lamellar assembly, as well as a more fluid-like hydrophobic chain region. Substitution of the maltosyl headgroup of βMal-C(12) by an isomaltosyl moiety leads to a significant decrease in bilayer spacing as well as a markedly altered pattern of inter-headgroup hydrogen bonding. The distinctive simulated structures of the two regioisomers provide insight into the difference of ~90 °C in their observed clearing temperatures. For all four simulated glycolipid systems, with the exception of the sn-2 chain of the branched maltoside, the alkyl chains are ordered and exhibit a distinct tilt, consistent with recent crystallographic analysis of a branched chain Guerbet glycoside. These insights into structure-property relationships from simulation provide an important molecular basis for future design of synthetic glycolipid materials.
Currently, there is no three-dimensional structure of D-specific dehalogenase (DehD) in the protein database. We modeled DehD using ab initio technique, performed molecular dynamics (MD) simulation and docking of D-2-chloropropionate (D-2CP), D-2-bromopropionate (D-2BP), monochloroacetate (MCA), monobromoacetate (MBA), 2,2-dichloropropionate (2,2-DCP), d,l-2,3-dichloropropionate (d,l-2,3-DCP), and 3-chloropropionate (3-CP) into the DehD active site. The sequences of DehD and D-2-haloacid dehalogenase (HadD) from Pseudomonas putida AJ1 have 15% sequence similarity. The model had 80% of the amino acid residues in the most favored region when compared to the crystal structure of DehI from Pseudomonas putida PP3. Docking analysis revealed that Arg107, Arg134 and Tyr135 interacted with D-2CP, and Glu20 activated the water molecule for hydrolytic dehalogenation. Single residue substitutions at 25-30 °C showed that polar residues of DehD were stable when substituted with nonpolar residues and showed a decrease in activity within the same temperature range. The molecular dynamics simulation of DehD and its variants showed that in R134A variant, Arg107 interacted with D-2CP, while in Y135A, Gln221 and Arg231 interacted with D-2CP. It is our emphatic belief that the new model will be useful for the rational design of DehDs with enhanced potentials.
The molecular dynamics of a synthetic branched chain glycolipid, 2-decyl-tetradecyl-β-d-maltoside (C14-10G2), in the dry assemblage of smectic and columnar liquid crystal phases has been studied by dielectric spectroscopy as a function of frequency and temperature during the cooling process. Strong relaxation modes were observed corresponding to the tilted smectic and columnar phases, respectively. At low frequency (∼900 Hz to 1 kHz) in the smectic phase, Process I* was observed due to the tilted sugar bilayer structure. The process continued in the columnar phase (Process I) with an abrupt dynamic change due to phase transition in the frequency range of ∼1.3 kHz to 22 kHz. An additional process (Process II) was observed in the columnar phase with a broader relaxation in the frequency range of ∼10 Hz to 1 kHz. A bias field dependence study was performed in the columnar phase and we found that the relaxation strength rapidly decreased with increased applied dc bias field. This relaxation originates from a collective motion of polar groups within the columns. The results of dielectric spectroscopy were supported by a molecular dynamics simulation study to identify the origin of the relaxation processes, which could be related to the chirality and hydrogen bonds of the sugar lipid.
Fruit bromelain is a cysteine protease accumulated in pineapple fruits. This proteolytic enzyme has received high demand for industrial and therapeutic applications. In this study, fruit bromelain sequences QIM61759, QIM61760 and QIM61761 were retrieved from the National Center for Biotechnology Information (NCBI) Genbank Database. The tertiary structure of fruit bromelain QIM61759, QIM61760 and QIM61761 was generated by using MODELLER. The result revealed that the local stereochemical quality of the generated models was improved by using multiple templates during modelling process. Moreover, by comparing with the available papain model, structural analysis provides an insight on how pro-peptide functions as a scaffold in fruit bromelain folding and contributing to inactivation of mature protein. The structural analysis also disclosed the similarities and differences between these models. Lastly, thermal stability of fruit bromelain was studied. Molecular dynamics simulation of fruit bromelain structures at several selected temperatures demonstrated how fruit bromelain responds to elevation of temperature.
Cell-penetrating peptides are used in the delivery of peptides and biologics, with some cell-penetrating peptides found to be more efficient than others. The exact mechanism of how they interact with the cell membrane and penetrate it, however, remains unclear. This study attempts to investigate the difference in free energy profiles of three cell-penetrating peptides (TAT, CPP1 and CPP9) with a model lipid bilayer (DOPC) using molecular dynamics pulling simulations with umbrella sampling. Potential mean force (PMF) and free energy barrier between the peptides and DOPC are determined using WHAM analysis and MM-PBSA analysis, respectively. CPP9 is found to have the smallest PMF value, followed by CPP1 and TAT, consistent with the experimental data. YDEGE peptide, however, does not give the highest PMF value, although it is a non-cell-permeable peptide. YDEGE is also found to form water pores, alongside with TAT and CPP9, suggesting that it is difficult to distinguish true water pore formation from artefacts arising from pulling simulations. On the contrary, free energy analysis of the peptide-DOPC complex at the lipid-water interface with MM-PBSA provides results consistent with experimental data with CPP9 having the least interaction with DOPC and lowest free energy barrier, followed by CPP1, TAT and YDEGE. These findings suggest that peptide-lipid interaction at the lipid-water interface has a direct correlation with the penetration efficiency of peptides across the lipid bilayer.
Polymorphisms of the ADIPOR2 gene are frequently linked to a higher risk of developing diseases including obesity, type 2 diabetes and cardiovascular diseases. Though mutations of the ADIPOR2 gene are detrimental, there is a lack of comprehensive in silico analyses of the functional and structural impacts at the protein level. Considering the involvement of ADIPOR2 in glucose uptake and fatty acid oxidation, an in silico functional analysis was conducted to explore the possible association between genetic mutations and phenotypic variations. A genomic analysis of 82 nonsynonymous SNPs in ADIPOR2 was initiated using SIFT followed by the SNAP2, nsSNPAnalyzer, PolyPhen-2, SNPs&GO, FATHMM and PROVEAN servers. A total of 10 mutations (R126W, L160Q, L195P, F201S, L235R, L235P, L256R, Y328H, E334K and Q349H) were predicted to have deleterious effects on the ADIPOR2 protein and were therefore selected for further analysis. Theoretical models of the variants were generated by comparative modeling via MODELLER 9.16. A protein structural analysis of these amino acid variants was performed using SNPeffect, I-Mutant, ConSurf, Swiss-PDB Viewer and NetSurfP to explore their solvent accessibility, molecular dynamics and energy minimization calculations. In addition, FTSite was used to predict the ligand binding sites, while NetGlycate, NetPhos2.0, UbPerd and SUMOplot were used to predict post-translational modification sites. All of the variants showed increased free energy, though F201S exhibited the highest energy increase. The root mean square deviation values of the modeled mutants strongly indicated likely pathogenicity. Remarkably, three binding sites were detected on ADIPOR2, and two mutations at positions 328 and 201 were found in the first and second binding pockets, respectively. Interestingly, no mutations were found at the post-translational modification sites. These genetic variants can provide a better understanding of the wide range of disease susceptibility associated with ADIPOR2 and aid the development of new molecular diagnostic markers for these diseases. The findings may also facilitate the development of novel therapeutic elements for associated diseases.
Human dihydrofolate reductase (DHFR) is a conserved enzyme that is central to folate metabolism and is widely targeted in pathogenic diseases as well as cancers. Although studies have reported the fact that genetic mutations in DHFR leads to a rare autosomal recessive inborn error of folate metabolism and drug resistance, there is a lack of an extensive study on how the deleterious non-synonymous SNPs (nsSNPs) disrupt its phenotypic effects. In this study, we aim at discovering the structural and functional consequences of nsSNPs in DHFR by employing a combined computational approach consisting of ten recently developed in silico tools for identification of damaging nsSNPs and molecular dynamics (MD) simulation for getting deeper insights into the magnitudes of damaging effects. Our study revealed the presence of 12 most deleterious nsSNPs affecting the native phenotypic effects, with three (R71T, G118D, Y122D) identified in the co-factor and ligand binding active sites. MD simulations also suggested that these three SNPs particularly Y122D, alter the overall structural flexibility and dynamics of the native DHFR protein which can provide more understandings into the crucial roles of these mutants in influencing the loss of DHFR function.
Less sedimentation and convection in a microgravity environment has become a well-suited condition for growing high quality protein crystals. Thermostable T1 lipase derived from bacterium Geobacilluszalihae has been crystallized using the counter diffusion method under space and earth conditions. Preliminary study using YASARA molecular modeling structure program for both structures showed differences in number of hydrogen bond, ionic interaction, and conformation. The space-grown crystal structure contains more hydrogen bonds as compared with the earth-grown crystal structure. A molecular dynamics simulation study was used to provide insight on the fluctuations and conformational changes of both T1 lipase structures. The analysis of root mean square deviation (RMSD), radius of gyration, and root mean square fluctuation (RMSF) showed that space-grown structure is more stable than the earth-grown structure. Space-structure also showed more hydrogen bonds and ion interactions compared to the earth-grown structure. Further analysis also revealed that the space-grown structure has long-lived interactions, hence it is considered as the more stable structure. This study provides the conformational dynamics of T1 lipase crystal structure grown in space and earth condition.
Antimicrobial drug resistance (AMR) is a severe global threat to public health. The increasing emergence of drug-resistant bacteria requires the discovery of novel antibacterial agents. Quinoline derivatives have previously been reported to exhibit antimalarial, antiviral, antitumor, antiulcer, antioxidant and, most interestingly, antibacterial properties. In this study, we evaluated the binding affinity of three newly designed hydroxyquinolines derived from sulfanilamide (1), 4-amino benzoic acid (2) and sulfanilic acid (3) towards five bacterial protein targets (PDB ID: 1JIJ, 3VOB, 1ZI0, 6F86, 4CJN). The three derivatives were designed considering the amino acid residues identified at the active site of each protein involved in the binding of each co-crystallized ligand and drug-likeness properties. The ligands displayed binding energy values with the target proteins ranging from -2.17 to -8.45 kcal/mol. Compounds (1) and (3) showed the best binding scores towards 1ZI0/3VOB and 1JIJ/4CJN, respectively, which may serve as new antibiotic scaffolds. Our in silico results suggest that sulfanilamide (1) or sulfanilic acid (3) hydroxyquinoline derivatives have the potential to be developed as bacterial inhibitors, particularly MRSA inhibitors. But before that, it must go through the proper preclinical and clinical trials for further scientific validation. Further experimental studies are warranted to explore the antibacterial potential of these compounds through preclinical and clinical studies.Communicated by Ramaswamy H. Sarma.
Clinacanthus nutans is a medicinal plant recognised for its anticancer properties. We previously discovered that the C. nutans extract had the most potent inhibitory effect on MCF7 breast cancer cell and significantly induced apoptosis. However, there is a scarcity of studies demonstrating the molecular interactions of C. nutans-derived chemical compounds associated with apoptosis-related proteins. Therefore, the objective of this study was to determine the potential chemical compounds found in the C. nutans extract and examine their interactions with the targeted apoptotic proteins using molecular docking and molecular dynamic simulations. To address this objective, the compounds found in the SF2 extract of C. nutans were analysed using Gas Chromatography-Mass Spectrometry (GC-MS). The molecular interaction of the compounds with the targeted apoptotic proteins were determined using molecular docking and molecular dynamic simulations. GC-MS analysis revealed a total of 32 compounds in the SF2 extract. Molecular docking analysis showed that compound β-amyrenol had the highest binding affinity for MDM2-P53 (-7.26 kcal/mol), BCL2 (-11.14 kcal/mol), MCL1-BAX (-6.42 kcal/mol), MCL1-BID (-6.91 kcal/mol), and caspase-9 (-12.54 kcal/mol), whereas campesterol had the highest binding affinity for caspase-8 (-10.11 kcal/mol) and caspase-3 (-10.14 kcal/mol). These selected compounds were subjected to molecular dynamic simulation at 310 K for 100 ns. The results showed that the selected protein-ligand conformation complexes were stable, compact, and did not alter much when compared to the protein references. The findings indicate that β-amyrenol and campesterol are potentially significant compounds that might provide insight into the molecular interactions of the compounds with the apoptosis-related proteins.Communicated by Ramaswamy H. Sarma.
A histidine acid phosphatase (HAP) (PhySc) with 99.50% protein sequence similarity with PHO5 from Saccharomyces cerevisiae was expressed functionally with the molecular mass of ∼110 kDa through co-expression along with the set of molecular chaperones dnaK, dnaJ, GroESL. The purified HAP illustrated the optimum activity of 28.75 ± 0.39 U/mg at pH 5.5 and 40 ˚C. The Km and Kcat values towards calcium phytate were 0.608 ± 0.09 mM and 650.89 ± 3.6 s- 1. The half-lives (T1/2) at 55 and 60 ˚C were 2.75 min and 55 s, respectively. The circular dichroism (CD) demonstrated that PhySc includes 30.5, 28.1, 21.3, and 20.1% of random coils, α-Helix, β-Turns, and β-Sheet, respectively. The Tm recorded by CD for PhySc was 56.5 ± 0.34˚C. The molecular docking illustrated that His59 and Asp322 act as catalytic residues in the PhySc. MD simulation showed that PhySc at 40 ˚C has higher structural stability over those of the temperatures 60 and 80 ˚C that support the thermodynamic in vitro investigations. Secondary structure content results obtained from MD simulation indicated that PhySc consists of 34.03, 33.09, 17.5, 12.31, and 3.05% of coil, helix, turn, sheet, and helix310, respectively, which is almost consistent with the experimental results.
The infection produced by the SARS-CoV-2 virus remains a significant health crisis worldwide. The lack of specific medications for COVID-19 necessitates a concerted effort to find the much-desired therapies for this condition. The main protease (Mpro) of SARS-CoV-2 is a promising target, vital for virus replication and transcription. In this study, fifty pyrazole derivatives were tested for their pharmacokinetics and drugability, resulting in eight hit compounds. Subsequent molecular docking simulations on SARS-CoV-2 main protease afforded two lead compounds with strong affinity at the active site. Additionally, the molecular dynamics (MD) simulations of lead compounds (17 and 39), along with binding free energy calculations, were accomplished to validate the stability of the docked complexes and the binding poses achieved in docking experiments. Based on these findings, compound 17 and 39, with their favorable projected pharmacokinetics and pharmacological characteristics, are the proposed potential antiviral candidates which require further investigation to be used as anti-SARS-CoV-2 medication.
Inactivation of smoothened protein (SMO) by the antagonists in SHH-driven cancer types is essential for inhibition of cancer progression. This article presents molecular dynamics (MD) trajectories of water solution of three protein-ligand complexes smoothened-β-sitosterol (SMO-BST), smoothened-sonidegib (SMO-SNG) and smoothened-cholesterol (SMO-CLR) using CHARMM36 and SPC/E water model combination. Additionally, the work presents the topologies and trajectories of GROMACS files that were employed to analyse the protein-ligand interaction types (PyContact) and binding energy calculation (g_mmpbsa). The data demonstrated that equilibrated models of SMO-SNG and SMO-CLR complexes showed crucial residues that almost similar for interaction and contribution energy as previously reported in laboratory setup (in vitro). Initial simulations confirmed the role of ARG451 and TRP535 in the dynamic regulation of SMO. These data then were used as a reference for understanding the molecular dynamics of SMO-BST complex and thus predicted its mechanism of action.
This paper focus to examine the best molecular interaction between Polyamide Thin Film Composite (PA TFC) layers with different properties of the support membrane. The support membrane of Nylon 66 (N66) and Polyvinylidene fluoride (PVDF) was chosen to represent the hydrophilic and hydrophobic model respectively in the Molecular Dynamic (MD) simulation. The Condensed-Phase Optimized Molecular Potential for Atomistic Simulation Studies (COMPASS) force field was used with the total simulation runs were set 1000 picoseconds run production ensembles. The temperature and pressure set for both ensembles were 298 K and 1 atm respectively. The validity of our model densities data was check and calculated where the deviation must be less than 6%. The comparison between hydrophobic and hydrophilic of the support membrane data was examined by the distance and magnitude of intensity of the Radial Distribution Function (RDF's) trends.
Recent breakthroughs in G protein-coupled receptor (GPCR) crystallography and the subsequent increase in number of solved GPCR structures has allowed for the unprecedented opportunity to utilize their experimental structures for structure-based drug discovery applications. As virtual screening represents one of the primary computational methods used for the discovery of novel leads, the GPCR-Bench dataset was created to facilitate comparison among various virtual screening protocols. In this study, we have benchmarked the performance of Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) in improving virtual screening enrichment in comparison to docking with Glide, using the entire GPCR-Bench dataset of 24 GPCR targets and 254,646 actives and decoys. Reranking the top 10% of the docked dataset using MM/PBSA resulted in improvements for six targets at EF1% and nine targets at EF5%, with the gains in enrichment being more pronounced at the EF1% level. We additionally assessed the utility of rescoring the top ten poses from docking and the ability of short MD simulations to refine the binding poses prior to MM/PBSA calculations. There was no clear trend of the benefit observed in both cases, suggesting that utilizing a single energy minimized structure for MM/PBSA calculations may be the most computationally efficient approach in virtual screening. Overall, the performance of MM/PBSA rescoring in improving virtual screening enrichment obtained from docking of the GPCR-Bench dataset was found to be relatively modest and target-specific, highlighting the need for validation of MM/PBSA-based protocols prior to prospective use.
This study aims at investigating the distortion of poly(dimethylsiloxane) (PDMS) nanostructures in a soft lithography demolding process using molecular dynamics simulation. Experimental results show that after peeling, PDMS nanopillars became 10-60% longer in height than the mold size. Molecular dynamics simulations have been employed to plot the stress-strain curve of the nanopillars when subjected to uniaxial stress. Three force fields (COMPASS, CVFF, and PCFF) were used for modeling. The demolding process in soft lithography and nanoimprint lithography causes significant deformation in replication. The experimental results show clear signs of elongation after demolding. Molecular dynamics simulations are employed to study the stress-strain relationship of the PDMS nanopillars. The results from the simulation show that a PDMS nanopillar at temperature T = 300 K under tensile stress shows characteristics of flexible plastic under tensile stress and has a lower Young's modulus, ultimate tensile stress, and Poisson's ratio.
This study presents the initial structural model of L-haloacid dehalogenase (DehLBHS1) from Bacillus megaterium BHS1, an alkalotolerant bacterium known for its ability to degrade halogenated environmental pollutants. The model provides insights into the structural features of DehLBHS1 and expands our understanding of the enzymatic mechanisms involved in the degradation of these hazardous pollutants. Key amino acid residues (Arg40, Phe59, Asn118, Asn176, and Trp178) in DehLBHS1 were identified to play critical roles in catalysis and molecular recognition of haloalkanoic acid, essential for efficient binding and transformation of haloalkanoic acid molecules. DehLBHS1 was modeled using I-TASSER, yielding a best TM-score of 0.986 and an RMSD of 0.53 Å. Validation of the model using PROCHECK revealed that 89.2% of the residues were located in the most favored region, providing confidence in its structural accuracy. Molecular docking simulations showed that the non-simulated DehLBHS1 preferred 2,2DCP over other substrates, forming one hydrogen bond with Arg40 and exhibiting a minimum energy of -2.5 kJ/mol. The simulated DehLBHS1 exhibited a minimum energy of -4.3 kJ/mol and formed four hydrogen bonds with Arg40, Asn176, Asp9, and Tyr11, further confirming the preference for 2,2DCP. Molecular dynamics simulations supported this preference, based on various metrics, including RMSD, RMSF, gyration, hydrogen bonding, and molecular distance. MM-PBSA calculations showed that the DehLBHS1-2,2-DCP complex had a markedly lower binding energy (-21.363 ± 1.26 kcal/mol) than the DehLBHS1-3CP complex (-14.327 ± 1.738 kcal/mol). This finding has important implications for the substrate specificity and catalytic function of DehLBHS1, particularly in the bioremediation of 2,2-DCP in contaminated alkaline environments. These results provide a detailed view of the molecular interactions between the enzyme and its substrate and may aid in the development of more efficient biocatalytic strategies for the degradation of halogenated compounds.Communicated by Ramaswamy H. Sarma.
Genomic surveillance is crucial for tracking emergence and spread of novel variants of pathogens, such as SARS-CoV-2, to inform public health interventions and to enforce control measures. However, in some settings especially in low- and middle- income counties, where sequencing platforms are limited, only certain patients get to be selected for sequencing surveillance. Here, we show that patients with multiple comorbidities potentially harbour SARS-CoV-2 with higher mutation rates and thus deserve more attention for genomic surveillance. The relationship between the patient comorbidities, and type of amino acid mutations was assessed. Correlation analysis showed that there was a significant tendency for mutations to occur within the ORF1a region for patients with higher number of comorbidities. Frequency analysis of the amino acid substitution within ORF1a showed that nsp3 P822L of the PLpro protease was one of the highest occurring mutations. Using molecular dynamics, we simulated that the P822L mutation in PLpro represents a system with lower Root Mean Square Deviation (RMSD) fluctuations, and consistent Radius of gyration (Rg), Solvent Accessible Surface Area (SASA) values-indicate a much stabler protein than the wildtype. The outcome of this study will help determine the relationship between the clinical status of a patient and the mutations of the infecting SARS-CoV-2 virus.