Search

Cyclic peptides represent a promising class of drug candidates. A significant obstacle limiting their development as therapeutics is the lack of an ability to predict their membrane permeability. We use molecular dynamics simulations to assess the ability of a set of widely used parameters in describing the membrane permeability of a set of model cyclic peptides; the parameters include polar surface area (PSA), the number of hydrogen bonds, and transfer free energy between an aqueous phase and a membrane mimicking phase. These parameters were found to generally correlate with the membrane permeability of the set of cyclic peptides. We propose two new descriptors, the charge reweighted PSA and the non-polar surface area to PSA ratio; both show enhanced correlation with membrane permeability. This inspired us to explore crosslinking of the peptide to reduce the accessible surface area of the backbone polar atoms, and we find that this can indeed result in reductions in the accessible PSA. This gives reason to speculate that crosslinking may result in increased permeability, thus suggesting a new scaffold for the development of cyclic peptides as potential therapeutics. [ABSTRACT FROM AUTHOR]

The N‐terminal transactivation domain (TAD) of p53 is a disordered region with multiple phosphorylation sites. Phosphorylation at Thr18 is crucial for the release of p53 from its negative regulator, MDM2. In stressed cells, CK1δ is responsible for phosphorylating Thr18, but requires Ser15 to be phosphorylated. To understand the mechanistic underpinnings of this sequential phosphorylation, molecular modeling and molecular dynamics simulation studies of these phosphorylation events were carried out. Our models suggest that a positively charged region on CK1δ near the adenosine triphosphate (ATP) binding pocket, which is conserved across species, sequesters the negatively charged pSer15, thereby constraining the positioning of the rest of the peptide, such that the side chain of Thr18 is positioned close to the γ‐phosphate of ATP. Furthermore, our studies show that the phosphorylated p53 TAD1 (p53pSer15) peptide binds more strongly to CK1δ than does p53. p53 adopts a helical structure when bound to CK1δ, which is lost upon phosphorylation at Ser15, thus gaining higher flexibility and ability to morph into the binding site. We propose that upon phosphorylation at Ser15 the p53 TAD1 peptide binds to CK1δ through an electrostatically driven induced fit mechanism resulting in a flanking fuzzy complex. [ABSTRACT FROM AUTHOR]

Background/Aim: Host defense peptides (HDPs) have the potential to provide a novel solution to antimicrobial resistance (AMR) in view of their unique and broad-spectrum antimicrobial activities. We had recently developed a novel hybrid HDP based on LL-37 and human beta-defensin-2, named CaD23, which was shown to exhibit good in vivo antimicrobial efficacy against Staphylococcus aureus in a bacterial keratitis murine model. This study aimed to examine the potential CaD23-antibiotic synergism and the secondary structure and underlying mechanism of action of CaD23. Methods: Peptide-antibiotic interaction was evaluated against S. aureus , methicillin-resistant S. aureus (MRSA), and Pseudomonas aeruginosa using established checkerboard and time-kill assays. Fractional inhibitory concentration index (FICI) was calculated and interpreted as synergistic (FIC<0.5), additive (FIC between 0.5–1.0), indifferent (FIC between >1.0 and ≤4), or antagonistic (FIC>4). SYTOX green uptake assay was performed to determine the membrane-permeabilising action of CaD23. Molecular dynamics (MD) simulations were performed to evaluate the interaction of CaD23 with bacterial and mammalian mimetic membranes. Circular dichroism (CD) spectroscopy was also performed to examine the secondary structures of CaD23. Results: CaD23-amikacin and CaD23-levofloxacin combination treatment exhibited a strong additive effect against S. aureus SH1000 (FICI = 0.60–0.69) and MRSA43300 (FICI = 0.56–0.60) but an indifferent effect against P. aeruginosa (FIC = 1.03–1.15). CaD23 (at 25 μg/ml; 2xMIC) completely killed S. aureus within 30 min. When used at sub-MIC concentration (3.1 μg/ml; 0.25xMIC), it was able to expedite the antimicrobial action of amikacin against S. aureus by 50%. The rapid antimicrobial action of CaD23 was attributed to the underlying membrane-permeabilising mechanism of action, evidenced by the SYTOX green uptake assay and MD simulations studies. MD simulations revealed that cationicity, alpha-helicity, amphiphilicity and hydrophobicity (related to the Trp residue at C-terminal) play important roles in the antimicrobial action of CaD23. The secondary structures of CaD23 observed in MD simulations were validated by CD spectroscopy. Conclusion: CaD23 is a novel alpha-helical, membrane-active synthetic HDP that can enhance and expedite the antimicrobial action of antibiotics against Gram-positive bacteria when used in combination. MD simulations serves as a powerful tool in revealing the peptide secondary structure, dissecting the mechanism of action, and guiding the design and optimisation of HDPs. [ABSTRACT FROM AUTHOR]

