Molecular Dynamics Simulation and Molecular Docking Studies of Angiotensin Converting Enzyme with Inhibitor Lisinopril and Amyloid Beta Peptide

The Contribution of the Renin-Angiotensin System to Alzheimer's Disease

The renin-angiotensin system (RAS) is becoming increasingly recognised as a biochemical pathway relevant to the development and progression of Alzheimer's disease (AD). RAS involvement in AD was initially linked to AD via numerous genetic association studies and more recent Genome-Wide Association Studies (GWAS), and in some cases in relation to classical hallmarks of AD pathology. Since these initial findings, which will be summarised here, several complementary areas of research are converging in support of what has been proposed as the Angiotensin Hypothesis for Alzheimer's disease. This hypothesis proposes how the RAS and disease-associated changes to the normal balance between opposing regulatory pathways within RAS warrant careful consideration in the pathogenesis of AD and its pathology. We discuss some of these in relation to RAS-targeting therapeutics, originally developed for the treatment of cardiovascular conditions, and how they might be repurposed as interventions for AD.

Anti-Melanogenic Effects of Takifugu flavidus Muscle Hydrolysate in B16F10 Melanoma Cells and Zebrafish

Abnormal melanogenesis can lead to hyperpigmentation. Tyrosinase (TYR), a key rate-limiting enzyme in melanin production, is an important therapeutic target for these disorders. We investigated the TYR inhibitory activity of hydrolysates extracted from the muscle tissue of Takifugu flavidus (TFMH). We used computer-aided virtual screening to identify a novel peptide that potently inhibited melanin synthesis, simulated its binding mode to TYR, and evaluated functional efficacy in vitro and in vivo. TFMH inhibited the diphenolase activities of mTYR, reducing TYR substrate binding activity and effectively inhibiting melanin synthesis. TFMH indirectly reduced cAMP response element-binding protein phosphorylation in vitro by downregulating melanocortin 1 receptor expression, thereby inhibiting expression of the microphthalmia-associated transcription factor, further decreasing TYR, tyrosinase related protein 1, and dopachrome tautomerase expression and ultimately impeding melanin synthesis. In zebrafish, TFMH significantly reduced black spot formation. TFMH (200 μg/mL) decreased zebrafish TYR activity by 43% and melanin content by 52%. Molecular dynamics simulations over 100 ns revealed that the FGFRSP (T-6) peptide stably binds mushroom TYR via hydrogen bonds and ionic interactions. T-6 (400 μmol/L) reduced melanin content in B16F10 melanoma cells by 71% and TYR activity by 79%. In zebrafish, T-6 (200 μmol/L) inhibited melanin production by 64%. TFMH and T-6 exhibit good potential for the development of natural skin-whitening cosmetic products.

An investigation of binding interactions of tumor-targeted peptide conjugated polyphenols with the kinase domain of ephrin B4 and B2 receptors

Recent studies have shown that Ephrin receptors may be upregulated in several types of cancers including breast, ovarian and endometrial cancers, making them a target for drug design. In this work, we have utilized a target-hopping approach to design new natural product-peptide conjugates and examined their interactions with the kinase-binding domain of EphB4 and EphB2 receptors. The peptide sequences were generated through point mutations of the known EphB4 antagonist peptide TNYLFSPNGPIA. Their anticancer properties and secondary structures were analyzed computationally. Conjugates of most optimum of peptides were then designed by binding the N-terminal of the peptides with the free carboxyl group of the polyphenols sinapate, gallate and coumarate, which are known for their inherent anticancer properties. To investigate if these conjugates have a potential to bind to the kinase domain, we carried out docking studies and MMGBSA free energy calculations of the trajectories based on the molecular dynamics simulations, with both the apo and the ATP bound kinase domains of both receptors. In most cases binding interactions occurred within the catalytic loop region, while in some cases the conjugates were found to spread out across the N-lobe and the DFG motif region. The conjugates were further tested for prediction of pharmacokinetic properties using ADME studies. Our results indicated that the conjugates were lipophilic and MDCK permeable with no CYP interactions. These findings provide an insight into the molecular interactions of these peptides and conjugates with the kinase domain of the EphB4 and EphB2 receptor. As a proof of concept, we synthesized and carried out SPR analysis with two of the conjugates (gallate-TNYLFSPNGPIA and sinapate-TNYLFSPNGPIA). Results indicated that the conjugates showed higher binding with the EphB4 receptor and minimal binding to EphB2 receptor. Sinapate-TNYLFSPNGPIA showed inhibitory activity against EphB4. These studies reveal that some of the conjugates may be developed for further investigation into in vitro and in vivo studies and potential development as therapeutics.

Proteomics analysis to investigate the potential mechanism of theacrine against UV-induced skin photodamage

AIM

This study aimed to explore the underlying mechanism of theacrine treatment of UV-induced skin photodamage.

MATERIALS AND METHODS

Tandem Mass Tag (TMT) relative quantitative proteomics analysis was used to characterize the proteins and pathways associated with the ability of theacrine to combat photodamage in mouse skin by modeling UV irradiation of the backs of ICR mice.

RESULTS

Apoptosis-related proteins and signaling pathways play a key role in the ability of theacrine to protect against skin photodamage, according to proteomic and bioinformatics analysis; molecular docking and Western blotting further revealed that theacrine was associated with apoptosis-related proteins (p53, Bcl-2, Bax, caspase-3, and cleaved-caspase-3) with strong binding affinity, which can significantly reduce skin cell apoptosis induced by UV exposure.

CONCLUSION

The findings revealed that theacrine can reduce UVB-induced epidermal damage by controlling the apoptosis signaling pathway, implying that theacrine could be a useful anti-UVB damage agent.

