Capparis cartilaginea decne (capparaceae): isolation of flavonoids by high-speed countercurrent chromatography and their anti-inflammatory evaluation

Introduction: Capparis cartilaginea Decne. (CC) originates from the dry regions of Asia and the Mediterranean basin. In traditional medicine, tea of CC leaves is commonly used to treat inflammatory conditions such as rheumatism, arthritis, and gout. Due to the limited studies on the phytochemistry and biological activity of CC compared to other members of the Capparaceae family, this work aims to: 1) Identify the chemical composition of CC extract and 2) Investigate the potential anti-inflammatory effect of CC extract, tea and the isolated compounds. Methods: To guarantee aim 1, high-speed countercurrent chromatography (HSCC) method; Nuclear Magnetic Resonance (NMR) and High-Performance Liquid Chromatography coupled to Electrospray Ionisation and Quadrupole Time-of-Flight Mass Spectrometry (HPLC-ESIQTOF-MS/MS) were employed for this purpose. To guarantee aim 2, we studied the effect of the isolated flavonoids on matrix metalloproteinases (MMPs) -9 and -2 in murine macrophages. Molecular docking was initially performed to assess the binding affinity of the isolated flavonoids to the active site of MMP-9. Results and discussion: In silico model was a powerful tool to predict the compounds that could strongly bind and inhibit MMPs. CC extract and tea have shown to possess a significant antioxidant and anti-inflammatory effect, which can partially explain their traditional medicinal use.

Plant Antimicrobial Peptides: State of the Art, In Silico Prediction and Perspectives in the Omics Era

Even before the perception or interaction with pathogens, plants rely on constitutively guardian molecules, often specific to tissue or stage, with further expression after contact with the pathogen. These guardians include small molecules as antimicrobial peptides (AMPs), generally cysteine-rich, functioning to prevent pathogen establishment. Some of these AMPs are shared among eukaryotes (eg, defensins and cyclotides), others are plant specific (eg, snakins), while some are specific to certain plant families (such as heveins). When compared with other organisms, plants tend to present a higher amount of AMP isoforms due to gene duplications or polyploidy, an occurrence possibly also associated with the sessile habit of plants, which prevents them from evading biotic and environmental stresses. Therefore, plants arise as a rich resource for new AMPs. As these molecules are difficult to retrieve from databases using simple sequence alignments, a description of their characteristics and in silico (bioinformatics) approaches used to retrieve them is provided, considering resources and databases available. The possibilities and applications based on tools versus database approaches are considerable and have been so far underestimated.

Recent Advances and Computational Approaches in Peptide Drug Discovery

BACKGROUND

Drug design and development is a vast field that requires huge investment along with a long duration for providing approval to suitable drug candidates. With the advancement in the field of genomics, the information about druggable targets is being updated at a fast rate which is helpful in finding a cure for various diseases.

METHODS

There are certain biochemicals as well as physiological advantages of using peptide-based therapeutics. Additionally, the limitations of peptide-based drugs can be overcome by modulating the properties of peptide molecules through various biomolecular engineering techniques. Recent advances in computational approaches have been helpful in studying the effect of peptide drugs on the biomolecular targets. Receptor - ligand-based molecular docking studies have made it easy to screen compatible inhibitors against a target.Furthermore, there are simulation tools available to evaluate stability of complexes at the molecular level. Machine learning methods have added a new edge by enabling accurate prediction of therapeutic peptides.

RESULTS

Peptide-based drugs are expected to take over many popular drugs in the near future due to their biosafety, lower off-target binding chances and multifunctional properties.

CONCLUSION

This article summarises the latest developments in the field of peptide-based therapeutics related to their usage, tools, and databases.

Improving Protein-Peptide Docking Results via Pose-Clustering and Rescoring with a Combined Knowledge-Based and MM-GBSA Scoring Function

