Search

Since the beginning of the COVID-19 pandemic, extensive drug repurposing efforts have sought to identify small-molecule antivirals with various mechanisms of action. Here, we aim to review research progress on small-molecule viral entry and fusion inhibitors that directly bind to the SARS-CoV-2 Spike protein. Early in the pandemic, numerous small molecules were identified in drug repurposing screens and reported to be effective in in vitro SARS-CoV-2 viral entry or fusion inhibitors. However, given minimal experimental information regarding the exact location of small-molecule binding sites on Spike, it was unclear what the specific mechanism of action was or where the exact binding sites were on Spike for some inhibitor candidates. The work of countless researchers has yielded great progress, with the identification of many viral entry inhibitors that target elements on the S1 receptor-binding domain (RBD) or N-terminal domain (NTD) and disrupt the S1 receptor-binding function. In this review, we will also focus on highlighting fusion inhibitors that target inhibition of the S2 fusion function, either by disrupting the formation of the postfusion S2 conformation or alternatively by stabilizing structural elements of the prefusion S2 conformation to prevent conformational changes associated with S2 function. We highlight experimentally validated binding sites on the S1/S2 interface and on the S2 subunit. While most substitutions to the Spike protein to date in variants of concern (VOCs) have been localized to the S1 subunit, the S2 subunit sequence is more conserved, with only a few observed substitutions in proximity to S2 binding sites. Several recent small molecules targeting S2 have been shown to have robust activity over recent VOC mutant strains and/or greater broad-spectrum antiviral activity for other more distantly related coronaviruses.

Clofazimine and Arbidol have both been reported to be effective in vitro SARS-CoV-2 fusion inhibitors. Both are promising drugs that have been repurposed for the treatment of COVID-19 and have been used in several previous and ongoing clinical trials. Small-molecule bindings to expressed constructs of the trimeric S2 segment of Spike and the full-length SARS-CoV-2 Spike protein were measured using a Surface Plasmon Resonance (SPR) binding assay. We demonstrate that Clofazimine, Toremifene, Arbidol and its derivatives bind to the S2 segment of the Spike protein. Clofazimine provided the most reliable and highest-quality SPR data for binding with S2 over the conditions explored. A molecular docking approach was used to identify the most favorable binding sites on the S2 segment in the prefusion conformation, highlighting two possible small-molecule binding sites for fusion inhibitors. Results related to molecular docking and modeling of the structure-activity relationship (SAR) of a newly reported series of Clofazimine derivatives support the proposed Clofazimine binding site on the S2 segment. When the proposed Clofazimine binding site is superimposed with other experimentally determined coronavirus structures in structure-sequence alignments, the changes in sequence and structure may rationalize the broad-spectrum antiviral activity of Clofazimine in closely related coronaviruses such as SARS-CoV, MERS, hCoV-229E, and hCoV-OC43.

Catecholamine-stimulated β 2 -adrenergic receptor (β 2 AR) signaling via the canonical G s -adenylyl cyclase-cAMP-PKA pathway regulates numerous physiological functions, including the therapeutic effects of exogenous β-agonists in the treatment of airway disease. β 2 AR signaling is tightly regulated by GRKs and β-arrestins, which together promote β 2 AR desensitization and internalization as well as downstream signaling, often antithetical to the canonical pathway. Thus, the ability to bias β 2 AR signaling toward the G s pathway while avoiding β-arrestin-mediated effects may provide a strategy to improve the functional consequences of β 2 AR activation. Since attempts to develop G s -biased agonists and allosteric modulators for the β 2 AR have been largely unsuccessful, here we screened small molecule libraries for allosteric modulators that selectively inhibit β-arrestin recruitment to the receptor. This screen identified several compounds that met this profile, and, of these, a difluorophenyl quinazoline (DFPQ) derivative was found to be a selective negative allosteric modulator of β-arrestin recruitment to the β 2 AR while having no effect on β 2 AR coupling to G s . DFPQ effectively inhibits agonist-promoted phosphorylation and internalization of the β 2 AR and protects against the functional desensitization of β-agonist mediated regulation in cell and tissue models. The effects of DFPQ were also specific to the β 2 AR with minimal effects on the β 1 AR. Modeling, mutagenesis, and medicinal chemistry studies support DFPQ derivatives binding to an intracellular membrane-facing region of the β 2 AR, including residues within transmembrane domains 3 and 4 and intracellular loop 2. DFPQ thus represents a class of biased allosteric modulators that targets an allosteric site of the β 2 AR.

