Protein Cavity Research Articles

Monitoring Dynamic Conformations of a Single Fluorescent Molecule Inside a Protein Cavity.

Fluorescence nanoscopy and single-molecule methods are entering the realm of structural biology, breaking new ground for dynamic structural measurements at room temperature and liquid environments. Here, single-molecule localization microscopy, polarization-dependent single-molecule excitation, and protein engineering are combined to determine the orientation of a fluorophore forming hydrogen bonds inside a protein cavity. The observed conformations are in good agreement with molecular dynamics simulations, enabling a new, more realistic interplay between experiments and simulations to identify stable conformations and the key interactions involved. Furthermore, jumps between conformations can be monitored with a precision of 3° and a time resolution of a few seconds, confirming the potential of this methodology for retrieving dynamic structural information of nanoscopic biological systems under physiologically compatible conditions.

Investigation into the effect of phenylalanine gating on anaerobic haem breakdown using the energy landscape approach.

We have recently demonstrated a novel anaerobic NADH-dependent haem breakdown reaction, which is carried out by a range of haemoproteins. The Yersinia enterocolitica protein, HemS, is the focus of further research presented in the current paper. Using conventional experimental methods, bioinformatics, and energy landscape theory (ELT), we provide new insight into the mechanism of the novel breakdown process. Of particular interest is the behavior of a double phenylalanine gate, which opens and closes according to the relative situations of haem and NADH within the protein pocket. This behavior suggests that the double phe-gate fulfills a regulatory role within the pocket, controlling the access of NADH to haem. Additionally, stopped-flow spectroscopy results provide kinetic comparisons between the wild-type and the selected mutants. We also present a fully resolved crystal structure for the F104AF199A HemS monomer, including its extensive loop, the first such structure to be completely resolved for HemS or any of its close homologues. The energy landscapes approach provided key information regarding the gating strategy employed by HemS, compensating for current limitations with conventional biophysical and molecular dynamics approaches. We propose that ELT become more widely used in the field, particularly in the investigation of the dynamics and interactions of weak-binding ligands, and for gating features, within protein cavities.

Investigation of the Structure and Functional Activity of the YqeK Protein in Streptococcus pyogenes with High Efficiency in Hydrolyzing Ap4A

Streptococcus pyogenes is an important zoonotic Gram-positive bacterium that appears in chains, without spores or flagella, and belongs to the beta-hemolytic streptococci. It can be transmitted through droplets or contact, with the preferred antibiotics being penicillin, erythromycin, or cephalosporins. However, the misuse of these drugs has led to antibiotic resistance, posing a significant threat to both human and animal health. Studying resistance genes encoding proteins is crucial for mitigating the emergence of resistant strains and improving treatment outcomes. Interestingly, a dinucleotide known as diadenosine tetraphosphate (Ap4A) exists in Streptococcus pyogenes; its accumulation in response to various stress signals can inhibit bacterial pathogenicity and enhance antibiotic susceptibility. Our research focuses on the Sp-yqeK protein, which we have identified as a hydrolase that symmetrically cleaves Ap4A. The Sp-yqeK protein effectively cleaves Ap4A, producing adenosine diphosphate (ADP) molecules. Results indicate that this enzyme exhibits optimal activity at pH 7.0 and a temperature of 45 °C. Furthermore, we determined the crystal structure of the Sp-yqeK, Mg2+, and ADP complex at a resolution of 2.0 Å, providing insights into the interactions crucial for catalytic efficiency between Sp-yqeK and ADP. This complex reveals unique folding characteristics of the HD domain superfamily proteins, accommodating both ADP and Mg2+. These components are securely embedded into the polar cavity of the yqeK protein through conserved residues (His29, Lys62, His91, His117, Asp135, Leu172, Phe180, and Thr183), highlighting the residues responsible for Ap4A hydrolysis and Mg2+ binding. Our research offers a deeper understanding of the hydrolysis mechanism of Ap4A and the specificity of Sp-yqeK, providing structural insights that may support future studies on antibiotic resistance in Streptococcus pyogenes and other Gram-positive bacteria.

