Molecular Dynamic Simulation of Space and Earth-Grown Crystal Structures of Thermostable T1 Lipase Geobacillus zalihae Revealed a Better Structure

CHAPERONg: A tool for automated GROMACS-based molecular dynamics simulations and trajectory analyses

Molecular dynamics (MD) simulation is a powerful computational tool used in biomolecular studies to investigate the dynamics, energetics, and interactions of a wide range of biological systems at the atomic level. GROMACS is a widely used free and open-source biomolecular MD simulation software recognized for its efficiency, accuracy, and extensive range of simulation options. However, the complexity of setting up, running, and analyzing MD simulations for diverse systems often poses a significant challenge, requiring considerable time, effort, and expertise. Here, we introduce CHAPERONg, a tool that automates the GROMACS MD simulation pipelines for protein and protein-ligand systems. CHAPERONg also integrates seamlessly with GROMACS modules and third-party tools to provide comprehensive analyses of MD simulation trajectories, offering up to 20 post-simulation processing and trajectory analyses. It also streamlines and automates established pipelines for conducting and analyzing biased MD simulations via the steered MD-umbrella sampling workflow. Thus, CHAPERONg makes MD simulations more accessible to beginner GROMACS users whilst empowering experts to focus on data interpretation and other less programmable aspects of MD simulation workflows. CHAPERONg is written in Bash and Python, and the source code is freely available at https://github.com/abeebyekeen/CHAPERONg. Detailed documentation and tutorials are available online at dedicated web pages accessible via https://abeebyekeen.com/chaperong-online.

Hollow Nanomaterials in Advanced Drug Delivery Systems: From Single- to Multiple Shells

Hollow-structured nanomaterials (HSNMs) have attracted increased interest in biomedical fields, owing to their excellent potential as drug delivery systems (DDSs) for clinical applications. Among HSNMs, hollow multi-shelled structures (HoMSs) exhibit properties such as high loading capacity, sequential drug release, and multi-functionalized modification and represent a new class of nanoplatforms for clinical applications. The remarkable properties of HoMS-based DDS can simultaneously satisfy and enhance DDSs for delivering small molecular drugs (e.g., antibiotics, chemotherapy drugs, and imaging agents) and macromolecular drugs (e.g., protein/peptide- and nucleic acid-based drugs). First, the latest research advances in delivering small molecular drugs are summarized and highlight the inherent advantages of HoMS-based DDSs for small molecular drug targeting, combining continuous therapeutic drug delivery and theranostics to optimize the clinical benefit. Meanwhile, the macromolecular drugs DDSs are in the initial development stage and currently offer limited delivery modes. There is a growing need to analyze the deficiency of other HSNMs and integrate the advantages of HSNMs, providing solutions for the safe, stable, and cascade delivery of macromolecular drugs to meet vast treatment requirements. Therefore, the latest advances in HoMS-based DDSs are comprehensively reviewed, mainly focusing on the characteristics, research progress by drug category, and future research prospects.

Knotting terminal ends of mutant T1 lipase with disulfide bond improved structure rigidity and stability

Lipase biocatalysts offer unique properties which are often impaired by low thermal and methanol stability. In this study, the rational design was employed to engineer a disulfide bond in the protein structure of Geobacillus zalihae T1 lipase in order to improve its stability. The selection of targeted disulfide bond sites was based on analysis of protein spatial configuration and change of Gibbs free energy. Two mutation points (S2C and A384C) were generated to rigidify the N-terminal and C-terminal regions of T1 lipase. The results showed the mutant 2DC lipase improved methanol stability from 35 to 40% (v/v) after 30 min of pre-incubation. Enhancement in thermostability for the mutant 2DC lipase at 70 °C and 75 °C showed higher half-life at 70 °C and 75 °C for 30 min and 52 min, respectively. The mutant 2DC lipase maintained the same optimum temperature (70 °C) as T1 lipase, while thermally induced unfolding showed the mutant maintained higher rigidity. The kcat/Km values demonstrated a relatively small difference between the T1 lipase (WT) and 2DC lipase (mutant). The kcat/Km (s^-1 mM^-1) of the T1 and 2DC showed values of 13,043 ± 224 and 13,047 ± 312, respectively. X-ray diffraction of 2DC lipase crystal structure with a resolution of 2.04 Å revealed that the introduced single disulfide bond did not lower initial structural interactions within the residues. Enhanced methanol and thermal stability are suggested to be strongly related to the newly disulfide bridge formation and the enhanced compactness and rigidity of the mutant structure. KEY POINTS: • Protein engineering via rational design revealed relative improved enzymatic performance. • The presence of disulfide bond impacts on the rigidity and structural function of proteins. • X-ray crystallography reveals structural changes accompanying protein modification.