Angiotensin I-Converting Enzyme Inhibitory Activity of Two Peptides Derived from In Vitro Digestion Products of Pork Sausage with Partial Substitution of NaCl by KCl

This study aimed to identify angiotensin I-converting enzyme (ACE) from in vitro digestion products of pork sausage with partial substitution of NaCl by KCl (PSRK). Peptides from in vitro digestion products of PSRK were identified through liquid chromatography with tandem mass spectrometry analysis coupled with de novo sequencing. Subsequently, the ACE inhibitory peptides LIVGFPAYGH and IVGFPAYGH were screened based on PeptideRanker, in silico absorption, molecular docking, and the determination of ACE inhibitory activity. In addition, the ACE inhibitory peptides LIVGFPAYGH and IVGFPAYGH were mixed-type inhibitors; these peptides' ACE inhibitory activities were expressed as the 50% inhibitory concentration (IC₅₀) values in vitro, which were 196.16 and 150.88 μM, respectively. After 2 h of incubation, LIVGFPAYGH and IVGFPAYGH could be transported through Caco-2 cell monolayers with paracellular passive diffusion. Furthermore, LIVGFPAYGH and IVGFPAYGH significantly increased the levels of ACE2 and nitric oxide while decreasing the levels of ACE, angiotensin II, and endothelin-1 in Ang I-treated human umbilical vein endothelial cells, indicating the ACE inhibitory effect of LIVGFPAYGH and IVGFPAYGH. In summary, LIVGFPAYGH and IVGFPAYGH from PSRK can be used as functional foods with antihypertensive activity.

Structure based virtual screening, molecular dynamic simulation to identify the oxadiazole derivatives as inhibitors of Enterococcus D-Ala-D-Ser ligase for combating vancomycin resistance

Vancomycin resistance in enterococci mainly arises due to alteration in terminal peptidoglycan dipeptide. A comprehensive structural analysis for substrate specificity of dipeptide modifying d-Alanine: d-Serine ligase (Ddls) is essential to screen its inhibitors for combating vancomycin resistance. In this study modeled 3D structure of EgDdls from E. gallinarum was used for structure based virtual screening (SBVS) of oxadiazole derivatives. Initially, fifteen oxadiazole derivatives were identified as inhibitors at the active site of EgDdls from PubChem database. Further, four EgDdls inhibitors were evaluated using pharmacokinetic profile and molecular docking. The results of molecular docking showed that oxadiazole inhibitors could bind preferentially at ATP binding pocket with the lowest binding energy. Further, molecular dynamics simulation results showed stable behavior of EgDdls in complex with screened inhibitors. The residues Phe172, Lys174, Glu217, Phe292, and Asn302 of EgDdls were mainly involved in interactions with screened inhibitors. Furthermore, MM-PBSA calculation showed electrostatic and van der Waals interactions mainly contribute to overall binding energy. The PCA analysis showed motion of central domain and omega loop of EgDdls. This is involved in the formation of native dipeptide and stabilized after binding of 2-(1-(Ethylsulfonyl) piperidin-4-yl)-5-(furan-2-yl)-1,3,4-oxadiazole, which could be reason for the inhibition of EgDdls. Hence, in this study we have screened inhibitors of EgDdls which could be useful to alleviate the vancomycin resistance problem in enterococci, involved in hospital-acquired infections, especially urinary tract infections (UTI).

Exploration of interaction between angiotensin I-converting enzyme (ACE) and the inhibitory peptide from Wakame (Undaria pinnatifida)

The interaction between angiotensin I-converting enzyme (ACE) and the inhibitory peptide KNFL from Wakame was explored using isothermal titration calorimetry, multiple spectroscopic techniques and molecular dynamics simulations, and an inhibition model was established based on free energy binding theory. The experiments revealed that the binding of KNFL to ACE was a spontaneous exothermic process driven by enthalpy and entropy and occurred via multiple binding sites to form stable complexes. The complexes may be formed through multiple steps of inducing fit and conformational selection. The peptide KNFL had a fluorescence quenching effect on ACE and its addition not only affected the microenvironment around the ACE Trp and Tyr residues, but also increased the diameter and altered the conformation of ACE. This study should prove useful for improving our understanding of the mechanism of ACE inhibitory peptides.

Design and Study of In Silico Binding Dynamics of Certain Isoxazole Bearing Leads Against Aβ-42 and BACE-1 Loop in Protein Fibrillation

Aims:

Design isoxazole bearing leads as dual inhibitors against Amyloid β and BACE-1 loop in protein fibrillation.

Background:

Protein fibrillation is one of the key reasons for several diseases, namely Alzheimer’s, Parkinson’s, and many others. One of the key strategies of preventing protein fibrillation is destabilizing the protein fibrils themselves or inhibiting the amyloid fibril-forming pathway in the initial stage.

Introduction:

Attempts have been taken to design newer leads to inhibit protein fibrillation by targeting the β-amyloidogenesis pathway in the brain. To exploit interfenestration between Amyloid β -42 protein and BACE-1 (β-site amyloid precursor protein cleaving enzyme) for amyloidogenesis, studies are undertaken to design dual inhibitors against the same.

Method:

In vitro binding interactions were found using docking, de novo ligand design, and MD simulation study.

Results:

Three compounds bearing an isoxazole heterocyclic nucleus were designed which could successfully bind to the hydrophobic raft and salt bridge residues Asp 23-Lys-26 of Amyloid β, destabilizing the growing fibril. Additionally, one of our candidate compounds exhibited force of interaction with Thr232 at the S3 pocket of BACE-1, interacted with key residue Asp228, Tyr71, and Thr72 of the β-hairpin flap and hydrogen bonding with Gly11 at loop 10s.