Background and Purpose: β-Adrenoceptor agonists relieve airflow obstruction by activating β 2 -adrenoceptors, which are G protein-coupled receptors (GPCRs) expressed on human airway smooth muscle (HASM) cells. The currently available β-adrenoceptor agonists are balanced agonists, however, and signal through both the stimulatory G protein (G s )- and β-arrestin-mediated pathways. While G s signalling is beneficial and promotes HASM relaxation, β-arrestin activation is associated with reduced G s efficacy. In this context, biased ligands that selectively promote β 2 -adrenoceptor coupling to G s signalling represent a promising strategy to treat asthma. Here, we examined several β-adrenoceptor agonists to identify G s -biased ligands devoid of β-arrestin-mediated effects., Experimental Approach: G s -biased ligands for the β 2 -adrenoceptor were identified by high-throughput screening and then evaluated for G s interaction, G i interaction, cAMP production, β-arrestin interaction, GPCR kinase (GRK) phosphorylation of the receptor, receptor trafficking, ERK activation, and functional desensitization of the β 2 -adrenoceptor., Key Results: We identified ractopamine, dobutamine, and higenamine as G s -biased agonists that activate the G s /cAMP pathway upon β 2 -adrenoceptor stimulation while showing minimal G i or β-arrestin interaction. Furthermore, these compounds did not induce any receptor trafficking and had reduced GRK5-mediated phosphorylation of the β 2 -adrenoceptor. Finally, we observed minimal physiological desensitization of the β 2 -adrenoceptor in primary HASM cells upon treatment with biased agonists., Conclusion and Implications: Our work demonstrates that G s -biased signalling through the β 2 -adrenoceptor may prove to be an effective strategy to promote HASM relaxation in the treatment of asthma. Such biased compounds may also be useful in identifying the molecular mechanisms that determine biased signalling and in design of safer drugs., (© 2022 British Pharmacological Society.)

Human parathyroid hormone (PTH) is an 84-amino acid peptide that contains two methionine (Met) residues located at positions 8 and 18. It has long been recognized that Met residues in PTH are subject to oxidation to become Met sulfoxide, resulting in a decreased biological function of the peptide. However, the mechanism of the lost biological function of PTH oxidation remains elusive. To characterize whether the shift from the hydrophobic nature of the native Met residue to the hydrophilic nature of Met sulfoxide plays a role in the reduction of biological activity upon PTH oxidation, we conducted in silico and in vitro site-directed mutagenesis of Met-8 and Met-18 to the hydrophilic residue asparagine (Asn) or to the hydrophobic residue leucine (Leu) and compared the behavior of these mutated peptides with that of PTH oxidized at Met-8 and/or Met-18. Our results showed that the biological activity of the Asn-8 and Asn-8/Asn-18 mutants was significantly reduced, similar to Met-8 sulfoxide and Met-8/Met-18 sulfoxide analogues, while the functions of Asn-18, Leu-8, Leu-8/Leu-18 mutants, or Met-18 sulfoxide analogues were similar to wild-type PTH. This is rationalized from molecular modeling and immunoprecipitation assay, demonstrating disruption of hydrophobic interactions between Met-8 and Met-18 of PTH and type-1 PTH receptor (PTHR1) upon mutation or oxidation. Thus, these novel findings support the notion that the loss of biological function of PTH upon oxidation of Met-8 is due, at least in part, to the conversion from a hydrophobic to a hydrophilic residue that disrupts direct hydrophobic interaction between PTH and PTHR1.

Background: Changes in the expression and activity of the AKT oncogene play an important role in psychiatric disease. We present translational data assessing the role of AKT in psychiatric symptoms., Methods: (1) We assessed the protein activity of an AKT3 mutant harboring a PH domain mutation (Q60H) detected in a patient with schizophrenia, the corresponding AKT1 mutant (Q61H), and wild-type AKT1 and AKT3 transduced in AKT-null mouse fibroblasts and modeled the Q61H mutation onto the crystal structure of the Akt1 PH domain. (2) We analyzed the results of earlier genome-wide association studies to determine the distribution of schizophrenia-associated single-nucleotide polymorphisms (SNPs) in the AKT3 gene. (3) We analyzed the psychiatric adverse events (AEs) of patients treated with M2698 (p70S6K/AKT1/AKT3 inhibitor) and with other PI3K/AKT/mTOR pathway inhibitors., Results: (1) Proteins encoded by AKT3 (AKT3Q60H) and AKT1 (AKT1Q61H) mutants had lower kinase activity than those encoded by wild-type AKT3 and AKT1, respectively. Molecular modeling of the AKT1-Q61H mutant suggested conformational changes that may reduce the binding of D3-phosphorylated phosphoinositides to the PH domain. (2) We identified multiple SNPs in the AKT3 gene that were strongly associated with schizophrenia (p < 0.5 × 10 -8 ). (3) Psychiatric AEs, mostly insomnia, anxiety, and depression, were noted in 29% of patients treated with M2698. In randomized studies, their incidence was higher in PI3K/AKT/mTOR inhibitor arms compared with placebo arms. All psychiatric AEs were reversible., Conclusions: Our data elucidate the incidence and mechanisms of psychiatric AEs in patients treated with PI3K/AKT/mTOR inhibitors and emphasize the need for careful monitoring., (© 2022. The Author(s).)

