Search

The human cytochrome P450 (CYP) 2C9 and 2C19 enzymes are two highly similar isoforms with key roles in drug metabolism. They are anchored to the endoplasmic reticulum membrane by their N-terminal transmembrane helix and interactions of their cytoplasmic globular domain with the membrane. However, their crystal structures were determined after N-terminal truncation and mutating residues in the globular domain that contact the membrane. Therefore, the CYP-membrane interactions are not structurally well-characterized and their dynamics and the influence of membrane interactions on CYP function are not well understood. We describe herein the modeling and simulation of CYP 2C9 and CYP 2C19 in a phospholipid bilayer. The simulations revealed that, despite high sequence conservation, the small sequence and structural differences between the two isoforms altered the interactions and orientations of the CYPs in the membrane bilayer. We identified residues (including K72, P73, and I99 in CYP 2C9 and E72, R73, and H99 in CYP 2C19) at the protein-membrane interface that contribute not only to the differing orientations adopted by the two isoforms in the membrane, but also to their differing substrate specificities by affecting the substrate access tunnels. Our findings provide a mechanistic interpretation of experimentally observed effects of mutagenesis on substrate selectivity.

ERRATUM Rudi Tong, Rebecca C. Wade, Neil J. Bruce "Comparative electrostatic analysis of adenylyl cyclase for isoform dependent regulation properties" Proteins: Structure, Function and Bioinformatics 2019 87(7), pp.619-620https://doi.org/10.1002/prot.25167 ORIGINAL PAPER: Rudi Tong, Rebecca C. Wade, Neil J. Bruce "Comparative electrostatic analysis of adenylyl cyclase for isoform dependent regulation properties" Proteins: Structure, Function and Bioinformatics 201684(12), pp.1844-1858https://doi.org/10.1002/prot.25167

Long-term potentiation and depression of synaptic activity in response to stimuli is a key factor in reinforcement learning. Strengthening of the corticostriatal synapses depends on the second messenger cAMP, whose synthesis is catalysed by the enzyme adenylyl cyclase 5 (AC5), which is itself regulated by the stimulatory G alpha(olf) and inhibitory G alpha(i) proteins. AC isoforms have been suggested to act as coincidence detectors, promoting cellular responses only when convergent regulatory signals occur close in time. However, the mechanism for this is currently unclear, and seems to lie in their diverse regulation patterns. Despite attempts to isolate the ternary complex, it is not known if G alpha(olf) and G alpha(i) can bind to AC5 simultaneously, nor what activity the complex would have. Using protein structure-based molecular dynamics simulations, we show that this complex is stable and inactive. These simulations, along with Brownian dynamics simulations to estimate protein association rates constants, constrain a kinetic model that shows that the presence of this ternary inactive complex is crucial for AC5's ability to detect coincident signals, producing a synergistic increase in cAMP. These results reveal some of the prerequisites for corticostriatal synaptic plasticity, and explain recent experimental data on cAMP concentrations following receptor activation. Moreover, they provide insights into the regulatory mechanisms that control signal processing by different AC isoforms., Author summary Adenylyl cyclases (ACs) are enzymes that can translate extracellular signals into the intracellular molecule cAMP, which is thus a 2nd messenger of extracellular events. The brain expresses nine membrane-bound AC variants, and AC5 is the dominant form in the striatum. The striatum is the input stage of the basal ganglia, a brain structure involved in reward learning, i.e. the learning of behaviors that lead to rewarding stimuli (such as food, water, sugar, etc). During reward learning, cAMP production is crucial for strengthening the synapses from cortical neurons onto the striatal principal neurons, and its formation is dependent on several neuromodulatory systems such as dopamine and acetylcholine. It is, however, not understood how AC5 is activated by transient (subsecond) changes in the neuromodulatory signals. Here we combine several computational tools, from molecular dynamics and Brownian dynamics simulations to bioinformatics approaches, to inform and constrain a kinetic model of the AC5-dependent signaling system. We use this model to show how the specific molecular properties of AC5 can detect particular combinations of co-occuring transient changes in the neuromodulatory signals which thus result in a supralinear/synergistic cAMP production. Our results also provide insights into the computational capabilities of the different AC isoforms.

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromo-lecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

The balance and the overall concentration of intracellular deoxyribonucleoside triphosphates (dNTPs) are important determinants of faithful DNA replication. Despite the established fact that changes in dNTP pools negatively influence DNA replication fidelity, it is not clear why certain dNTP pool alterations are more mutagenic than others. As intracellular dNTP pools are mainly controlled by ribonucleotide reductase (RNR), and given the limited number of eukaryotic RNR mutations characterized so far, we screened for RNR1 mutations causing mutator phenotypes in Saccharomyces cerevisiae. We identified 24 rnr1 mutant alleles resulting in diverse mutator phenotypes linked in most cases to imbalanced dNTPs. Among the identified rnr1 alleles the strongest mutators presented a dNTP imbalance in which three out of the four dNTPs were elevated (dCTP, dTTP and dGTP), particularly if dGTP levels were highly increased. These rnr1 alleles caused growth defects/lethality in DNA replication fidelity-compromised backgrounds, and caused strong mutator phenotypes even in the presence of functional DNA polymerases and mismatch repair. In summary, this study pinpoints key residues that contribute to allosteric regulation of RNR’s overall activity or substrate specificity. We propose a model that distinguishes between different dNTP pool alterations and provides a mechanistic explanation why certain dNTP imbalances are particularly detrimental.

The past few years have seen increasing recognition of the importance of understanding molecular binding kinetics. This has led to the development of myriad computational methods for studying the kinetics of binding processes and predicting their associated rate constants that show varying ranges of application, degrees of accuracy, and computational requirements. In order to help researchers decide which method might be suitable for their projects, we have developed KBbox, a web server that guides users in choosing the methods they should consider on the basis of the information they wish to obtain, the data they currently have available, and the computational resources to which they have access. KBbox provides information on the toolbox of available methods, their associated software tools, an expanding list of curated examples of published applications, and tutorials explaining how to apply some of the methods. It has been designed to allow the easy addition of new methods, tools, and examples as they are developed and published. KBbox is available at https://kbbox.h-its.org/ .