Heterozygous HNF1A gene mutations can cause maturity onset diabetes of the young 3 (MODY3), characterized by insulin secretion defects. However, specific mechanisms of MODY3 in humans remain unclear due to lack of access to diseased human pancreatic cells. Here, we utilize MODY3 patient-derived human induced pluripotent stem cells (hiPSCs) to study the effect(s) of a causal HNF1A+/H126D mutation on pancreatic function. Molecular dynamics simulations predict that the H126D mutation could compromise DNA binding and gene target transcription. Genome-wide RNA-Seq and ChIP-Seq analyses on MODY3 hiPSC-derived endocrine progenitors reveal numerous HNF1A gene targets affected by the mutation. We find decreased glucose transporter GLUT2 expression, which is associated with reduced glucose uptake and ATP production in the MODY3 hiPSC-derived β-like cells. Overall, our findings reveal the importance of HNF1A in regulating GLUT2 and several genes involved in insulin secretion that can account for the insulin secretory defect clinically observed in MODY3 patients. Heterozygous HNF1A mutations can give rise to maturity onset diabetes of the young 3 (MODY3), characterized by insulin secretion defects. Here the authors show that MODY3-related HNF1A mutation in patient hiPSCderived pancreatic cells decreases glucose transporter GLUT2 expression due to compromised DNA binding. [ABSTRACT FROM AUTHOR]

Barnacles employ a protein-based cement to firmly attach to immersed substrates. The cement proteins (CPs) have previously been identified and sequenced. However, the molecular mechanisms of adhesion are not well understood, in particular, because the three-dimensional molecular structure of CPs remained unknown to date. Here, we conducted multi-dimensional nuclear magnetic resonance (NMR) studies and molecular dynamics (MD) simulations of recombinant Megabalanus rosa Cement Protein 20 (rMrCP20). Our NMR results show that rMrCP20 contains three main folded domain regions intervened by two dynamic loops, resulting in multiple protein conformations that exist in equilibrium. We found that 12 out of 32 Cys in the sequence engage in disulfide bonds that stabilize the β-sheet domains owing to their placement at the extremities of β-strands. Another feature unveiled by NMR is the location of basic residues in turn regions that are exposed to the solvent, playing an important role for intermolecular contact with negatively charged surfaces. MD simulations highlight a highly stable and conserved β-motif (β7-β8), which may function as nuclei for amyloid-like nanofibrils previously observed in the cured adhesive cement. To the best of our knowledge, this is the first report describing the tertiary structure of an extracellular biological adhesive protein at the molecular level. This article is part of the theme issue 'Transdisciplinary approaches to the study of adhesion and adhesives in biological systems'. [ABSTRACT FROM AUTHOR]

The ability of DENV2 to display different morphologies (hence different antigenic properties) complicates vaccine and therapeutics development. Previous studies showed most strains of laboratory adapted DENV2 particles changed from smooth to “bumpy” surfaced morphology when the temperature is switched from 29°C at 37°C. Here we identified five envelope (E) protein residues different between two alternative passage history DENV2 NGC strains exhibiting smooth or bumpy surface morphologies. Several mutations performed on the smooth DENV2 infectious clone destabilized the surface, as observed by cryoEM. Molecular dynamics simulations demonstrated how chemically subtle substitution at various positions destabilized dimeric interactions between E proteins. In contrast, three out of four DENV2 clinical isolates showed a smooth surface morphology at 37°C, and only at high fever temperature (40°C) did they become “bumpy”. These results imply vaccines should contain particles representing both morphologies. For prophylactic and therapeutic treatments, this study also informs on which types of antibodies should be used at different stages of an infection, i.e., those that bind to monomeric E proteins on the bumpy surface or across multiple E proteins on the smooth surfaced virus. [ABSTRACT FROM AUTHOR]