Umifenovir (Arbidol) has been reported to exhibit some degree of efficacy in multiple clinical trials for the treatment of COVID-19 as a monotherapy. It has also demonstrated synergistic inhibition of SARS-CoV-2 with other direct-acting antivirals such as Remdesivir. A computational approach was used to identify the most favorable binding site to the SARS-CoV-2 Spike S2 segment and to perform virtual screening. Compounds selected from modeling were evaluated in a live SARS-CoV-2 infection assay. An Arbidol (ARB) derivative with substitutions at both the C-4 and C-6 positions was found to exhibit a modest improvement in activity and solubility properties in comparison to ARB. However, all of the derivatives were found to only be partial inhibitors, rather than full inhibitors in a virus-induced cytopathic effect-based assay. The binding mode is also corroborated by parallel modeling of a series of oleanolic acid trisaccharide saponin fusion inhibitors shown to bind to the S2 segment. Recently determined experimental structures of the Spike protein allowed atomic resolution modeling of fusion inhibitor binding as a function of pH, and the implications for the molecular mechanism of direct-acting fusion inhibitors targeting the S2 segment are discussed.

The 2019 emergence of, SARS-CoV-2 has tragically taken an immense toll on human life and far reaching impacts on society. There is a need to identify effective antivirals with diverse mechanisms of action in order to accelerate preclinical development. This study focused on five of the most established drug target proteins for direct acting small molecule antivirals: Nsp5 Main Protease, Nsp12 RNA-dependent RNA polymerase, Nsp13 Helicase, Nsp16 2'-O methyltransferase and the S2 subunit of the Spike protein. A workflow of solvent mapping and free energy calculations was used to identify and characterize favorable small-molecule binding sites for an aromatic pharmacophore (benzene). After identifying the most favorable sites, calculated ligand efficiencies were compared utilizing computational fragment screening. The most favorable sites overall were located on Nsp12 and Nsp16, whereas the most favorable sites for Nsp13 and S2 Spike had comparatively lower ligand efficiencies relative to Nsp12 and Nsp16. Utilizing fragment screening on numerous possible sites on Nsp13 helicase, we identified a favorable allosteric site on the N-terminal zinc binding domain (ZBD) that may be amenable to virtual or biophysical fragment screening efforts. Recent structural studies of the Nsp12:Nsp13 replication-transcription complex experimentally corroborates ligand binding at this site, which is revealed to be a functional Nsp8:Nsp13 protein-protein interaction site in the complex. Detailed structural analysis of Nsp13 ZBD conformations show the role of induced-fit flexibility in this ligand binding site and identify which conformational states are associated with efficient ligand binding. We hope that this map of over 200 possible small-molecule binding sites for these drug targets may be of use for ongoing discovery, design, and drug repurposing efforts. This information may be used to prioritize screening efforts or aid in the process of deciphering how a screening hit may bind to a specific target protein., Competing Interests: The authors have declared that no competing interests exist.

Phosphorylation of specific residues in the activation loops of AGC kinase group (protein kinase A, G, and C families) is required for activity of most of these kinases, including the catalytic subunit of PKA (PKAc). Although many phosphorylated AGC kinases are sensitive to phosphatase-mediated dephosphorylation, the PKAc activation loop uniquely resists dephosphorylation, rendering it "constitutively" phosphorylated in cells. Previous biophysical experiments and structural modeling have suggested that the N-terminal myristoylation signal and the C-terminal F XX F motif in PKAc regulate its thermal stability and catalysis. Here, using site-directed mutagenesis, molecular modeling, and in cell-free and cell-based systems, we demonstrate that substitutions of either the PKAc myristoylation signal or the F XX F motif only modestly reduce phosphorylation and fail to affect PKAc function in cells. However, we observed that these two sites cooperate with an N-terminal F XX W motif to cooperatively establish phosphatase resistance of PKAc while not affecting kinase-dependent phosphorylation of the activation loop. We noted that this tripartite cooperative mechanism of phosphatase resistance is functionally relevant, as demonstrated by changes in morphology, adhesion, and migration of human airway smooth muscle cells transfected with PKAc variants containing amino acid substitutions in these three sites. These findings establish that three allosteric sites located at the PKAc N and C termini coordinately regulate the phosphatase sensitivity of this enzyme. This cooperative mechanism of phosphatase resistance of AGC kinase opens new perspectives toward therapeutic manipulation of kinase signaling in disease., (© 2020 Chan et al.)