We describe a transferable multiresolution computational approach to build and simulate complexes of two proteins—cytochrome P450 (CYP) and CYP reductase (CPR)—in a membrane bilayer using Brownian dynamics (BD) and all-atom molecular dynamics (MD) simulations. Our benchmarks showed that MD simulations of these systems could be carried out efficiently with up to 180 nodes (4320 cores) using NAMD version 2.12. Our results provide a basis for defining the ensemble of electron transfer-competent arrangements of CYP-CPR-membrane complexes and for understanding differences in the interactions with CPR of different CYPs, which have implications for CYP-mediated drug metabolism and the exploitation of CYPs as drug targets. This work was carried out in the DYNATHOR (DYNAmics of THe complex of cytOchrome P450 and cytochrome P450 Reductase in a phospholipid bilayer) project at HLRS.

Flavonoids represent a potential source of new antitrypanosomatidic leads. Starting from a library of natural products, we combined target-based screening on pteridine reductase 1 with phenotypic screening on Trypanosoma brucei for hit identification. Flavonols were identified as hits, and a library of 16 derivatives was synthesized. Twelve compounds showed EC50 values against T. brucei below 10 μM. Four X-ray crystal structures and docking studies explained the observed structure-activity relationships. Compound 2 (3,6-dihydroxy-2-(3-hydroxyphenyl)-4H-chromen-4-one) was selected for pharmacokinetic studies. Encapsulation of compound 2 in PLGA nanoparticles or cyclodextrins resulted in lower in vitro toxicity when compared to the free compound. Combination studies with methotrexate revealed that compound 13 (3-hydroxy-6-methoxy-2-(4-methoxyphenyl)-4H-chromen-4-one) has the highest synergistic effect at concentration of 1.3 μM, 11.7-fold dose reduction index and no toxicity toward host cells. Our results provide the basis for further chemical modifications aimed at identifying novel antitrypanosomatidic agents showing higher potency toward PTR1 and increased metabolic stability.

Signal recognition particle (SRP) is a universally conserved protein-RNA complex that mediates co-translational protein translocation and membrane insertion by targeting translating ribosomes to membrane translocons. The existence of parallel co- and post-translational transport pathways, however, raises the question of the cellular substrate pool of SRP and the molecular basis of substrate selection. Here we determine the binding sites of bacterial SRP within the nascent proteome of Escherichia coli at amino acid resolution, by sequencing messenger RNA footprints of ribosome-nascent-chain complexes associated with SRP. SRP, on the basis of its strong preference for hydrophobic transmembrane domains (TMDs), constitutes a compartment-specific targeting factor for nascent inner membrane proteins (IMPs) that efficiently excludes signal-sequence-containing precursors of periplasmic and outer membrane proteins. SRP associates with hydrophobic TMDs enriched in consecutive stretches of hydrophobic and bulky aromatic amino acids immediately on their emergence from the ribosomal exit tunnel. By contrast with current models, N-terminal TMDs are frequently skipped and TMDs internal to the polypeptide sequence are selectively recognized. Furthermore, SRP binds several TMDs in many multi-spanning membrane proteins, suggesting cycles of SRP-mediated membrane targeting. SRP-mediated targeting is not accompanied by a transient slowdown of translation and is not influenced by the ribosome-associated chaperone trigger factor (TF), which has a distinct substrate pool and acts at different stages during translation. Overall, our proteome-wide data set of SRP-binding sites reveals the underlying principles of pathway decisions for nascent chains in bacteria, with SRP acting as the dominant triaging factor, sufficient to separate IMPs from substrates of the SecA-SecB post-translational translocation and TF-assisted cytosolic protein folding pathways.

Sequence-specific cleavage of collagen by mammalian collagenase plays a pivotal role in cell function. Collagenases are matrix metalloproteinases that cleave the peptide bond at a specific position on fibrillar collagen. The collagenase Hemopexin-like (HPX) domain has been proposed to be responsible for substrate recognition, but the mechanism by which collagenases identify the cleavage site on fibrillar collagen is not clearly understood. In this study, Brownian dynamics simulations coupled with atomic-detail and coarse-grained molecular dynamics simulations were performed to dock matrix metalloproteinase-1 (MMP-1) on a collagen IIIα1 triple helical peptide. We find that the HPX domain recognizes the collagen triple helix at a conserved R-X11-R motif C-terminal to the cleavage site to which the HPX domain of collagen is guided electrostatically. The binding of the HPX domain between the two arginine residues is energetically stabilized by hydrophobic contacts with collagen. From the simulations and analysis of the sequences and structural flexibility of collagen and collagenase, a mechanistic scheme by which MMP-1 can recognize and bind collagen for proteolysis is proposed.

The addition of three N-terminal histidines to β-lactamase inhibitor protein was shown experimentally to increase its binding potency to an Au(111) surface substantially but the binding mechanism was not resolved. Here, we propose a complete adsorption mechanism for this fusion protein by means of a multi-scale simulation approach and free energy calculations. We find that adsorption is a three-step process: (i) recognition of the surface predominantly by the histidine fusion peptide and formation of an encounter complex facilitated by a reduced dielectric screening of water in the interfacial region, (ii) adsorption of the protein on the surface and adoption of a specific binding orientation, and (iii) adaptation of the protein structure on the metal surface accompanied by induced fit. We anticipate that the mechanistic features of protein adsorption to an Au(111) surface revealed here can be extended to other inorganic surfaces and proteins and will therefore aid the design of specific protein-surface interactions.