• CPPs can partake in multiple cellular uptake mechanisms. • MD simulations can be used to gain atomistic insight into CPP uptake mechanisms. • Biomolecular MD simulations are challenging owing to complex energy landscapes. • Enhanced sampling methods can increase sampling of the energy landscape. Cell-penetrating peptides (CPPs) offer an exciting approach to tackle the pharmacokinetic challenges associated with the delivery of large, polar molecules to intracellular targets. Since the discovery of the first CPPs in the early 1990s, vast amounts of research have been undertaken to characterise their cellular uptake mechanisms. Despite this, the precise mechanisms of cellular internalisation of many CPPs remain elusive owing to inconsistent experimental results. Molecular dynamics (MD) simulations provide an approach to probe specific aspects of the internalisation process and many published CPP studies have incorporated simulation data. This review provides a critical evaluation of the current approaches that are being used to simulate CPPs interacting with artificial lipid bilayers. This review discusses the simulation techniques currently being used to study cell-penetrating peptides, along with the mechanistic insights gained and challenges associated with such studies. [ABSTRACT FROM AUTHOR]

R248 in the DNA binding domain (DBD) of p53 interacts directly with the minor groove of DNA. Earlier nuclear magnetic resonance (NMR) studies indicated that the R248Q mutation resulted in conformation changes in parts of DBD far from the mutation site. However, how information propagates from the mutation site to the rest of the DBD is still not well understood. We performed a series of all-atom molecular dynamics (MD) simulations to dissect sterics and charge effects of R248 on p53-DBD conformation: (i) wild-type p53 DBD; (ii) p53 DBD with an electrically neutral arginine side-chain; (iii) p53 DBD with R248A; (iv) p53 DBD with R248W; and (v) p53 DBD with R248Q. Our results agree well with experimental observations of global conformational changes induced by the R248Q mutation. Our simulations suggest that both charge- and sterics are important in the dynamics of the loop (L3) where the mutation resides. We show that helix 2 (H2) dynamics is altered as a result of a change in the hydrogen bonding partner of D281. In turn, neighboring L1 dynamics is altered: in mutants, L1 predominantly adopts the recessed conformation and is unable to interact with the major groove of DNA. We focused our attention the R248Q mutant that is commonly found in a wide range of cancer and observed changes at the zinc-binding pocket that might account for the dominant negative effects of R248Q. Furthermore, in our simulations, the S6/S7 turn was more frequently solvent exposed in R248Q, suggesting that there is a greater tendency of R248Q to partially unfold and possibly lead to an increased aggregation propensity. Finally, based on the observations made in our simulations, we propose strategies for the rescue of R248Q mutants. [ABSTRACT FROM AUTHOR]

eIF4E is frequently over-expressed in different cancers and causes increased translation of oncogenic proteins via deregulated cap-dependent translation. Inhibitors of the eIF4E:eIF4G interactions represent an approach that would normalize cap-dependent translation. Stapled peptides represent an emerging class of therapeutics that can target protein: protein interactions. We present here molecular dynamics simulations for a set of rationally designed stapled peptides in solution and in complex with eIF4E, supported with biophysical and crystallographic data. Clustering of the simulated structures revealed the favoured conformational states of the stapled peptides in their bound or free forms in solution. Identifying these populations has allowed us to design peptides with improved affinities by introducing mutations into the peptide sequence to alter their conformational distributions. These studies emphasise the effects that engineered mutations have on the conformations of free and bound peptides, and illustrate that both states must be considered in efforts to attain high affinity binding. [ABSTRACT FROM AUTHOR]

