Inhibition of matrix metalloproteinases by HIV-1 integrase strand transfer inhibitors

More than fifteen million women with the human immunodeficiency virus type-1 (HIV-1) infection are of childbearing age world-wide. Due to improved and affordable access to antiretroviral therapy (ART), the number of in utero antiretroviral drug (ARV)-exposed children has exceeded a million and continues to grow. While most recommended ART taken during pregnancy suppresses mother to child viral transmission, the knowledge of drug safety linked to fetal neurodevelopment remains an area of active investigation. For example, few studies have suggested that ARV use can be associated with neural tube defects (NTDs) and most notably with the integrase strand transfer inhibitor (INSTI) dolutegravir (DTG). After risk benefit assessments, the World Health Organization (WHO) made recommendations for DTG usage as a first and second-line preferred treatment for infected populations including pregnant women and those of childbearing age. Nonetheless, long-term safety concerns remain for fetal health. This has led to a number of recent studies underscoring the need for biomarkers to elucidate potential mechanisms underlying long-term neurodevelopmental adverse events. With this goal in mind, we now report the inhibition of matrix metalloproteinases (MMPs) activities by INSTIs as an ARV class effect. Balanced MMPs activities play a crucial role in fetal neurodevelopment. Inhibition of MMPs activities by INSTIs during neurodevelopment could be a potential mechanism for adverse events. Thus, comprehensive molecular docking testing of the INSTIs, DTG, bictegravir (BIC), and cabotegravir (CAB), against twenty-three human MMPs showed broad-spectrum inhibition. With a metal chelating chemical property, each of the INSTI were shown to bind Zn++ at the MMP's catalytic domain leading to MMP inhibition but to variable binding energies. These results were validated in myeloid cell culture experiments demonstrating MMP-2 and 9 inhibitions by DTG, BIC and CAB and even at higher degree than doxycycline (DOX). Altogether, these data provide a potential mechanism for how INSTIs could affect fetal neurodevelopment.

HIV-Tat Exacerbates the Actions of Atazanavir, Efavirenz, and Ritonavir on Cardiac Ryanodine Receptor (RyR2)

The incidence of sudden cardiac death (SCD) in people living with HIV infection (PLWH), especially those with inadequate viral suppression, is high and the reasons for this remain incompletely characterized. The timely opening and closing of type 2 ryanodine receptor (RyR2) is critical for ensuring rhythmic cardiac contraction-relaxation cycles, and the disruption of these processes can elicit Ca²⁺ waves, ventricular arrhythmias, and SCD. Herein, we show that the HIV protein Tat (HIV-Tat: 0-52 ng/mL) and therapeutic levels of the antiretroviral drugs atazanavir (ATV: 0-25,344 ng/mL), efavirenz (EFV: 0-11,376 ng/mL), and ritonavir (RTV: 0-25,956 ng/mL) bind to and modulate the opening and closing of RyR2. Abacavir (0-14,315 ng/mL), bictegravir (0-22,469 ng/mL), Rilpivirine (0-14,360 ng/mL), and tenofovir disoproxil fumarate (0-18,321 ng/mL) did not alter [³H]ryanodine binding to RyR2. Pretreating RyR2 with low HIV-Tat (14 ng/mL) potentiated the abilities of ATV and RTV to bind to open RyR2 and enhanced their ability to bind to EFV to close RyR2. In silico molecular docking using a Schrodinger Prime protein-protein docking algorithm identified three thermodynamically favored interacting sites for HIV-Tat on RyR2. The most favored site resides between amino acids (AA) 1702-1963; the second favored site resides between AA 467-1465, and the third site resides between AA 201-1816. Collectively, these new data show that HIV-Tat, ATV, EFV, and RTV can bind to and modulate the activity of RyR2 and that HIV-Tat can exacerbate the actions of ATV, EFV, and RTV on RyR2. Whether the modulation of RyR2 by these agents increases the risk of arrhythmias and SCD remains to be explored.

