Search

Atomic-scale simulations have progressed tremendously over the past decade, largely thanks to the availability of machine-learning interatomic potentials. These potentials combine the accuracy of electronic structure calculations with the ability to reach extensive length and time scales. The i-PI package facilitates integrating the latest developments in this field with advanced modeling techniques thanks to a modular software architecture based on inter-process communication through a socket interface. The choice of Python for implementation facilitates rapid prototyping but can add computational overhead. In this new release, we carefully benchmarked and optimized i-PI for several common simulation scenarios, making such overhead negligible when i-PI is used to model systems up to tens of thousands of atoms using widely adopted machine learning interatomic potentials, such as Behler-Parinello, DeePMD, and MACE neural networks. We also present the implementation of several new features, including an efficient algorithm to model bosonic and fermionic exchange, a framework for uncertainty quantification to be used in conjunction with machine-learning potentials, a communication infrastructure that allows for deeper integration with electronic-driven simulations, and an approach to simulate coupled photon-nuclear dynamics in optical or plasmonic cavities., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

β-Glucocerebrosidase (GBA/GCase) mutations leading to misfolded protein cause Gaucher's disease and are a major genetic risk factor for Parkinson's disease and dementia with Lewy bodies. The identification of small molecule pharmacological chaperones that can stabilize the misfolded protein and increase delivery of degradation-prone mutant GCase to the lysosome is a strategy under active investigation. Here, we describe the first use of fragment-based drug discovery (FBDD) to identify pharmacological chaperones of GCase. The fragment hits were identified by using X-ray crystallography and biophysical techniques. This work led to the discovery of a series of compounds that bind GCase with nM potency and positively modulate GCase activity in cells.

Lithium ortho-thiophosphate (Li 3 PS 4 ) has emerged as a promising candidate for solid-state electrolyte batteries, thanks to its highly conductive phases, cheap components, and large electrochemical stability range. Nonetheless, the microscopic mechanisms of Li-ion transport in Li 3 PS 4 are far from being fully understood, the role of PS 4 dynamics in charge transport still being controversial. In this work, we build machine learning potentials targeting state-of-the-art DFT references (PBEsol, r 2 SCAN, and PBE0) to tackle this problem in all known phases of Li 3 PS 4 (α, β, and γ), for large system sizes and time scales. We discuss the physical origin of the observed superionic behavior of Li 3 PS 4 : the activation of PS 4 flipping drives a structural transition to a highly conductive phase, characterized by an increase in Li-site availability and by a drastic reduction in the activation energy of Li-ion diffusion. We also rule out any paddle-wheel effects of PS 4 tetrahedra in the superionic phases-previously claimed to enhance Li-ion diffusion-due to the orders-of-magnitude difference between the rate of PS 4 flips and Li-ion hops at all temperatures below melting. We finally elucidate the role of interionic dynamical correlations in charge transport, by highlighting the failure of the Nernst-Einstein approximation to estimate the electrical conductivity. Our results show a strong dependence on the target DFT reference, with PBE0 yielding the best quantitative agreement with experimental measurements not only for the electronic band gap but also for the electrical conductivity of β- and α-Li 3 PS 4 ., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

Fragment-based ligand discovery was successfully applied to histone deacetylase HDAC2. In addition to the anticipated hydroxamic acid- and benzamide-based fragment screening hits, a low affinity (∼1 mM) α-amino-amide zinc binding fragment was identified, as well as fragments binding to other regions of the catalytic site. This alternative zinc-binding fragment was further optimized, guided by the structural information from protein-ligand complex X-ray structures, into a sub-μM, brain penetrant, HDAC2 inhibitor ( 17 ) capable of modulating histone acetylation levels in vivo ., Competing Interests: The authors declare no competing financial interest., (© 2022 American Chemical Society.)

Ubiquitin specific protease USP15 is a deubiquitinating enzyme reported to regulate several biological and cellular processes, including TGF-β signaling, regulation of immune response, neuro-inflammation and mRNA splicing. Here we study the USP15 D1D2 catalytic domain and present the crystal structure in its catalytically-competent conformation. We compare this apo-structure to a previous misaligned state in the same crystal lattice. In both structures, mitoxantrone, an FDA approved antineoplastic drug and a weak inhibitor of USP15 is bound, indicating that it is not responsible for inducing a switch in the conformation of active site cysteine in the USP15 D1D2 structure. Instead, mitoxantrone contributes to crystal packing, by forming a stack of 12 mitoxantrone molecules. We believe this reflects how mitoxantrone can be responsible for e.g. nuclear condensate partitioning. We conclude that USP15 can switch between active and inactive states in the absence of ubiquitin, and that this is independent of mitoxantrone binding. These insights can be important for future drug discovery targeting USP15., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