Atomisticsimulations of a set of stapled alpha helical peptides derived from the BH3 helixof MCL-1 (Stewart et al. (2010) Nat Chem Biol 6: 595-601) complexed to a fragment (residues 172-320) of MCL-1 revealed that the highest affinity is achieved when the staples engage the surface of MCL-1 as has also been demonstrated for p53-MDM2 (Joseph et al. (2010) Cell Cycle 9: 4560-4568; Baeket al. (2012) J Am Chem Soc 134: 103-106). Affinity is also modulated by the abilityof the staples to pre-organize the peptides as helices. Molecular dynamics simulations of these stapled BH3 peptides were carried out followed by determination of the energies of interactions using MM/GBSA methods. These show that the location of the staple is a key determinant of a good binding stapled peptide from a bad binder. The good binder derives binding affinity from interactions between the hydrophobic staple and a hydrophobic patch on MCL-1. The position of the staple was varied, guiding the design of new stapled peptides with higher affinities. [ABSTRACT FROM AUTHOR]

Proteins of the major histocompatibility complex (MHC) class I, or human leukocyte antigen (HLA) in humans interact with endogenous peptides and present them to T cell receptors (TCR), which in turn tune the immune system to recognize and discriminate between self and foreign (non-self) peptides. Of especial importance are peptides derived from tumor-associated antigens. T cells recognizing these peptides are found in cancer patients, but not in cancer-free individuals. What stimulates this recognition, which is vital for the success of checkpoint based therapy? A peptide derived from the protein p53 (residues 161–169 or p161) was reported to show this behavior. T cells recognizing this unmodified peptide could be further stimulated in vitro to create effective cancer killing CTLs (cytotoxic T lymphocytes). We hypothesize that the underlying difference may arise from post-translational glycosylation of p161 in normal individuals, likely masking it against recognition by TCR. Defects in glycosylation in cancer cells may allow the presentation of the native peptide. We investigate the structural consequences of such peptide glycosylation by investigating the associated structural dynamics. [ABSTRACT FROM AUTHOR]

S100B(ββ) proteins are a family of multifunctional proteins that are present in several tissues and regulate a wide variety of cellular processes. Their altered expression levels have been associated with several human diseases, such as cancer, inflammatory disorders and neurodegenerative conditions, and hence are of interest as a therapeutic target and a biomarker. Small molecule inhibitors of S100B(ββ) have achieved limited success. Guided by the wealth of available experimental structures of S100B(ββ) in complex with diverse peptides from various protein interacting partners, we combine comparative structural analysis and molecular dynamics simulations to design a series of peptides and their analogues (stapled) as S100B(ββ) binders. The stapled peptides were subject to in silico mutagenesis experiments, resulting in optimized analogues that are predicted to bind to S100B(ββ) with high affinity, and were also modified with imaging agents to serve as diagnostic tools. These stapled peptides can serve as theranostics, which can be used to not only diagnose the levels of S100B(ββ) but also to disrupt the interactions of S100B(ββ) with partner proteins which drive disease progression, thus serving as novel therapeutics. [ABSTRACT FROM AUTHOR]

Intrinsically disordered regions (IDRs) in proteins can regulate their activity by facilitating protein-protein interactions (PPIs) as exemplified in the recruitment of the eukaryotic translation initiation factor 4E (eIF4E) protein by the protein eIF4G. Deregulation of this PPI module is central to a broad spectrum of cancer related malignancies and its targeted inhibition through bioactive peptides is a promising strategy for therapeutic intervention. We employed molecular dynamics simulations coupled with biophysical assays to rationally develop peptide derivatives from the intrinsically disordered eIF4G scaffold by incorporating non-natural amino acids that facilitates disorder-to-order transition. The conformational heterogeneity of these peptides and the degree of structural reorganization required to adopt the optimum mode of interaction with eIF4E underscores their differential binding affinities. The presence of a pre-structured local helical element in the ensemble of structures was instrumental in the efficient docking of the peptides on to the protein surface. The formation of Y4: P38 hydrogen-bond interaction between the peptide and eIF4E is a rate limiting event in the efficient recognition of the protein since it occurs through the disordered region of the peptide. These insights were exploited to further design features into the peptide to propagate bound-state conformations in solution which resulted in the generation of a potent eIF4E binder. The study illustrates the molecular basis of eIF4E recognition by a disordered epitope from eIF4G and its modulation to generate peptides that can potentially attenuate translation initiation in oncology. Unlabelled Image • Intrinsically disordered region is involved in eIF4G: eIF4E protein interaction. • Order inducing Aib peptides were derived from the disordered epitope of eIF4G. • Pre-structured helical element in the peptides was critical for recognizing eIF4E. • Propagation of bound state conformation in solution resulted in potent eIF4E binder. • Developed a potential lead compound to attenuate translation initiation in oncology. [ABSTRACT FROM AUTHOR]

