Rational peptide design for targeting cancer cell invasion.

Cell invasion is an important process in cancer progression and recurrence. Invasion and implantation of cancer cells from their original place to other tissues, by disabling vital organs, challenges the treatment of cancer patients. Given the importance of the matter, many molecular treatments have been developed to inhibit cancer cell invasion. Because of their low production cost and ease of production, peptides are valuable therapeutic molecules for inhibiting cancer cell invasion. In recent years, advances in the field of computational biology have facilitated the design of anti-cancer peptides. In our investigation, using computational biology approaches such as evolutionary analysis, residue scanning, protein-peptide interaction analysis, molecular dynamics, and free energy analysis, our team designed a peptide library with about 100 000 candidates based on A6 (acetyl-KPSSPPEE-amino) sequence which is an anti-invasion peptide. During computational studies, two of the designed peptides that give the highest scores and showed the greatest sequence similarity to A6 were entered into the experimental analysis workflow for further analysis. In experimental analysis steps, the anti-metastatic potency and other therapeutic effects of designed peptides were evaluated using MTT assay, RT-qPCR, zymography analysis, and invasion assay. Our study disclosed that the IK1 (acetyl-RPSFPPEE-amino) peptide, like A6, has great potency to inhibit the invasion of cancer cells.

In silico approaches to investigate enzyme immobilization: a comprehensive systematic review.

Enzymes are popular catalysts with many applications, especially in industry. Biocatalyst usage on a large scale is facing some limitations, such as low operational stability, low recyclability, and high enzyme cost. Enzyme immobilization is a beneficial strategy to solve these problems. Bioinformatics tools can often correctly predict immobilization outcomes, resulting in a cost-effective experimental phase with the least time consumed. This study provides an overview of in silico methods predicting immobilization processes via a comprehensive systematic review of published articles till 11 December 2022. It also mentions the strengths and weaknesses of the processes and explains the computational analyses in each method that are required for immobilization assessment. In this regard, Web of Science and Scopus databases were screened to gain relevant publications. After screening the gathered documents (n = 3873), 60 articles were selected for the review. The selected papers have applied in silico procedures including only molecular dynamics (MD) simulations (n = 20), parallel tempering Monte Carlo (PTMC) and MD simulations (n = 3), MD and docking (n = 1), density functional theory (DFT) and MD (n = 1), only docking (n = 11), metal ion binding site prediction (MIB) server and docking (n = 2), docking and DFT (n = 1), docking and analysis of enzyme surfaces (n = 1), only DFT (n = 1), only MIB server (n = 2), analysis of an enzyme structure and surface (n = 12), rational design of immobilized derivatives (RDID) software (n = 3), and dissipative particle dynamics (DPD; n = 2). In most included studies (n = 51), enzyme immobilization was investigated experimentally in addition to in silico evaluation.

Design, construction and in vivo functional assessment of a hinge truncated sFLT01.

Gene therapy for the treatment of ocular neovascularization has reached clinical trial phases. The AAV2-sFLT01 construct was already evaluated in a phase 1 open-label trial administered intravitreally to patients with advanced neovascular age-related macular degeneration. SFLT01 protein functions by binding to VEGF and PlGF molecules and inhibiting their activities simultaneously. It consists of human VEGFR1/Flt-1 (hVEGFR1), a polyglycine linker, and the Fc region of human IgG1. The IgG1 upper hinge region of the sFLT01 molecule makes it vulnerable to radical attacks and prone to causing immune reactions. This study pursued two goals: (i) minimizing the immunogenicity and vulnerability of the molecule by designing a truncated molecule called htsFLT01 (hinge truncated sFLT01) that lacked the IgG1 upper hinge and lacked 2 amino acids from the core hinge region; and (ii) investigating the structural and functional properties of the aforesaid chimeric molecule at different levels (in silico, in vitro, and in vivo). Molecular dynamics simulations and molecular mechanics energies combined with Poisson-Boltzmann and surface area continuum solvation calculations revealed comparable free energy of binding and binding affinity for sFLT01 and htsFLT01 to their cognate ligands. Conditioned media from human retinal pigment epithelial (hRPE) cells that expressed htsFLT01 significantly reduced tube formation in HUVECs. The AAV2-htsFLT01 virus suppressed vascular development in the eyes of newborn mice. The htsFLT01 gene construct is a novel anti-angiogenic tool with promising improvements compared to existing treatments.