We present a novel crystallographic screening methodology (MiniFrags) that employs high-concentration aqueous soaks with a chemically diverse and ultra-low-molecular-weight library (heavy atom count 5-7) to identify ligand-binding hot and warm spots on proteins. We propose that MiniFrag screening represents a highly effective method for guiding optimisation of fragment-derived lead compounds or chemical tools and that the high screening hit rates reflect enhanced sampling of chemical space., (Copyright © 2019 Astex Pharmaceuticals. Published by Elsevier Ltd.. All rights reserved.)

Historically chloroquine was used to treat the most deadly form of malaria, caused by the parasite Plasmodium falciparum. The selective pressure of chloroquine therapy led to the rapid emergence of chloroquine resistant parasites. Resistance has been attributed to the Plasmodium falciparum Chloroquine Resistance Transporter (PfCRT), an integral membrane protein of unknown structure. A PfCRT structure would provide new insights into how the protein confers chloroquine resistance and thereby also yield novel opportunities for developing anti-malarial therapies. Although PfCRT is an attractive target for characterisation and structure determination, very little work has been published on its expression and purification. Here we present a medium throughput protocol, employing Sf9 insect cells, for testing the expression, stability and purification yield of rationally designed PfCRT mutant constructs and constructs of a PfCRT orthologue from Neospora caninum (NcCRT). We have identified a conserved cysteine residue in PfCRT that results in elevated protein stability when mutated. Combining this mutation with the insertion of T4-lysozyme into a specific surface loop further augments PfCRT protein yield and thermostability. Screening also identified an NcCRT construct with an elevated purification yield. Furthermore it was possible to purify both PfCRT and NcCRT constructs at milligram-scales, with high purities and with size exclusion chromatography profiles that were consistent with monodispersed, homogeneous protein., (Crown Copyright © 2017. Published by Elsevier Inc. All rights reserved.)

Background: Homeless adults frequently use emergency departments (EDs), yet previous studies investigating ED utilization by the older segment received little attention. This study sought to characterize older homeless adults who utilized local urban EDs., Methods: ED encounters at three hospitals in Hamilton (Ont.) were analyzed, and demographic and clinical characteristics of the older homeless (age > 50) vs. younger counterparts (age ≤ 50) were compared during a 24-month period., Results: Of all adults, 1,330 were homeless, of whom 66% were above age 50. Older homeless adults sought less acute care within 30 days from an index visit compared with their younger counterparts. Non-acute illnesses constituted only 18% of triaged cases. Older homeless women with access to a primary care physician (PCP) were 3.3 times more likely to return to ED within 30 days, whereas older homeless men (irrespective of PCP access) were less likely to return to ED., Conclusions: Despite high homeless patient acuity, a lesser number of ED visits with increasing age remains concerning because of previously reported high morbidity and mortality rates. Access to primary care may not be enough to reduce ED utilization. Further research is needed to evaluate acute care interventions and their effectiveness in ED, and to identify homeless patients requiring more targeted services., Competing Interests: DISCLOSURES The authors declare that no conflicts of interest exist.

The members of the NSD subfamily of lysine methyl transferases are compelling oncology targets due to the recent characterization of gain-of-function mutations and translocations in several hematological cancers. To date, these proteins have proven intractable to small molecule inhibition. Here, we present initial efforts to identify inhibitors of MMSET (aka NSD2 or WHSC1) using solution phase and crystal structural methods. On the basis of 2D NMR experiments comparing NSD1 and MMSET structural mobility, we designed an MMSET construct with five point mutations in the N-terminal helix of its SET domain for crystallization experiments and elucidated the structure of the mutant MMSET SET domain at 2.1 Å resolution. Both NSD1 and MMSET crystal systems proved resistant to soaking or cocrystallography with inhibitors. However, use of the close homologue SETD2 as a structural surrogate supported the design and characterization of N-alkyl sinefungin derivatives, which showed low micromolar inhibition against both SETD2 and MMSET.

Herein we describe the application of fragment-based drug design to bacterial DNA ligase. X-ray crystallography was used to guide structure-based optimization of a fragment-screening hit to give novel, nanomolar, AMP-competitive inhibitors. The lead compound 13 showed antibacterial activity across a range of pathogens. Data to demonstrate mode of action was provided using a strain of S. aureus, engineered to overexpress DNA ligase.

Here, we describe the identification of a clinical candidate via structure-based optimization of a ligand efficient pyrazole-benzimidazole fragment. Aurora kinases play a key role in the regulation of mitosis and in recent years have become attractive targets for the treatment of cancer. X-ray crystallographic structures were generated using a novel soakable form of Aurora A and were used to drive the optimization toward potent (IC(50) approximately 3 nM) dual Aurora A/Aurora B inhibitors. These compounds inhibited growth and survival of HCT116 cells and produced the polyploid cellular phenotype typically associated with Aurora B kinase inhibition. Optimization of cellular activity and physicochemical properties ultimately led to the identification of compound 16 (AT9283). In addition to Aurora A and Aurora B, compound 16 was also found to inhibit a number of other kinases including JAK2 and Abl (T315I). This compound demonstrated in vivo efficacy in mouse xenograft models and is currently under evaluation in phase I clinical trials.