In Silico Strategies for Characterizing Inner Cavities of Lipid-Binding Proteins.

Cavities in proteins perform diverse functions such as substrate binding, enzyme catalysis, passage for transportation of small molecules, and protein oligomerization. Often, the physical properties of these cavities are closely linked to the protein function; such as the hydrophobic lipid-binding cavities in lipid-binding proteins (LBPs) that protect lipid substrates from the larger aqueous milieu. Therefore, the characterization of protein cavities can provide valuable insights into protein structure-function relationships, hinting toward their mechanism of action while aiding in the identification of ligand binding sites that are essential for drug discovery approaches. Several algorithms have historically been designed to identify and characterize the different types of cavities in protein structures. We summarize these algorithms and provide a step-by-step guide for locating and characterizing internal cavities in proteins using CICLOP by using ATP-binding cassette transporter A1 (ABCA1) as an example.

Computational Insights for Interactions between nsP2 and nsP3 of CHIKV and Hormones through DFT Computations and Molecular Dynamics Simulations.

The non-structural protein (nsP2 & nsP3) of the Chikungunya virus (CHIKV) is responsible for the transmission of viral infection. The main role of non-structural proteins are involved in the transcription process at an early stage of the infection. In this work, authors have studied the impact of nsP2 and nsP3 of CHIKV on hormones present in the human body using a computational approach. The ten hormones of chemical properties such as 4-Androsterone-2,17-dione, aldosterone, androsterone, corticosterone, cortisol, cortisone, estradiol, estrone, progesterone and testosterone were taken as a potency. From the molecular docking, the binding energy of the complexes is estimated, and cortisone was found to be the highest negative binding energy (-6.57 kcal/mol) with the nsP2 and corticosterone with the nsP3 (-6.47 kcal/mol). This is based on the interactions between hormones and nsP2/nsP3, which are types of noncovalent intermolecular interactions categorized into three types: electrostatic interactions, van der Waals (vdW) interactions, and hydrogen-bonding (H-bonding) interactions. To validate the docking results, additional molecular dynamics simulations and MM-GBSA methods were performed. The change in enthalpy, entropy, and free energy were calculated using MM-GBSA methods. The nsP2 and nsP3 of CHIKV interact strongly with the cortisone and corticosterone with free energy changes of -20.55 & -36.08 kcal/mol, respectively. Methods: The crystal structures of 3TKR and 3GPO proteins of nsP2 and nsP3 were extracted from the RCSB Protein Data Bank. Initially, unnecessary atoms like extra cations or anions and missing explicit hydrogen atoms were removed and added from the native domain of nsP2 and nsP3. The alignment of coordinated in the native domain was performed using Chimera and Notepad++ tools. The molecular docking of protein and ligand was performed usingAutoDock tool; it is essential for the prediction of the orientation of the ligand into the cavity of the target protein based on binding affinity. Based on thermodynamic parameters, MD Simulations were employed to calculate the change in binding free energies of various complexes followed by a change in enthalpy and entropy with time. According to MD production, the CPPTAJ and PTRAJ programs were used to analyse the trajectories, such as dynamic stability (RMSD), residual fluctuation (RMSF), compatibility, and hydrogen bonds of the newly formed complexes. After that, the Density Functional Theory (DFT) were used to calculate the electronic properties of selected molecules by Gaussian 16 on applying the B3LYP method with the 6-311G (d, p) basis set.

Lapachol, a natural food component, interacts with human serum albumin: Insights of its impact on the pharmacokinetics of clinically used drugs