Reconstructing the cell-cell interaction network among mouse immune cells.

Intercellular interactions and cell-cell communication are critical to regulating cell functions, especially in normal immune cells and immunotherapies. Ligand-receptor pairs mediating these cell-cell interactions can be identified using diverse experimental and computational approaches. Here, we reconstructed the intercellular interaction network between Mus musculus immune cells using publicly available receptor-ligand interaction databases and gene expression data from the immunological genome project. This reconstructed network accounts for 50,317 unique interactions between 16 cell types between 731 receptor-ligand pairs. Analysis of this network shows that cells of hematopoietic lineages use fewer communication pathways for interacting with each other, while nonhematopoietic stromal cells use the most network communications. We further observe that the WNT, BMP, and LAMININ pathways are the most significant contributors to the overall number of cell-cell interactions among the various pathways in the reconstructed communication network. This resource will enable the systematic analysis of normal and pathologic immune cell interactions, along with the study of emerging immunotherapies.

Exome sequencing reveals novel rare variants in Iranian familial multiple sclerosis: The importance of POLD2 in the disease pathogenesis.

The prevalence of familial multiple sclerosis (FMS) is increasing worldwide which endorses the heritability of the disease. Given that many genome variations are ethnicity-specific and consanguineous marriage could affect genetic diseases, hereditary disease gene analysis among FMS patients from Iran, a country with high rates of parental consanguinity, could be highly effective in finding mutations underlying disease pathogenesis. To examine rare genetic mutations, we selected three Iranian FMS cases with ≥3 MS patients in more than one generation and performed whole exome sequencing. We identified a homozygous rare missense variant in POLD2 (p. Arg141Cys; rs372336011). Molecular dynamics analysis showed reduced polar dehydration energy and conformational changes in POLD2 mutant. Further, we found a heterozygote rare missense variant in NBFP1 (p. Gly487Asp; rs778806175). Our study revealed the possible role of novel rare variants in FMS. Molecular dynamic simulation provided the initial evidence of the structural changes behind POLD2 mutant.

RepCOOL: computational drug repositioning via integrating heterogeneous biological networks.

BACKGROUND: It often takes more than 10 years and costs more than 1 billion dollars to develop a new drug for a particular disease and bring it to the market. Drug repositioning can significantly reduce costs and time in drug development. Recently, computational drug repositioning attracted a considerable amount of attention among researchers, and a plethora of computational drug repositioning methods have been proposed. This methodology has widely been used in order to address various medical challenges, including cancer treatment. The most common cancers are lung and breast cancers. Thus, suggesting FDA-approved drugs via drug repositioning for breast cancer would help us to circumvent the approval process and subsequently save money as well as time. METHODS: In this study, we propose a novel network-based method, named RepCOOL, for drug repositioning. RepCOOL integrates various heterogeneous biological networks to suggest new drug candidates for a given disease. RESULTS: The proposed method showed a promising performance on benchmark datasets via rigorous cross-validation. The final drug repositioning model has been built based on a random forest classifier after examining various machine learning algorithms. Finally, in a case study, four FDA approved drugs were suggested for breast cancer stage II. CONCLUSION: Results show the potency of the proposed method in detecting true drug-disease relationships. RepCOOL suggested four new drugs for breast cancer stage II namely Doxorubicin, Paclitaxel, Trastuzumab, and Tamoxifen.

Multifunctional Acidocin 4356 Combats Pseudomonas aeruginosa through Membrane Perturbation and Virulence Attenuation: Experimental Results Confirm Molecular Dynamics Simulation.