Nonsense-mediated mRNA decay (NMD) is a quality control mechanism that recognizes post-transcriptionally abnormal transcripts and mediates their degradation. The master regulator of NMD is UPF1, an enzyme with intrinsic ATPase and helicase activities. The cancer genomic sequencing data has identified frequently mutated residues in the CH-domain and ATP-binding site of UPF1. In silico screening of UPF1 stability change as a function over 41 cancer mutations has identified five variants with significant effects: K164R, R253W, T499M, E637K, and E833K. To explore the effects of these mutations on the associated energy landscape of UPF1, molecular dynamics simulations (MDS) were performed. MDS identified stable H-bonds between residues S152, S203, S205, Q230/R703, and UPF2/AMPPNP, and suggest that phosphorylation of Serine residues may control UPF1-UPF2 binding. Moreover, the alleles K164R and R253W in the CH-domain improved UPF1-UPF2 binding. In addition, E637K and E833K alleles exhibited improved UPF1-AMPPNP binding compared to the T499M variant; the lower binding is predicted from hindrance caused by the side-chain of T499M to the docking of the tri-phosphate moiety (AMPPNP) into the substrate site. The dynamics of wild-type/mutant systems highlights the flexible nature of the ATP-binding region in UPF1. These insights can facilitate the development of drug discovery strategies for manipulating NMD signaling in cell systems using chemical tools. [ABSTRACT FROM AUTHOR]

The transcription factor p53 regulates cellular integrity in response to stress. p53 is mutated in more than half of cancerous cells, with a majority of the mutations localized to the DNA binding domain (DBD). In order to map the structural and dynamical features of the DBD, we carried out multiple copy molecular dynamics simulations (totaling 0.8 μs). Simulations show the loop 1 to be the most dynamic element among the DNA-contacting loops (loops 1-3). Loop 1 occupies two major conformational states: extended and recessed; the former but not the latter displays correlations in atomic fluctuations with those of loop 2 (~24 Å apart). Since loop 1 binds to the major groove whereas loop 2 binds to the minor groove of DNA, our results begin to provide some insight into the possible mechanism underpinning the cooperative nature of DBD binding to DNA. We propose (1) a novel mechanism underlying the dynamics of loop 1 and the possible tread-milling of p53 on DNA and (2) possible mutations on loop 1 residues to restore the transcriptional activity of an oncogenic mutation at a distant site. [ABSTRACT FROM AUTHOR]

The Bcl2 family of proteins is capable of switching the apoptotic machinery by directly controlling the release of apoptotic factors from the mitochondrial outer membrane. They have ‘pro’ and ‘anti’-apoptotic subgroups of proteins which antagonize each other’s function; however a detailed atomistic understanding of their mechanisms based on the dynamical events, particularly in the membrane, is lacking. Using molecular dynamics simulations totaling 1.6µs we outline the major differences between the conformational dynamics in water and in membrane. Using implicit models of solvent and membrane, the simulated results reveal a picture that is in agreement with the ‘hit-and run’ concept which states that BH3-only peptides displace the tail (which acts as a pseudo substrate of the protein itself) from its binding pocket; this helps the membrane association of the protein after which the BH3 peptide becomes free. From simulations, Bcl-xL appears to be auto-inhibited by its C-terminal tail that embeds into and covers the hydrophobic binding pocket. However the tail is unable to energetically compete with BH3-peptides in water. In contrast, in the membrane, neither the tail nor the BH3-peptides are stable in the binding pocket and appear to be easily dissociated off as the pocket expands in response to the hydrophobic environment. This renders the binding pocket large and open, thus receptive to interactions with other protein partners. Principal components of the motions are dramatically different in the aqueous and in the membrane environments and provide clues regarding the conformational transitions that Bcl-xL undergoes in the membrane, in agreement with the biochemical data. [ABSTRACT FROM AUTHOR]