The diarylisoxazole molecular scaffold is found in several NSAIDs, especially those with high selectivity for COX-1. Here, we have determined the structural basis for COX-1 binding to two diarylisoxazoles: mofezolac, which is polar and ionizable, and 3-(5-chlorofuran-2-yl)-5-methyl-4-phenylisoxazole (P6) that has very low polarity. X-ray analysis of the crystal structures of COX-1 bound to mofezolac and 3-(5-chlorofuran-2-yl)-5-methyl-4-phenylisoxazole allowed the identification of specific binding determinants within the enzyme active site, relevant to generate structure/activity relationships for diarylisoxazole NSAIDs., (Copyright © 2017 Elsevier Masson SAS. All rights reserved.)

Following insights from recent crystal structures of the muscarinic acetylcholine receptor, binding modes of Positive Allosteric Modulators (PAMs) were predicted under the assumption that PAMs should bind to the extracellular surface of the active state. A series of well-characterized PAMs for adenosine (A 1 R, A 2A R, A 3 R) and muscarinic acetylcholine (M 1 R, M 5 R) receptors were modeled using both rigid and flexible receptor CHARMM-based molecular docking. Studies of adenosine receptors investigated the molecular basis of the probe-dependence of PAM activity by modeling in complex with specific agonist radioligands. Consensus binding modes map common pharmacophore features of several chemical series to specific binding interactions. These models provide a rationalization of how PAM binding slows agonist radioligand dissociation kinetics. M 1 R PAMs were predicted to bind in the analogous M 2 R PAM LY2119620 binding site. The M 5 R NAM (ML-375) was predicted to bind in the PAM (ML-380) binding site with a unique induced-fit receptor conformation. © 2017 Wiley Periodicals, Inc., (© 2017 Wiley Periodicals, Inc.)

The Akt protein kinase, also known as protein kinase B, plays key roles in insulin receptor signalling and regulates cell growth, survival and metabolism. Recently, we described a mechanism to enhance Akt phosphorylation that restricts access of cellular phosphatases to the Akt activation loop (Thr(308) in Akt1 or protein kinase B isoform alpha) in an ATP-dependent manner. In the present paper, we describe a distinct mechanism to control Thr(308) dephosphorylation and thus Akt deactivation that depends on intramolecular interactions of Akt C-terminal sequences with its kinase domain. Modifications of amino acids surrounding the Akt1 C-terminal mTORC2 (mammalian target of rapamycin complex 2) phosphorylation site (Ser(473)) increased phosphatase resistance of the phosphorylated activation loop (pThr(308)) and amplified Akt phosphorylation. Furthermore, the phosphatase-resistant Akt was refractory to ceramide-dependent dephosphorylation and amplified insulin-dependent Thr(308) phosphorylation in a regulated fashion. Collectively, these results suggest that the Akt C-terminal hydrophobic groove is a target for the development of agents that enhance Akt phosphorylation by insulin., (© 2015 Authors; published by Portland Press Limited.)

Correctly predicting off-targets for a given molecular structure, which would have the ability to bind a large range of ligands, is both particularly difficult and important if they share no significant sequence or fold similarity with the respective molecular target ("distant off-targets"). A novel approach for identification of off-targets by direct superposition of protein binding pocket surfaces is presented and applied to a set of well-studied and highly relevant drug targets, including representative kinases and nuclear hormone receptors. The entire Protein Data Bank is searched for similar binding pockets and convincing distant off-target candidates were identified that share no significant sequence or fold similarity with the respective target structure. These putative target off-target pairs are further supported by the existence of compounds that bind strongly to both with high topological similarity, and in some cases, literature examples of individual compounds that bind to both. Also, our results clearly show that it is possible for binding pockets to exhibit a striking surface similarity, while the respective off-target shares neither significant sequence nor significant fold similarity with the respective molecular target ("distant off-target").

Poly(ADP-ribose) polymerases (PARPs) are a family of enzymes that use NAD(+) as a substrate to synthesize polymers of ADP-ribose (PAR) as post-translational modifications of proteins. PARPs have important cellular roles that include preserving genomic integrity, telomere maintenance, transcriptional regulation, and cell fate determination. The diverse biological roles of PARPs have made them attractive therapeutic targets, which have fueled the pursuit of small molecule PARP inhibitors. The design of PARP inhibitors has matured over the past several years resulting in several lead candidates in clinical trials. PARP inhibitors are mainly used in clinical trials to treat cancer, particularly as sensitizing agents in combination with traditional chemotherapy to reduce side effects. An exciting aspect of PARP inhibitors is that they are also used to selectivity kill tumors with deficiencies in DNA repair proteins (e.g., BRCA1/2) through an approach termed "synthetic lethality." In the midst of the tremendous efforts that have brought PARP inhibitors to the forefront of modern chemotherapy, most clinically used PARP inhibitors bind to conserved regions that permits cross-selectivity with other PARPs containing homologous catalytic domains. Thus, the differences between therapeutic effects and adverse effects stemming from pan-PARP inhibition compared to selective inhibition are not well understood. In this review, we discuss current literature that has found ways to gain selectivity for one PARP over another. We furthermore provide insights into targeting other domains that make up PARPs, and how new classes of drugs that target these domains could provide a high degree of selectivity by affecting specific cellular functions. A clear understanding of the inhibition profiles of PARP inhibitors will not only enhance our understanding of the biology of individual PARPs, but may provide improved therapeutic options for patients.