Targeting SARS-CoV-2 Receptor Binding Domain with Stapled Peptides: An In Silico Study

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved into a pandemic of unprecedented scale. This coronavirus enters cells by the interaction of the receptor binding domain (RBD) with the human angiotensin-converting enzyme 2 receptor (hACE2). In this study, we employed a rational structure-based design to propose 22-mer stapled peptides using the structure of the hACE2 α1 helix as a template. These peptides were designed to retain the α-helical character of the natural structure, to enhance binding affinity, and to display a better solubility profile compared to other designed peptides available in the literature. We employed different docking strategies (PATCHDOCK and ZDOCK) followed by a double-step refinement process (FIBERDOCK) to rank our peptides, followed by stability analysis/evaluation of the interaction profile of the best docking predictions using a 500 ns molecular dynamics (MD) simulation, and a further binding affinity analysis by molecular mechanics with generalized Born and surface area (MM/GBSA) method. Our most promising stapled peptides presented a stable profile and could retain important interactions with the RBD in the presence of the E484K RBD mutation. We predict that these peptides can bind to the viral RBD with similar potency to the control NYBSP-4 (a 30-mer experimentally proven peptide inhibitor). Furthermore, our study provides valuable information for the rational design of double-stapled peptide as inhibitors of SARS-CoV-2 infection.

Small molecule binding to inhibitor of nuclear factor kappa-B kinase subunit beta in an ATP non-competitive manner

Inhibitor of nuclear factor kappa-B kinase subunit beta (IKKβ) is a key regulator of the cannonical NF-κB pathway. IKKβ has been validated as a drug target for pathological conditions, which include chronic inflammatory diseases and cancer. Pharmacological studies revealed that chronic administration of ATP-competitive IKKβ inhibitors resulted in unexpected toxicity. We previously reported the discovery of 13-197 as a non-toxic IKKβ inhibitor that reduced tumor growth. Here, we show that 13-197 inhibits IKKβ in a ATP non-competitive manner and an allosteric pocket at the interface of the kinase and ubiquitin like domains was identified as the potential binding site.

Targeting 14-3-3ε activates apoptotic signaling to prevent cutaneous squamous cell carcinoma

More than a million cases of cutaneous squamous cell carcinoma are diagnosed in the USA each year, and its incidence is increasing. Most of these malignancies arise from premalignant lesions, providing an opportunity for intervention before malignant progression. We previously documented how cytoplasmic mislocalization of CDC25A in premalignant and malignant skin cancers confers resistance to apoptotic cell death via a mechanism that depends on its interaction with 14-3-3ε. From these data, we hypothesized that 14-3-3ε overexpression drives skin tumor development and progression, such that targeting 14-3-3ε may be a useful strategy for skin cancer treatment. Like CDC25A, 14-3-3ε was overexpressed and mislocalized to the cytoplasm of both benign and malignant human skin cancer. Skin-targeted deletion of the 14-3-3ε gene reduced skin tumor development by 75% and blocked malignant progression. 14-3-3ε suppressed apoptosis through activation of Akt, leading to inhibition of BCL2 associated agonist of cell death and upregulation of Survivin. Using virtual tetrapeptide libraries, we developed a novel peptide that specifically blocked 14-3-3ε heterodimerization and thereby prevented its interaction with CDC25A. The peptide reduced prosurvival signaling, killed skin cancer cells and reduced skin tumor growth in xenograft. Normal skin keratinocytes were unaffected by inhibition or deletion of 14-3-3ε. Thus, targeting of 14-3-3ε dimerization is a promising strategy for the treatment of premalignant skin lesions.

Peptide truncation leads to a twist and an unusual increase in affinity for casitas B-lineage lymphoma tyrosine kinase binding domain

We describe truncation and SAR studies to identify a pentapeptide that binds Cbl tyrosine kinase binding domain with a higher affinity than the parental peptide. The pentapeptide has an alternative binding mode that allows occupancy of a previously uncharacterized groove. A peptide library was used to map the binding site and define the interface landscape. Our results suggest that the pentapeptide is an ideal starting point for the development of inhibitors against Cbl driven diseases.

Hexapeptide fragment of carcinoembryonic antigen which acts as an agonist of heterogeneous ribonucleoprotein M

Colorectal cancers with metastatic potential secrete the glycoprotein carcinoembryonic antigen (CEA). CEA has been implicated in colorectal cancer metastasis by inducing Kupffer cells to produce inflammatory cytokines which, in turn, make the hepatic micro-environment ideal for tumor cell implantation. CEA binds to the heterogeneous ribonucleoprotein M (hnRNP M) which acts as a cell surface receptor in Kupffer cells. The amino acid sequence in CEA, which binds the hnRNP M receptor, is Tyr-Pro-Glu-Leu-Pro-Lys. In this study, the structure of Ac-Tyr-Pro-Glu-Leu-Pro-Lys-NH₂ (YPELPK) was investigated using electronic circular dichroism, vibrational circular dichroism, and molecular dynamics simulations. The binding of the peptide to hnRNP M was also investigated using molecular docking calculations. The biological activity of YPELPK was studied using differentiated human THP-1 cells, which express hnRNP M on their surface and secrete IL-6 when stimulated by CEA. YPELPK forms a stable polyproline-II helix and stimulates IL-6 production of THP-1 cells at micromolar concentrations.