A longstanding awareness in generating resistance to common antimicrobial therapies by Gram-negative bacteria has made them a major threat to global health. The application of antimicrobial peptides as a therapeutic agent would be a great opportunity to combat bacterial diseases. Here, we introduce a new antimicrobial peptide (∼8.3 kDa) from probiotic strain Lactobacillus acidophilus ATCC 4356, designated acidocin 4356 (ACD). This multifunctional peptide exerts its anti-infective ability against Pseudomonas aeruginosa through an inhibitory action on virulence factors, bacterial killing, and biofilm degradation. Reliable performance over tough physiological conditions and low hemolytic activity confirmed a new hope for the therapeutic setting. Antibacterial kinetic studies using flow cytometry technique showed that the ACD activity is related to the change in permeability of the membrane. The results obtained from molecular dynamic (MD) simulation were perfectly suited to the experimental data of ACD behavior. The structure-function relationship of this natural compound, along with the results of transmission electron microscopy analysis and MD simulation, confirmed the ability of the ACD aimed at enhancing bacterial membrane perturbation. The peptide was effective in the treatment of P. aeruginosa infection in mouse model. The results support the therapeutic potential of ACD for the treatment of Pseudomonas infections.IMPORTANCE Multidrug-resistant bacteria are a major threat to global health, and the Pseudomonas bacterium with the ability to form biofilms is considered one of the main causative agents of nosocomial infections. Traditional antibiotics have failed because of increased resistance. Thus, finding new biocompatible antibacterial drugs is essential. Antimicrobial peptides are produced by various organisms as a natural defense mechanism against pathogens, inspiring the possible design of the next generation of antibiotics. In this study, a new antimicrobial peptide was isolated from Lactobacillus acidophilus ATCC 4356, counteracting both biofilm and planktonic cells of Pseudomonas aeruginosa A detailed investigation was then conducted concerning the functional mechanism of this peptide by using fluorescence techniques, electron microscopy, and in silico methods. The antibacterial and antibiofilm properties of this peptide may be important in the treatment of Pseudomonas infections.

177Lu-labeled cyclic RGD peptide as an imaging and targeted radionuclide therapeutic agent in non-small cell lung cancer: Biological evaluation and preclinical study.

Non-small cell lung carcinoma (NSCLC) is among the most lethal lung cancers responsible for 80-85% of death. αvß3 integrin receptor subtype has been identified as a lung cancer biomarker since its expression correlates with tumor progression and metastasis. The extracellular domain of the receptor forms a binding site for RGD-based sequences. Therefore, specific targeting of αvß3 integrin receptors by these short peptides can be an excellent candidate for cancer imaging and therapy. In this research, the radiolabeling of DOTA-E(cRGDfK)2 with 177Lu was efficiently implemented. The Log P value, in vivo, in vitro, metabolic stability, cellular uptake and specific binding of the radiopeptide was determined. The tumor targeting capacity and the therapeutic potential of the radiotracer was studied in A549 tumor-bearing mice. Imaging studies at different time intervals were performed by SPECT/CT. Radiochemical purity of more than 99% and Log P of -3.878 was obtained for 177Lu-labelled peptide. Radiotracer showed favorable in vivo, in vitro and metabolic stability. The radiopeptide dissociation constant (Kd) was 15.07 nM. Radiopeptide specific binding was more than 95%. Biodistribution studies showed high accumulation of the radiopeptide in tumor and rapid excretion by urinary route. Maximum tumor uptake was at 4 h post-injection. Following administration of this radiopeptide to mice, not only tumor growth was suppressed, but significant tumor shrinkage was also observed. In conclusion, this radiopeptide can be employed for staging, follow-up imaging and as peptide receptor radionuclide therapeutic agent allowing efficient therapy for NSCLC and other cancers overexpressing αvß3 integrin receptors.

Production and characterization of virus-like particles of grapevine fanleaf virus presenting L2 epitope of human papillomavirus minor capsid protein.