DPP-4 inhibition mediated antidiabetic potential of phytoconstituents of an aqueous fruit extract of Withania coagulans (Stocks) Dunal: in-silico, in-vitro and in-vivo assessments

The DPP-4 inhibition is an interesting target for the development of antidiabetic agents which promotes the longevity of GPL-1(Glucagon-like peptide 1). The current study was intended to assess DPP-4(Dipeptidyl Peptidase-4) inhibition mediated antidiabetic effect of phytocompounds of an aqueous fruit extract of Withania coagulans (Stocks) Dunal by in-vitro, in-silico and in-vivo approaches. The phytoconstituents screening was executed by LCMS (Liquid Chromatography with tandem mass spectrometry). The in-vitro and in-vivo, DPP-4 assays were performed by using available kits. The in-vitro DPP-4 activity was inhibited up to 68.3% by the test extract. Accordingly, in-silico determinations of molecular docking, molecular dynamics and pharmacokinetics were performed between the target enzyme DPP-4 and leading phytocompounds. The molecular dynamics authenticated the molecular docking data by crucial parameters of cytosolic milieu by the potential energy, RSMD (Root Mean Square Deviation), RSMF (Root Mean Square Fluctuation), system density, NVT (Number of particles at fixed volume, ensemble) and NPT (Number of particles at fixed pressure, ensemble). Accordingly, ADMET predictions assessed the druggability profile. Subsequently, the course of the test extract and the sitagliptin (positive control), instigated significant (p ≤ 0.001) ameliorations in HOMA indices and the equal of antioxidants in nicotinamide-streptozotocin induced type 2 diabetic animal model. Compassionately, the histopathology represented increased pancreatic cellular mass which caused in restoration of histoarchitectures. It has been concluded that phytoconstituents in W. coagulans aqueous fruit extract can regulate DPP-4, resulting in improved glucose homeostasis and enhanced endocrinal pancreatic cellular mass.Communicated by Ramaswamy H. Sarma.

In silico discovery of novel inhibitors from Northern African natural products database against main protease (Mpro) of SARS-CoV-2

The recent outbreak of COVID-19 (Coronavirus Disease 2019), caused by a novel SARS-CoV-2 virus, has led to public health emergencies worldwide where time is as important as equipment to save lives. Antimalarial drugs such as hydroxychloroquine and chloroquine derivatives are used in emergencies but they are not suitable for patients with high blood pressure, diabetes and heart problems. Since there are no approved drugs for this disease, science is challenged to find vaccines and new drugs. Therefore, as part of our Silico drug design strategy, we identified drug-like compounds that inhibit replication of the main protease (Mpro) of SARS-CoV-2 based on receptor-based virtual database screening, molecular docking, molecular dynamics, and drug-similarity profiling from the NANPDB natural products database available at North African. The two resulting hit compounds named 5- Chloro-Omega-hydroxy-1-O-methylemodin and cystodion E showed the highest binding energy with Mpro of SARS-CoV-2 and strong inhibitory activity compared with the previously published N3 inhibitor. The complexes of these two compounds were validated by molecular dynamics analysis (RMSD, RMSF, Rg, total number of hydrogen bonds and secondary structure fractions of the protein in the complex) as the best method to evaluate the biological stability of the system. Therefore, these molecules deserve more attention in drug development compared to COVID-19. HighlightsA large database of natural compounds was screened against nCoV-2's Mpro.Molecular docking and Molecular dynamics were used as powerful methods.Two compounds were found are very attractive to inhibit Mpro of nCoV-2.ADME-Tox profiling is evaluated the active compounds are not cancerogenic.Communicated by Ramaswamy H. Sarma.