Conclusion:

Protein flexibility dynamics of the Aβ-42 protein revealed that there is a considerable conformational change of the same with or without ligand binding. The lower RMSF of the bound region and reprogramming residual contacts within the Aβ-42 protein suggested successful binding of the ligand with the protein, lowering the access for further β-β dimerization.

Gaussian Accelerated Molecular Dynamics Simulations Investigation on the Mechanism of Angiotensin-Converting Enzyme (ACE) C-Domain Inhibition by Dipeptides

Angiotensin-converting enzyme (ACE)-inhibitory peptides extracted from food proteins can lower blood pressure by inhibiting ACE activity. A recent study showed that the inhibitory activity of IY (Ile-Tyr, a dipeptide derived from soybean protein) against ACE was much higher than that of LL (Leu-Leu), although they had similar hydrophobic and predicted activity values. It was difficult to reveal the deep molecular mechanism underlying this phenomenon by traditional experimental methods. The Apo and two complex systems (i.e., ACE–LL and ACE–IY) were therefore subjected to 1 μs long Gaussian accelerated molecular dynamics (GaMD) simulations. The results showed that the binding of IY can cause obvious contraction of the active site of ACE, mainly manifested by a significant lateral shift of α13, α14, and α15. In addition, hinge 2 and hinge 3 were more stable in the ACE–IY system, while these phenomena were not present in the ACE–LL system. Moreover, the α10 of the IY-bound ACE kept an inward state during the simulation progress, which facilitated the ACE to remain closed. However, for the LL-bound ACE, the α10 switched between two outward states. To sum up, our study provides detailed insights into inhibitor-induced conformational changes in ACE that may help in the design of specific inhibitors targeting ACE for the treatment of hypertension.

A Novel Angiotensin-I-Converting Enzyme (ACE) Inhibitory Peptide from Takifugu flavidus

Alcalase, neutral protease, and pepsin were used to hydrolyze the skin of Takifugu flavidus. The T. flavidus hydrolysates (TFHs) with the maximum degree of hydrolysis (DH) and angiotensin-I-converting enzyme (ACE)-inhibitory activity were selected and then ultra-filtered to obtain fractions with components of different molecular weights (MWs) (<1, 1-3, 3-10, 10-50, and >50 kDa). The components with MWs < 1 kDa showed the strongest ACE-inhibitory activity with a half-maximal inhibitory concentration (IC50) of 0.58 mg/mL. Purification and identification using semi-preparative liquid chromatography, Sephadex G-15 gel chromatography, RP-HPLC, and LC-MS/MS yielded one new potential ACE-inhibitory peptide, PPLLFAAL (non-competitive suppression mode; IC50 of 28 μmmol·L-1). Molecular docking and molecular dynamics simulations indicated that the peptides should bind well to ACE and interact with amino acid residues and the zinc ion at the ACE active site. Furthermore, a short-term assay of antihypertensive activity in spontaneously hypertensive rats (SHRs) revealed that PPLLFAAL could significantly decrease the systolic blood pressure (SBP) and diastolic blood pressure (DBP) of SHRs after intravenous administration. These results suggested that PPLLFAAL may have potential applications in functional foods or pharmaceuticals as an antihypertensive agent.

Heptapeptide Isolated from Isochrysis zhanjiangensis Exhibited Anti-Photoaging Potential via MAPK/AP-1/MMP Pathway and Anti-Apoptosis in UVB-Irradiated HaCaT Cells

Marine microalgae can be used as sustainable protein sources in many fields with positive effects on human and animal health. DAPTMGY is a heptapeptide isolated from Isochrysis zhanjiangensis which is a microalga. In this study, we evaluated its anti-photoaging properties and mechanism of action in human immortalized keratinocytes cells (HaCaT). The results showed that DAPTMGY scavenged reactive oxygen species (ROS) and increase the level of endogenous antioxidants. In addition, through the exploration of its mechanism, it was determined that DAPIMGY exerted anti-photoaging effects. Specifically, the heptapeptide inhibits UVB-induced apoptosis through down-regulation of p53, caspase-8, caspase-3 and Bax and up-regulation of Bcl-2. Thus, DAPTMGY, isolated from I. zhanjiangensis, exhibits protective effects against UVB-induced damage.

A Multi-target Drug Designing for BTK, MMP9, Proteasome and TAK1 for the Clinical Treatment of Mantle Cell Lymphoma

BACKGROUND

Mantle cell lymphoma (MCL) is a type of non-Hodgkin lymphoma characterized by the mutation and overexpression of the cyclin D1 protein by the reciprocal chromosomal translocation t(11;14)(q13:q32).

AIM

The present study aims to identify potential inhibition of MMP9, Proteasome, BTK, and TAK1 and determine the most suitable and effective protein target for the MCL.

METHODOLOGY

Nine known inhibitors for MMP9, 24 for proteasome, 15 for BTK and 14 for TAK1 were screened. SB-3CT (PubChem ID: 9883002), oprozomib (PubChem ID: 25067547), zanubrutinib (PubChem ID: 135565884) and TAK1 inhibitor (PubChem ID: 66760355) were recognized as drugs with high binding capacity with their respective protein receptors. 41, 72, 102 and 3 virtual screened compounds were obtained after the similarity search with compound (PubChem ID:102173753), PubChem compound SCHEMBL15569297 (PubChem ID:72374403), PubChem compound SCHEMBL17075298 (PubChem ID:136970120) and compound CID: 71814473 with best virtual screened compounds.

RESULT