Protein-peptide docking, which predicts the complex structure between a protein and a peptide, is a valuable computational tool in peptide therapeutics development and the mechanistic investigation of peptides involved in cellular processes. Although current peptide docking approaches are often able to sample near-native peptide binding modes, correctly identifying those near-native modes from decoys is still challenging because of the extremely high complexity of the peptide binding energy landscape. In this study, we have developed an efficient postdocking rescoring protocol using a combined scoring function of knowledge-based ITScorePP potentials and physics-based MM-GBSA energies. Tested on five benchmark/docking test sets, our postdocking strategy showed an overall significantly better performance in binding mode prediction and score-rmsd correlation than original docking approaches. Specifically, our postdocking protocol outperformed original docking approaches with success rates of 15.8 versus 10.5% for pepATTRACT on the Global_57 benchmark, 5.3 versus 5.3% for CABS-dock on the Global_57 benchmark, 17.0 versus 11.3% for FlexPepDock on the LEADS-PEP data set, 40.3 versus 33.9% for HPEPDOCK on the Local_62 benchmark, and 64.2 versus 52.8% for HPEPDOCK on the LEADS-PEP data set when the top prediction was considered. These results demonstrated the efficacy and robustness of our postdocking protocol.

Highly Flexible Ligand Docking: Benchmarking of the DockThor Program on the LEADS-PEP Protein-Peptide Data Set

Protein-peptide interactions play a crucial role in many cellular and biological functions, which justify the increasing interest in the development of peptide-based drugs. However, predicting experimental binding modes and affinities in protein-peptide docking remains a great challenge for most docking programs due to some particularities of this class of ligands, such as the high degree of flexibility. In this paper, we present the performance of the DockThor program on the LEADS-PEP data set, a benchmarking set composed of 53 diverse protein-peptide complexes with peptides ranging from 3 to 12 residues and with up to 51 rotatable bonds. The DockThor performance for pose prediction on redocking studies was compared with some state-of-the-art docking programs that were also evaluated on the LEADS-PEP data set, AutoDock, AutoDock Vina, Surflex, GOLD, Glide, rDock, and DINC, as well as with the task-specific docking protocol HPepDock. Our results indicate that DockThor could dock 40% of the cases with an overall backbone RMSD below 2.5 Å when the top-scored docking pose was considered, exhibiting similar results to Glide and outperforming other protein-ligand docking programs, whereas rDock and HPepDock achieved superior results. Assessing the docking poses closest to the crystal structure (i.e., best-RMSD pose), DockThor achieved a success rate of 60% in pose prediction. Due to the great overall performance of handling peptidic compounds, the DockThor program can be considered as suitable for docking highly flexible and challenging ligands, with up to 40 rotatable bonds. DockThor is freely available as a virtual screening Web server at https://www.dockthor.lncc.br/ .

Hierarchical Flexible Peptide Docking by Conformer Generation and Ensemble Docking of Peptides

Given the importance of peptide-mediated protein interactions in cellular processes, protein-peptide docking has received increasing attention. Here, we have developed a Hierarchical flexible Peptide Docking approach through fast generation and ensemble docking of peptide conformations, which is referred to as HPepDock. Tested on the LEADS-PEP benchmark data set of 53 diverse complexes with peptides of 3-12 residues, HPepDock performed significantly better than the 11 docking protocols of five small-molecule docking programs (DOCK, AutoDock, AutoDock Vina, Surflex, and GOLD) in predicting near-native binding conformations. HPepDock was also evaluated on the 19 bound/unbound and 10 unbound/unbound protein-peptide complexes of the Glide SP-PEP benchmark and showed an overall better performance than Glide SP-PEP+MM-GBSA and FlexPepDock in both bound and unbound docking. HPepDock is computationally efficient, and the average running time for docking a peptide is ∼15 min with the range from about 1 min for short peptides to around 40 min for long peptides.

Computer-aided design of amino acid-based therapeutics: a review