Improved resolution crystal structure of Acanthamoeba actophorin reveals structural plasticity not induced by microgravity

Actophorin, a protein that severs actin filaments isolated from the amoeba Acanthamoeba castellanii, was employed as a test case for crystallization under microgravity. Crystals of purified actophorin were grown under microgravity conditions aboard the International Space Station (ISS) utilizing an interactive crystallization setup between the ISS crew and ground-based experimenters. Crystals grew in conditions similar to those grown on earth. The structure was solved by molecular replacement at a resolution of 1.65 Å. Surprisingly, the structure reveals conformational changes in a remote β-turn region that were previously associated with actophorin phosphorylated at the terminal residue Ser1. Although crystallization under microgravity did not yield a higher resolution than crystals grown under typical laboratory conditions, the conformation of actophorin obtained from solving the structure suggests greater flexibility in the actophorin β-turn than previously appreciated and may be beneficial for the binding of actophorin to actin filaments.

Structure elucidation and docking analysis of 5M mutant of T1 lipase Geobacillus zalihae

5M mutant lipase was derived through cumulative mutagenesis of amino acid residues (D43E/T118N/E226D/E250L/N304E) of T1 lipase from Geobacillus zalihae. A previous study revealed that cumulative mutations in 5M mutant lipase resulted in decreased thermostability compared to wild-type T1 lipase. Multiple amino acids substitution might cause structural destabilization due to negative cooperation. Hence, the three-dimensional structure of 5M mutant lipase was elucidated to determine the evolution in structural elements caused by amino acids substitution. A suitable crystal for X-ray diffraction was obtained from an optimized formulation containing 0.5 M sodium cacodylate trihydrate, 0.4 M sodium citrate tribasic pH 6.4 and 0.2 M sodium chloride with 2.5 mg/mL protein concentration. The three-dimensional structure of 5M mutant lipase was solved at 2.64 Å with two molecules per asymmetric unit. The detailed analysis of the structure revealed that there was a decrease in the number of molecular interactions, including hydrogen bonds and ion interactions, which are important in maintaining the stability of lipase. This study facilitates understanding of and highlights the importance of hydrogen bonds and ion interactions towards protein stability. Substrate specificity and docking analysis on the open structure of 5M mutant lipase revealed changes in substrate preference. The molecular dynamics simulation of 5M-substrates complexes validated the substrate preference of 5M lipase towards long-chain p-nitrophenyl-esters.

Exploring the potential of novel phenolic compounds as potential therapeutic candidates against SARS-CoV-2, using quantum chemistry, molecular docking and dynamic studies

In the current study, the interaction of SARS-CoV-2 protein (A and B chains of nsp13) with different recently synthesized phenolic compounds (Sreenivasulu et al., Synthetic Communications, 2020, 112-122) has been studied. The interactions have been investigated by using molecular docking, quantum chemical and molecular dynamics simulations methods. The molecular structures of all the ligands are studied quantum chemically in terms of their optimized structures, 3-D orbital distributions, global chemical descriptors, molecular electrostatic potential plots and HOMO-LUMO orbital energies. All the ligands show reasonably good binding affinities with nsp-13 protein. The ligand L2 shows to have better binding affinities to Chain A and Chain B of nsp13 protein, which are -6.7 and -6.4 kcal/mol. The study of intermolecular interactions indicates that L2 shows different hydrophobic and hydrogen bond interactions with both chains. Furthermore, molecular dynamic simulations of the nsp13-L2 complex are obtained over a time scale of 60 ns, which indicates its stability and flexibility behavior as assessed in terms of its RMSD and RMSF graphs. The ADMET analysis also shows no violation of Lipinski rule (RO5) by studied phenolic compounds. We believe that the current findings will be further confirmed by in vitro and in vivo studies of these recent phenolic compounds for their potential as inhibitors for SARS-Co-V-2 virus.

Effect and mechanism analysis of different linkers on efficient catalysis of subunit-fused nitrile hydratase