Exploiting the P-1 pocket of BRCT domains toward a structure guided inhibitor design

Breast cancer gene 1 carboxy terminus (BRCT) domains are found in a number of proteins that are important for DNA damage response (DDR). The BRCT domains bind phosphorylated proteins and these protein-protein interactions are essential for DDR and DNA repair. High affinity domain specific inhibitors are needed to facilitate the dissection of the protein-protein interactions in the DDR signaling. The BRCT domains of BRCA1 bind phosphorylated protein through a pSXXF consensus recognition motif. We identified a hydrophobic pocket at the P-1 position of the pSXXF binding site. Here we conducted a structure-guided synthesis of peptide analogs with hydrophobic functional groups at the P-1 position. Evaluation of these led to the identification of a peptide mimic 15 with a inhibitory constant (K(i)) of 40 nM for BRCT(BRCA1). Analysis of the TopBP1 and MDC1 BRCT domains suggests a similar approach is viable to design high affinity inhibitors.

Beyond the frog: the evolution of homology models of human IKKβ

Nuclear Factor κ B is implicated in tumor progression and chronic inflammatory diseases and is regulated by IκB kinase β (IKKβ). The crystal structure of IKKβ has been recently solved for Xenopus laevis. Homology models of human IKKβ have been developed prior to and after the crystal structure was solved. Here, we compare four models of human IKKβ and evaluate their performance in both broad and focused library docking studies.

Single-molecule atomic force microscopy force spectroscopy study of Aβ-40 interactions

Misfolding and aggregation of amyloid β-40 (Aβ-40) peptide play key roles in the development of Alzheimer's disease (AD). However, very little is known about the molecular mechanisms underlying these molecular processes. We developed a novel experimental approach that can directly probe aggregation-prone states of proteins and their interactions. In this approach, the proteins are anchored to the surface of the atomic force microscopy substrate (mica) and the probe, and the interaction between anchored molecules is measured in the approach-retraction cycles. We used dynamic force spectroscopy (DFS) to measure the stability of transiently formed dimers. One of the major findings from DFS analysis of α-synuclein (α-Syn) is that dimeric complexes formed by misfolded α-Syn protein are very stable and dissociate over a range of seconds. This differs markedly from the dynamics of monomers, which occurs on a microsecond to nanosecond time scale. Here we applied the same approach to quantitatively characterize interactions of Aβ-40 peptides over a broad range of pH values. These studies showed that misfolded dimers are characterized by lifetimes in the range of seconds. This value depends on pH and varies between 2.7 s for pH 2.7 and 0.1 s for pH 7, indicating that the aggregation properties of Aβ-40 are modulated by the environmental conditions. The analysis of the contour lengths revealed the existence of various pathways for dimer dissociation, suggesting that dimers with different conformations are formed. These structural variations result in different aggregation pathways, leading to different types of oligomers and higher-order aggregates, including fibrils.

Calculation of weakly polar interaction energies in polypeptides using density functional and local Møller-Plesset perturbation theory

The interaction energies of ubiquitous weakly polar interactions in proteins are comparable with those of hydrogen bonds, consequently, they stabilize local, secondary, and tertiary structures. However, the most widely-used density functionals fail to describe the weakly polar interactions. Thus, it is important to find and test functionals which adequately describe and quantify the energetics of such interactions. For this purpose, interaction energies in the hydrophobic core of rubredoxin (PDB id: 1rb9) and in the S22 subset of the JSCH-2005 benchmark database were computed with the BHandHLYP and PWPW91 functionals and with the pseudospectral implementation of the local MP2 (PS-LMP2) method. The cc-pVDZ, cc-pVTZ(-f), cc-pVTZ, cc-pVQZ(-g), aug-cc-pVDZ, aug-cc-VTZ(-f), and aug-cc-pVTZ basis sets were used for the calculations. In the S22 subset the PS-LMP2 results were extrapolated to the complete basis set limit. Furthermore, the a posteriori counterpoise method of Boys and Bernardi was used to correct the basis set superposition errors in the calculation of interaction energies. Calculations using the BHandHLYP functional, both for the various weakly polar interactions in rubredoxin and for the dispersion interactions in the S22 subset, were in good agreement with those using the coupled cluster (CCSD(T)) and the resolution of identity MP2 (RIMP2) methods and clearly outperformed both the PWPW91 functional and the PS-LMP2 method. The results for the S22 hydrogen bonded subset, obtained with PWPW91 calculations, were closest to those of the reference high level calculations. For the "mixed" (hydrogen bonded and dispersive) interactions in the S22 subset, results obtained with the BHandHLYP and PS-LMP2 calculations agreed well with the reference calculations.

