Studies on molecular mechanism between ACE and inhibitory peptides in different bioactivities by 3D-QSAR and MD simulations

2-Phenylcyclopropylmethylamine (PCPMA) Derivatives as D3R-Selective Ligands for 3D-QSAR, Docking and Molecular Dynamics Simulation Studies

Dopamine D3 receptor (D3R) is a key receptor for regulating motor, cognitive, and other functions. In this study, 50 2-phenylcyclopropylmethylamine (PCPMA) derivatives with good selectivity for D3R were investigated using a three-dimensional quantitative structure-activity relationship (3D-QSAR) method. The CoMFA and CoMSIA model results showed good predictive ability, as evidenced by high r2 and q2 values. 3D-QSAR results showed that steric, electrostatic, and hydrophobic fields played important roles in the binding of PCPMAs to D3R. Based on above results, four novel PCPMAs were designed, which were predicted to have a stronger affinity with D3R. Molecular docking combined with 300 ns molecular dynamics simulations were performed to reveal the mode of interaction between D3R and PCPMAs. Additionally, a combination of free energy calculations and energy decomposition results indicated strong interaction between the ligands and residues in the binding pocket of the D3 receptor. This work provides suggestions for exploring more selective D3R ligands, and this theoretical framework also lays the foundation for future experimental investigations to evaluate the pharmacological characteristics and binding affinities of novel derivatives.

The novel angiotensin-I-converting enzyme inhibitory peptides from Scomber japonicus muscle protein hydrolysates: QSAR-based screening, molecular docking, kinetic and stability studies

Food-derived angiotensin-converting enzyme-inhibitory (ACE-I) peptides have attracted extensive attention. Herein, the ACE-I peptides from Scomber japonicus muscle hydrolysates were screened, and their mechanisms of action and inhibition stability were explored. The quantitative structure-activity relationship (QSAR) model based on 5z-scale metrics was developed to rapidly screen for ACE-I peptides. Two novel potential ACE-I peptides (LTPFT, PLITT) were predicted through this model coupled with in silico screening, of which PLITT had the highest activity (IC50: 48.73 ± 7.59 μM). PLITT inhibited ACE activity with a mixture of non-competitive and competitive mechanisms, and this inhibition mainly contributed to the hydrogen bonding based on molecular docking study. PLITT is stable under high temperatures, pH, glucose, and NaCl. The zinc ions (Zn2+) and copper ions (Cu2+) enhanced ACE-I activity. The study suggests that the QSAR model is effective in rapidly screening for ACE-I inhibitors, and PLITT can be supplemented in foods to lower blood pressure.

ACE inhibitory peptides from enzymatic hydrolysate of fermented black sesame seed: Random forest-based optimization, screening, and molecular docking analysis

In this study, black sesame seeds were fermented by Lactobacillus Plantarum NCU116 and then hydrolyzed using acid protease to improve Angiotensin-I-converting enzyme (ACE) inhibitory activity. The random forest-particle swarm optimization (RF-PSO) model was applied to predict the ACE inhibitory activity during the hydrolysis process based on the experimental data. After separating by adsorption chromatography, gel filtration chromatography, and reversed phased-high performance liquid chromatography and then screening in silico method, eight peptides were identified from fermented black sesame seed hydrolysates as ITAPHW, SLPNYHPSPR, QYLPR, IRPNGL, YHNAPIL, LSYPR, GFAGDDAPRA, and LDPNPRSF with IC50 values of 51.69 μM, 146.67 μM, 655.02 μM, 752.60 μM, 1.02 mM, 2.01 mM, 1.97 mM, and 3.43 mM, respectively. ITAPHW and SLPNYHPSPR exhibited high antioxidant activity and inhibited the ACE activity in a non-competitive pattern. Molecular docking revealed that the strong ACE inhibition of ITAPHW and SLPNYHPSPR is probably attributed to the interaction with Zn2+ of ACE.

Screening and molecular dynamics simulation of ACE inhibitory tripeptides derived from milk fermented with Lactobacillus delbrueckii QS306