MMP9 inhibitors show commendable affinity and good interaction profile of compound holding PubChem ID:102173753 over the most effective established inhibitor SB-3CT. The pharmacophore study of the best virtual screened compound reveals its high efficacy based on various interactions. The virtual screened compound's better affinity with the target MMP9 protein was deduced using toxicity and integration profile studies.

CONCLUSION

Based on the ADMET profile, the compound (PubChem ID: 102173753) could be a potent drug for MCL treatment. Similar to the established SB-3CT, the compound was non-toxic with LD50 values for both the compounds lying in the same range.

Structural insights and inhibition mechanism of TMPRSS2 by experimentally known inhibitors Camostat mesylate, Nafamostat and Bromhexine hydrochloride to control SARS-coronavirus-2: A molecular modeling approach

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) has been responsible for the cause of global pandemic Covid-19 and to date, there is no effective treatment available. The spike ‘S’ protein of SARS-CoV-2 and ACE2 of the host cell are being targeted to design new drugs to control Covid-19. Similarly, a transmembrane serine protease, TMPRSS2 of the host cell plays a significant role in the proteolytic cleavage of viral ‘S’ protein helpful for the priming of ACE2 receptors and viral entry into human cells. However, three-dimensional structural information and the inhibition mechanism of TMPRSS2 is yet to be explored experimentally. Hence, we have used a molecular dynamics (MD) simulated homology model of TMPRSS2 to study the inhibition mechanism of experimentally known inhibitors Camostat mesylate, Nafamostat and Bromhexine hydrochloride (BHH) using molecular modeling techniques. Prior to docking, all three inhibitors were geometry optimized by semi-empirical quantum chemical RM1 method. Molecular docking analysis revealed that Camostat mesylate and its structural analogue Nafamostat interact strongly with residues His296 and Ser441 present in the catalytic triad of TMPRSS2, whereas BHH binds with Ala386 along with other residues. Comparative molecular dynamics simulations revealed the stable behavior of all the docked complexes. MM-PBSA calculations also revealed the stronger binding of Camostat mesylate to TMPRSS2 active site residues as compared to Nafamostat and BHH. Thus, this structural information could be useful to understand the mechanistic approach of TMPRSS2 inhibition, which may be helpful to design new lead compounds to prevent the entry of SARS-Coronavirus 2 in human cells.

Renin-Angiotensin System and Alzheimer's Disease Pathophysiology: From the Potential Interactions to Therapeutic Perspectives

New roles of the Renin-Angiotensin System (RAS), apart from fluid homeostasis and Blood Pressure (BP) regulation, are being progressively unveiled, since the discoveries of RAS alternative axes and local RAS in different tissues, including the brain. Brain RAS is reported to interact with pathophysiological mechanisms of many neurological and psychiatric diseases, including Alzheimer's Disease (AD). Even though AD is the most common cause of dementia worldwide, its pathophysiology is far from elucidated. Currently, no treatment can halt the disease course. Successive failures of amyloid-targeting drugs have challenged the amyloid hypothesis and increased the interest in the inflammatory and vascular aspects of AD. RAS compounds, both centrally and peripherally, potentially interact with neuroinflammation and cerebrovascular regulation. This narrative review discusses the AD pathophysiology and its possible interaction with RAS, looking forward to potential therapeutic approaches. RAS molecules affect BP, cerebral blood flow, neuroinflammation, and oxidative stress. Angiotensin (Ang) II, via angiotensin type 1 receptors may promote brain tissue damage, while Ang-(1-7) seems to elicit neuroprotection. Several studies dosed RAS molecules in AD patients' biological material, with heterogeneous results. The link between AD and clinical conditions related to classical RAS axis overactivation (hypertension, heart failure, and chronic kidney disease) supports the hypothesized role of this system in AD. Additionally, RAStargeting drugs as Angiotensin Converting Enzyme inhibitors (ACEis) and Angiotensin Receptor Blockers (ARBs) seem to exert beneficial effects on AD. Results of randomized controlled trials testing ACEi or ARBs in AD are awaited to elucidate whether AD-RAS interaction has implications on AD therapeutics.

Proteome-wide mapping and reverse vaccinology-based B and T cell multi-epitope subunit vaccine designing for immune response reinforcement against Porphyromonas gingivalis

Porphyromonas gingivalis, a prominent pathogen responsible for acute periodontal diseases, is widely studied by the scientific community for its successful evasion of the host immune system. P. gingivalis is associated with rheumatoid arthritis, dementia, and Alzheimer's. The pathogen successfully survives itself against the heavy load of conventional antibiotics because of its ability to evade the host immune system. Subtractive proteomics and reverse vaccinology approaches were employed in order to prioritize the best proteins for vaccine designing. Three vaccine candidates with Uniprot ID: Q7MWZ2 (histidine Kinase), Q7MVL1 (Fe (2+) transporter), and Q7MWZ2 (Capsular polysaccharide transport protein) were identified for vaccine designing. These proteins are antigenic and essential for pathogen survival. A wide range of immunoinformatics tools was applied for the prediction of epitopes, B, and T cells, for the vaccine candidate proteins. Molecular docking of the predicted epitopes against the MHC molecules were carried out. In-silico vaccine was constructed using carefully evaluated epitopes and consequently modeled for docking with human Toll-like receptor 2. Chain C of Pam3CSK4 (PDB ID; 2Z7X) was linked to the vaccine as an adjuvant to boost immune response towards the vaccine. For stability evaluation of the vaccine-TLR-2 docked complex, Molecular Dynamics simulations were performed. The reverse-translated nucleotide sequence cloned in Eschericia coli to attain the maximal expression of the vaccine protein. The maximal expression was ensured by CAI score of 0.96. The current vaccine requires future experimental validation to confirm its effectiveness. The vaccine developed will be helpful to protect against P. gingivalis associated infections.Communicated by Ramaswamy H. Sarma.