During the last two decades, the pharmaceutical industry has progressed from detecting small molecules to designing biologic-based therapeutics. Amino acid-based drugs are a group of biologic-based therapeutics that can effectively combat the diseases caused by drug resistance or molecular deficiency. Computational techniques play a key role to design and develop the amino acid-based therapeutics such as proteins, peptides and peptidomimetics. In this study, it was attempted to discuss the various elements for computational design of amino acid-based therapeutics. Protein design seeks to identify the properties of amino acid sequences that fold to predetermined structures with desirable structural and functional characteristics. Peptide drugs occupy a middle space between proteins and small molecules and it is hoped that they can target "undruggable" intracellular protein-protein interactions. Peptidomimetics, the compounds that mimic the biologic characteristics of peptides, present refined pharmacokinetic properties compared to the original peptides. Here, the elaborated techniques that are developed to characterize the amino acid sequences consistent with a specific structure and allow protein design are discussed. Moreover, the key principles and recent advances in currently introduced computational techniques for rational peptide design are spotlighted. The most advanced computational techniques developed to design novel peptidomimetics are also summarized.

Structure-based Methods for Binding Mode and Binding Affinity Prediction for Peptide-MHC Complexes

Understanding the mechanisms involved in the activation of an immune response is essential to many fields in human health, including vaccine development and personalized cancer immunotherapy. A central step in the activation of the adaptive immune response is the recognition, by T-cell lymphocytes, of peptides displayed by a special type of receptor known as Major Histocompatibility Complex (MHC). Considering the key role of MHC receptors in T-cell activation, the computational prediction of peptide binding to MHC has been an important goal for many immunological applications. Sequence- based methods have become the gold standard for peptide-MHC binding affinity prediction, but structure-based methods are expected to provide more general predictions (i.e., predictions applicable to all types of MHC receptors). In addition, structural modeling of peptide-MHC complexes has the potential to uncover yet unknown drivers of T-cell activation, thus allowing for the development of better and safer therapies. In this review, we discuss the use of computational methods for the structural modeling of peptide-MHC complexes (i.e., binding mode prediction) and for the structure-based prediction of binding affinity.

Structure-Based Design of Non-natural Macrocyclic Peptides That Inhibit Protein-Protein Interactions

Macrocyclic peptides can interfere with challenging biomolecular targets including protein–protein interactions. Whereas there are various approaches that facilitate the identification of peptide-derived ligands, their evolution into higher affinity binders remains a major hurdle. We report a virtual screen based on molecular docking that allows the affinity maturation of macrocyclic peptides taking non-natural amino acids into consideration. These macrocycles bear large and flexible substituents that usually complicate the use of docking approaches. A virtual library containing more than 1400 structures was screened against the target focusing on docking poses with the core structure resembling a known bioactive conformation. Based on this screen, a macrocyclic peptide 22 involving two non-natural amino acids was evolved showing increased target affinity and biological activity. Predicted binding modes were verified by X-ray crystallography. The presented workflow allows the screening of large macrocyclic peptides with diverse modifications thereby expanding the accessible chemical space and reducing synthetic efforts.

LEADS-PEP: A Benchmark Data Set for Assessment of Peptide Docking Performance

With increasing interest in peptide-based therapeutics also the application of computational approaches such as peptide docking has gained more and more attention. In order to assess the suitability of docking programs for peptide placement and to support the development of peptide-specific docking tools, an independently constructed benchmark data set is urgently needed. Here we present the LEADS-PEP benchmark data set for assessing peptide docking performance. Using a rational and unbiased workflow, 53 protein-peptide complexes with peptide lengths ranging from 3 to 12 residues were selected. The data set is publicly accessible at www.leads-x.org . In a second step we evaluated several small molecule docking programs for their potential to reproduce peptide conformations as present in LEADS-PEP. While most tested programs were capable to generate native-like binding modes of small peptides, only Surflex-Dock and AutoDock Vina performed reasonably well for peptides consisting of more than five residues. Rescoring of docking poses with scoring functions ChemPLP, ChemScore, and ASP further increased the number of top-ranked near-native conformations. Our results suggest that small molecule docking programs are a good and fast alternative to specialized peptide docking programs.