Molecular docking of small-molecules is an important procedure for computer-aided drug design. Modeling receptor side chain flexibility is often important or even crucial, as it allows the receptor to adopt new conformations as induced by ligand binding. However, the accurate and efficient incorporation of receptor side chain flexibility has proven to be a challenge due to the huge computational complexity required to adequately address this problem. Here we describe a new docking approach with a very fast, graph-based optimization algorithm for assignment of the near-optimal set of residue rotamers. We extensively validate our approach using the 40 DUD target benchmarks commonly used to assess virtual screening performance and demonstrate a large improvement using the developed side chain optimization over rigid receptor docking (average ROC AUC of 0.693 vs. 0.623). Compared to numerous benchmarks, the overall performance is better than nearly all other commonly used procedures. Furthermore, we provide a detailed analysis of the level of receptor flexibility observed in docking results for different classes of residues and elucidate potential avenues for further improvement., (Copyright © 2013 Wiley Periodicals, Inc.)

AGC kinases, including the three Akt (protein kinase B) isoforms, protein kinase A (PKA) and all protein kinase C (PKC) isoforms, require activation loop phosphorylation (threonine 308 in Akt1) as well as phosphorylation of a C-terminal residue (serine 473 in Akt1) for catalytic activity and phosphorylation of downstream targets. Conversely, phosphatases reverse these phosphorylations. Virtually all cellular processes are affected by AGC kinases, a circumstance that has led to intense scrutiny of the molecular mechanisms that regulate phosphorylation of these kinases. Here, we review a new layer of control of phosphorylation in Akt, PKA and PKC pointing to ATP binding pocket occupancy as a means to decelerate dephosphorylation of these and, potentially, other kinases. This additional level of kinase regulation opens the door to search for new functional motifs for the rational design of non- ATP-competitive kinase inhibitors that discriminate within and between protein kinase families.

Small- and intermediate-conductance Ca(2+)-activated potassium channels, activated by Ca(2+)-bound calmodulin, have an important role in regulating membrane excitability. These channels are also linked to clinical abnormalities. A tremendous amount of effort has been devoted to developing small molecule compounds targeting these channels. However, these compounds often suffer from low potency and lack of selectivity, hindering their potential for clinical use. A key contributing factor is the lack of knowledge of the binding site(s) for these compounds. Here we demonstrate by X-ray crystallography that the binding pocket for the compounds of the 1-ethyl-2-benzimidazolinone (1-EBIO) class is located at the calmodulin-channel interface. We show that, based on structure data and molecular docking, mutations of the channel can effectively change the potency of these compounds. Our results provide insight into the molecular nature of the binding pocket and its contribution to the potency and selectivity of the compounds of the 1-EBIO class.

Phosphorylation of a threonine residue (T308 in Akt1) in the activation loop of Akt kinases is a prerequisite for deregulated Akt activity frequently observed in neoplasia. Akt phosphorylation in vivo is balanced by the opposite activities of kinases and phosphatases. Here we describe that targeting Akt kinase to the cell membrane markedly reduced sensitivity of phosphorylated Akt to dephosphorylation by protein phosphatase 2A. This effect was amplified by occupancy of the ATP binding pocket by either ATP or ATP-competitive inhibitors. Mutational analysis revealed that R273 in Akt1 and the corresponding R274 in Akt2 are essential for shielding T308 in the activation loop against dephosphorylation. Thus, occupancy of the nucleotide binding pocket of Akt kinases enables intramolecular interactions that restrict phosphatase access and sustain Akt phosphorylation. This mechanism provides an explanation for the "paradoxical" Akt hyperphosphorylation induced by ATP-competitive inhibitor, A-443654. The lack of phosphatase resistance further contributes insight into the mechanism by which the human Akt2 R274H missense mutation may cause autosomal-dominant diabetes mellitus.