Deciphering the Inhibition Mechanism of under Trial Hsp90 Inhibitors and Their Analogues: A Comparative Molecular Dynamics Simulation

Heat shock protein 90 (Hsp90) performs functions in cellular activities together with other signaling pathways. Hsp90 is evolutionarily conserved and universally articulated as a human cancer-causing agent involved in lung cancer and breast cancer followed by colon and rectum cancers. It has emerged as an effective drug candidate, and inhibition may affect several signaling pathways associated with cancer spread. Therefore, in-silico approaches, molecular docking, molecular dynamics simulation, and binding free energy calculations were applied to create insights into the inhibition mechanism against Hsp90 to identify new cancer therapeutic drugs. Top-docked Hsp90-inhibitor complexes with their analogues were selected as the best complexes based on the GOLD fitness score and orientation. The significant interaction of Hsp90 inhibitors and their analogues were observed to be bound with active site residues as well as residing within the same cavity region. System stability factors RMSD, RMSF, beta-factor, and radius of gyration were analyzed for top-docked complexes and ensure strong binding interaction between inhibitors and the Hsp90 cavity. Cavity bound inhibitors were found to retain consistent hydrogen bonding during the simulation. The radial distribution function (RDF) illustrated that interacting active site residues drive the binding and stability of the inhibitors. Similarly, the axial frequency distribution, which is an indigenously developed analytical tool, produced noteworthy knowledge of the hydrogen-bonding pattern. Results yielded new insights into the design of cancer therapeutic drugs against Hsp90. This finding suggests that under trial Hsp90 inhibitors MPC-3100 could be a potential starting point into the development of potential anticancer agents with the possibility of future directions for the improvement of early existing Hsp90 inhibitors CNF-2024 and SNX-5422 as an anticancer agent.

Spontaneous Hinge-Bending Motions of Angiotensin I Converting Enzyme: Role in Activation and Inhibition

The inhibition of human angiotensin I converting enzyme (ACE) has been regarded as a promising approach for the treatment of hypertension. Despite research attempts over many years, our understanding the mechanisms of activation and inhibition of ACE is still far from complete. Here, we present results of all atom molecular dynamics simulations of ACE with and without ligands. Two types of inhibitors, competitive and mixed non-competitive, were used to model the ligand bound forms. In the absence of a ligand the simulation showed spontaneous large hinge-bending motions of multiple conversions between the closed and open states of ACE, while the ligand bound forms were stable in the closed state. Our simulation results imply that the equilibrium between pre-existing backbone conformations shifts in the presence of a ligand. The hinge-bending motion of ACE is considered as an essential to the enzyme function. A mechanistic model of activation and the inhibition may provide valuable information for novel inhibitors of ACE.

Structural and molecular bases of angiotensin-converting enzyme inhibition by bovine casein-derived peptides: an in silico molecular dynamics approach

The angiotensin-converting enzyme (ACE) plays a key role in blood pressure regulation process, and its inhibition is one of the main drug targets for the treatment of hypertension. Though various peptides from milk proteins are well-known for their ACE-inhibitory capacity, research devoted to understand the molecular bases of such property remain scarce, specifically for such peptides. Therefore, in this work, computational molecular docking and molecular dynamics calculations were performed to enlighten the intermolecular interactions involved in ACE inhibition by six different casein-derived peptides (FFVAPFPEVFGK, FALPQYLK, ALNEINQFYQK, YLGYLEQLLR, HQGLPQEVLNENLLR and NAVPITPTLNR). Two top ranked docking poses for each peptide (one with N- and the other C-terminal peptide extremity oriented towards the ACE active site) were selected for dynamic simulations (50 ns; GROMOS53A6 force field), and the results were correlated to in vitro ACE inhibition capacity. Two molecular features appeared to be essential for peptides to present high ACE inhibition capacity in vitro: i) to interact with the S1 active site residues (Ala354, Glu384, and Tyr523) by hydrogen bonds; ii) to interact with Zn²⁺ coordinated residues (His383, His387, and Glu411) by short-lenght hydrogen bonds, as observed in the cases of ALNEINQFYQK (IACE = 80.7%), NAVPITPTLNR (IACE = 80.7%), and FALPQYLK (IACE = 79.0%). Regardless of the temporal stability of these strong interactions, they promoted some disruption of Zn²⁺ tetrahedral coordination during the molecular dynamics trajectories, and were pointed as the main reason for the greatest ACE inhibition by these peptides. On the other hand, peptides with intermediate inhibition capacity (50% < I_ACE < 45%) interacted mainly by weaker interactions (e.g.: electrostatic and hydrophobic) with the Zn²⁺ coordinated residues, and were not able to change significantly its tetrahedral coordination structure. These findings may: i) assist the discrimination in silico of "good" and "bad" ACE-inhibitory peptides from other food sources, and/or ii) aid in designing de novo new molecules with ACE-inhibitory capacity. Communicated by Ramaswamy Sarma.

Antihypertensive agents in Alzheimer's disease: beyond vascular protection