Quantum chemical quantification of weakly polar interaction energies in the TC5b miniprotein

The tertiary structure of the TC5b miniprotein is stabilized by inter-residue interactions of the Trp-cage, which is composed of a Tyr and several Pro residues surrounding a central Trp residue. The interactions include Ar-Ar (aromatic side-chain-aromatic side-chain), Ar-NH (aromatic side-chain-backbone amide), and CH-pi (aromatic side-chain-aliphatic hydrogen) interactions. In the present work, the strength of the weakly polar interactions found in the TC5b miniprotein was quantified using all of the available 38 NMR structures (1L2Y) from the Protein Data Bank with DFT quantum chemical calculations at the BHandHLYP/cc-pVTZ level of theory and molecular fragmentation with capping of the partial structures. The energies of interaction between the individual residues of the Trp-cage range between -5.85+/-1.41 and -21.30+/-0.88 kcal mol(-1), leading to a significant total structural stabilization energy of -52.13+/-2.56 kcal mol(-1) of which about 50% is from the weakly polar interactions. Furthermore, the strengths of the individual weakly polar interactions are between -2.32+/-0.17 and -2.93+/-0.12 kcal mol(-1) for the CH-pi interactions, between -2.48+/-0.97 and -3.09+/-1.02 kcal mol(-1) for the Ar-NH interaction and -2.74+/-1.06 kcal mol(-1) for the Ar-Ar interaction.

Evaluation of methods to cap molecular fragments in calculating energies of interaction in avian pancreatic polypeptide

The accuracy of the determination of the energy of interaction between Phe20 and the Pro5-Thr6-Tyr7-Pro8 complex inside the hydrophobic core of avian pancreatic polypeptide was investigated using three capping strategies for molecular fractionation with conjugated caps and DFT quantum chemical calculations at the BHandHLYP/cc-pVTZ level of theory. The most accurate determination resulted from acetylation of the alpha-amino group combined with methyl amidation of the alpha-carbonyl group with relative deviations less than 10%. Combinations of hydrogenation of the alpha-amino group with the replacement of the alpha-carbonyl group with a hydrogen and the hydrogenation of the alpha-amino group with methylation of the alpha-carbonyl group were less accurate, leading to relative deviations up to 35%. Choice of capping methods depends on the structural features of the polypeptide system, the desired accuracy and the available computational resources.

The Role of Aromatic Residues in Stabilizing the Secondary and Tertiary Structure of Avian Pancreatic Polypeptide

Avian Pancreatic Polypeptide is a 36 residue protein that exhibits a tertiary fold. Results of previous experimental and computational studies indicate that the structure of aPP is stabilized more by non-bonded interactions than by the hydrophobic effect. Aromatic residues are known to participate in a variety of long range non-bonded interactions, with both backbone atoms and the atoms of other side-chains, which could be responsible, in part, for the stability of both the local secondary structure and the tertiary fold. The effect of these aromatic interactions on the stability of aPP was calculated using BHandHLYP/cc-pVTZ. Aromatic residues were shown to participate in multiple hydrogen bonded and weakly polar interactions in the secondary structure. The energies of the weakly polar interactions are comparable with those of hydrogen bonds. Aromatic residues were also shown to participate in multiple weakly polar interactions across the tertiary fold, again with energies similar to those of hydrogen bonds.

Aromatic-backbone interactions in model α-helical peptides

The effects on helical stability of weak polar interactions between aromatic side-chains and the peptide backbone were examined. alpha-Helical model peptides, hexa-Ala, with sequential Tyr replacement, were investigated computationally to obtain the geometries and energetics of the interactions. Geometries were obtained with the B3LYP/6-31G* level of theory. Interaction energies were calculated using BHandHLYP/cc-pVTZ and an improved method to correct for basis set superposition error when fragmentation caused steric clashes. Both i, i + 1 and i, i - 4 interactions were observed when Tyr was in position i = 5. The position of the aromatic residue in the amino acid sequence was crucial in facilitating aromatic-backbone interactions. The distance between the center of the aromatic ring of Tyr and the individual interacting backbone atoms ranged from 3.65 to 5.50 A. The interactions have energies of the same order as hydrogen bonds and, thus, could have a significant impact on the stability of the helix.