Fermentation optimization experiments are ideally performed at small scale to reduce time, cost and resource requirements. Currently microwell plates (MWPs) are under investigation for this purpose as the format is ideally suited to automated high-throughput experimentation. In order to translate an optimized small-scale fermentation process to laboratory and pilot scale stirred-tank reactors (STRs) it is necessary to characterize key engineering parameters at both scales given the differences in geometry and the mechanisms of aeration and agitation. In this study oxygen mass transfer coefficients are determined in three MWP formats and in 7.5 L and 75 L STRs. k(L)a values were determined in cell-free media using the dynamic gassing-out technique over a range of agitation conditions. Previously optimized culture conditions at the MWP scale were then scaled up to the larger STR scales on the basis of matched k(L)a values. The accurate reproduction of MWP (3 mL) E. coli BL21 (DE3) culture kinetics at the two larger scales was shown in terms of cell growth, protein expression, and substrate utilization for k(L)a values that provided effective mixing and gas-liquid distribution at each scale. This work suggests that k(L)a provides a useful initial scale-up criterion for MWP culture conditions which enabled a 15,000-fold scale translation in this particular case. This work complements our earlier studies on the application of DoE techniques to MWP fermentation optimization and in so doing provides a generic framework for the generation of large quantities of soluble protein in a rapid and cost-effective manner.

We describe the structure-guided optimization of the molecular fragments 2-amino-3-benzyloxypyridine 1 (IC(50) 1.3 mM) and 3-(2-(4-pyridyl)ethyl)indole 2 (IC(50) 35 microM) identified using X-ray crystallographic screening of p38alpha MAP kinase. Using two separate case studies, the article focuses on the key compounds synthesized, the structure-activity relationships and the binding mode observations made during this optimization process, resulting in two potent lead series that demonstrate significant increases in activity. We describe the process of compound elaboration either through the growing out from fragments into adjacent pockets or through the conjoining of overlapping fragments and demonstrate that we have exploited the mobile conserved activation loop, consisting in part of Asp168-Phe169-Gly170 (DFG), to generate significant improvements in potency and kinase selectivity.

Protein tyrosine phosphatases regulate signal transduction pathways involving tyrosine phosphorylation and have been implicated in the development of cancer, diabetes, rheumatoid arthritis and hypertension. Increasing evidence suggests that the cellular redox state is involved in regulating tyrosine phosphatase activity through the reversible oxidization of the catalytic cysteine to sulphenic acid (Cys-SOH). But how further oxidation to the irreversible sulphinic (Cys-SO2H) and sulphonic (Cys-SO3H) forms is prevented remains unclear. Here we report the crystal structures of the regulatory sulphenic and irreversible sulphinic and sulphonic acids of protein tyrosine phosphatase 1B (PTP1B), an important enzyme in the negative regulation of the insulin receptor and a therapeutic target in type II diabetes and obesity. We also identify a sulphenyl-amide species that is formed through oxidation of its catalytic cysteine. Formation of the sulphenyl-amide causes large changes in the PTP1B active site, which are reversible by reduction with the cellular reducing agent glutathione. The sulphenyl-amide is a protective intermediate in the oxidative inhibition of PTP1B. In addition, it may facilitate reactivation of PTP1B by biological thiols and signal a unique state of the protein.

Fasciclin I is an insect neural cell adhesion molecule consisting of four FAS1 domains, homologs of which are present in many bacterial, plant, and animal proteins. The crystal structure of FAS1 domains 3 and 4 of Drosophila fasciclin I reveals a novel domain fold, consisting of a seven-stranded beta wedge and a number of alpha helices. The two domains are arranged in a linear fashion and interact through a substantial polar interface. Missense mutations in the FAS1 domains of the human protein betaig-h3 cause corneal dystrophies. Many mutations alter highly conserved core residues, but the two most common mutations, affecting Arg-124 and Arg-555, map to exposed alpha-helical regions, suggesting reduced protein solubility as the disease mechanism.