[Display omitted] Carbon nanotubes (CNTs) represent a novel platform for cellular delivery of therapeutic peptides. Chemically-functionalized CNTs may enhance peptide uptake by improving their membrane targeting properties. Using coarse-grained (CG) molecular dynamics (MD) simulations, we investigate membrane interactions of a peptide conjugated to pristine and chemically-modified CNTs. As proof of principle, we focus on their interactions with PM2, an amphipathic stapled peptide that inhibits the E3 ubiquitin ligase HDM2 from negatively regulating the p53 tumor suppressor. CNT interaction with both simple planar lipid bilayers as well as spherical lipid vesicles was studied, the latter as a surrogate for curved cellular membranes. Membrane permeation was rapid and spontaneous for both pristine and oxidized CNTs when unconjugated. This was slowed upon addition of a noncovalently attached peptide surface "sheath", which may be an effective way to slow CNT entry and avert membrane rupture. The CNT conjugates were observed to "desheath" their peptide layer at the bilayer interface upon insertion, leaving their cargo behind in the outer leaflet. This suggests that a synergy may exist to optimize CNT safety whilst enhancing the delivery efficiency of "hitchhiking" therapeutic molecules. [ABSTRACT FROM AUTHOR]

Naturally occurring linear polyamines are known to enable bacteria to be resistant to cationic membrane active peptides. To understand this protective mechanism, molecular dynamics simulations are employed to probe their effect on a model bacterial outer membrane. Being protonated at physiological pH, the amine groups of the polyamine engage in favorable electrostatic interactions with the negatively charged phosphate groups of the membrane. Additionally, the amine groups form large number of hydrogen bonds with the phosphate groups. At high concentrations, these hydrogen bonds and the electrostatic network can non-covalently crosslink the lipid A molecules, resulting in stabilization of the outer membrane against membrane active antibiotics such as colistin and polymyxin B. Moreover, large polyamine molecules (e.g., spermidine) have a stronger stabilization effect than small polyamine molecules (e.g., ethylene diamine). The atomistic insights provide useful guidance for the design of next generation membrane active amine-rich antibiotics, especially to tackle the growing threat of multi-drug resistance of Gram negative bacteria. Unlabelled Image • Polymyxin B disrupt salt bridges of Gram negative outer membrane • Linear polyamines can protect the outer membrane against the action of polymyxin B. • MD simulations revealed the protective mechanism of polyamine through electrostatic networks. • Larger polyamines display more significant protective effect. [ABSTRACT FROM AUTHOR]

Summary Ligand binding pockets in proteins contain water molecules, which play important roles in modulating protein-ligand interactions. Available crystallographic data for the 5′ mRNA cap-binding pocket of the translation initiation factor protein eIF4E shows several structurally conserved waters, which also persist in molecular dynamics simulations. These waters engage an intricate hydrogen-bond network between the cap and protein. Two crystallographic waters in the cleft of the pocket show a high degree of conservation and bridge two residues, which are part of an evolutionarily conserved scaffold. This appears to be a preformed recognition module for the cap with the two structural waters facilitating an efficient interaction. This is also recapitulated in a new crystal structure of the apo protein. These findings open new windows for the design and screening of compounds targeting eIF4E. [ABSTRACT FROM AUTHOR]

Targeting protein–protein interacting sites for potential therapeutic applications is a challenge in the development of inhibitors, and this becomes more difficult when these interfaces are relatively planar, as in the eukaryotic translation initiation factor 4E (eIF4E) protein. eIF4E is an oncogene that is overexpressed in numerous forms of cancer, making it a prime target as a therapeutic molecule. We report here the presence of a cryptic pocket at the protein-binding interface of eIF4E, which opens transiently during molecular dynamics simulations of the protein in solvent water and is observed to be stable when solvent water is mixed with benzene molecules. This pocket can also be seen in the ensemble of structures available from the solution-state conformations of eIF4E. The accessibility of the pocket is gated by the side-chain transitions of an evolutionarily conserved tryptophan residue. It is found to be feasible for accommodating clusters of benzene molecules, which signify the plasticity and ligandability of the pocket. We also observe that the newly formed cavity provides a favorable binding environment for interaction of a well-recognized small molecule inhibitor of eIF4E. The occurrence of this transiently accessible cavity highlights the existence of a more pronounced binding groove in a region that has traditionally been considered to be planar. Together, the data suggest that an alternate binding cavity exists on eIF4E and could be exploited for the rational design and development of a new class of lead compounds against the protein. [ABSTRACT FROM AUTHOR]