Introduction: Midlife hypertension has been consistently linked with increased risk of cognitive decline and Alzheimer's disease (AD). Observational studies and randomized trials show that the use of antihypertensive therapy is associated with a lesser incidence or prevalence of cognitive impairment and dementia. However, whether antihypertensive agents specifically target the pathological process of AD remains elusive.Areas covered: This review of literature provides an update on the clinical and preclinical arguments supporting anti-AD properties of antihypertensive drugs. The authors focused on validated all classes of antihypertensive treatments such as angiotensin-converting enzyme inhibitors (ACEi), angiotensin receptor blockers (ARB), calcium channel blockers (CCB), β-blockers, diuretics, neprilysin inhibitors, and other agents. Three main mechanisms can be advocated: action on the concurrent vascular pathology, action on the vascular component of Alzheimer's pathophysiology, and action on nonvascular targets.Expert opinion: In 2019, while there is no doubt that hypertension should be treated in primary prevention of vascular disease and in secondary prevention of stroke and mixed dementia, the place of antihypertensive agents in the secondary prevention of 'pure' AD remains an outstanding question.

Novel ACE inhibitory tripeptides from ovotransferrin using bioinformatics and peptidomics approaches

Food-derived ACE inhibitory peptides have recently attracted increased attention. This work focused on a more efficient in silico method to find ACE inhibitory peptides from ovotransferrin. In this work, ovotransferrin was digested into peptides by virtual enzymolysis. Subsequently, in vitro ACE inhibitory activity of potential tripeptides was conducted following the peptide score, toxicity, and water solubility prediction. Both pharmacophore study and flexible docking were applied to analyze ACE inhibition mechanism of tripeptides. Our results demonstrated that EWL was a potent ACE inhibitory tripeptide with IC₅₀ value of 380 ± 10 μM. Besides, pharmacophore and flexible docking showed that the pi interaction and hydrogen bond were the key interactions in ACE-EWL complex. It appears that the in vitro ACE inhibitory activity of tripeptide EWL was consistent with its molecular modeling.

Dipeptidyl peptidase-4 inhibitors can inhibit angiotensin converting enzyme

Angiotensin-1 converting enzyme inhibitors (ACEIs) improve insulin sensitivity. Inhibitors of dipeptidyl peptidase-4 (DPP-4) are anti-diabetic drugs with several cardio-renal effects. Both ACE and DPP-4 share common features. Thus, we tested if they could be inhibited by one inhibitor. First, in silico screening was used to investigate the ability of different DPP-4 inhibitors or ACEIs to interact with DPP-4 and ACE. The results of screening were then extrapolated into animal study. Fifty Sprague Dawley rats were randomly assigned into 5 groups treated with vehicle, captopril, enalapril, linagliptin or sitagliptin. Both low and high doses of each drug were tested. Baseline blood samples and samples at days 1, 8, 10, 14 were used to measure plasma DPP-4 and ACE activities and angiotensin II levels. Active glucagon-like peptide-1 (GLP-1) levels were measured after oral glucose challenge. All tested DPP-4 inhibitors could interact with ACE at a relatively reasonable binding energy while most of the ACEIs only interacted with DPP-4 at a predicted high inhibition constant. In rats, high dose of sitagliptin was able to inhibit ACE activity and reduce angiotensin II levels while linagliptin had only a mild effect. ACEIs did not significantly affect DPP-4 activity or prevent GLP-1 degradation. It seems that some DPP-4 inhibitors could inhibit ACE and this could partially explain the cardio-renal effects of these drugs. Further studies are required to determine if such inhibition could take place in clinical settings.

The potential peptides against angiotensin-I converting enzyme through a virtual tripeptide-constructing library

Peptides derived from food proteins are promising bioactive source for inhibiting Angiotensin-I converting enzyme (ACE) activity. Bioactive peptides (BP) have received much attention, particularly from the pharmaceutical industry. As they not only own potent properties but also possess less side-effects than synthetic drugs. In this work, an 8000 possible tripeptides library was constructed to predict the potential ACE inhibitory peptides by using in silico tools. GOLD molecular docking was then applied to determine the binding mode of action between ACE and each of tripeptide from this in-house library. The first 662 high-ranking tripeptides by ChemScore were chosen to create association rules of tripeptides-ACE complexes. An orientation pattern of amino acid in the binding tunnel of ACE has been examined by frequency analysis. The association rules (confident values over 90%) illustrated that hydrophobic factor has been displayed as main components in the ACE tripeptides inhibitor from four factors in equation, hydrophobic, aromatic, polar, charged. According to in silico output, five tripeptides were chosen to test in vitro study of ACE-inhibitory activity. The half-maximal inhibitory concentration (IC₅₀) of these selective five peptides, WCW, IWW, WWW, WWI and WLW for inhibiting ACE were 49.50 ± 3.88 μM, 489.14 ± 8.84 μM, 536.02 ± 38.57 μM, 752.91 ± 41.89 μM and 1783 ± 0.113 μM, respectively. Molecular dynamics simulations approach was applied to study the interaction of WCW (Trp-Cys-Trp) within ACE pocket site. This ligand was stabilized by strong hydrogen bonding interactions with ACE active site, Tyr523-Trp'1 (99.76%) and His353-Trp'1 (95.68%). Our computational protocol could be considered as a new tool for identifying active peptide against ACE from hydrolysated peptides of natural sources.

Production of antioxidant and ACE-inhibitory peptides from Kluyveromyces marxianus protein hydrolysates: Purification and molecular docking