Receptor tyrosine kinases of the Axl family are activated by Gas6, the product of growth arrest-specific gene 6. Gas6-Axl signaling is implicated in cell survival, adhesion, and migration. The receptor-binding site of Gas6 is located within a C-terminal pair of laminin G-like (LG) domains that do not resemble any other receptor tyrosine kinase ligand. We report the crystal structure at 2.2-A resolution of a Gas6 fragment spanning both LG domains (Gas6-LG). The structure reveals a V-shaped arrangement of LG domains strengthened by an interdomain calcium-binding site. LG2 of Gas6-LG contains two unusual features: an alpha-helix cradled by one edge of the LG beta-sandwich and a conspicuous patch of surface-exposed hydrophobic residues. Mutagenesis of some residues in this patch reduces Gas6-LG binding to the extracellular domain of Axl as well as Axl activation in glioblastoma cells, identifying a component of the receptor-binding site of Gas6.

Creatine kinase (CK) catalyses the reversible transfer of the phosphate moiety from phosphocreatine (PCr) to ADP, generating creatine and ATP. The crystal structure of a cytosolic brain-type creatine kinase is reported at 2.3 A. The biological dimer sits on a crystallographic twofold axis. The N-terminal residues of both subunits come very close to the crystallographic twofold at the dimer interface. The electron density observed is consistent with two alternative conformations for the N-termini, as previously found for chicken brain-type creatine kinase.

Endostatin is a fragment of the C-terminal domain NC1 of collagen XVIII that inhibits angiogenesis and tumor growth. We report the characterization of a collagen XV endostatin analogue and its parent NC1 domain, obtained by recombinant expression in mammalian cells. Both NC1 domains contain a trimerization domain, a hinge region that is more sensitive to proteolysis in collagen XVIII and the endostatin domain. Unlike endostatin-XVIII, endostatin-XV does not bind zinc or heparin, which is explained by the crystal structure of endostatin-XV. The collagen XV and XVIII fragments inhibited chorioallantoic membrane angiogenesis induced by basic fibroblast growth factor (FGF-2) or vascular endothelial growth factor (VEGF), but there are striking differences depending on which cytokine is used and whether free endostatins or NC1 domains are applied. The collagen XV and XVIII fragments showed a similar binding repertoire for extracellular matrix proteins. Differences were found in the immunohistological localization in vessel walls and basement membrane zones. Together, these data indentify endostatin-XV as an angiogenesis inhibitor, which differs from endostatin-XVIII in several important functional details., (Copyright 2000 Academic Press.)

Laminin G domain-like (LG) modules of approximately 180-200 residues are found in a number of extracellular and receptor proteins and often are present in tandem arrays. LG modules are implicated in interactions with cellular receptors (integrins, alpha-dystroglycan), sulfated carbohydrates and other extracellular ligands. The recently determined crystal structures of LG modules of the laminin alpha2 chain reveal a compact beta sandwich fold and identify a novel calcium binding site. Binding epitopes for heparin, sulfatides and alpha-dystroglycan have been mapped by site-directed mutagenesis and show considerable overlap. The epitopes are located in surface loops around the calcium site, which in other proteins (agrin, neurexins) are modified by alternative splicing. Efficient ligand binding often requires LG modules to be present in tandem. The close proximity of the N- and C-termini in the LG module, as well as a unique link region between laminin LG3 and LG4, impose certain constraints on the arrangement of LG tandems. Further modifications may be introduced by proteolytic processing of laminin G domains, which is known to occur in the alpha2, alpha3 and alpha4 chains.

The laminins are large heterotrimeric glycoproteins with fundamental roles in basement membrane architecture and function. The C-terminus of the laminin alpha chain contains a tandem of five laminin G-like (LG) domains. We report the 2.0 A crystal structure of the laminin alpha2 LG4-LG5 domain pair, which harbours binding sites for heparin and the cell surface receptor alpha-dystroglycan, and is 41% identical to the laminin alpha1 E3 fragment. LG4 and LG5 are arranged in a V-shaped fashion related by a 110 degrees rotation about an axis passing near the domain termini. An extended N-terminal segment is disulfide bonded to LG5 and stabilizes the domain pair. Two calcium ions, one each in LG4 and LG5, are located 65 A apart at the tips of the domains opposite the polypeptide termini. An extensive basic surface region between the calcium sites is proposed to bind alpha-dystroglycan and heparin. The LG4-LG5 structure was used to construct a model of the laminin LG1-LG5 tandem and interpret missense mutations underlying protein S deficiency.

Laminin G-like (LG) modules in the extracellular matrix glycoproteins laminin, perlecan, and agrin mediate the binding to heparin and the cell surface receptor alpha-dystroglycan (alpha-DG). These interactions are crucial to basement membrane assembly, as well as muscle and nerve cell function. The crystal structure of the laminin alpha 2 chain LG5 module reveals a 14-stranded beta sandwich. A calcium ion is bound to one edge of the sandwich by conserved acidic residues and is surrounded by residues implicated in heparin and alpha-DG binding. A calcium-coordinated sulfate ion is suggested to mimic the binding of anionic oligosaccharides. The structure demonstrates a conserved function of the LG module in calcium-dependent lectin-like alpha-DG binding.