Peptides in milk fermented with Lactobacillus delbrueckii QS306 before and after ultrahigh pressure treatment were identified using proteomics. Subsequently, 16 stable tripeptides were screened out based on activity score prediction, PeptideCutter analysis, and hydrophobicity calculations. Among them, WRP, WSR, and YRP showed the best angiotensin-converting enzyme (ACE) inhibitory activity, and their semi-inhibitory concentrations were 46.707, 300.121, and 89.555 μM, respectively. WRP and WSR were competitive inhibitors, whereas YRP was non-competitive. Gastrointestinal simulation revealed that WRP and YRP had better gastrointestinal stability. The values of RMSD, ΔGbind, ΔGpol, and RSMF obtained from molecular dynamics simulation indicated that the interaction of WRP and ACE was stable. Thus, Lactobacillus delbrueckii QS306-fermented milk can serve as an important source of ACE inhibitory peptides both before and after ultrahigh pressure treatment. The strategy of in silico screening, activity evaluation, and molecular dynamics simulation adopted in this study can be applied to the large-scale screening of novel peptides with high ACE inhibitory activity.

In silico study to recognize novel angiotensin-converting-enzyme-I inhibitors by 2D-QSAR and constraint-based molecular simulations

Cardiovascular diseases (CVD) such as heart failure, stroke, and hypertension affect 64.3 million people worldwide and are responsible for 30% of all deaths. Primary inhibition of the angiotensin-converting enzyme (ACE) is significant in the management of CVD. In the present study, the genetic algorithm-multiple linear regressions (GA-MLR) method is used to generate highly predictive and statistically significant (R2 = 0.70-0.75, Q2LOO=0.67-0.73, Q2LMO=0.66-0.72, CCCex=0.70-0.78) quantitative structure-activity relationships (QSAR) models conferring to OECD requirements using a dataset of 255 structurally diverse and experimentally validated ACE inhibitors. The models contain simply illustratable Padel, Estate, and PyDescriptors that correlate structural scaffold requisite for ACE inhibition. Also, constraint-based molecular docking reveals an interaction profile between ligands and enzymes which is then correlated with the essential structural features associated with the QSAR models. The QSAR-based virtual screening was utilized to find novel lead molecules from a designed database of 102 thiadiazole derivatives. The Applicability domain (AD), Molecular Docking, Molecular dynamics, and ADMET analysis suggest two compound D24 and D40 are inflexibly linked to the protein binding site and follows drug-likeness properties.Communicated by Ramaswamy H. Sarma.

Screening of Oral Potential Angiotensin-Converting Enzyme Inhibitory Peptides from Zizyphus jujuba Proteins Based on Gastrointestinal Digestion In Vivo

Plant proteins are a good source of active peptides, which can exert physiological effects on the body. Predicting the possible activity of plant proteins and obtaining active peptides with oral potential are challenging. In this study, the potential activity of peptides from Zizyphus jujuba proteins after in silico simulated gastrointestinal digestion was predicted using the BIOPEP-UWM™ database. The ACE-inhibitory activity needs to be further investigated. The actual peptides in mouse intestines after the oral administration of Zizyphus jujuba protein were collected and analyzed, 113 Zizyphus jujuba peptides were identified, and 3D-QSAR models of the ACE-inhibitory activity were created and validated using a training set (34 peptides) and a test set (12 peptides). Three peptides, RLPHV, TVKPGL and KALVAP, were screened using the 3D-QSAR model and were found to bind to the active sites of the ACE enzyme, and their IC50 values were determined. Their values were 6.01, 3.81, and 17.06 μM, respectively. The in vitro digestion stabilities of the RLPHV, TVKPGL, and KALVAP peptides were 82%, 90%, and 78%. This article provides an integrated method for studying bioactive peptides derived from plant proteins.

Screening, Characterization, and Mechanistic Evaluation of Angiotensin Converting Enzyme Inhibitory Peptides Derived from Milk Fermented with Lactobacillus delbrueckii QS306 with and without Ultrahigh-Pressure Treatment