Lapachol (LAP), a natural 1,4-naphthoquinone used in popular medicine in South America, is an antioxidant and antimicrobial compound in teas and infusions and used as a food additive; however, its interactive profile with the main protein carrier of compounds in the human bloodstream (human serum albumin, HSA) was not still characterized. Additionally, the impact of LAP in binding clinically drugs to albumin is still unknown. Thus, the present work describes the interaction HSA:LAP using different biophysical techniques, i.e., 1H saturation-transfer difference nuclear magnetic resonance (1H STD-NMR), isothermal titration calorimetry (ITC), steady-state and time-resolved fluorescence measurements combined with molecular docking calculations. LAP interacts with subdomain region IIA (site I), mainly driven by enthalpy effects, while subdomain region IB (site III) was identified as the second binding site, mainly driven by entropy effects. The binding is spontaneous, strong (binding constant average, Kaverage ≈ 4.45 × 105 M−1), and there is a positive cooperativity in the presence of ibuprofen, with the LAP structure fully buried into the protein cavities. Overall, LAP might impact the residence time (pharmacokinetic profile) of drugs that bind to subdomains regions IIA and IB of albumin, e.g., warfarin, phenylbutazone, diflunisal, naproxen, camptothecin, doxorubicin, daunorubicin, suramin, and tyrosine kinase inhibitors.

A comprehensive analysis of the spectroscopic and drug-like properties of dimethyl 5-hydroxybenzene-1,3-dicarboxylate: insights from DFT and MD simulations

ABSTRACT An in-depth investigation was undertaken to explore the ground-state molecular and electronic structure of Dimethyl 5-hydroxybenzene-1,3-dicarboxylate (D5HD). The potential energy scan of D5HD at various dihedral angles identified the lowest energy conformer, which was re-optimised at the higher basis set. The recorded FT-Raman and FT-IR spectra were compared with corresponding theoretical spectra, demonstrating a high degree of coherence. A comprehensive Natural Bond Orbital analysis was done for D5HD to gain insights for the overall electronic distribution within the molecule. Additionally, nonlinear optical (NLO) properties along with various DFT-based descriptors were calculated. Molecular docking studies against human acetylcholinesterase enzyme (PDB 2JF0), D5HD displayed good binding affinity. The study revealed that hydrogen bonded as well as van der Waals interactions exist between D5HD and different residues in the binding cavity of the target protein. The 100 ns MD simulation confirms the conformational stability of protein (2JF0)–ligand (D5HD) complex. The binding energy and residue decomposition analysis performed using MM-GBSA approach showed positive results with binding energy of −17.76 kcal/mol. Furthermore, D5HD effectively passed all pharmacokinetic filters, establishing its potential as a promising candidate in the quest for novel acetylcholinesterase inhibitors. Highlights Systematic exploration of the structure and spectroscopic properties of D5HD. Molecular docking between D5HD and acetylcholinesterase (PDB 2JF0) reveals the establishment of hydrogen bonds and Van der Waals interactions. D5HD successfully complies with all relevant pharmacokinetic filters. MD simulation analysis on 2JF0–D5HD complex confirms the conformational stability. Energy and residue decomposition based on MMGBSA approach.

Comprehensive Phytochemical Profiling, GC-MS Analysis, Molecular Docking and Antiproliferative Activity of Ethanol Fraction of <i>Tabernaemontana coronaria</i> and <i>Thunbergia alata</i>