BACKGROUND: Virus-like particle (VLP) platform represents a promising approach for the generation of efficient and immunogenic subunit vaccines. Here, the feasibility of using grapevine fanleaf virus (GFLV) VLPs as a new carrier for the presentation of human papillomavirus (HPV) L2 epitope was studied. To achieve this goal, a model of the HPV L2 epitope secondary structure was predicted and its insertion within 5 external loops in the GFLV capsid protein (CP) was evaluated. RESULTS: The epitope sequence was genetically inserted in the αB-αB" domain C of the GFLV CP, which was then over-expressed in Pichia pastoris and Escherichia coli. The highest expression yield was obtained in E. coli. Using this system, VLP formation requires a denaturation-refolding step, whereas VLPs with lower production yield were directly formed using P. pastoris, as confirmed by electron microscopy and immunostaining electron microscopy. Since the GFLV L2 VLPs were found to interact with the HPV L2 antibody under native conditions in capillary electrophoresis and in ELISA, it can be assumed that the inserted epitope is located at the VLP surface with its proper ternary structure. CONCLUSIONS: The results demonstrate that GFLV VLPs constitute a potential scaffold for surface display of the epitope of interest.

Improving the thermal stability of azoreductase from Halomonas elongata by introducing a disulfide bond via site-directed mutagenesis.

Azoreductases mainly reduce azo dyes, the largest class of colorants, to colorless aromatic amines. AzoH, a new azoreductase from the halophilic bacterium, Halomonas elongata, has been recently cloned and expressed in Escherichia coli. The aim of this study was to improve thermal stability of this enzyme by introducing new disulfide bonds. Since X-ray crystallography was not available, homology modeling and molecular dynamics was used to construct the enzyme three-dimensional structure. Potential disulfide bonds for increasing thermal stability were found using DIScover online software. Appropriate mutations (L49C/D108C) to form a disulfide bond were introduced by the Quik-Change method. Mutant protein expressed in E. coli showed increased thermal stability at 50 °C (increased half-life from 12.6 Min in AzoH to 26.66 Min in a mutated enzyme). The mutated enzyme could also tolerate 5% (w/v) NaCl and retained 30% of original activity after 24 H incubation, whereas the wild-type enzyme was completely inactivated. According to circular dichroism studies, the secondary structure was not altered by this mutation; however, a blue shift in intrinsic florescent graph revealed changes in the tertiary structure. This is the first study to improve thermal stability and salt tolerance of a halophilic azoreductase.

The inhibitory effects of biomimetically designed peptides on α-synuclein aggregation.

Parkinson's disease is characterized by accumulation of inclusion bodies in dopaminergic neurons, where insoluble and fibrillar α-synuclein makes up the major component of these inclusion bodies. So far, several strategies have been applied in order to suppress α-synuclein aggregation and toxicity in Parkinson's disease. In the present study, a new database has been established by segmentation of all the proteins deposited in protein Data Bank. The database data base was searched for the sequences which adopt ß structure and are identical or very similar to the regions of α-synuclein which are involved in aggregation. The adjacent ß strands of the found sequences were chosen as the peptide inhibitors of α-synuclein aggregation. Two of the predicted peptides, namely KISVRV and GQTYVLPG, were experimentally proved to be efficient in suppressing aggregation of α-synuclein in vitro. Moreover, KISVRV exhibited the ability to disrupt oligomers of α-syn which are assumed to be the pathogenic species in Parkinson's disease.

Increase of Bacillus badius Phenylalanine dehydrogenase specificity towards phenylalanine substrate by site-directed mutagenesis.

Phenylalanine dehydrogenase (PheDH) is a key enzyme in medical diagnostic for determining the amount of phenylalanine to detect phenylketonuria (PKU) disease. However, determination of phenylalanine can be usually disturbed in presence of tyrosine in blood samples. Position N145 of B.sphaericus PheDH, has been previously showed a crucial role in substrate binding, which corresponded by position V144 in B. badius PheDH. In this study, the PheDH of B. badius due to reasonable activity was cloned and subjected to site-directed mutagenesis at mentioned position, followed by kinetic and structural studies to find more exclusive mutants. The results showed that the V144L mutant considerably increases specificity toward phenylalanine and decreases toward l-tyrosine, while in V144N mutant, the specificity reduces toward phenylalanine and increases toward tyrosine. Moreover, concerning the mutated V144D, significantly reduced kcat and also decreased km value for phenylalanine relative to that of wild type. The Phe/Tyr specificity constant in V144L increased more than 4-fold compared to wild type, makes it to be a suitable candidate for more specific identification of PKU. Finally, docking and molecular dynamic simulation on wild type and mutants clarified the structural basis behind more specificity of V144L mutant for phenylalanine substrate.