Background Cytochrome P450 sterol 14α-demethylase (CYP51) is an essential enzyme for sterol biosynthesis and a target for anti-parasitic drug design. However, the design of parasite-specific drugs that inhibit parasitic CYP51 without severe side effects remains challenging. The active site of CYP51 is situated in the interior of the protein. Here, we characterize the potential ligand egress routes and mechanisms in Trypanosoma brucei and human CYP51 enzymes. Methods We performed Random Acceleration Molecular Dynamics simulations of the egress of four different ligands from the active site of models of soluble and membrane-bound T. brucei CYP51 and of soluble human CYP51. Results In the simulations, tunnel 2 f, which leads to the membrane, was found to be the predominant ligand egress tunnel for all the ligands studied. Tunnels S, 1 and W, which lead to the cytosol, were also used in T. brucei CYP51, whereas tunnel 1 was the only other tunnel used significantly in human CYP51. The common tunnels found previously in other CYPs were barely used. The ligand egress times were shorter for human than T. brucei CYP51, suggesting lower barriers to ligand passage. Two gating residues, F105 and M460, in T. brucei CYP51 that modulate the opening of tunnels 2 f and S were identified. Conclusions Although the main egress tunnel was the same, differences in the tunnel-lining residues, ligand passage and tunnel usage were found between T. brucei and human CYP51s. General Significance The results provide a basis for the design of selective anti-parasitic agents targeting the ligand tunnels.

A growing consensus is emerging that optimizing the drug–target affinity alone under equilibrium conditions does not necessarily translate into higher potency in vivo and that instead binding kinetic parameters should be optimized to ensure better efficacy. Therefore, in silico methods are needed to predict the kinetic parameters and the mechanistic determinants of drug–protein binding. Here we demonstrate the application of COMparative BINding Energy (COMBINE) analysis to derive quantitative structure–kinetics relationships (QSKRs) for the dissociation rate constants (koff) of inhibitors of heat shock protein 90 (HSP90) and HIV-1 protease. We derived protein-specific scoring functions by correlating koff rate constants with a subset of weighted interaction energy components determined from the energy-minimized structures of drug–protein complexes. As the QSKRs derived for these sets of chemically diverse compounds have good predictive ability and provide insights into important drug–protein interactions for optimizing koff, COMBINE analysis offers a promising approach for binding kinetics-guided lead optimization.

Background: S100A1 is an EF-hand calcium (Ca 2+ ) sensor protein with predominant expression in the mammalian heart. Although its outstanding influence on cardiac performance is well described, only few target molecules are known so far, such as the sarcoplasmic reticulum Ca 2+ ATPase (SERCA2a). Objective and Results: The discrepancy between the known interactors and the diversity of S100A1`s functions brings up the question which key targets are undetected so far and can, in addition to the known ones, explain the importance of S100A1 for the heart function. S100A1 knock-out (SKO -/- ) mouse heart extracts and lysates of isolated SKO -/- cardiomyocytes were subjected to HPLC S100A1 affinity-purification (AP) columns coupled to tandem mass spectrometry (MS) followed by network modeling of the interactome, to identify novel and to confirm previously known targets in the heart. Human recombinant S100A1 protein AP columns were used to trap EGTA- and Ca 2+ -dependent interactions. Subsequent MS and computational analyses revealed more than 90 novel cardiac interactors for S100A1 in a strict Ca 2+ -dependent manner, while previously known targets were confirmed. We showed that helix 4 of S100A1 (S100A1ct) alone can improve main parameters for the heart performance and might thereby function as a therapeutic biologic. Investigating the molecular mechanisms of S100A1`s interactions with its target molecules might help to understand this effect and unveil relevant regions for other therapeutic effects. In order to uncover the mechanisms of the various interactions, docking experiments were performed to gain an understanding, which areas of S100A1 and S100A1ct might be relevant for the biological functions on selected targets. Thus we detected a region on SERCA2a where S100A1 as well as S100A1ct might compete with the SERCA2a inhibitor phospholamban and by that may improve the SERCA2a activity. These in silico findings will be combined with experimental investigations of the mode of interaction by Cross-link MS experiments. Conclusion: Combining target screening via AP-MS/MS with computer-supported epitope mapping provides a powerful strategy to understand the molecular action of known therapies and by that helps to develop new treatments against heart failure.

Drug-target residence time (τ), one of the main determinants of drug efficacy, remains highly challenging to predict computationally and, therefore, is usually not considered in the early stages of drug design. Here, we present an efficient computational method, τ-random acceleration molecular dynamics (τRAMD), for the ranking of drug candidates by their residence time and obtaining insights into ligand-target dissociation mechanisms. We assessed τRAMD on a data set of 70 diverse drug-like ligands of the N-terminal domain of HSP90α, a pharmaceutically important target with a highly flexible binding site, obtaining computed relative residence times with an accuracy of about 2.3τ for 78% of the compounds and less than 2.0τ within congeneric series. Analysis of dissociation trajectories reveals features that affect ligand unbinding rates, including transient polar interactions and steric hindrance. These results suggest that τRAMD will be widely applicable as a computationally efficient aid to improving drug residence times during lead optimization.

Prolonged drug residence times may result in longer‐lasting drug efficacy, improved pharmacodynamic properties, and “kinetic selectivity” over off‐targets with high drug dissociation rates. However, few strategies have been elaborated to rationally modulate drug residence time and thereby to integrate this key property into the drug development process. Herein, we show that the interaction between a halogen moiety on an inhibitor and an aromatic residue in the target protein can significantly increase inhibitor residence time. By using the interaction of the serine/threonine kinase haspin with 5‐iodotubercidin (5‐iTU) derivatives as a model for an archetypal active‐state (type I) kinase–inhibitor binding mode, we demonstrate that inhibitor residence times markedly increase with the size and polarizability of the halogen atom. The halogen–aromatic π interactions in the haspin–inhibitor complexes were characterized by means of kinetic, thermodynamic, and structural measurements along with binding‐energy calculations.