Computational Modeling of Small Molecule Ligand Binding Interactions and Affinities

Understanding and controlling biological phenomena via structure-based drug screening efforts often critically rely on accurate description of protein-ligand interactions. However, most of the currently available computational techniques are affected by severe deficiencies in both protein and ligand conformational sampling as well as in the scoring of the obtained docking solutions. To overcome these limitations, we have recently developed MedusaDock, a novel docking methodology, which simultaneously models ligand and receptor flexibility. Coupled with MedusaScore, a physical force field-based scoring function that accounts for the protein-ligand interaction energy, MedusaDock, has reported the highest success rate in the CSAR 2011 exercise. Here, we present a standard computational protocol to evaluate the binding properties of the two enantiomers of the non-selective β-blocker propanolol in the β2 adrenergic receptor's binding site. We describe details of our protocol, which have been successfully applied to several other targets.

Rational, computer-enabled peptide drug design: principles, methods, applications and future directions

Peptides provide promising templates for developing drugs to occupy a middle space between small molecules and antibodies and for targeting 'undruggable' intracellular protein-protein interactions. Importantly, rational or in cerebro design, especially when coupled with validated in silico tools, can be used to efficiently explore chemical space and identify islands of 'drug-like' peptides to satisfy diverse drug discovery program objectives. Here, we consider the underlying principles of and recent advances in rational, computer-enabled peptide drug design. In particular, we consider the impact of basic physicochemical properties, potency and ADME/Tox opportunities and challenges, and recently developed computational tools for enabling rational peptide drug design. Key principles and practices are spotlighted by recent case studies. We close with a hypothetical future case study.

Modeling of protein-peptide interactions using the CABS-dock web server for binding site search and flexible docking

Protein-peptide interactions play essential functional roles in living organisms and their structural characterization is a hot subject of current experimental and theoretical research. Computational modeling of the structure of protein-peptide interactions is usually divided into two stages: prediction of the binding site at a protein receptor surface, and then docking (and modeling) the peptide structure into the known binding site. This paper presents a comprehensive CABS-dock method for the simultaneous search of binding sites and flexible protein-peptide docking, available as a user's friendly web server. We present example CABS-dock results obtained in the default CABS-dock mode and using its advanced options that enable the user to increase the range of flexibility for chosen receptor fragments or to exclude user-selected binding modes from docking search. Furthermore, we demonstrate a strategy to improve CABS-dock performance by assessing the quality of models with classical molecular dynamics. Finally, we discuss the promising extensions and applications of the CABS-dock method and provide a tutorial appendix for the convenient analysis and visualization of CABS-dock results. The CABS-dock web server is freely available at http://biocomp.chem.uw.edu.pl/CABSdock/.

Docking small peptides remains a great challenge: an assessment using AutoDock Vina

There is a growing interest in the mechanisms and the prediction of how flexible peptides bind proteins, often in a highly selective and conserved manner. While both existing small-molecule docking methods and custom protocols can be used, even short peptides make difficult targets owing to their high torsional flexibility. Any benchmarking should therefore start with those. We compiled a meta-data set of 47 complexes with peptides up to five residues, based on 11 related studies from the past decade. Although their highly varying strategies and constraints preclude direct, quantitative comparisons, we still provide a comprehensive overview of the reported results, using a simple yet stringent measure: the quality of the top-scoring peptide pose. Using the entire data set, this is augmented by our own benchmark of AutoDock Vina, a freely available, fast and widely used docking tool. It particularly addresses non-expert users and was therefore implemented in a highly integrated manner. Guidelines addressing important issues such as the amount of sampling required for result reproducibility are so far lacking. Using peptide docking as an example, this is the first study to address these issues in detail. Finally, to encourage further, standardized benchmarking efforts, the compiled data set is made available in an accessible, transparent and extendable manner.