Objective: Current study involves the phytochemical examination, GC MS testing of ethanol fraction of leaves of two plants Tabernaemontana coronaria and Thunbergia alata. Methods: Both ethanol fractions of selected plants were subjected to In vitro antiproliferative activity by employing MTT assay on A549 cell lines. Zebra Fish fin model and zebra fish embryo tests were employed to assess the fin regeneration and effect on angiogenesis respectively. Results: The phytochemical screening discovered existence of terpenoids, proteins, carbohydrates, phenols, tannins, saponins, flavonoids, glycosides, and alkaloids in both T. coronaria and T. alata. The GC MS profile of the ethanol portion of T. coronaria leaves identified 16 components, while the ethanol fraction of T. alata leaves had 14 components. The molecular docking experiments showed that compound 1 and compound 4 had favorable docking energies of -8.7 kcal.mol-1and -8.2 kcal.mol-1, correspondingly, in the site of JNK-1 kinase. Compound 4 established hydrogen bond interactions with Ser34 and Asp169 in the catalytic and DFG motif regions of the JNK-1, respectively. Compounds 3 and 7, with docking energies of -6.4 and -7.9 kcal.mol-1, correspondingly, also resided in active motif of JNK-1. Compound 2 had docking energy of -5.4 kcal.mol-1 and was well placed in the protein cavity. Regarding the binding of compounds in the KAS III, compound 4 had an excellent docking energy of -8.0kcal.mol-1, and compound 2 had docking energy -5.9 kcal.mol-1. Both compounds were well placed in the active pocket of KAS III macromolecule and established hydrogen bond interactions with Asn260 and Arg262. Compound 4 also established hydrophobic contacts with Arg46 and Arg223. Conclusion: The study states that T. coronaria and T. alata treatment strongly inhibited A549 cells viability, and cell volume expansion, which result in cell proliferation. Likewise, a noteworthy decrease in fin regeneration and reduction in percentage vessel growth was observed in zebra fish and embryo assays.

A thiopyridine-bound mirror-image copper center in an artificial non-heme metalloenzyme

Artificial metalloenzymes, in which a metal complex and protein matrix are combined, have been synthesized to catalyze stereoselective reactions using the chiral environment provided by the protein cavity. Artificial metalloenzymes can be engineered by the chemical modification and mutagenesis of the protein matrix. We developed artificial non-heme metalloenzymes using a cupin superfamily protein (TM1459) with a 4-His tetrad-metal-binding motif. The Cu-bound H52A/C106D mutant with 3-His triad showed a S-enantioselective Michael addition of nitromethane to α,β-unsaturated ketone, 2-aza-chalcone 1. In this study, we demonstrated a chemical modification near the copper-binding site of this mutant to reverse its enantioselectivity. For chemical modification, the amino acid on the Si-face of the binding state of 1 to the copper center was replaced with Cys, followed by reaction with 4,4′-dithiopyridine (4-PDS) to form S-(pyridin-4-ylthio)cysteine (Cys-4py). Cu-bound I49C-4py/H52A/C106D showed reversal of the enantioselectivity from S-form to R-form (ee = 71%, (R)). The effect of steric hindrance of the amino acids at position 49 on enantioselectivity was investigated using I49X/H52A/C106D mutants (X = A, C, I, F, and W). Additionally, chemical modification with 2,2′-dithiopyridine (2-PDS) produced I49-2py/H52A/C106D, which showed lower R-enantioselectivity than I49-4py/H52A/C106D. Among the mutants, the 4py-modification on the Si-face was the most effective in reversing the enantioselectivity. By tuning the Re-face side, the H54A mutation introduced into the I49C-4py/H52A/C106D increased the R-enantioselectivity (ee = 88%, (R)). X-ray crystallography revealed a coordinated structure with ligation of thiopyridine in Cu-bound I49C-4py/H52A/H54A/C106D.

Building shape-focused pharmacophore models for effective docking screening