Kinetoplastid parasites—trypanosomes and leishmanias—infect millions of humans and cause economically devastating diseases of livestock, and the few existing drugs have serious deficiencies. Benzoxaborole-based compounds are very promising potential novel anti-trypanosomal therapies, with candidates already in human and animal clinical trials. We investigated the mechanism of action of several benzoxaboroles, including AN7973, an early candidate for veterinary trypanosomosis. In all kinetoplastids, transcription is polycistronic. Individual mRNA 5'-ends are created by trans splicing of a short leader sequence, with coupled polyadenylation of the preceding mRNA. Treatment of Trypanosoma brucei with AN7973 inhibited trans splicing within 1h, as judged by loss of the Y-structure splicing intermediate, reduced levels of mRNA, and accumulation of peri-nuclear granules. Methylation of the spliced leader precursor RNA was not affected, but more prolonged AN7973 treatment caused an increase in S-adenosyl methionine and methylated lysine. Together, the results indicate that mRNA processing is a primary target of AN7973. Polyadenylation is required for kinetoplastid trans splicing, and the EC50 for AN7973 in T. brucei was increased three-fold by over-expression of the T. brucei cleavage and polyadenylation factor CPSF3, identifying CPSF3 as a potential molecular target. Molecular modeling results suggested that inhibition of CPSF3 by AN7973 is feasible. Our results thus chemically validate mRNA processing as a viable drug target in trypanosomes. Several other benzoxaboroles showed metabolomic and splicing effects that were similar to those of AN7973, identifying splicing inhibition as a common mode of action and suggesting that it might be linked to subsequent changes in methylated metabolites. Granule formation, splicing inhibition and resistance after CPSF3 expression did not, however, always correlate and prolonged selection of trypanosomes in AN7973 resulted in only 1.5-fold resistance. It is therefore possible that the modes of action of oxaboroles that target trypanosome mRNA processing might extend beyond CPSF3 inhibition., Author summary Trypanosomes and leishmanias infect millions of humans and cause economically devastating diseases of livestock; the few existing drugs have serious deficiencies. Trypanosomosis of cattle caused by Trypanosoma congolense and Trypanosoma vivax, is a serious problem in Africa because cattle are used not only for food but also for traction, and new drugs are needed. A single injection of the benzoxaborole compound AN7973 cured T. congolense infection in cattle and goats. Although slightly lower effectiveness against T. vivax precludes development of AN7973 as a commercially viable treatment against cattle trypanosomosis, it could still have potential for diseases caused by other trypanosomes. We used a large range of methods to find out how AN7973 kills trypanosomes, and compared it with other benzoxaboroles. AN7973 and some of the other compounds had effects on parasite metabolism that resembled those previously seen for a benzoxaborole that is being tested for human sleeping sickness. The most rapid effect of AN7973, however, was on processing of trypanosome mRNA. As a consequence, amounts of mRNA decreased and synthesis of proteins stopped. We conclude that AN7973 and some other benzoxaboroles kill trypanosomes by stopping gene expression.

Prolonged drug residence times may result in longer lasting drug efficacy, improved pharmacodynamic properties and “kinetic selectivity” over off-targets with fast drug dissociation rates. However, few strategies have been elaborated to rationally modulate drug residence time and thereby to integrate this key property into the drug development process. Here, we show that the interaction between a halogen moiety on an inhibitor and an aromatic residue in the target protein can significantly increase inhibitor residence time. By using the interaction of the serine/threonine kinase haspin with 5-iodotubercidin (5-iTU) derivatives as a model for an archetypal active state (type I) kinase-inhibitor binding mode, we demonstrate that inhibitor residence times markedly increase with the size and polarizability of the halogen atom. This key interaction is dependent on the interactions with an aromatic residue in the gate keeper position and we observe this interaction in other kinases with an aromatic gate keeper residue. We provide a detailed mechanistic characterization of the halogen-aromatic π interactions in the haspin-inhibitor complexes by means of kinetic, thermodynamic, and structural measurements along with binding energy calculations. Since halogens are frequently used in drugs and aromatic residues are often present in the binding sites of proteins, our results provide a compelling rationale for introducing aromatic-halogen interactions to prolong drug-target residence times.

Structure-based drug design has often been restricted by the rather static picture of protein–ligand complexes presented by crystal structures, despite the widely accepted importance of protein flexibility in biomolecular recognition. Here we report a detailed experimental and computational study of the drug target, human heat shock protein 90, to explore the contribution of protein dynamics to the binding thermodynamics and kinetics of drug-like compounds. We observe that their binding properties depend on whether the protein has a loop or a helical conformation in the binding site of the ligand-bound state. Compounds bound to the helical conformation display slow association and dissociation rates, high-affinity and high cellular efficacy, and predominantly entropically driven binding. An important entropic contribution comes from the greater flexibility of the helical relative to the loop conformation in the ligand-bound state. This unusual mechanism suggests increasing target flexibility in the bound state by ligand design as a new strategy for drug discovery., An understanding of the dynamics of drug binding and unbinding processes is important for drug discovery. Here, the authors give insights into the binding mechanism of small drug-like molecules to human Hsp90 by combining thermodynamics and kinetics studies as well as molecular dynamics simulations.

The conformational properties of unbound multi-Cys2 His2 (mC2H2) zinc finger proteins, in which zinc finger domains are connected by flexible linkers, are studied by a multiscale approach. Three methods on different length scales are utilized. First, atomic detail molecular dynamics simulations of one zinc finger and its adjacent flexible linker confirmed that the zinc finger is more rigid than the flexible linker. Second, the end-to-end distance distributions of mC2H2 zinc finger proteins are computed using an efficient atomistic pivoting algorithm, which only takes excluded volume interactions into consideration. The end-to-end distance distribution gradually changes its profile, from left-tailed to right-tailed, as the number of zinc fingers increases. This is explained by using a worm-like chain model. For proteins of a few zinc fingers, an effective bending constraint favors an extended conformation. Only for proteins containing more than nine zinc fingers, is a somewhat compacted conformation preferred. Third, a mesoscale model is modified to study both the local and the global conformational properties of multi-C2H2 zinc finger proteins. Simulations of the CCCTC-binding factor (CTCF), an important mC2H2 zinc finger protein for genome spatial organization, are presented.