Evaluating the accuracy of protein design using native secondary sub-structures.

BACKGROUND: According to structure-dependent function of proteins, two main challenging problems called Protein Structure Prediction (PSP) and Inverse Protein Folding (IPF) are investigated. In spite of IPF essential applications, it has not been investigated as much as PSP problem. In fact, the ultimate goal of IPF problem or protein design is to create proteins with enhanced properties or even novel functions. One of the major computational challenges in protein design is its large sequence space, namely searching through all plausible sequences is impossible. Inasmuch as, protein secondary structure represents an appropriate primary scaffold of the protein conformation, undoubtedly studying the Protein Secondary Structure Inverse Folding (PSSIF) problem is a quantum leap forward in protein design, as it can reduce the search space. In this paper, a novel genetic algorithm which uses native secondary sub-structures is proposed to solve PSSIF problem. In essence, evolutionary information can lead the algorithm to design appropriate amino acid sequences respective to the target secondary structures. Furthermore, they can be folded to tertiary structures almost similar to their reference 3D structures. RESULTS: The proposed algorithm called GAPSSIF benefits from evolutionary information obtained by solved proteins in the PDB. Therefore, we construct a repository of protein secondary sub-structures to accelerate convergence of the algorithm. The secondary structure of designed sequences by GAPSSIF is comparable with those obtained by Evolver and EvoDesign. Although we do not explicitly consider tertiary structure features through the algorithm, the structural similarity of native and designed sequences declares acceptable values. CONCLUSIONS: Using the evolutionary information of native structures can significantly improve the quality of designed sequences. In fact, the combination of this information and effective features such as solvent accessibility and torsion angles leads IPF problem to an efficient solution. GAPSSIF can be downloaded at http://bioinformatics.aut.ac.ir/GAPSSIF/ .

Appraisal of role of the polyanionic inducer length on amyloid formation by 412-residue 1N4R Tau protein: A comparative study.

In many neurodegenerative diseases, formation of protein fibrillar aggregates has been observed as a major pathological change. Neurofibrillary tangles, mainly composed of fibrils formed by the microtubule-associated protein; Tau, are a hallmark of a group of neurodegenerative diseases such as Alzheimer's disease. Tau belongs to the class of natively unfolded proteins and partially folds into an ordered ß-structure during aggregation. Polyanionic cofactors such as heparin are commonly used as inducer of Tau aggregation in vitro. The role of heparin in nucleation and elongation steps during Tau fibril formation is not fully understood. In the current study, aggregation kinetics as well as structure of Tau amyloid fibrils, by using the 1N4R isoform, have been reproducibly determined in the presence of heparin and the shorter molecule; enoxaparin. The kinetic studies demonstrated that heparin (not enoxaparin) efficiently accelerates Tau amyloid formation and revealed, mechanistically, that the molecular weight of the inducer is important in accelerating amyloidogenesis. The kinetic parameter values of Tau amyloid aggregation, especially, the amyloid aggregation extent, were relatively different in the presence of heparin and enoxaparin, at various stoichiometries of the inducers binding. Also, based on the results, obtained from CD, FTIR, AFM and XRD studies, it may be suggested that the inducer length plays a critical role mainly in the nucleation process, so that it determines that oligomers lie on or off the pathway of Tau fibrillization. The biochemical results herein suggest that the chemical environment of the extracellular matrix as well as localization of distinct glycosaminoglycans may influence deposition behavior of Tau amyloidosis.

Inhibition of guinea pig aldehyde oxidase activity by different flavonoid compounds: An in vitro study.

