Issue #112: PlantP450Dock: an Automated Molecular Docking Pipeline of Plant Cytochrome P450s
Protein Design Digest #112: PlantP450Dock: an Automated Molecular Docking Pipeline of Plant Cytochro…

Building something in Protein Design?
I love collaborating on new challenges. Let's build together.
Subscribe to Protein Design Digest
Daily curated signals from arXiv, PubMed, and BioRxiv.
Signal of the Day
PlantP450Dock: an Automated Molecular Docking Pipeline of Plant Cytochrome P450s
Cytochrome P450 enzymes (CYPs) are the primary drivers of chemical diversity in plant secondary metabolism, yet fewer than 10% of plant P450s have been functionally characterized. Computational docking offers a scalable approach to prioritize candidates for experimental validation, but existing workflows are ill-suited for plant P450s due to the absence of the heme cofactor in AlphaFold-predicted structures and the lack of objective criteria for flexible residue selection. Here we present PlantP450Dock, an automated pipeline that integrates heme implantation, molecular dynamics-based conformational sampling, data-driven flexible residue selection, and semi-flexible docking into a single streamlined workflow. The heme cofactor is transferred from a crystallographic reference template to the AlphaFold model via a local coordinate transformation algorithm, yielding a positional deviation of less than 0.2 Å relative to the experimentally determined structure of CYP73A33 (PDB: 6VBY). A 100 ns molecular dynamics simulation confirmed stable Fe–S coordination geometry throughout (2.61 ± 0.08 Å), and a singular value decomposition-based heme plane filtering strategy objectively identified active-site flexible residues without operator input. Cross-family validation across four phylogenetically distinct P450s belonging to the CYP73, CYP711, CYP706, and CYP701 families produced catalytically competent binding poses with substrate-to-iron distances of 2.8–4.4 Å without any enzyme-specific parameter adjustment. PlantP450Dock will be made freely accessible as a web server, providing the community with a standardized and reproducible computational framework to accelerate the functional annotation of the largely uncharacterized plant P450 superfamily.
Why this matters: Provides actionable mutations to enhance catalytic efficiency or thermostability.
Also Worth Reading
Exploring the mechanism of saffron in treating viral myocarditis using network pharmacology and molecular docking.
Viral myocarditis (VM) is a cardiovascular disorder that can lead to heart failure and cardiogenic shock. Saffron, a traditional Chinese medicinal herb, has shown therapeutic potential against VM in numerous studies. However, the mechanisms through which saffron exerts its effects on VM remain poorly understood. Thus, this study aimed to elucidate the active compounds, molecular targets, and signaling pathways involved in saffron’s therapeutic action against VM by employing network pharmacology and molecular docking approaches. The active compounds and corresponding targets of saffron were retrieved from the Traditional Chinese Medicine Systems Pharmacology database. VM-associated targets were sourced from the GeneCards database. Overlapping targets between saffron and VM were then identified. Protein-protein interaction networks were established and analyzed utilizing the STRING platform and Cytoscape software to determine core targets. Furthermore, gene ontology and Kyoto encyclopedia of genes and genomes enrichment analyses were carried out utilizing Bioconductor in R to explore the potential biological activities and signaling pathways through which saffron may act against VM. Finally, molecular docking and model visualization were carried out using AutoDock Tools and PyMOL open-source software. From the database, we identified 4 active compounds in saffron with potential effects against VM: crocetin, isorhamnetin, kaempferol, and quercetin. A total of 60 corresponding targets were observed, with TNF, IL-6, IL-1β, CXCL8, and JUN emerging as core targets. Kyoto encyclopedia of genes and genomes enrichment analysis revealed 155 regulatory signaling pathways, among which the TNF, AGE-RAGE, and IL-17 signaling pathways, lipid metabolism, and atherosclerosis were the most prominent. Molecular docking results indicated that quercetin showed the strongest binding affinity toward IL-1β and CXCL8. The therapeutic effect of saffron against VM is not driven by a single factor, but rather involves multiple active compounds, targets, and signaling pathways.
Unraveling the anti-neuroinflammatory mechanisms of Cervus cucumis polypeptide injection in Alzheimer’s disease: insights from network pharmacology, molecular docking, molecular dynamics simulation, and experimental validation.
Objective Alzheimer’s disease (AD) is a progressive neurodegenerative disorder with increasing global prevalence, in which neuroinflammation serves as a critical pathological driver exacerbating cognitive decline. While current therapies offer limited symptomatic relief, multi-target strategies are urgently needed. Cervus cucumis polypeptide injection (CCPI), a traditional Chinese medicine (TCM) formulation, has demonstrated anti-inflammatory properties; however, its mechanisms of action against AD remain unclear. This study aimed to elucidate the anti-AD potential mechanisms of CCPI using an integrated approach combining network pharmacology, molecular docking, molecular dynamics (MD) simulation, and experimental validation. Methods Active components and corresponding targets of CCPI were retrieved from the TCMSP database, while AD-related targets were collected from Genecards, OMIM, and DrugBank. Potential therapeutic targets were identified by intersecting drug and disease targets, followed by protein-protein interaction (PPI) network construction, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Molecular docking and MD simulations were performed to evaluate interactions between potential active components and key targets. In vitro experiments were conducted on Aβ 25-35 -induced BV2 microglial cells to assess cell viability (CCK-8 assay), inflammatory cytokine levels (ELISA), and protein expression (Western blot) related to the neuroinflammation pathway and microglial polarization. Results A total of 28 active components and 50 common targets of CCPI for AD treatment were identified. Linoleic acid (LA) was determined to be a potential active component, with IL-6 as the key target based on PPI network topology. Molecular docking and MD simulation confirmed a stable binding affinity between LA and IL-6. KEGG analysis revealed significant enrichment in the HIF-1 signaling pathway, particularly the IL-6/STAT3/VEGF signaling pathway. In vitro , CCPI treatment significantly enhanced cell viability and attenuated the pro-inflammatory response, as evidenced by reduced levels of IL-6, IL-1β, and TNF-α, decreased the expression of the pro-inflammatory marker iNOS. Concurrently, it elevated the expression of the anti-inflammatory/repair-associated marker CD206. Western blot analysis further verified that CCPI suppressed IL-6/STAT3 activation while upregulating VEGF expression. Additionally, LA alone significantly reduced IL-6 levels and STAT3 phosphorylation, decreased the expression of iNOS, and increased the expression of CD206, with therapeutic efficacy comparable to CCPI. Conclusion CCPI exerts neuroprotective effects in AD models by regulating the IL-6/STAT3/VEGF pathway, downregulating the expression of the inflammation-related iNOS protein, upregulating the expression of the CD206 protein associated with anti-inflammatory and reparative functions, remodeling the functional state of microglia, inhibiting their pro-inflammatory responses, and enhancing their reparative functions. Its potential active component, LA, likely mediates this effect by stably binding to and inhibiting IL-6, thus suppressing the downstream STAT3 phosphorylation that drives inflammatory activation.
A multimodal approach integrating spectroscopy, deep learning guided molecular docking, and molecular dynamics simulation for predictive assessment of pioglitazone to albumin binding for formulation development.
Binding affinity is a critical parameter that can influence the state of the drug in vivo and help to define the formulation strategy. The current study implements a multimodal approach to analyse the binding affinity between human serum albumin (HSA) and pioglitazone. Ultraviolet (UV) absorbance and fluorescence spectrometry analyses were performed on different combinations of HSA and pioglitazone complexes, and the absorbance and fluorescence intensities were mapped to calculate the binding constant. DynamicBind, a distinct deep-learning artificial intelligence tool, was implemented to perform in silico docking studies using a non-conventional approach. Furthermore, molecular dynamics simulation was also performed to generate root mean square deviation, radius of gyration, and root mean square fluctuation values, followed by principal component analysis, probability distribution function, and free energy landscape analysis. The simulation output was analysed to interpret the binding affinity and associated conformation of the protein-active pharmaceutical ingredient (API) complex. The binding constant calculated through UV analysis was 1.1 × 10 4 M -1 . Fluorescence spectroscopic analysis derived a value of 1.7 × 10 5 M -1 . At the same time, DynamicBind predicted the cLDDT score for the top predicted model to be 0.634, and a binding affinity value of greater than 5, indicating a relatively moderate binding between pioglitazone and HSA. The results from molecular dynamics simulations further complemented our earlier observations, indicating non-covalent binding interactions and a stable protein-API complex, which is desirable for developing a formulation using HSA as a carrier polymer. This orthogonal approach also provided critical information on the fate of the API and possible considerations that needed to be made during the design of the formulation process, highlighting the need for similar approaches that could provide multifaceted advantages and help in optimising R&D costs and timelines.
Research & AI Updates
- UK-based open data resource ‘behind £12bn productivity gain’ - Research Professional News — UK-based open data resource ‘behind £12bn productivity gain’ Research Professional News.
- Google DeepMind and Edison Are Building the AI Scientist - Genetic Engineering and Biotechnology News — Google DeepMind and Edison Are Building the AI Scientist Genetic Engineering and Biotechnology News.
- OpenBind releases first open dataset and AI model for drug discovery - Technology Org — OpenBind releases first open dataset and AI model for drug discovery Technology Org.
- Clarence resident recognized for structural biological research - Bee Group Newspapers — Clarence resident recognized for structural biological research Bee Group Newspapers.
From the Industry
- FDA in Flux — May 2026 Newsletter | Mintz - Health Care Viewpoints - JD Supra — FDA in Flux — May 2026 Newsletter | Mintz - Health Care Viewpoints JD Supra.
- A migraine drugmaker heads to the public markets through reverse merger - BioPharma Dive — A migraine drugmaker heads to the public markets through reverse merger BioPharma Dive.
- Spine fusion protein from Bone Biologics now lasts 29 months - Stock Titan — Spine fusion protein from Bone Biologics now lasts 29 months Stock Titan.
- GenSight Biologics to Hold its Combined General Meeting on May 19, 2026 - BioSpace — GenSight Biologics to Hold its Combined General Meeting on May 19, 2026 BioSpace.
- Regeneron pens $2.3B Parabilis pact to develop new twist on ADCs - Fierce Biotech — Regeneron pens $2.3B Parabilis pact to develop new twist on ADCs Fierce Biotech.
- Is Stacking Biologics Safe if You Have Crohn’s Disease? - Everyday Health — Is Stacking Biologics Safe if You Have Crohn’s Disease? Everyday Health.
- Sen. David McCormick tours AI-powered biotech labs at Penn to promote AI and federal funding admit NSF turmoil - WHYY — Sen.
Quick Reads
Molecular docking and molecular dynamics simulations of compounds targeting GABAA receptor with potential relevance to anesthesia.
The gamma-aminobutyric acid type A (GABA A ) receptor is the primary mediator of inhibitory neurotransmission in the central nervous system and represents an important pharmacological target for sedative and anesthetic agents. Read more →
Hybrid Computational Framework Integrating Ensemble Learning, Molecular Docking, and Dynamics for Predicting Antimalarial Efficacy of Malaria Box Compounds.
The emergence of drug-resistant strains of Plasmodium falciparum continues to challenge global malaria control efforts, underscoring the urgent need for novel therapeutic strategies. Read more →
Accelerating Drug Discovery with HyperLab: An Easy-to-Use AI-Driven Platform
HyperLab, developed by HITS, is a web-based, AI-driven drug discovery platform designed to increase research efficiency for experimental drug discovery researchers. Read more →
Structural bias in machine learning-guided peptide design
ABSTRACT Machine learning continues to accelerate peptide and protein design through the rapid prediction and generation of sequences with desired characteristics. Read more →
An Interpretable and Robust Multi-Parameter Prioritization Framework for BACE1 Inhibitors Integrating Meta-Ensemble QSAR, Protein Language Model–Guided Residue Weighting, and Sensitivity-Validated Ranking
Alzheimer’s disease remains a major therapeutic challenge, and no β-secretase (BACE1) inhibitor has achieved clinical approval. Read more →
Luteolin as a Multitarget Agent against Colorectal Cancer Explored Through Integrative Network Pharmacology, Molecular Docking, and Dynamics Simulations.
Background Colorectal cancer (CRC) is a leading cause of cancer-related death worldwide. Read more →
Recent advances in the structural biology of photosystem II.
Photosystem II (PSII) is the membrane protein-pigment complex responsible for the light-driven oxidation of water to molecular oxygen, a reaction that enables the aerobic biosphere and powers biological carbon fixation. Read more →
AI and pharmacophore-based screening, molecular docking, and dynamic simulation for identification of CDK6 inhibitors to combat NSCLC.
The overactivated CDK6/Cyclin D3 complex is considered to be associated with poor prognosis in patients with NSCLC, thus developing CDK6/Cyclin D3 inhibitors is expected to provide new options for NSCLC patients. Read more →
Pipeline Tip
Index your BigWig files before visualization to save memory.
Resources & Tools
- Dataset: MGnify - Metagenomics resource for microbiome sequence data.
- Dataset: PDBbind - Binding affinity data with 3D structures of protein-ligand complexes.
- Tool: Boltz-1 - Open-source biomolecular structure prediction model. View all tools →
- Tool: ProteinSolver - Graph-based neural network for protein sequence design. View all tools →
- Event: Protein Design Hub (LinkedIn Group) (Ongoing)
- Event: Structural Biology Events (Open)
- Job: Are non-antibiotic drugs contributing to antimicrobial resistance? - Nature at Nature Careers
- Job: What the AI era doctor should know: a scoping review of proposed artificial intelligence competencies for medical education - Nature at Nature Careers
The protein structure is the language of life; design is its poetry. — Recep Adiyaman