In this study, the peptides in milk fermented with Lactobacillus delbrueckii QS306 with and without ultrahigh-pressure treatment were identified using UPLC-Q-exactive-HF-X-MS/MS. In total, 27 novel pentapeptides with potential angiotensin converting enzyme inhibitory (ACEI) activity were screened via bioinformatic analysis, and the activities of seven novel pentapeptides were assessed. Among them, HLPLP, PYPQR, and VAPFP exhibited better IC50 values. Stability assessment via in vitro simulation revealed that the three pentapeptides were significantly reduced (P < 0.05) during digestion but exhibited 85% activity after digestion. HLPLP was a competitive inhibitor, while PYPQR and VAPFP were noncompetitive inhibitors of ACE. Molecular docking indicated that the three peptides could stably bind to ACE. Molecular dynamics (MD) simulation and activity verification indicated that PYPQR and VAPFP had better stability and activity. This study demonstrated that novel ACEI pentapeptides in milk fermented with L. delbrueckii QS306 with and without ultrahigh-pressure treatment could be considered promising candidates for controlling hypertension.

Structure-Activity Relationship of Novel ACE Inhibitory Undecapeptides from Stropharia rugosoannulata by Molecular Interactions and Activity Analyses

Undecapeptide is the central peptide molecule in the peptide base material of Stropharia rugosoannulata, and angiotensin-converting enzyme (ACE) plays a crucial role in hypertension. To fully explore the interaction mechanism and ACE-inhibitory activity of long-chain peptides from Stropharia rugosoannulata, the binding conformations of twenty-seven undecapeptides with the ACE receptor were revealed by molecule docking. The undecapeptide GQEDYDRLRPL with better receptor binding capacity and higher secondary mass spectral abundance was screened. All amino acid residues except proline in GQEDYDRLRPL interacted with the ACE receptor. GQEDYDRLRPL interfered with the receptor's overall structure, with significant fluctuations in amino acid residues 340-355, including two residues in the receptor's active pockets. The binding constants of GQEDYDRLRPL to the ACE receptors were at the μM level, with a kinetic binding constant of 9.26 × 10-7 M, which is a strong binding, and a thermodynamic binding constant of 3.06 × 10-6 M. Intermolecular interaction were exothermic, enthalpy-driven, and specific binding reactions. GQEDYDRLRPL had an IC50 value of 164.41 μmol/L in vitro and superior antihypertensive effects at low-gavage administration in vivo. Obtaining information on the interaction mechanism of ACE-inhibitory undecapeptides from S. rugosoannulata with the ACE receptor will help to develop and utilize ACE inhibitors of natural origin.

Revealing the Sequence Characteristics and Molecular Mechanisms of ACE Inhibitory Peptides by Comprehensive Characterization of 160,000 Tetrapeptides

Chronic diseases, such as hypertension, cause great harm to human health. Conventional drugs have promising therapeutic effects, but also cause significant side effects. Food-sourced angiotensin-converting enzyme (ACE) inhibitory peptides are an excellent therapeutic alternative to pharmaceuticals, as they have fewer side effects. However, there is no systematic and effective screening method for ACE inhibitory peptides, and the lack of understanding of the sequence characteristics and molecular mechanism of these inhibitory peptides poses a major obstacle to the development of ACE inhibitory peptides. Through systematically calculating the binding effects of 160,000 tetrapeptides with ACE by molecular docking, we found that peptides with Tyr, Phe, His, Arg, and especially Trp were the characteristic amino acids of ACE inhibitory peptides. The tetrapeptides of WWNW, WRQF, WFRV, YYWK, WWDW, and WWTY rank in the top 10 peptides exhibiting significantly high ACE inhibiting behaviors, with IC₅₀ values between 19.98 ± 8.19 μM and 36.76 ± 1.32 μM. Salt bridges, π-π stacking, π-cations, and hydrogen bonds contributed to the high binding characteristics of the inhibitors and ACE. Introducing eight Trp into rabbit skeletal muscle protein (no Trp in wide sequence) endowed the protein with a more than 90% ACE inhibition rate, further suggesting that meat with a high content of Trp could have potential utility in hypertension regulation. This study provides a clear direction for the development and screening of ACE inhibitory peptides.

Molecular modeling study of pyrrolidine derivatives as novel myeloid cell leukemia-1 inhibitors through combined 3D-QSAR, molecular docking, ADME/Tox and MD simulation techniques