Aldehyde oxidase (AO), a cytosolic molybdenum-containing hydroxylase, is predominantly active in liver and other tissues of mammalian species and involved in the metabolism of extensive range of aldehydes and nitrogen-containing compounds. A wide range of natural components including polyphenols are able to interfere with AO-catalyzed reactions. Polyphenols and flavonoids are one of the extensive secondary plant metabolites ubiquitously present in plants considered an important part of the human diet. The aim of the present study was to investigate inhibitory effect of selected phenolic compounds from three subclasses of aurone, flavanone and phenolic lactone compounds on the activity of AO, spectrophotometrically. AO enzyme was partially purified from liver of guinea pig. Then, inhibitory effects of 10 flavonoid compounds including 8 derivatives of 2-benzylidenebenzofuran-3(2H)-ones, as well as naringenin and ellagic acid on the activity of aldehyde oxidase were assessed compared with the specific inhibitor of AO, menadione. Among the phenolic compounds with inhibitory effects on the enzyme, ellagic acid (IC50=14.47 µM) was the most potent agent with higher inhibitory action than menadione (IC50=31.84 µM). The mechanisms by which flavonoid compounds inhibit AO activity have been also determined. The inhibitory process of the assessed compounds occurs via either a non-competitive or mixed mode. Although flavonoid compounds extensively present in the nature, mainly in dietary regimen, aurones with promising biological properties are not widely distributed in nature, so synthesis of aurone derivatives is of great importance. Additionally, aurones seem to provide a promising scaffold in medicinal chemistry for the skeleton of new developing drugs, so the results of the current study can be valuable in order to better understanding drug-food as well as drug-drug interaction and also appears to be worthwhile in drug development strategies.

MemBuilder: a web-based graphical interface to build heterogeneously mixed membrane bilayers for the GROMACS biomolecular simulation program.

MOTIVATION: Molecular dynamics (MD) simulations have had a profound impact on studies of membrane proteins during past two decades, but the accuracy of MD simulations of membranes is limited by the quality of membrane models and the applied force fields. Membrane models used in MD simulations mostly contain one kind of lipid molecule. This is far from reality, for biological membranes always contain more than one kind of lipid molecule. Moreover, the lipid composition and their distribution are functionally important. As a result, there is a necessity to prepare more realistic lipid membranes containing different types of lipids at physiological concentrations. RESULTS: To automate and simplify the building process of heterogeneous lipid bilayers as well as providing molecular topologies for included lipids based on both united and all-atom force fields, we provided MemBuilder as a web-based graphical user interface. AVAILABILITY AND IMPLEMENTATION: MemBuilder is a free web server available from www.membuilder.org.

NETAL: a new graph-based method for global alignment of protein-protein interaction networks.

MOTIVATION: The interactions among proteins and the resulting networks of such interactions have a central role in cell biology. Aligning these networks gives us important information, such as conserved complexes and evolutionary relationships. Although there have been several publications on the global alignment of protein networks; however, none of proposed methods are able to produce a highly conserved and meaningful alignment. Moreover, time complexity of current algorithms makes them impossible to use for multiple alignment of several large networks together. RESULTS: We present a novel algorithm for the global alignment of protein-protein interaction networks. It uses a greedy method, based on the alignment scoring matrix, which is derived from both biological and topological information of input networks to find the best global network alignment. NETAL outperforms other global alignment methods in terms of several measurements, such as Edge Correctness, Largest Common Connected Subgraphs and the number of common Gene Ontology terms between aligned proteins. As the running time of NETAL is much less than other available methods, NETAL can be easily expanded to multiple alignment algorithm. Furthermore, NETAL overpowers all other existing algorithms in term of performance so that the short running time of NETAL allowed us to implement it as the first server for global alignment of protein-protein interaction networks. AVAILABILITY: Binaries supported on linux are freely available for download at http://www.bioinf.cs.ipm.ir/software/netal. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Kallistatin as an inhibitory protein against colorectal cancer cells through binding to LRP6.