The performances of several two-step scoring approaches for molecular docking were assessed for their ability to predict binding geometries and free energies. Two new scoring functions designed for "step 2 discrimination" were proposed and compared to our CHARMM implementation of the linear interaction energy (LIE) approach using the Generalized-Born with Molecular Volume (GBMV) implicit solvation model. A scoring function S1 was proposed by considering only "interacting" ligand atoms as the "effective size" of the ligand and extended to an empirical regression-based pair potential S2. The S1 and S2 scoring schemes were trained and 5-fold cross-validated on a diverse set of 259 protein-ligand complexes from the Ligand Protein Database (LPDB). The regression-based parameters for S1 and S2 also demonstrated reasonable transferability in the CSARdock 2010 benchmark using a new data set (NRC HiQ) of diverse protein-ligand complexes. The ability of the scoring functions to accurately predict ligand geometry was evaluated by calculating the discriminative power (DP) of the scoring functions to identify native poses. The parameters for the LIE scoring function with the optimal discriminative power (DP) for geometry (step 1 discrimination) were found to be very similar to the best-fit parameters for binding free energy over a large number of protein-ligand complexes (step 2 discrimination). Reasonable performance of the scoring functions in enrichment of active compounds in four different protein target classes established that the parameters for S1 and S2 provided reasonable accuracy and transferability. Additional analysis was performed to definitively separate scoring function performance from molecular weight effects. This analysis included the prediction of ligand binding efficiencies for a subset of the CSARdock NRC HiQ data set where the number of ligand heavy atoms ranged from 17 to 35. This range of ligand heavy atoms is where improved accuracy of predicted ligand efficiencies is most relevant to real-world drug design efforts.

The arenavirus genome encodes for a Z-protein, which contains a RING domain that coordinates two zinc ions, and has been identified as having several functional roles at various stages of the virus life cycle. Z-protein binds to multiple host proteins and has been directly implicated in the promotion of viral budding, repression of mRNA translation, and apoptosis of infected cells. Using homology models of the Z-protein from Lassa strain arenavirus, replica exchange molecular dynamics (MD) was used to refine the structures, which were then subsequently clustered. Population-weighted ensembles of low-energy cluster representatives were predicted based upon optimal agreement of the chemical shifts computed with the SPARTA program with the experimental NMR chemical shifts. A member of the refined ensemble was identified to be a potential binder of budding factor Tsg101 based on its correspondence to the structure of the HIV-1 Gag late domain when bound to Tsg101. Members of these ensembles were docked against the crystal structure of human eIF4E translation initiation factor. Two plausible binding modes emerged based upon their agreement with experimental observation, favorable interaction energies and stability during MD trajectories. Mutations to Z are proposed that would either inhibit both binding mechanisms or selectively inhibit only one mode. The C-terminal domain conformation of the most populated member of the representative ensemble shielded protein-binding recognition motifs for Tsg101 and eIF4E and represents the most populated state free in solution. We propose that C-terminal flexibility is key for mediating the different functional states of the Z-protein., ((c) 2010 Wiley-Liss, Inc.)

Orthogonal aminoacyl-tRNA synthetase/tRNA pairs from archaea have been evolved to facilitate site specific in vivo incorporation of unnatural amino acids into proteins in Escherichia coli. Using this approach, unnatural amino acids have been successfully incorporated with high translational efficiency and fidelity. In this study, CHARMM-based molecular docking and free energy calculations were used to evaluate rational design of specific protein-ligand interactions for aminoacyl-tRNA synthetases. A series of novel unnatural amino acid ligands were docked into the p-benzoyl-L-phenylalanine tRNA synthetase, which revealed that the binding pocket of the enzyme does not provide sufficient space for significantly larger ligands. Specific binding site residues were mutated to alanine to create additional space to accommodate larger target ligands, and then mutations were introduced to improve binding free energy. This approach was used to redesign binding sites for several different target ligands, which were then tested against the standard 20 amino acids to verify target specificity. Only the synthetase designed to bind Man-alpha-O-Tyr was predicted to be sufficiently selective for the target ligand and also thermodynamically stable. Our study suggests that extensive redesign of the tRNA synthatase binding pocket for large bulky ligands may be quite thermodynamically unfavorable.

Incorporating receptor flexibility into molecular docking should improve results for flexible proteins. However, the incorporation of explicit all-atom flexibility with molecular dynamics for the entire protein chain may also introduce significant error and "noise" that could decrease docking accuracy and deteriorate the ability of a scoring function to rank native-like poses. We address this apparent paradox by comparing the success of several flexible receptor models in cross-docking and multiple receptor ensemble docking for p38α mitogen-activated protein (MAP) kinase. Explicit all-atom receptor flexibility has been incorporated into a CHARMM-based molecular docking method (CDOCKER) using both molecular dynamics (MD) and torsion angle molecular dynamics (TAMD) for the refinement of predicted protein-ligand binding geometries. These flexible receptor models have been evaluated, and the accuracy and efficiency of TAMD sampling is directly compared to MD sampling. Several flexible receptor models are compared, encompassing flexible side chains, flexible loops, multiple flexible backbone segments, and treatment of the entire chain as flexible. We find that although including side chain and some backbone flexibility is required for improved docking accuracy as expected, docking accuracy also diminishes as additional and unnecessary receptor flexibility is included into the conformational search space. Ensemble docking results demonstrate that including protein flexibility leads to to improved agreement with binding data for 227 active compounds. This comparison also demonstrates that a flexible receptor model enriches high affinity compound identification without significantly increasing the number of false positives from low affinity compounds.