Subangstrom accuracy in pHLA-I modeling by Rosetta FlexPepDock refinement protocol

Flexible peptides binding to human leukocyte antigen (HLA) play a key role in mediating human immune responses and are also involved in idiosyncratic adverse drug reactions according to recent research. However, the structural determinations of pHLA complexes remain challenging under the present conditions. In this paper, the performance of a new peptide docking method, namely FlexPepDock, was systematically investigated by a benchmark of 30 crystallized structures of peptide-HLA class I complexes. The docking results showed that the near-native pHLA-I models with peptide bb-RMSD less than 2 Å were ranked in the top 1 model for 100% (70/70) docking cases, and the subangstrom models with peptide bb-RMSD less than 1 Å were ranked in the top 5 lowest-energy models for 65.7% (46/70) docking cases. Furthermore, 10 out of 70 docking cases ranked the subangstrom all-atom models in the top 5 lowest-energy models. The results showed that the FlexPepDock can generate high-quality models of pHLA-I complexes and can be widely applied to pHLA-I modeling and mechanism research of peptide-mediated immune responses.

PaFlexPepDock: parallel ab-initio docking of peptides onto their receptors with full flexibility based on Rosetta

Structural information related to protein–peptide complexes can be very useful for novel drug discovery and design. The computational docking of protein and peptide can supplement the structural information available on protein–peptide interactions explored by experimental ways. Protein–peptide docking of this paper can be described as three processes that occur in parallel: ab-initio peptide folding, peptide docking with its receptor, and refinement of some flexible areas of the receptor as the peptide is approaching. Several existing methods have been used to sample the degrees of freedom in the three processes, which are usually triggered in an organized sequential scheme. In this paper, we proposed a parallel approach that combines all the three processes during the docking of a folding peptide with a flexible receptor. This approach mimics the actual protein–peptide docking process in parallel way, and is expected to deliver better performance than sequential approaches. We used 22 unbound protein–peptide docking examples to evaluate our method. Our analysis of the results showed that the explicit refinement of the flexible areas of the receptor facilitated more accurate modeling of the interfaces of the complexes, while combining all of the moves in parallel helped the constructing of energy funnels for predictions.

Optimized atomic statistical potentials: assessment of protein interfaces and loops

MOTIVATION

Statistical potentials have been widely used for modeling whole proteins and their parts (e.g. sidechains and loops) as well as interactions between proteins, nucleic acids and small molecules. Here, we formulate the statistical potentials entirely within a statistical framework, avoiding questionable statistical mechanical assumptions and approximations, including a definition of the reference state.

RESULTS

We derive a general Bayesian framework for inferring statistically optimized atomic potentials (SOAP) in which the reference state is replaced with data-driven 'recovery' functions. Moreover, we restrain the relative orientation between two covalent bonds instead of a simple distance between two atoms, in an effort to capture orientation-dependent interactions such as hydrogen bonds. To demonstrate this general approach, we computed statistical potentials for protein-protein docking (SOAP-PP) and loop modeling (SOAP-Loop). For docking, a near-native model is within the top 10 scoring models in 40% of the PatchDock benchmark cases, compared with 23 and 27% for the state-of-the-art ZDOCK and FireDock scoring functions, respectively. Similarly, for modeling 12-residue loops in the PLOP benchmark, the average main-chain root mean square deviation of the best scored conformations by SOAP-Loop is 1.5 Å, close to the average root mean square deviation of the best sampled conformations (1.2 Å) and significantly better than that selected by Rosetta (2.1 Å), DFIRE (2.3 Å), DOPE (2.5 Å) and PLOP scoring functions (3.0 Å). Our Bayesian framework may also result in more accurate statistical potentials for additional modeling applications, thus affording better leverage of the experimentally determined protein structures.

AVAILABILITY AND IMPLEMENTATION

SOAP-PP and SOAP-Loop are available as part of MODELLER (http://salilab.org/modeller).