Kluyveromyces marxianus protein hydrolysates were prepared by two different sonicated-enzymatic (trypsin and chymotrypsin) hydrolysis treatments to obtain antioxidant and ACE-inhibitory peptides. Trypsin and chymotrypsin hydrolysates obtained by 5 h, exhibited the highest antioxidant and ACE-inhibitory activities. After fractionation using ultra-filtration and reverse phase high performance liquid chromatography (RP-HPLC) techniques, two new peptides were identified. One fragment (LL-9, MW = 1180 Da) with the amino acid sequence of Leu-Pro-Glu-Ser-Val-His-Leu-Asp-Lys showed significant ACE inhibitory activity (IC₅₀ = 22.88 μM) while another peptide fragment (VL-9, MW = 1118 Da) with the amino acid sequence of Val-Leu-Ser-Thr-Ser-Phe-Pro-Pro-Lys showed the highest antioxidant and ACE inhibitory properties (IC₅₀ = 15.20 μM, 5568 μM TE/mg protein). The molecular docking studies revealed that the ACE inhibitory activities of VL-9 is due to interaction with the S2 (His513, His353, Glu281) and S’1 (Glu162) pockets of ACE and LL-9 can fit perfectly into the S1 (Thr345) and S2 (Tyr520, Lys511, Gln281) pockets of ACE.

Virtual screening following rational drug design based approach for introducing new anti amyloid beta aggregation agent

Amyloid β (Aβ) sheets aggregations is the main reason of Alzheimer disease. The interacting areas between monomers are residue number 38 to 42. Inhibition of interaction between Aβ molecules prevents plaque formation. In the present study, we have performed a high-throughput virtual screening among ZINC database and top 1000 hits were checked again regarding binding affinity by AutoDock software. Top 4 successive second step screening hits was considered for drug design purpose against aggregation site of Aβ molecules. The toxicity and pharmacological properties of new designed ligands was assessed by PROTOX and FAFdrugs3 webservers. Several steps of modifications performed in the structures of hit#1 and hit#2 and finally new designed ligand based on hit 1, 1-RD-3 (3-[(Z)-6-Hydroxy-4-{[5-(2-methoxyethyl)-6-methyltetrahydro-2H-pyran-2-yl]methyl}-1-methyl-3-hexenyloxy]tetrahydro-2Hpyran- 4-ol) and a designed ligand based on hit 2, 2-RD-2 (6-(Hydroxymethyl)-4-{5-hydroxy-6-methyl-4-[(3- methylcyclohexyl)methyl]tetrahydro-2H-pyran-2-yloxy}tetrahydro-2H-pyran-2,3,5-triol) could successfully pass pharmacological filters. The LD50 of 37000 mg/kg for 1-RD-3 and 2000 mg/kg for 2-RD-2 indicates that the designed ligands can be considered as new candidates for anti Aβ aggregation to treat Alzheimer's disease. Interestingly, after performing several modification steps still a considerable binding affinity of -9.3 kcal/mol for 1-RD-3 and -9.8 kcal/mol for 2-RD-2 still remained. Theoretically, the new designed molecules can reduce the deposition of Aβ in the cerebral cortex and as the results the Alzheimer symptoms could be decreased.

Molecular modeling approach to explore the role of cathepsin B from Hordeum vulgare in the degradation of Aβ peptides

The pathological hallmark of Alzheimer's disease is the accumulation of Aβ peptides in human brains. These Aβ peptides can be degraded by several enzymes such as hACE, hECE, hIDE and cathepsin B. Out of which cathepsin B also belongs to the papain super family and has been found in human brains, it has a role in Aβ peptide degradation through limited proteolysis. The Aβ concentrations are maintained properly by its production and clearance via receptor-mediated cellular uptake and direct enzymatic degradation. However, the reduced production of Aβ degrading enzymes as well as their Aβ degrading activity in human brains initiate the process of accumulation of Aβ peptides. So it becomes essential to investigate the molecular interactions involved in the process of Aβ degradation in detail at the atomic level. Hence, homology modeling, molecular docking and molecular dynamics simulation techniques have been used to explore the possible role of cathepsin B from Hordeum vulgare in the degradation of amyloid beta (Aβ) peptides. The homology model of cathepsin B from Hordeum vulgare shows good similarity with human cathepsin B. Molecular docking and MD simulation results revealed that the active site residues Cys32, HIS112, HIS113 are involved in the catalytic activity of cathepsin B. The sulfhydryl group of the Cys32 residue of cathepsin B from Hordeum vulgare cleaves the Aβ peptide from the carboxylic end of Glu11. Hence, this structural study might be helpful in designing alternative strategies for the treatment of AD.

Structural Characterizations of the Fas Receptor and the Fas-Associated Protein with Death Domain Interactions

The Fas receptor is a representative death receptor, and the Fas-associated protein with death domain (FADD) is a crucial adapter protein needed to support the Fas receptor's activity. The Fas-FADD interactions constitute an important signaling pathway that ultimately induces apoptosis or programmed cell death in biological systems. The interactions responsible for this cell-death process are governed by the binding process of the Fas ligand to the Fas, followed by the caspase cascade activation. Using a computational approach, the present communication explores certain essential structural aspects of the Fas-FADD death domains and their interfacial interactions.

Conformational transition pathway in the inhibitor binding process of human monoacylglycerol lipase

Human monoacylglycerol lipase (MGL) catalyzes the hydrolysis of 2-arachidonoylglycerol to arachidonic and glycerol, which plays a pivotal role in the normal biological processes of brain. Co-crystal structure of the MGL in complex with its inhibitor, compound 1, shows that the helix α4 undergoes large-scale conformational changes in response to the compound 1 binding compared to the apo MGL. However, the detailed conformational transition pathway of the helix α4 in the inhibitor binding process of MGL has remained unclear. Here, conventional molecular dynamics (MD) and nudged elastic band (NEB) simulations were performed to explore the conformational transition pathway of the helix α4. Conventional MD simulations unveiled that the compound 1 induced the closed conformation of the active site of MGL, reduced the conformational flexibility of the helix α4, and elicited the large-scale conformational rearrangement of the helix α4, leading to the complete folding of the helix α4. Moreover, NEB simulations revealed that the conformational transition pathway of helix α4 underwent an almost 180° counter-clockwise rotation of the helix α4. Our computational results advance the structural and mechanistic understanding of the inhibitory mechanism.