Transthyretin (TTR)-containing amyloid fibrils are deposited in cardiac tissue as a natural consequence of aging. A large number of inherited mutations lead to amyloid diseases by accelerating TTR deposition in other organs. Amyloid formation is preceded by a disruption of the quaternary structure of TTR and conformational changes in the monomer. To study conformational changes preceding the formation of amyloid, we performed molecular dynamics simulations of the wild-type monomer, amyloidogenic variants (V30M, L55P, V122I) and a protective variant (T119M) at neutral and low pH. At low pH, the D strand dissociated from the beta-sheet to expose the A strand, consistent with experimental studies. In amyloidogenic variants and in the wild-type at low pH, there was a conformational change in the beta-sheets into alpha-sheet via peptide bond flips that was not observed at neutral pH in the wild-type monomer. The same residues participated in conversion in each amyloidogenic variant simulation, originating in the G strand between residues 106 and 109, with accelerated conversion at low pH. The T119M protective variant changed the local conformation of the H strand and suppressed the conversion observed in amyloidogenic variants.

The self-assembly of beta(2)-microglobulin into fibrils leads to dialysis-related amyloidosis. pH-mediated partial unfolding is required for the formation of the amyloidogenic intermediate that then self-assembles into amyloid fibrils. Two partially folded intermediates of beta(2)-microglobulin have been identified experimentally and linked to the formation of fibrils of distinct morphology, yet it remains difficult to characterize these partially unfolded states at high resolution using experimental approaches. Consequently, we have performed molecular dynamics simulations at neutral and low pH to determine the structures of these partially unfolded amyloidogenic intermediates. In the low-pH simulations, we observed the formation of alpha-sheet structure, which was first proposed by Pauling and Corey. Multiple simulations were performed, and two distinct intermediate state ensembles were identified that may account for the different fibril morphologies. The predominant early unfolding intermediate was nativelike in structure, in agreement with previous NMR studies. The late unfolding intermediate was significantly disordered, but it maintained an extended elongated structure, with hydrophobic clusters and residual alpha-extended chain strands in specific regions of the sequence that map to amyloidogenic peptides. We propose that the formation of alpha-sheet facilitates self-assembly into partially unfolded prefibrillar amyloidogenic intermediates.

Several neurodegenerative diseases are linked to expanded repeats of glutamine residues, which lead to the formation of amyloid fibrils and neuronal death. The length of the repeats correlates with the onset of Huntington's disease, such that healthy individuals have <38 residues and individuals with >38 repeats exhibit symptoms. Because it is difficult to obtain atomic-resolution structural information for poly(l-glutamine) (polyQ) in aqueous solution experimentally, we performed molecular dynamics simulations to investigate the conformational behavior of this homopolymer. In simulations of 20-, 40-, and 80-mer polyQ, we observed the formation of the "alpha-extended chain" conformation, which is characterized by alternating residues in the alpha(L) and alpha(R) conformations to yield a sheet. The structural transition from disordered random-coil conformations to the alpha-extended chain conformation exhibits modest length and temperature dependence, in agreement with the experimental observation that aggregation depends on length and temperature. We propose that fibril formation in polyQ may occur through an alpha-sheet structure, which was proposed by Pauling and Corey. Also, we propose an atomic-resolution model of how the inhibitory peptide QBP1 (polyQ-binding peptide 1) may bind to polyQ in an alpha-extended chain conformation to inhibit fibril formation.

The correct treatment of van der Waals and electrostatic nonbonded interactions in molecular force fields is essential for performing realistic molecular dynamics (MD) simulations of solvated polypeptides. The most computationally tractable treatment of nonbonded interactions in MD utilizes a spherical distance cutoff (typically, 8-12 A) to reduce the number of pairwise interactions. In this work, we assess three spherical atom-based cutoff approaches for use with all-atom explicit solvent MD: abrupt truncation, a CHARMM-style electrostatic shift truncation, and our own force-shifted truncation. The chosen system for this study is an end-capped 17-residue alanine-based alpha-helical peptide, selected because of its use in previous computational and experimental studies. We compare the time-averaged helical content calculated from these MD trajectories with experiment. We also examine the effect of varying the cutoff treatment and distance on energy conservation. We find that the abrupt truncation approach is pathological in its inability to conserve energy. The CHARMM-style shift truncation performs quite well but suffers from energetic instability. On the other hand, the force-shifted spherical cutoff method conserves energy, correctly predicts the experimental helical content, and shows convergence in simulation statistics as the cutoff is increased. This work demonstrates that by using proper and rigorous techniques, it is possible to correctly model polypeptide dynamics in solution with a spherical cutoff. The inherent computational advantage of spherical cutoffs over Ewald summation (and related) techniques is essential in accessing longer MD time scales.