A series of pyrrolidine derivatives have been used to study the main structural requirements for designing novel Mcl-1 inhibitors. For this purpose, three models CoMSIA, CoMFA and HQSAR were generated using QSAR molecular modeling techniques. The statistical results of the CoMFA (Q2 = 0.689; R = 0.999; R2pred = 0.986), CoMSIA (Q2 = 0.614; R2 = 0.923; R2pred = 0.815) and HQSAR (Q2= 0.603; R2 = 0.662; R2pred = 0.743) models showed good stability and predictability. The results of the models were presented as contours and colored fragments indicating the favorable and unfavorable contribution to the inhibitory activity of Mcl-1. Based on the obtained results, four new compounds were designed with more potent predicted pIC50 inhibitory activity. The ADME/Tox results and the pharmacokinetic properties revealed that these four compounds are orally bioavailable and show good permeability. In addition the four compounds showing non-inhibitors of CYP3A4 and CYP2D6 with the exception of Pred03. At the level of toxicity profile, the compounds Pred01, Pred02 and Pred03 showed interesting results and showed no AMES toxicity, no hERG inhibition and no skin sensitization. Molecular docking results were used to uncover the mode of interaction between the ligand and key residues of protein binding site. Molecular docking results were supported by molecular simulation and binding free energy estimation (MMPBSA). These results demonstrate the stability of the analyzed compounds in the target protein binding site during a 100 ns trajectory. Finally, all these results create a strong lead to develop promising new Pyrrolidine-based inhibitors against Mcl-1.Communicated by Ramaswamy H. Sarma.

Angiotensin converting enzyme inhibitors from medicinal plants: a molecular docking and dynamic simulation approach

Angiotensin converting enzyme (ACE) is a key enzyme and mediator in the aetiology of high blood pressure (HBP) and hypertension. As one of the leading cause of untimely death worldwide, there is a lot of research and studies on the management and treatment of hypertension. The usage of medicinal plants in the management of hypertension as alternative to synthetic allopathic drugs is a common practice in folkloric and traditional medicine. Therefore, this study was aimed to investigate the ACE inhibitory activity of some medicinal plants which are commonly used in the treatment of HBP in southwestern part of Nigeria using extensive in-silico approach. Compounds identified in the plants through GC-MS technique, together with Lisinopril were docked against ACE protein. It was observed that only 40 of the compounds had binding affinity ≥ - 6.8 kcal/mol which was demonstrated by the standard drug (lisinopril). Interaction between the compounds and ACE was via conventional hydrogen, carbon hydrogen, alkyl, pi-alkyl, pi-carbon, and Van Der Wall bonds among others. Most of these compounds exhibited drug like properties, without violating majority of the physicochemical descriptors and Lipinski rule of 5. The ADMET evaluation revealed that only 2 compounds (cyclopentadecanone and oxacycloheptadecan-2-one) which were identified in Bacopa florinbunda plant were predicted non-toxic and thus were subjected to molecular dynamics and simulation with ACE. From the molecular dynamics and mechanics analysis, both cyclopentadecanone and oxacycloheptadecan-2-one showed high stability and inhibitory potentials when bound to ACE. Oxacycloheptadecan-2-one was more stable than lisinopril and cyclopentadecanone in the ligand-ACE complex; we therefore suggested its experimental and clinical validation as drug candidates for the treatment of hypertension.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s40203-022-00135-z.

Bioactive Peptides: An Understanding from Current Screening Methodology

Bioactive peptides with high potency against numerous human disorders have been regarded as a promising therapy in disease control. These peptides could be released from various dietary protein sources through hydrolysis processing using physical conditions, chemical agents, microbial fermentation, or enzymatic digestions. Considering the diversity of the original proteins and the complexity of the multiple structural peptides that existed in the hydrolysis mixture, the screening of bioactive peptides will be a challenge task. Well-organized and well-designed methods are necessarily required to enhance the efficiency of studying the potential peptides. This article, hence, provides an overview of bioactive peptides with an emphasis on the current strategy used for screening and characterization methods. Moreover, the understanding of the biological activities of peptides, mechanism inhibitions, and the interaction of the complex of peptide–enzyme is commonly evaluated using specific in vitro assays and molecular docking analysis.

A QSAR-ICE-SSD Model Prediction of the PNECs for Per- and Polyfluoroalkyl Substances and Their Ecological Risks in an Area of Electroplating Factories