The simulation of diffusional association (SDA) Brownian dynamics software package has been widely used in the study of biomacromolecular association. Initially developed to calculate bimolecular protein–protein association rate constants, it has since been extended to study electron transfer rates, to predict the structures of biomacromolecular complexes, to investigate the adsorption of proteins to inorganic surfaces, and to simulate the dynamics of large systems containing many biomacromolecular solutes, allowing the study of concentration‐dependent effects. These extensions have led to a number of divergent versions of the software. In this article, we report the development of the latest version of the software (SDA 7). This release was developed to consolidate the existing codes into a single framework, while improving the parallelization of the code to better exploit modern multicore shared memory computer architectures. It is built using a modular object‐oriented programming scheme, to allow for easy maintenance and extension of the software, and includes new features, such as adding flexible solute representations. We discuss a number of application examples, which describe some of the methods available in the release, and provide benchmarking data to demonstrate the parallel performance. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Kar1 interacts with the C-terminal (CT) region of Sfi1 to tether the yeast centrosome’s bridge to the nuclear envelope; Sfi1-CT and C-terminal Cdc31-binding regions form the antiparallel Sfi1 overlap in the bridge., The yeast spindle pole body (SPB) is the functional equivalent of the mammalian centrosome. The half bridge is a SPB substructure on the nuclear envelope (NE), playing a key role in SPB duplication. Its cytoplasmic components are the membrane-anchored Kar1, the yeast centrin Cdc31, and the Cdc31-binding protein Sfi1. In G1, the half bridge expands into the bridge through Sfi1 C-terminal (Sfi1-CT) dimerization, the licensing step for SPB duplication. We exploited photo-activated localization microscopy (PALM) to show that Kar1 localizes in the bridge center. Binding assays revealed direct interaction between Kar1 and C-terminal Sfi1 fragments. kar1Δ cells whose viability was maintained by the dominant CDC31-16 showed an arched bridge, indicating Kar1’s function in tethering Sfi1 to the NE. Cdc31-16 enhanced Cdc31–Cdc31 interactions between Sfi1–Cdc31 layers, as suggested by binding free energy calculations. In our model, Kar1 binding is restricted to Sfi1-CT and Sfi1 C-terminal centrin-binding repeats, and centrin and Kar1 provide cross-links, while Sfi1-CT stabilizes the bridge and ensures timely SPB separation.

Human thymidylate synthase (hTS), a target for antiproliferative drugs, is an obligate homodimer. Single-point mutations to alanine at the monomer–monomer interface may enable the identification of specific residues that delineate sites for drugs aimed at perturbing the protein–protein interactions critical for activity. We computationally identified putative hotspot residues at the interface and designed mutants to perturb the intersubunit interaction. Dimer dissociation constants measured by a FRET-based assay range from 60 nM for wild-type hTS up to about 1 mM for single-point mutants and agree with computational predictions of the effects of these mutations. Mutations that are remote from the active site retain full or partial activity, although the substrate KM values were generally higher and the dimer was less stable. The lower dimer stability of the mutants can facilitate access to the dimer interface by small molecules and thereby aid the design of inhibitors that bind at the dimer interface.Graphical abstract

Sterol 14α-demethylase (cytochrome P450 family 51 (CYP51)) is an essential enzyme occurring in all biological kingdoms. In eukaryotes, it is located in the membrane of the endoplasmic reticulum. Selective inhibitors of trypanosomal CYP51s that do not affect the human CYP51 have been discovered in vitro and found to cure acute and chronic mouse Chagas disease without severe side effects in vivo. Crystal structures indicate that CYP51 may be more rigid than most CYPs, and it has been proposed that this property may facilitate antiparasitic drug design. Therefore, to investigate the dynamics of trypanosomal CYP51, we built a model of membrane-bound Trypanosoma brucei CYP51 and then performed molecular dynamics simulations of T. brucei CYP51 in membrane-bound and soluble forms. We compared the dynamics of T. brucei CYP51 with those of human CYP51, CYP2C9, and CYP2E1. In the simulations, the CYP51s display low mobility in the buried active site although overall mobility is similar in all the CYPs studied. The simulations suggest that in CYP51, pathway 2f serves as the major ligand access tunnel, and both pathways 2f (leading to membrane) and S (leading to solvent) can serve as ligand egress tunnels. Compared with the other CYPs, the residues at the entrance of the ligand access tunnels in CYP51 have higher mobility that may be necessary to facilitate the passage of its large sterol ligands. The water (W) tunnel is accessible to solvent during most of the simulations of CYP51, but its width is affected by the conformations of the heme's two propionate groups. These differ from those observed in the other CYPs studied because of differences in their hydrogen-bonding network. Our simulations give insights into the dynamics of CYP51 that complement the available experimental data and have implications for drug design against CYP51 enzymes.

Profilin is an actin monomer binding protein that provides ATP-actin for incorporation into actin filaments. In contrast to higher eukaryotic cells with their large filamentous actin structures, apicomplexan parasites typically contain only short and highly dynamic microfilaments. In apicomplexans, profilin appears to be the main monomer-sequestering protein. Compared to classical profilins, apicomplexan profilins contain an additional arm-like β-hairpin motif, which we show here to be critically involved in actin binding. Through comparative analysis using two profilin mutants, we reveal this motif to be implicated in gliding motility of Plasmodium berghei sporozoites, the rapidly migrating forms of a rodent malaria parasite transmitted by mosquitoes. Force measurements on migrating sporozoites and molecular dynamics simulations indicate that the interaction between actin and profilin fine-tunes gliding motility. Our data suggest that evolutionary pressure to achieve efficient high-speed gliding has resulted in a unique profilin-actin interface in these parasites., Author summary The malaria parasite Plasmodium has two invasive forms that migrate across different tissue barriers, the ookinete and the very rapidly migrating sporozoite. Previous work has shown that the motility of these and related parasites (e.g. Toxoplasma gondii) depends on a highly dynamic actin cytoskeleton and retrograde flow of surface adhesins. These unusual actin dynamics are due to the divergent structure of protozoan actins and the actions of actin-binding proteins, which can have non-canonical functions in these parasites. Profilin is one of the most important and most investigated actin-binding proteins, which binds ADP-actin and catalyzes ADP-ATP exchange to then promote actin polymerization. Parasite profilins bind monomeric actin and contain an additional domain compared to canonical profilins. Here we show that this additional domain of profilin is critical for actin binding and rapid sporozoite motility but has little impact on the slower ookinete. Sporozoites of a parasite line carrying mutations in this domain cannot translate force production and retrograde flow into optimal parasite motility. Using molecular dynamics simulations, we find that differences between mutant parasites in their capacity to migrate can be traced back to a single hydrogen bond at the actin-profilin interface.