The homotetramer of transthyretin (TTR) dissociates into a monomeric amyloidogenic intermediate that self-assembles into amyloid fibrils at low pH. We have performed molecular dynamics simulations of monomeric TTR at neutral and low pH at physiological (310 K) and very elevated temperature (498 K). In the low-pH simulations at both temperatures, one of the two beta-sheets (strands CBEF) becomes disrupted, and alpha-sheet structure forms in the other sheet (strands DAGH). alpha-sheet is formed by alternating alphaL and alphaR residues, and it was first proposed by Pauling and Corey. Overall, the simulations are in agreement with the available experimental observations, including solid-state NMR results for a TTR-peptide amyloid. In addition, they provide a unique explanation for the results of hydrogen exchange experiments of the amyloidogenic intermediate-results that are difficult to explain with beta-structure. We propose that alpha-sheet may represent a key pathological conformation during amyloidogenesis., (Copyright 2004 Elsevier Ltd.)

Transthyretin, beta(2)-microglobulin, lysozyme, and the prion protein are four of the best-characterized proteins implicated in amyloid disease. Upon partial acid denaturation, these proteins undergo conformational change into an amyloidogenic intermediate that can self-assemble into amyloid fibrils. Many experiments have shown that pH-mediated changes in structure are required for the formation of the amyloidogeneic intermediate, but it has proved impossible to characterize these conformational changes at high resolution using experimental means. To probe these conformational changes at atomic resolution, we have performed molecular dynamics simulations of these proteins at neutral and low pH. In low-pH simulations of all four proteins, we observe the formation of alpha-pleated sheet secondary structure, which was first proposed by L. Pauling and R. B. Corey [(1951) Proc. Natl. Acad. Sci. USA 37, 251-256]. In all beta-sheet proteins, transthyretin and beta(2)-microglobulin, alpha-pleated sheet structure formed over the strands that are highly protected in hydrogen-exchange experiments probing amyloidogenic conditions. In lysozyme and the prion protein, alpha-sheets formed in the specific regions of the protein implicated in the amyloidogenic conversion. We propose that the formation of alpha-pleated sheet structure may be a common conformational transition in amyloidosis.

We have determined consensus protein-fold classifications on the basis of three classification methods, SCOP, CATH, and Dali. These classifications make use of different methods of defining and categorizing protein folds that lead to different views of protein-fold space. Pairwise comparisons of domains on the basis of their fold classifications show that much of the disagreement between the classification systems is due to differing domain definitions rather than assigning the same domain to different folds. However, there are significant differences in the fold assignments between the three systems. These remaining differences can be explained primarily in terms of the breadth of the fold classifications. Many structures may be defined as having one fold in one system, whereas far fewer are defined as having the analogous fold in another system. By comparing these folds for a nonredundant set of proteins, the consensus method breaks up broad fold classifications and combines restrictive fold classifications into metafolds, creating, in effect, an averaged view of fold space. This averaged view requires that the structural similarities between proteins having the same metafold be recognized by multiple classification systems. Thus, the consensus map is useful for researchers looking for fold similarities that are relatively independent of the method used to compare proteins. The 30 most populated metafolds, representing the folds of about half of a nonredundant subset of the PDB, are presented here. The full list of metafolds is presented on the Web.

Neutrophil gelatinase associated lipocalin (NGAL), a constituent of neutrophil granules, is a member of the lipocalin family of binding proteins. NGAL can also be highly induced in epithelial cells in both inflammatory and neoplastic colorectal disease. NGAL is proposed to mediate inflammatory responses by sequestering neutrophil chemoattractants, particularly N-formylated tripeptides and possibly leukotriene B(4) and platelet activating factor. The crystal structures of NGAL display a typical lipocalin fold, albeit with an unusually large and atypically polar binding site, or calyx. The fold of NGAL is most similar to the epididymal retinoic acid-binding protein, another lipocalin, though the overall architecture of the calyces are very different. The crystal structures also reveal either sulfate ions or an adventitiously copurified fatty acid bound in the binding site. Neither ligand is displaced by added N-formylated tripeptides. The size, shape, and character of the NGAL calyx, as well as the low relative affinity for N-formylated tripeptides, suggest that neither the copurified fatty acid nor any of the proposed ligands are likely to be the preferred ligand of this protein. Comparisons between the crystal structures and the recently reported solution structure of NGAL reveal significant differences, in terms of both the details of the structure and the overall flexibility of the fold.

We present a new method for the determination of bilayer structure based on a combination of computational studies and laboratory experiments. From molecular dynamics simulations, the volumes of submolecular fragments of saturated and unsaturated phosphatidylcholines in the liquid crystalline state have been extracted with a precision not available experimentally. Constancy of component volumes, both among different lipids and as a function of membrane position for a given lipid, have been examined. The component volumes were then incorporated into the liquid crystallographic method described by Wiener and White (1992. Biophys. J. 61:434-447, and references therein) for determining the structure of a fluid-phase dioleoylphosphatidylcholine bilayer from x-ray and neutron diffraction experiments.