Protein fusion using a linker plays an important role for protein evolution. However, designing suitable linkers for protein evolution is yet challenging and under-explored. To further clarify the regular pattern of suitable type of linker for fusion proteins, one nitrile hydratase (NHase) was used as a target protein and subunit fusion strategy was carried out to improve its efficient catalysis. Subunit-fused variants with three different types of linkers were constructed and characterized. All variants exhibited higher stability than that of the wild type. The longer the linker was, the higher stability NHase showed, however, too long linker affected NHase activity and expression. Among the three types of linkers, the α-helical linker seemed more suitable for NHase than flexible or rigid linkers. Though it is not clear how the linkers affecting the activity, structure analysis indicated that the stability improvement is dependent on the additional salt bridge, H-bond, and the subunit interface area increasing due to the linker insertion, among which the additional salt bridge and interface area were more important factors. The results described here may be useful for redesigning other enzymes through subunit fusion.

Ion-Pair Interaction and Hydrogen Bonds as Main Features of Protein Thermostability in Mutated T1 Recombinant Lipase Originating from Geobacillus zalihae

A comparative structure analysis between space- and an Earth-grown T1 recombinant lipase from Geobacillus zalihae had shown changes in the formation of hydrogen bonds and ion-pair interactions. Using the space-grown T1 lipase validated structure having incorporated said interactions, the recombinant T1 lipase was re-engineered to determine the changes brought by these interactions to the structure and stability of lipase. To understand the effects of mutation on T1 recombinant lipase, five mutants were developed from the structure of space-grown T1 lipase and biochemically characterized. The results demonstrate an increase in melting temperature up to 77.4 °C and 76.0 °C in E226D and D43E, respectively. Moreover, the mutated lipases D43E and E226D had additional hydrogen bonds and ion-pair interactions in their structures due to the improvement of stability, as observed in a longer half-life and an increased melting temperature. The biophysical study revealed differences in β-Sheet percentage between less stable (T118N) and other mutants. As a conclusion, the comparative analysis of the tertiary structure and specific residues associated with ion-pair interactions and hydrogen bonds could be significant in revealing the thermostability of an enzyme with industrial importance.

Changes of Thermostability, Organic Solvent, and pH Stability in Geobacillus zalihae HT1 and Its Mutant by Calcium Ion

Thermostable T1 lipase from Geobacillus zalihae has been crystallized using counter-diffusion method under space and Earth conditions. The comparison of the three-dimensional structures from both crystallized proteins show differences in the formation of hydrogen bond and ion interactions. Hydrogen bond and ion interaction are important in the stabilization of protein structure towards extreme temperature and organic solvents. In this study, the differences of hydrogen bond interactions at position Asp43, Thr118, Glu250, and Asn304 and ion interaction at position Glu226 was chosen to imitate space-grown crystal structure, and the impact of these combined interactions in T1 lipase-mutated structure was studied. Using space-grown T1 lipase structure as a reference, subsequent simultaneous mutation D43E, T118N, E226D, E250L, and N304E was performed on recombinant wild-type T1 lipase (wt-HT1) to generate a quintuple mutant term as 5M mutant lipase. This mutant lipase shared similar characteristics to its wild-type in terms of optimal pH and temperature. The stability of mutant 5M lipase improved significantly in acidic and alkaline pH as compared to wt-HT1. 5M lipase was highly stable in organic solvents such as dimethyl sulfoxide (DMSO), methanol, and n-hexane compared to wt-HT1. Both wild-type and mutant lipases were found highly activated in calcium as compared to other metal ions due to the presence of calcium-binding site for thermostability. The presence of calcium prolonged the half-life of mutant 5M and wt-HT1, and at the same time increased their melting temperature (T_m). The melting temperature of 5M and wt-HT1 lipases increased at 8.4 and 12.1 °C, respectively, in the presence of calcium as compared to those without. Calcium enhanced the stability of mutant 5M in 25% (v/v) DMSO, n-hexane, and n-heptane. The lipase activity of wt-HT1 also increased in 25% (v/v) ethanol, methanol, acetonitrile, n-hexane, and n-heptane in the presence of calcium. The current study showed that the accumulation of amino acid substitutions D43E, T118N, E226D, E250L, and N304E produced highly stable T1 mutant when hydrolyzing oil in selected organic solvents such as DMSO, n-hexane, and n-heptane. It is also believed that calcium ion plays important role in regulating lipase thermostability.