Amyloid cascade hypothesis: Pathogenesis and therapeutic strategies in Alzheimer's disease

Alzheimer's disease is an irreversible, progressive neurodegenerative disorder. Various therapeutic approaches are being used to improve the cholinergic neurotransmission, but their role in AD pathogenesis is still unknown. Although, an increase in tau protein concentration in CSF has been described in AD, but several issues remains unclear. Extensive and accurate analysis of CSF could be helpful to define presence of tau proteins in physiological conditions, or released during the progression of neurodegenerative disease. The amyloid cascade hypothesis postulates that the neurodegeneration in AD caused by abnormal accumulation of amyloid beta (Aβ) plaques in various areas of the brain. The amyloid hypothesis has continued to gain support over the last two decades, particularly from genetic studies. Therefore, current research progress in several areas of therapies shall provide an effective treatment to cure this devastating disease. This review critically evaluates general biochemical and physiological functions of Aβ directed therapeutics and their relevance.

Evaluation of a possible role of Stigmatella aurantiaca ACE in Aβ peptide degradation: a molecular modeling approach

Amyloid-β (Aβ)-degrading enzymes are known to degrade Aβ peptides, a causative agent of Alzheimer's disease. These enzymes are responsible for maintaining Aβ concentration. However, loss of such enzymes or their Aβ-degrading activity because of certain genetic as well as nongenetic reasons initiates the accumulation of Aβ peptides in the human brain. Considering the limitations of the human enzymes in clearing Aβ peptide, the search for microbial enzymes that could cleave Aβ is necessary. Hence, we built a three-dimensional model of angiotensin-converting enzyme (ACE) from Stigmatella aurantiaca using homology modeling technique. Molecular docking and molecular dynamics simulation techniques were used to outline the possible cleavage mechanism of Aβ peptide. These findings suggest that catalytic residue Glu 434 of the model could play a crucial role to degrade Aβ peptide between Asp 7 and Ser 8. Thus, ACE from S. aurantiaca might cleave Aβ peptides similar to human ACE and could be used to design new therapeutic strategies against Alzheimer's disease.

The impact of molecular dynamics on drug design: applications for the characterization of ligand-macromolecule complexes

Among all tools available to design new drugs, molecular dynamics (MD) simulations have become an essential technique. Initially developed to investigate molecular models with a limited number of atoms, computers now enable investigations of large macromolecular systems with a simulation time reaching the microsecond range. The reviewed articles cover four years of research to give an overview on the actual impact of MD on the current medicinal chemistry landscape with a particular emphasis on studies of ligand-protein interactions. With a special focus on studies combining computational approaches with data gained from other techniques, this review shows how deeply embedded MD simulations are in drug design strategies and articulates what the future of this technique could be.

Molecular basis of laccase bound to lignin: insight from comparative studies on the interaction of Trametes versicolor laccase with various lignin model compounds

Binding orientation of lignin model compounds in laccase.

Insight into molecular interactions of Aβ peptide and gelatinase from Enterococcus faecalis: a molecular modeling approach

Alzheimer's disease is characterized by the presence of extracellular deposition of amyloid beta (Aβ) peptides.

AutoDock4(Zn): an improved AutoDock force field for small-molecule docking to zinc metalloproteins

Zinc is present in a wide variety of proteins and is important in the metabolism of most organisms. Zinc metalloenzymes are therapeutically relevant targets in diseases such as cancer, heart disease, bacterial infection, and Alzheimer's disease. In most cases a drug molecule targeting such enzymes establishes an interaction that coordinates with the zinc ion. Thus, accurate prediction of the interaction of ligands with zinc is an important aspect of computational docking and virtual screening against zinc containing proteins. We have extended the AutoDock force field to include a specialized potential describing the interactions of zinc-coordinating ligands. This potential describes both the energetic and geometric components of the interaction. The new force field, named AutoDock4Zn, was calibrated on a data set of 292 crystal complexes containing zinc. Redocking experiments show that the force field provides significant improvement in performance in both free energy of binding estimation as well as in root-mean-square deviation from the crystal structure pose. The new force field has been implemented in AutoDock without modification to the source code.

Homology modeling, molecular docking and MD simulation studies to investigate role of cysteine protease from Xanthomonas campestris in degradation of Aβ peptide

Cysteine protease is known to degrade amyloid beta peptide which is a causative agent of Alzheimer's disease. This cleavage mechanism has not been studied in detail at the atomic level. Hence, a three-dimensional structure of cysteine protease from Xanthomonas campestris was constructed by homology modeling using Geno3D, SWISS-MODEL, and MODELLER 9v7. All the predicted models were analyzed by PROCHECK and PROSA. Three-dimensional model of cysteine protease built by MODELLER 9v7 shows similarity with human cathepsin B crystal structure. This model was then used further for docking and simulation studies. The molecular docking study revealed that Cys17, His87, and Gln88 residues of cysteine protease form an active site pocket similar to human cathepsin B. Then the docked complex was refined by molecular dynamic simulation to confirm its stable behavior over the entire simulation period. The molecular docking and MD simulation studies showed that the sulfhydryl hydrogen atom of Cys17 of cysteine protease interacts with carboxylic oxygen of Lys16 of Aβ peptide indicating the cleavage site. Thus, the cysteine protease model from X. campestris having similarity with human cathepsin B crystal structure may be used as an alternate approach to cleave Aβ peptide a causative agent of Alzheimer's disease.