Per- and polyfluoroalkyl substances (PFASs) are a class of highly fluorinated aliphatic compounds that are persistent and bioaccumulate, posing a potential threat to the aquatic environment. The electroplating industry is considered to be an important source of PFASs. Due to emerging PFASs and many alternatives, the acute toxicity data for PFASs and their alternatives are relatively limited. In this study, a QSAR–ICE–SSD composite model was constructed by combining quantitative structure-activity relationship (QSAR), interspecies correlation estimation (ICE), and species sensitivity distribution (SSD) models in order to obtain the predicted no-effect concentrations (PNECs) of selected PFASs. The PNECs for the selected PFASs ranged from 0.254 to 6.27 mg/L. The ΣPFAS concentrations ranged from 177 to 983 ng/L in a river close to an electroplating industry in Shenzhen. The ecological risks associated with PFASs in the river were below 2.97 × 10⁻⁴.

Accurate prediction for adsorption rate of peptides with high ACE-inhibitory activity from sericin hydrolysate on thiophene hypercross-linked polymer using CoMSIA in 3D-QSAR model

Efficient screening of angiotensin converting enzyme inhibitory (ACE-I) peptides from agricultural or edible sources attract increasing attention. However, their purification process from the complex natural system is still semi-empirical or even uncontrollable, which has seriously reduced their screening efficiency. Herein, inspired by the prediction of ACE-I activity, 3D-QSAR was proposed to predict the adsorption performance of peptides from sericin hydrolysate (SH) having high ACE-I activity on porous hypercross-linked polymers according to their molecular structures. Thiophene hypercross-linked polymer (T-HCP) possessing better screening capacity for ACE-I peptides was chosen as our research object in this work. The sequence and relative adsorption rate of 101 peptides in SH were analyzed by LC-MS and was used as a database to construct the relationship of peptide's chemical structure and adsorption performance on T-HCP by Comparative molecular similarity indices analysis (CoMSIA) from 3D-QSAR. Optimum CoMSIA revealed that enhanced interaction of hydrophobicity and H-bond between T-HCP and the peptide was conducive to increase the adsorption performance of di- to hexa- peptides. Based on these relationships, the adsorption capability of 24 designed peptides with distinguished hydrophobic and H-bond fields was predicted on T-HCP by using optimum CoMSIA and the results of half of these were verified, which showed high consistency with their predicted adsorption rate. Interestingly, these peptides having higher adsorption capacities on T-HCP also possessed higher ACE-I activity. This can be attributed to the high concentration of aromatic surface with π-π interaction and weak-polar CSC group with H-bond interaction on T-HCP material, which is ideal for the selective adsorption of peptides with higher ACE-I activity from SH. This study provides important theoretical guidance for the industrial screening of bio-functional peptides from complex protein mixtures.

Structure characteristics of flavonoids for heterocyclic aromatic amines inhibition using quantitative structure-activity relationship modeling

The objective of this study was to investigate the structure characteristics of flavonoids that act as inhibitors for heterocyclic aromatic amines (HAAs) formation. Five quantitative structure-activity relationship models for predicting the inhibitory rates of HAAs (norharman, harman, PhIP, MeIQx, and 4,8-DiMeIQx) were established using selected chemometric parameters (R² : 0.591-0.920), and indicated that the hydrophobicity, hydroxyl groups, and topological structure of flavonoids played important roles in the inhibition of HAAs formation. The 5,7-dihydroxyls in meta-position of the A-ring and the 4'-hydroxyl in the B-ring of flavonoids were critical for the inhibitory effects of HAAs, whereas the introduction of 3-hydroxyl and 3-O-glucoside in the C-ring reduced the inhibitory effects. Catechin served as the most effective inhibitor of HAAs followed by luteolin and genistein. The study can bring us a broader idea for controlling the formation of HAAs according to the structure of flavonoids. PRACTICAL APPLICATIONS: Heterocyclic aromatic amines (HAAs) are a class of organic substances with carcinogenic and mutagenic effect formed during the heating process of meat products. The formation of HAAs can be inhibited by adding natural antioxidants such as flavonoids to the meat during pretreatment. This inhibition is influenced by the unique structure of flavonoids. Thus, there has been an increasing demand to exploit the effective HAAs inhibitors from flavonoids by structure characteristics. Our study showed that the inhibitory effect of flavonoids on the formation of HAAs was mainly depended on their hydrophobicity, hydroxyl groups, and topological structure using the multiple QSAR models. Thus, effective HAAs inhibitors can be explored from dietary flavonoids according their structure characteristics.