A considerable number of approved drugs show non-equilibrium binding characteristics, emphasizing the potential role of drug residence times for in vivo efficacy. Therefore, a detailed understanding of the kinetics of association and dissociation of a target-ligand complex might provide crucial insight into the molecular mechanism-of-action of a compound. This deeper understanding will help to improve decision making in drug discovery, thus leading to a better selection of interesting compounds to be profiled further. In this review, we highlight the contributions of the Kinetics for Drug Discovery (K4DD) Consortium, which targets major open questions related to binding kinetics in an industry-driven public-private partnership.

ARL13B encodes for the ADP-ribosylation factor-like 13B GTPase, which is required for normal cilia structure and Sonic hedgehog (Shh) signaling. Disruptions in cilia structure or function lead to a class of human disorders called ciliopathies. Joubert syndrome is characterized by a wide spectrum of symptoms, including a variable degree of intellectual disability, ataxia, and ocular abnormalities. Here we report a novel homozygous missense variant c.[223G>A] (p.(Gly75Arg) in the ARL13B gene, which was identified by whole-exome sequencing of a trio from a consanguineous family with multiple-affected individuals suffering from intellectual disability, ataxia, ocular defects, and epilepsy. The same variant was also identified in a second family. We saw a striking difference in the severity of ataxia between affected male and female individuals in both families. Both ARL13B and ARL13B-c.[223G>A] (p.(Gly75Arg) expression rescued the cilia length and Shh defects displayed by Arl13b hennin (null) cells, indicating that the variant did not disrupt either ARL13B function. In contrast, ARL13B-c.[223G>A] (p.(Gly75Arg) displayed a marked loss of ARL3 guanine nucleotide-exchange factor activity, with retention of its GTPase activities, highlighting the correlation between its loss of function as an ARL3 guanine nucleotide-exchange factor and Joubert syndrome.

Flavonoids have previously been identified as antiparasitic agents and pteridine reductase 1 (PTR1) inhibitors. Herein, we focus our attention on the chroman-4-one scaffold. Three chroman-4-one analogues (1-3) of previously published chromen-4-one derivatives were synthesized and biologically evaluated against parasitic enzymes (Trypanosoma brucei PTR1-TbPTR1 and Leishmania major-LmPTR1) and parasites (Trypanosoma brucei and Leishmania infantum). A crystal structure of TbPTR1 in complex with compound 1 and the first crystal structures of LmPTR1-flavanone complexes (compounds 1 and 3) were solved. The inhibitory activity of the chroman-4-one and chromen-4-one derivatives was explained by comparison of observed and predicted binding modes of the compounds. Compound 1 showed activity both against the targeted enzymes and the parasites with a selectivity index greater than 7 and a low toxicity. Our results provide a basis for further scaffold optimization and structure-based drug design aimed at the identification of potent anti-trypanosomatidic compounds targeting multiple PTR1 variants.

Background Multi-target approaches are necessary to properly analyze or modify the function of a biochemical pathway or a protein family. An example of such a problem is the repurposing of the known human anti-cancer drugs, antifolates, as selective anti-parasitic agents. This requires considering a set of experimentally validated protein targets in the folate pathway of major pathogenic trypanosomatid parasites and humans: ( i ) the primary parasite on-targets: pteridine reductase 1 (PTR1) (absent in humans) and bifunctional dihydrofolate reductase-thymidylate synthase (DHFR–TS), ( ii ) the primary off-targets: human DHFR and TS, and ( iii ) the secondary on-target: human folate receptor β , a folate/antifolate transporter. Methods We computationally compared the structural, dynamic and physico-chemical properties of the targets. We based our analysis on available inhibitory activity and crystallographic data, including a crystal structure of the bifunctional T. cruzi DHFR–TS with tetrahydrofolate bound determined in this work. Due to the low sequence and structural similarity of the targets analyzed, we employed a mapping of binding pockets based on the known common ligands, folate and methotrexate. Results Our analysis provides a set of practical strategies for the design of selective trypanosomatid folate pathway inhibitors, which are supported by enzyme inhibition measurements and crystallographic structures. Conclusions The ligand-based comparative computational mapping of protein binding pockets provides a basis for repurposing of anti-folates and the design of new anti-trypanosmatid agents. General significance Apart from the target-based discovery of selective compounds, our approach may be also applied for protein engineering or analyzing evolutionary relationships in protein families.

Hsp70 participates in a broad spectrum of protein folding processes extending from nascent chain folding to protein disaggregation. This versatility in function is achieved through a diverse family of J-protein cochaperones that select substrates for Hsp70. Substrate selection is further tuned by transient complexation between different classes of J-proteins, which expands the range of protein aggregates targeted by metazoan Hsp70 for disaggregation. We assessed the prevalence and evolutionary conservation of J-protein complexation and cooperation in disaggregation. We find the emergence of a eukaryote-specific signature for interclass complexation of canonical J-proteins. Consistently, complexes exist in yeast and human cells, but not in bacteria, and correlate with cooperative action in disaggregation in vitro. Signature alterations exclude some J-proteins from networking, which ensures correct J-protein pairing, functional network integrity and J-protein specialization. This fundamental change in J-protein biology during the prokaryote-to-eukaryote transition allows for increased fine-tuning and broadening of Hsp70 function in eukaryotes. DOI: http://dx.doi.org/10.7554/eLife.24560.001, eLife digest All cells must maintain their proteins in a correctly folded shape to survive. The task of sustaining a healthy set of proteins has increased with the rise of complex life from prokaryotes (such as bacteria) that form simple single-celled organisms to eukaryotes (such as yeast, plants and multicellular animals). As a result of organisms ageing or acquiring genetic mutations, or under stressful conditions such as high temperature, proteins can lose their normal shape and clump together to form “aggregates”. These aggregates are potentially toxic to cells and have been linked to many human diseases including neurodegeneration and cancer. Cells contain molecular machines that help break down aggregates and subsequently recycle or rescue trapped proteins. Some of these machines are based around a protein called Hsp70, which can perform a wide range of protein folding processes. So-called J-proteins help Hsp70 to select aggregates to be targeted for break down. It used to be thought that different classes of J-proteins interacted with Hsp70 separately. However, in 2015, researchers showed that in humans, two different classes of J-proteins can bind to each other to form a “complex”, which has distinct aggregate selection properties. Now, Nillegoda et al. – including several of the researchers involved in the 2015 study – have examined the evolutionary history of these J-protein complexes. This revealed that different classes (A and B) of J-proteins first cooperated after prokaryotes and eukaryotes diverged from each other. In particular, the molecular machinery that breaks down aggregates in yeast cells – but not the machinery found in bacteria – depends on complexes formed from the two classes of J-proteins. Further investigation revealed that in humans, J-proteins have structural features that ensure they pair up correctly to perform unique activities. Furthermore, Nillegoda et al. suggest that cooperation between J-proteins may have enabled organisms such as humans – which contain over 40 distinct J-proteins – to carry out further specialized protein-folding tasks that do not occur in prokaryotes. Overall, the findings presented by Nillegoda et al. reveal another important layer to protein quality control in eukaryotic cells. The next step is to understand the possible roles of different J-protein complexes play in J-protein associated cellular protein quality control processes such as preventing protein aggregation, refolding or recycling abnormal proteins. This knowledge could ultimately be used to develop treatments for diseases and disorders in which protein aggregates form. DOI: http://dx.doi.org/10.7554/eLife.24560.002

This Special Issue contains thirteen articles that provide a vivid snapshot of the state-of-the-art of molecular modeling in drug design, illustrating recent advances and critically discussing important challenges [...]

Both subunits of αβ-tubulin that comprise the core components of microtubules bind GTP. GTP binding to α-tubulin has a structural role, whereas β-tubulin binds and hydrolyses GTP to regulate microtubule dynamics. γ-tubulin, another member of the tubulin superfamily that seeds microtubule nucleation at microtubule-organizing centres, also binds GTP; however, the importance of this association remains elusive. To address the role of GTP binding to γ-tubulin, we systematically mutagenized the GTP contact residues in the yeast γ-tubulin Tub4. Tub4(GTP)-mutant proteins that exhibited greatly reduced GTP affinity still assembled into the small γ-tubulin complex. However, tub4(GTP) mutants were no longer viable, and had defects in interaction between γ-tubulin and αβ-tubulin, decreased microtubule nucleation and defects in microtubule organization. In vitro and in vivo data show that only γ-tubulin loaded with GTP nucleates microtubules. Our results suggest that GTP recruitment to γ-tubulin enhances its interaction with αβ-tubulin similarly to GTP recruitment to β-tubulin.

Newly synthesized polypeptides undergo various cotranslational maturation steps, including N-terminal enzymatic processing, chaperone-assisted folding and membrane targeting, but the spatial and temporal coordination of these steps is unclear. We show that Escherichia coli methionine aminopeptidase (MAP) associates with ribosomes through a charged loop that is crucial for nascent-chain processing and cell viability. MAP competes with peptide deformylase (PDF), the first enzyme to act on nascent chains, for binding sites at the ribosomal tunnel exit. PDF has extremely fast association and dissociation kinetics, which allows it to frequently sample ribosomes and ensure the processing of nascent chains after their emergence. Premature recruitment of the chaperone trigger factor, or polypeptide folding, negatively affect processing efficiency. Thus, the fast ribosome association kinetics of PDF and MAP are crucial for the temporal separation of nascent-chain processing from later maturation events, including chaperone recruitment and folding.

Phenolnaphthalein derivatives show potential for pharmacological activity as inhibitors of thymidylate synthase (TS) but difficulties in their synthesis and derivatization hinder their development. A deconstruction approach aimed at identifying a suitable new scaffold was proposed. A new scaffold was identified and two compound libraries based on this scaffold were designed. The carboxamide library (Library B) showed specific inhibition activity against Escherichia coli TS, whereas the sulfonamide library (Library C) showed a non-specific inhibition profile against hTS. N-(1,3-Dioxo-1H,3H-benzo[de]isochromen-5-yl)-sulfonamide derivatives, 1C and 9C, showed one order of magnitude improvement in inhibition constant against hTS with respect to the starting lead and represent potential compounds for further lead development.

Resistance to drugs targeting human thymidylate synthase (TS) poses a major challenge in the field of anti-cancer therapeutics. Overexpression of the TS protein has been implicated as one of the factors leading to the development of resistance. Therefore, repressing translation by targeting the TS mRNA could help to overcome this problem. In this study, we report that the compound Hoechst 33258 (HT) can reduce cellular TS protein levels without altering TS mRNA levels, suggesting that it modulates TS expression at the translation level. We have combined nuclear magnetic resonance, UV-visible and fluorescence spectroscopy methods with docking and molecular dynamics simulations to study the interaction of HT with a region in the TS mRNA. The interaction predominantly involves intercalation of HT at a CC mismatch in the region near the translational initiation site. Our results support the use of HT-like compounds to guide the design of therapeutic agents targeting TS mRNA.

The mammalian cytochrome P450 (CYP) enzymes play important roles in drug metabolism, steroid biosynthesis, and xenobiotic degradation. The active site of CYPs is buried in the protein and thus the ligands have to enter and exit the active site via ligand tunnels. Conformational changes of flexible parts of the protein usually accompany the entrance and exit of ligands. Comparison of the crystal structures of mammalian CYPs in closed, open, and partially open states reveals that the greatest conformational diversity associated with ligand tunnel opening is in the regions of the B-C and F-G loops. Some CYPs have been observed to adopt different open and closed conformations when bound to different ligands, suggesting that the ligand entrance and exit routes might differ according to the ligand properties. Mammalian CYPs are mostly membrane-bound enzymes, making them difficult to characterize structurally and dynamically. A range of molecular dynamics simulation techniques has been applied to investigate the dynamics and the ligand tunnels of these proteins both in the aqueous environment, and more recently, in lipid bilayers. These simulations not only reveal multiple tunnels through which ligands can pass but also show that different tunnels are preferred by different ligands and that the lipid bilayer can influence the protein dynamics and tunnel opening. The results indicate that not only the active site but also the ligand tunnels can contribute to the different substrate specificity profiles of the mammalian CYPs.

The ring-forming AAA+ protein ClpB cooperates with the DnaK chaperone system to refold aggregated proteins in Escherichia coli. The M domain, a ClpB-specific coiled-coil structure with two wings, motif 1 and motif 2, is essential to disaggregation, but the positioning and mechanistic role of M domains in ClpB hexamers remain unresolved. We show that M domains nestle at the ClpB ring surface, with both M-domain motifs contacting the first ATPase domain (AAA-1). Both wings contribute to maintaining a repressed ClpB activity state. Motif 2 docks intramolecularly to AAA-1 to regulate ClpB unfolding power, and motif 1 contacts a neighboring AAA-1 domain. Mutations that stabilize motif 2 docking repress ClpB, whereas destabilization leads to derepressed ClpB activity with greater unfolding power that is toxic in vivo. Our results underline the vital nature of tight ClpB activity control and elucidate a regulated M-domain toggle control mechanism.

There is a need for improved and generally applicable scoring functions for fragment-based approaches to ligand design. Here, we evaluate the performance of a computationally efficient model for inhibitory activity estimation, which is composed only of multipole electrostatic energy and dispersion energy terms that approximate long-range ab initio quantum mechanical interaction energies. We find that computed energies correlate well with inhibitory activity for a compound series with varying substituents targeting two subpockets of the binding site of Trypanosoma brucei pteridine reductase 1. For one subpocket, we find that the model is more predictive for inhibitory activity than the ab initio interaction energy calculated at the MP2 level. Furthermore, the model is found to outperform a commonly used empirical scoring method. Finally, we show that the results for the two subpockets can be combined, which suggests that this simple nonempirical scoring function could be applied in fragment–based drug design. Electronic supplementary material The online version of this article (doi:10.1007/s10822-017-0035-4) contains supplementary material, which is available to authorized users.

We present the results for CAPRI Round 30, the first joint CASP-CAPRI experiment, which brought together experts from the protein structure prediction and protein–protein docking communities. The Round comprised 25 targets from amongst those submitted for the CASP11 prediction experiment of 2014. The targets included mostly homodimers, a few homotetramers, and two heterodimers, and comprised protein chains that could readily be modeled using templates from the Protein Data Bank. On average 24 CAPRI groups and 7 CASP groups submitted docking predictions for each target, and 12 CAPRI groups per target participated in the CAPRI scoring experiment. In total more than 9500 models were assessed against the 3D structures of the corresponding target complexes. Results show that the prediction of homodimer assemblies by homology modeling techniques and docking calculations is quite successful for targets featuring large enough subunit interfaces to represent stable associations. Targets with ambiguous or inaccurate oligomeric state assignments, often featuring crystal contact-sized interfaces, represented a confounding factor. For those, a much poorer prediction performance was achieved, while nonetheless often providing helpful clues on the correct oligomeric state of the protein. The prediction performance was very poor for genuine tetrameric targets, where the inaccuracy of the homology-built subunit models and the smaller pair-wise interfaces severely limited the ability to derive the correct assembly mode. Our analysis also shows that docking procedures tend to perform better than standard homology modeling techniques and that highly accurate models of the protein components are not always required to identify their association modes with acceptable accuracy. We are most grateful to the PDBe at the European Bioinformatics Institute in Hinxton, UK, for hosting the CAPRI website. Our deepest thanks go to all the structural biologists and to the following structural genomics initiatives: Northeast Structural Genomics Consortium, Joint Center for Structural Genomics, NatPro PSI:Biology, New York Structural Genomics Research Center, Midwest Center for Structural Genomics, Structural Genomics Consortium, for contributing the targets for this joint CASP-CAPRI experiment. MFL acknowledges support from the FRABio FR3688 Research Federation “Structural & Functional Biochemistry of Biomolecular Assemblies.”

Computer simulations of molecular systems often make use of regular rectangular grids with equidistant spacing to store information on their molecular interaction fields, e.g., electrostatic potential. These grids provide an easy way to store the data as they do not require any particular specification of the structure of the data. However, such grids may easily become large, and the storage and memory demands may become so high that calculations become infeasible. To overcome this problem, we show here how the data structure DT-Grid can be adapted and applied to efficiently represent macromolecular interaction grids by exploiting the nonuniformity of information on the grid; at the same time, this data structure ensures fast random data access. We demonstrate use of the DT-Grid data structure for potential of mean force and Brownian dynamics simulations of protein-surface binding and macromolecular association with the SDA software. We further demonstrate that the DT-Grid structure enables systems of large size, such as a viral capsid, and high resolution grids to be handled that are beyond current computational feasibility.

Simulations of the long-time scale motions of a ligand binding pocket in a protein may open up new perspectives for the design of compounds with steric or chemical properties differing from those of known binders. However, slow motions of proteins are difficult to access using standard molecular dynamics (MD) simulations and are thus usually neglected in computational drug design. Here, we introduce two nonequilibrium MD approaches to identify conformational changes of a binding site and detect transient pockets associated with these motions. The methods proposed are based on the rotamerically induced perturbation (RIP) MD approach, which employs perturbation of side-chain torsional motion for initiating large-scale protein movement. The first approach, Langevin-RIP (L-RIP), entails a series of short Langevin MD simulations, each starting with perturbation of one of the side-chains lining the binding site of interest. L-RIP provides extensive sampling of conformational changes of the binding site. In less than 1 ns of MD simulation with L-RIP, we observed distortions of the α-helix in the ATP binding site of HSP90 and flipping of the DFG loop in Src kinase. In the second approach, RIPlig, a perturbation is applied to a pseudoligand placed in different parts of a binding pocket, which enables flexible regions of the binding site to be identified in a small number of 10 ps MD simulations. The methods were evaluated for four test proteins displaying different types and degrees of binding site flexibility. Both methods reveal all transient pocket regions in less than a total of 10 ns of simulations, even though many of these regions remained closed in 100 ns conventional MD. The proposed methods provide computationally efficient tools to explore binding site flexibility and can aid in the functional characterization of protein pockets, and the identification of transient pockets for ligand design.