The performance of molecular docking can be improved by comparing the shape similarity of the flexibly sampled poses against the target proteins’ inverted binding cavities. The effectiveness of these pseudo-ligands or negative image-based models in docking rescoring is boosted further by performing enrichment-driven optimization. Here, we introduce a novel shape-focused pharmacophore modeling algorithm O-LAP that generates a new class of cavity-filling models by clumping together overlapping atomic content via pairwise distance graph clustering. Top-ranked poses of flexibly docked active ligands were used as the modeling input and multiple alternative clustering settings were benchmark-tested thoroughly with five demanding drug targets using random training/test divisions. In docking rescoring, the O-LAP modeling typically improved massively on the default docking enrichment; furthermore, the results indicate that the clustered models work well in rigid docking. The C+ +/Qt5-based algorithm O-LAP is released under the GNU General Public License v3.0 via GitHub (https://github.com/jvlehtonen/overlap-toolkit).Scientific contributionThis study introduces O-LAP, a C++/Qt5-based graph clustering software for generating new type of shape-focused pharmacophore models. In the O-LAP modeling, the target protein cavity is filled with flexibly docked active ligands, the overlapping ligand atoms are clustered, and the shape/electrostatic potential of the resulting model is compared against the flexibly sampled molecular docking poses. The O-LAP modeling is shown to ensure high enrichment in both docking rescoring and rigid docking based on comprehensive benchmark-testing.Graphical

Chemogenetic Evolution of Diversified Photoenzymes for Enantioselective [2 + 2] Cycloadditions in Whole Cells.

Artificial photoenzymes with novel catalytic modes not found in nature are in high demand; yet, they also present significant challenges in the field of biocatalysis. In this study, a chemogenetic modification strategy is developed to facilitate the rapid diversification of photoenzymes. This strategy integrates site-specific chemical conjugation of various artificial photosensitizers into natural protein cavities and the iterative mutagenesis in cell lysates. Through rounds of directed evolution, prominent visible-light-activatable photoenzyme variants were developed, featuring a thioxanthone chromophore. They successfully enabled the enantioselective [2 + 2] photocycloaddition of 2-carboxamide indoles, a class of UV-sensitive substrates that are traditionally challenging for known photoenzymes. Furthermore, the versatility of this photoenzyme is demonstrated in enantioselective whole-cell photobiocatalysis, enabling the efficient synthesis of enantioenriched cyclobutane-fused indoline tetracycles. These findings significantly expand the photophysical properties of artificial photoenzymes, a critical factor in enhancing their potential for harnessing excited-state reactivity in stereoselective transformations.

Molecules targeting a novel homotrimer cavity of Spike protein attenuate replication of SARS-CoV-2

The SARS-CoV-2 Spike glycoprotein (S) utilizes a unique trimeric conformation to interact with the ACE2 receptor on host cells, making it a prime target for inhibitors that block viral entry. We have previously identified a novel proteinaceous cavity within the Spike protein homotrimer that could serve as a binding site for small molecules. However, it is not known whether these molecules would inhibit, stimulate, or have no effect on viral replication. To address this, we employed structural-based screening to identify small molecules that dock into the trimer cavity and assessed their impact on viral replication. Our findings show that a cohort of identified small molecules binding to the Spike trimer cavity effectively reduces the replication of various SARS-CoV-2 variants. These molecules exhibited inhibitory effects on B.1 (European original, D614G, EDB2) and B.1.617.2 (δ) variants, while showing moderate activity against the B.1.1.7 (α) variant. We further categorized these molecules into distinct groups based on their structural similarities. Our experiments demonstrated a dose-dependent viral replication inhibitory activity of these compounds, with some, like BCC0040453 exhibiting no adverse effects on cell viability even at high concentrations. Further investigation revealed that pre-incubating virions with compounds like BCC0031216 at different temperatures significantly inhibited viral replication, suggesting their specificity towards the S protein. Overall, our study highlights the inhibitory impact of a diverse set of chemical molecules on the biological activity of the Spike protein. These findings provide valuable insights into the role of the trimer cavity in the viral replication cycle and aid drug discovery programs aimed at targeting the coronavirus family.

PharmaCore: The Automatic Generation of 3D Structure-Based Pharmacophore Models from Protein/Ligand Complexes.

In this work, we present PharmaCore: a new, completely automatic workflow aimed at generating three-dimensional (3D) structure-based pharmacophore models toward any target of interest. The proposed approach relies on using cocrystallized ligands to create the input files for generating the pharmacophore hypotheses, integrating not only the three-dimensional structural information on the ligand but also data concerning the binding mode of these molecules put in the protein cavity. We developed a Python library that, starting from the specific UniProt ID of the protein under investigation as the only element that requires user intervention, subsequently collects and aligns the corresponding structures bearing a known ligand in a fully automated fashion, bringing them all into the same coordinate system. The protocol includes a final phase in which the aligned small molecules are used to produce the pharmacophore hypotheses directly onto the protein structure using a specific software, e.g., Phase (Schrödinger LLC). To validate the entire procedure and highlight the possible applications in the field of drug discovery and repositioning, we first generated pharmacophores for soluble epoxide hydrolase (sEH) and compared with already-published ones. Then, we reproduced the binding profile of a reported selective binder of ATAD2 bromodomain (AM879), testing it against a panel of 1741 pharmacophores related to 16 epigenetic proteins and automatically generated with PharmaCore, finally disclosing putative unprecedented off-targets. The computational predictions were successfully validated with AlphaScreen assays, highlighting the applicability of the proposed workflow in drug discovery and repositioning. Finally, the process was also validated on tankyrase 2 and SARS-CoV-2 MPro, confirming the robustness of PharmaCore.

Discovery of a Novel Chemo-Type for TAAR1 Agonism via Molecular Modeling.

The search for novel effective TAAR1 ligands continues to draw great attention due to the wide range of pharmacological applications related to TAAR1 targeting. Herein, molecular docking studies of known TAAR1 ligands, characterized by an oxazoline core, have been performed in order to identify novel promising chemo-types for the discovery of more active TAAR1 agonists. In particular, the oxazoline-based compound S18616 has been taken as a reference compound for the computational study, leading to the development of quite flat and conformationally locked ligands. The choice of a "Y-shape" conformation was suggested for the design of TAAR1 ligands, interacting with the protein cavity delimited by ASP103 and aromatic residues such as PHE186, PHE195, PHE268, and PHE267. The obtained results allowed us to preliminary in silico screen an in-house series of pyrimidinone-benzimidazoles (1a-10a) as a novel scaffold to target TAAR1. Combined ligand-based (LBCM) and structure based (SBCM) computational methods suggested the biological evaluation of compounds 1a-10a, leading to the identification of derivatives 1a-3a (hTAAR1 EC50 = 526.3-657.4 nM) as promising novel TAAR1 agonists.

VirtuousPocketome: a computational tool for screening protein–ligand complexes to identify similar binding sites

Protein residues within binding pockets play a critical role in determining the range of ligands that can interact with a protein, influencing its structure and function. Identifying structural similarities in proteins offers valuable insights into their function and activation mechanisms, aiding in predicting protein–ligand interactions, anticipating off-target effects, and facilitating the development of therapeutic agents. Numerous computational methods assessing global or local similarity in protein cavities have emerged, but their utilization is impeded by complexity, impractical automation for amino acid pattern searches, and an inability to evaluate the dynamics of scrutinized protein–ligand systems. Here, we present a general, automatic and unbiased computational pipeline, named VirtuousPocketome, aimed at screening huge databases of proteins for similar binding pockets starting from an interested protein–ligand complex. We demonstrate the pipeline's potential by exploring a recently-solved human bitter taste receptor, i.e. the TAS2R46, complexed with strychnine. We pinpointed 145 proteins sharing similar binding sites compared to the analysed bitter taste receptor and the enrichment analysis highlighted the related biological processes, molecular functions and cellular components. This work represents the foundation for future studies aimed at understanding the effective role of tastants outside the gustatory system: this could pave the way towards the rationalization of the diet as a supplement to standard pharmacological treatments and the design of novel tastants-inspired compounds to target other proteins involved in specific diseases or disorders. The proposed pipeline is publicly accessible, can be applied to any protein–ligand complex, and could be expanded to screen any database of protein structures.

Non-oxidized and oxidized flaxseed orbitides differently induce HepG2 cell apoptosis: involvement of cellular uptake and membrane death receptor DR4.

Flaxseed orbitides have health-promoting properties, particularly potent anti-cancer activity. However, flaxseed orbitides containing a methionine structure, such as [1-9-NαC]-linusorb B2 (CLB), are easily oxidized to sulfoxide ([1-9-NαC],[1-Rs,Ss-MetO]-linusorb-B2 (CLC)) and sulfone ([1-9-NαC], [1-MetO]-linusorb B2 (CLK)), with CLC having less anti-cancer ability than CLB. It is unclear why oxidized flaxseed orbitides are less effective against cancer than non-oxidized flaxseed orbitide. Non-oxidized ([1-9-NαC]-linusorb-B3 (CLA) and CLB) and oxidized (CLC and CLK) flaxseed orbitides were found to significantly upregulate the levels of pro-apoptotic proteins, including Bax/Bcl-2, CytoC, caspase-3, and caspase-8, in a dose-dependent manner, with non-oxidized flaxseed orbitides being more effective than oxidized flaxseed orbitides. Mechanically, the cellular absorption of non-oxidized flaxseed orbitides was higher than that of oxidized flaxseed orbitides. Moreover, the significant fluorescence quenching of DR4 protein by flaxseed orbitides (especially non-oxidized orbitides) indicated the formation of a DR4-orbitide complex. Molecular docking demonstrated that non-oxidized orbitides could easily dock into the active cavity of DR4 protein. Further blocking DR4 significantly reduced the ability of non-oxidized flaxseed orbitides to stimulate caspase-3 expression, whereas oxidized flaxseed orbitides retained this ability. Non-oxidized flaxseed orbitides are more effective against cancer than oxidized flaxseed orbitides due to higher cellular uptake and activation of the DR4-mediated death receptor signaling pathway. © 2024 Society of Chemical Industry.

Efficient Sampling of Cavity Hydration in Proteins with Nonequilibrium Grand Canonical Monte Carlo and Polarizable Force Fields.

Prediction of the hydration levels of protein cavities and active sites is important to both mechanistic analysis and ligand design. Due to the unique microscopic environment of these buried water molecules, a polarizable model is expected to be crucial for an accurate treatment of protein internal hydration in simulations. Here we adapt a nonequilibrium candidate Monte Carlo approach for conducting grand canonical Monte Carlo simulations with the Drude polarizable force field. The GPU implementation enables the efficient sampling of internal cavity hydration levels in biomolecular systems. We also develop an enhanced sampling approach referred to as B-walking, which satisfies detailed balance and readily combines with grand canonical integration to efficiently calculate quantitative binding free energies of water to protein cavities. Applications of these developments are illustrated in a solvent box and the polar ligand binding site in trypsin. Our simulation results show that including electronic polarization leads to a modest but clear improvement in the description of water position and occupancy compared to the crystal structure. The B-walking approach enhances the range of water sampling in different chemical potential windows and thus improves the accuracy of water binding free energy calculations.

Targeted Affinity Purification and Mechanism of Action of Angiotensin-Converting Enzyme (ACE) Inhibitory Peptides from Sea Cucumber Gonads.

Protein hydrolysates from sea cucumber (Apostichopus japonicus) gonads are rich in active materials with remarkable angiotensin-converting enzyme (ACE) inhibitory activity. Alcalase was used to hydrolyze sea cucumber gonads, and the hydrolysate was separated by the ultrafiltration membrane to produce a low-molecular-weight peptide component (less than 3 kDa) with good ACE inhibitory activity. The peptide component (less than 3 kDa) was isolated and purified using a combination method of ACE gel affinity chromatography and reverse high-performance liquid chromatography. The purified fractions were identified by liquid chromatography-tandem mass spectrometry (LC-MS/MS), and the resulting products were filtered using structure-based virtual screening (SBVS) to obtain 20 peptides. Of those, three noncompetitive inhibitory peptides (DDQIHIF with an IC50 value of 333.5 μmol·L-1, HDWWKER with an IC50 value of 583.6 μmol·L-1, and THDWWKER with an IC50 value of 1291.8 μmol·L-1) were further investigated based on their favorable pharmacochemical properties and ACE inhibitory activity. Molecular docking studies indicated that the three peptides were entirely enclosed within the ACE protein cavity, improving the overall stability of the complex through interaction forces with the ACE active site. The total free binding energies (ΔGtotal) for DDQIHIF, HDWWKER, and THDWWKER were -21.9 Kcal·mol-1, -71.6 Kcal·mol-1, and -69.1 Kcal·mol-1, respectively. Furthermore, a short-term assay of antihypertensive activity in spontaneously hypertensive rats (SHRs) revealed that HDWWKER could significantly decrease the systolic blood pressure (SBP) of SHRs after intravenous administration. The results showed that based on the better antihypertensive activity of the peptide in SHRs, the feasibility of targeted affinity purification and computer-aided drug discovery (CADD) for the efficient screening and preparation of ACE inhibitory peptide was verified, which provided a new idea of modern drug development method for clinical use.

Modeling the SDF-1/CXCR4 protein using advanced artificial intelligence and antagonist screening for Japanese anchovy.

SDF-1/CXCR4 chemokine signaling are indispensable for cell migration, especially the Primordial Germ Cell (PGC) migration towards the gonadal ridge during early development. We earlier found that this signaling is largely conserved in the Japanese anchovy (Engraulis japonicus, EJ), and a mere treatment of CXCR4 antagonist, AMD3100, leads to germ cell depletion and thereafter gonad sterilization. However, the effect of AMD3100 was limited. So, in this research, we scouted for CXCR4 antagonist with higher potency by employing advanced artificial intelligence deep learning-based computer simulations. Three potential candidates, AMD3465, WZ811, and LY2510924, were selected and in vivo validation was conducted using Japanese anchovy embryos. We found that seven transmembrane motif of EJ CXCR4a and EJ CXCR4b were extremely similar with human homolog while the CXCR4 chemokine receptor N terminal (PF12109, essential for SDF-1 binding) was missing in EJ CXCR4b. 3D protein analysis and cavity search predicted the cavity in EJ CXCR4a to be five times larger (6,307ų) than that in EJ CXCR4b (1,241ų). Docking analysis demonstrated lower binding energy of AMD3100 and AMD3465 to EJ CXCR4a (Vina score -9.6) and EJ CXCR4b (Vina score -8.8), respectively. Furthermore, we observed significant PGC mismigration in microinjected AMD3465 treated groups at 10, 100 and 1 × 105nM concentration in 48h post fertilized embryos. The other three antagonists showed various degrees of PGC dispersion, but no significant effect compared to their solvent control at tested concentrations was observed. Cumulatively, our results suggests that AMD3645 might be a better candidate for abnormal PGC migration in Japanese anchovy and warrants further investigation.

CRAFT: a web-integrated cavity prediction tool based on flow transfer algorithm

Numerous computational methods, including evolutionary-based, energy-based, and geometrical-based methods, are utilized to identify cavities inside proteins. Cavity information aids protein function annotation, drug design, poly-pharmacology, and allosteric site investigation. This article introduces “flow transfer algorithm” for rapid and effective identification of diverse protein cavities through multidimensional cavity scan. Initially, it identifies delimiter and susceptible tetrahedra to establish boundary regions and provide seed tetrahedra. Seed tetrahedron faces are precisely scanned using the maximum circle radius to transfer seed flow to neighboring tetrahedra. Seed flow continues until terminated by boundaries or forbidden faces, where a face is forbidden if the estimated maximum circle radius is less or equal to the user-defined maximum circle radius. After a seed scanning, tetrahedra involved in the flow are clustered to locate the cavity. The CRAFT web interface integrates this algorithm for protein cavity identification with enhanced user control. It supports proteins with cofactors, hydrogens, and ligands and provides comprehensive features such as 3D visualization, cavity physicochemical properties, percentage contribution graphs, and highlighted residues for each cavity. CRAFT can be accessed through its web interface at http://pitools.niper.ac.in/CRAFT, complemented by the command version available at https://github.com/PGlab-NIPER/CRAFT/.Scientific contributionFlow transfer algorithm is a novel geometric approach for accurate and reliable prediction of diverse protein cavities. This algorithm employs a distinct concept involving maximum circle radius within the 3D Delaunay triangulation to address diverse van der Waals radii while existing methods overlook atom specific van der Waals radii or rely on complex weighted geometric techniques.