Mitogen-activated protein kinases-interacting kinase 1 and 2 (Mnk1/2) activate the oncogene eukaryotic initiation factor 4E (eIF4E) by phosphorylation. High level of phosphorylated eIF4E is associated with various types of cancers. Inhibition of Mnk prevents eIF4E phosphorylation, making them potential therapeutic targets for cancer. Recently, we have designed and synthesized a series of novel imidazopyridine and imidazopyrazine derivatives that inhibit Mnk1/2 kinases with a potency in the nanomolar to micromolar range. In the current work we model the inhibition of Mnk kinase activity by these inhibitors using various computational approaches. Combining homology modeling, docking, molecular dynamics simulations, and free energy calculations, we find that all compounds bind similarly to the active sites of both kinases with their imidazopyridine and imidazopyrazine cores anchored to the hinge regions of the kinases through hydrogen bonds. In addition, hydrogen bond interactions between the inhibitors and the catalytic Lys78 (Mnk1), Lys113 (Mnk2) and Ser131 (Mnk1), Ser166 (Mnk2) appear to be important for the potency and stability of the bound conformations of the inhibitors. The computed binding free energies (ΔGPred) of these inhibitors are in accord with experimental bioactivity data (pIC50) with correlation coefficients (r²) of 0.70 and 0.68 for Mnk1 and Mnk2 respectively. van der Waals energies and entropic effects appear to dominate the binding free energy (ΔGPred) for each Mnk–inhibitor complex studied. The models suggest that the activities of these small molecule inhibitors arise from interactions with multiple residues in the active sites, particularly with the hydrophobic residues. [ABSTRACT FROM AUTHOR]

The mechanism of coupling of ion pumping in the membrane-bound AO sector with ATP synthesis in the A3B3 headpiece of the A1 sector in the A1AO ATP synthase is a puzzle. Previously, crosstalk between the stalk and nucleotide-binding subunits FMm and BMm of the Methanosarcina mazei Gö1 A-ATP synthase has been observed by nucleotide-dependent cross-link formation of both subunits inside the enzyme. The recently determined NMR solution structure of FMm depicts the protein as a two-domain structure, with a well-folded N-terminus having 78 residues and a flexible C-terminal part (residues 79-101), proposed to become structured after binding to its partner, BMm. Here, we detail the crucial interactions between subunits BMm and FMm by determining the NMR structure of the very C-terminus of FMm, consisting of 20 residues and hereafter termed FMm(81-101), and performing molecular dynamics simulations on the resulting structure. These data demonstrate that the flexibility of the C-terminus enables FMm to switch between an elongated and retracted state. Docking and MD in conjunction with previously conducted and published NMR results, biochemical cross-linking, and fluorescence spectroscopy data were used to reconstruct a model of a BMm-FMm assembly. The model of the BMm-FMm complex shows the detailed interactions of helices 1 and 2 of the C-terminal domain of BMm with the C-terminal residues of FMm. Movements of both helices of BMm accommodate the incoming C-terminus of FMm and connect the events of ion pumping and nucleotide binding in the A1AO ATP synthase. [ABSTRACT FROM AUTHOR]

Summary Despite being initially identified in the blood filtrate, LEKTI is a 15-domain Kazal-type inhibitor mostly known in the regulation of skin desquamation. In the current study, screening of serine proteases in blood coagulation cascade showed that LEKTI domain 4 has inhibitory activity toward only FXIa, whereas LEKTI domain 6 inhibits both FXIa and FXaB (bovine FXa). Nuclear magnetic resonance structural and dynamic experiments plus molecular dynamics simulation revealed that LEKTI domain 4 has enhanced backbone flexibility at the reactive-site loop. A model of the LEKTI-protease complex revealed that FXaB has a narrower S4 pocket compared with FXIa and hence prefers only small side-chain residues at the P4 position, such as Ala in LEKTI domain 6. Mutational studies combined with a molecular complex model suggest that both a more flexible reactive-site loop and a bulky residue at the P4 position make LEKTI domain 4 a weaker but highly selective inhibitor of FXIa. [ABSTRACT FROM AUTHOR]