Kallistatin (KL) is a member of the serine proteinase inhibitor (serpin) family regulating oxidative stress, vascular relaxation, inflammation, angiogenesis, cell proliferation, and invasion. The heparin-binding site of Kallistatin has an important role in the interaction with LRP6 leading to the blockade of the Wnt signaling pathway. In this study, we aimed to explore the structural basis of the Kallistatin-LRP6E1E4 complex using in silico approaches and evaluating the anti-proliferative, apoptotic, and cell cycle arrest activities of Kallistatin in colon cancer lines. The molecular docking showed Kallistatin could bind to the LRP6E3E4 much stronger than LRP6E1E2. The Kallistatin-LRP6E1E2 and Kallistatin-LRP6E3E4 complexes were stable during Molecular Dynamics (MD) simulation. The Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) showed that the Kallistatin-LRP6E3E4 has a higher binding affinity compared to Kallistatin-LRP6E1E2. Kallistatin induced higher cytotoxicity and apoptosis in HCT116 compared to the SW480 cell line. This protein-induced cell-cycle arrest in both cell lines at the G1 phase. The B-catenin, cyclin D1, and c-Myc expression levels were decreased in response to treatment with Kallistatin in both cell lines while the LRP6 expression level was decreased in the HCT116 cell line. Kallistatin has a greater effect on the HCT116 cell line compared to the SW480 cell line. Kallistatin can be used as a cytotoxic and apoptotic-inducing agent in colorectal cancer cell lines.

Molecular dynamics simulation as a promising approach for computational study of liquid crystal-based aptasensors.

As a potent computational methodology, molecular dynamics (MD) simulation provides advantageous knowledge about biological compounds from the molecular viewpoint. In particular, MD simulation gives exact information about aptamer strands, such as the short synthetic oligomers, their orientation, binding sites, folding-unfolding state, and conformational re-arrangement. Also, the effect of the different chemicals and biochemicals as the components of aptamer-based sensors (aptasensors) on the aptamer-target interaction can be investigated by MD simulation. Liquid crystals (LCs) as soft substances with characteristics of both solid anisotropy and liquid fluidity are new candidates for designing label-free aptasensors. To now, diverse aptasensors have been developed experimentally based on the optical anisotropy, fluidity, and long-range orientational order of LCs. Here, we represent a computational model of an LC-based aptasensor through a detailed MD simulation study. The different parameters are defined and studied to achieve a comprehensive understanding of the computational design of the LC-based aptasensor, including the density of LCs, their orientation angle, and lognormal distribution in the absence and presence of aptamer strands, both aptamer and target molecules with various concentrations, and interfering substance. As a case study, the tobramycin antibiotic is considered the target molecule for the computational model of the LC-based aptasensor.Communicated by Ramaswamy H. Sarma.

In silico design of potential Mcl-1 peptide-based inhibitors.

CONTEXT: BIM (Bcl-2 interacting mediator of apoptosis)-derived peptides that specifically target over-expressed Mcl-1 (myeloid cell leukemia-1) protein and induce apoptosis are potentially anti-cancer agents. Since the helicity of BIM-derived peptides has a crucial role in their functionality, a range of strategies have been used to increase the helicity including the introduction of unnatural residues and stapling methods that have some drawbacks such as the accumulation in the liver. To avoid these drawbacks, this study aimed to design a more helical peptide by utilizing bioinformatics algorithms and molecular dynamics simulations without exploiting unnatural residues and stapling methods. MM-PBSA results showed that the mutations of A4fE and A2eE in analogue 5 demonstrate a preference towards binding with Mcl-1. As evidenced by Circular dichroism results, the helicity increases from 18 to 34%, these findings could enhance the potential of analogue 5 as an anti-cancer agent targeting Mcl-1. The applied strategies in this research could shed light on the in silico peptide design. Moreover, analogue 5 as a drug candidate can be evaluated in vitro and in vivo studies. METHODS: The sequence of the lead peptide was determined using the ApInAPDB database and PRALINE program. Contact finder and PDBsum web server softwares were used to determine the contact involved amino acids in complex with Mcl-1. All identified salt bridge contributing residues were unaltered to preserve the binding affinity. After proposing novel analogues, their secondary structures were predicted by Cham finder web server software and GOR, Neural Network, and Chou-Fasman algorithms. Finally, molecular dynamics simulations run for 100 ns were done using the GROMACS, version 5.0.7, with the CHARMM36 force field. MM-PBSA was used to assess binding affinity specificity in targeting Mcl-1 and Bcl-xL (B-cell lymphoma extra-large).