The pattern of growth hormone action in the mouse brain

Growth hormone (GH) acts in numerous organs expressing the GH receptor (GHR), including the brain. However, the mechanisms behind the brain's permeability to GH and how this hormone accesses different brain regions remain unclear. It is well-known that an acute GH administration induces phosphorylation of the signal transducer and activator of transcription 5 (pSTAT5) in the mouse brain. Thus, the pattern of pSTAT5 immunoreactive cells was analyzed at different time points after intraperitoneal or intracerebroventricular GH injections. After a systemic GH injection, the first cells expressing pSTAT5 were those near circumventricular organs, such as arcuate nucleus neurons adjacent to the median eminence. Both systemic and central GH injections induced a medial-to-lateral pattern of pSTAT5 immunoreactivity over time, as GH-responsive cells were initially observed in periventricular areas and were progressively detected in lateral brain structures. Very few choroid plexus cells exhibited GH-induced pSTAT5. Additionally, Ghr mRNA was poorly expressed in the mouse choroid plexus. In contrast, some tanycytes lining the floor of the third ventricle expressed Ghr mRNA and exhibited GH-induced pSTAT5. The transport of radiolabeled GH into the hypothalamus did not differ between wild-type and dwarf Ghr knockout mice, indicating that GH transport into the mouse brain is GHR-independent. Also, single-photon emission computed tomography confirmed that radiolabeled GH rapidly reaches the ventral part of the tuberal hypothalamus. In conclusion, our study provides novel and valuable information about the pattern and mechanisms behind GH transport into the mouse brain.

MISATO: machine learning dataset of protein-ligand complexes for structure-based drug discovery

Large language models have greatly enhanced our ability to understand biology and chemistry, yet robust methods for structure-based drug discovery, quantum chemistry and structural biology are still sparse. Precise biomolecule-ligand interaction datasets are urgently needed for large language models. To address this, we present MISATO, a dataset that combines quantum mechanical properties of small molecules and associated molecular dynamics simulations of ~20,000 experimental protein-ligand complexes with extensive validation of experimental data. Starting from the existing experimental structures, semi-empirical quantum mechanics was used to systematically refine these structures. A large collection of molecular dynamics traces of protein-ligand complexes in explicit water is included, accumulating over 170 μs. We give examples of machine learning (ML) baseline models proving an improvement of accuracy by employing our data. An easy entry point for ML experts is provided to enable the next generation of drug discovery artificial intelligence models.

SARS-CoV-2 Mpro oligomerization as a potential target for therapy

The main protease (Mpro) of SARS-CoV-2 is critical in the virus's replication cycle, facilitating the maturation of polyproteins into functional units. Due to its conservation across taxa, Mpro is a promising target for broad-spectrum antiviral drugs. Targeting Mpro with small molecule inhibitors, such as nirmatrelvir combined with ritonavir (Paxlovid™), which the FDA has approved for post-exposure treatment and prophylaxis, can effectively interrupt the replication process of the virus. A key aspect of Mpro's function is its ability to form a functional dimer. However, the mechanics of dimerization and its influence on proteolytic activity remain less understood. In this study, we utilized biochemical, structural, and molecular modelling approaches to explore Mpro dimerization. We evaluated critical residues, specifically Arg4 and Arg298, that are essential for dimerization. Our results show that changes in the oligomerization state of Mpro directly affect its enzymatic activity and dimerization propensity. We discovered a synergistic relationship influencing dimer formation, involving both intra- and intermolecular interactions. These findings highlight the potential for developing allosteric inhibitors targeting Mpro, offering promising new directions for therapeutic strategies.

Correction: Miranda et al. Anion Binding by Fluorescent Ureido-Hexahomotrioxacalix[3]arene Receptors: An NMR, Absorption and Emission Spectroscopic Study. Molecules 2022, 27, 3247

In the original publication [...].

When catchers meet - a computational study on the dimerization of the Buckycatcher

We study the dimerization of the buckycatcher in gas phase and in toluene. We created an extensive library of 36 different complexes, which were characterized at semi-empirical and DFT levels. Semi-empirical geometries and dimerization energies compare well against reference data or Density Functional Theory calculations we performed. Born-Oppenheimer molecular dynamics was used to understand what happens when two molecules of the buckycatcher meet, allowing us to infer on the lack of kinetic barriers when dimers form. Thermodynamically, it is possible that room temperature solutions contain dimerized buckycatcher. Using a very simple exchange model, it is shown, however, that dimerization cannot compete thermodynamically against complexation with fullerenes, which accounts for experimental observations.

Silmitasertib (CX-4945), a Clinically Used CK2-Kinase Inhibitor with Additional Effects on GSK3β and DYRK1A Kinases: A Structural Perspective

A clinical casein kinase 2 inhibitor, CX-4945 (silmitasertib), shows significant affinity toward the DYRK1A and GSK3β kinases, involved in down syndrome phenotypes, Alzheimer's disease, circadian clock regulation, and diabetes. This off-target activity offers an opportunity for studying the effect of the DYRK1A/GSK3β kinase system in disease biology and possible line extension. Motivated by the dual inhibition of these kinases, we solved and analyzed the crystal structures of DYRK1A and GSK3β with CX-4945. We built a quantum-chemistry-based model to rationalize the compound affinity for CK2α, DYRK1A, and GSK3β kinases. Our calculations identified a key element for CK2α's subnanomolar affinity to CX-4945. The methodology is expandable to other kinase selectivity modeling. We show that the inhibitor limits DYRK1A- and GSK3β-mediated cyclin D1 phosphorylation and reduces kinase-mediated NFAT signaling in the cell. Given the CX-4945's clinical and pharmacological profile, this inhibitory activity makes it an interesting candidate with potential for application in additional disease areas.

A Buckycatcher in Solution-A Computational Perspective

In this work, we study the buckycatcher (C₆₀H₂₈) in solution using quantum chemical models. We investigate the conformational equilibria in several media and the effects that molecules of solvent might have in interconversion barriers between the different conformers. These are studied in a hypothetical gas phase, in the dielectric of a solvent, as well as with hybrid solvation. In the latter case, due to a disruption of π-stacking interactions, the transition states are destabilized. We also evaluate the complexation of the buckycatcher with solvent-like molecules. In most cases studied, there should be no adducts formed because the enthalpy driving force cannot overcome entropic penalties.

Acid Rain and Flue Gas: Quantum Chemical Hydrolysis of NO2

Despite decades of efforts, much is still unknown about the hydrolysis of nitrogen dioxide (NO₂ ), a reaction associated with the formation of acid rain. From the experimental point of view, quantitative analyses are hard, and without pH control the products decompose to some reagents. We resort to high-level quantum chemistry to compute Gibbs energies for a network of reactions relevant to the hydrolysis of NO₂ . With COSMO-RS solvation corrections we calculate temperature dependent thermodynamic data in liquid water. Using the computed reaction energies, we determine equilibrium concentrations for a gas-liquid system at controlled pH. For different temperatures and initial concentrations of the different species, we observe that nitrogen dioxide should be fully converted to nitric and nitrous acid. The thermodynamic data in this work can have a potential major impact for several industries with regards to the understanding of atmospheric chemistry and in the reduction of anthropomorphic pollution.

ULYSSES: An Efficient and Easy to Use Semiempirical Library for C+

We introduce ULYSSES, a user-friendly and robust C++ library for semiempirical quantum chemical calculations. In its current version, ULYSSES is equipped with a large set of different semiempirical models, most of which are based on the Neglect of Diatomic Differential Overlap (NDDO) approximation. Empirical corrections for dispersion and hydrogen bonding are available for most methods, so that higher quality is achieved in the calculation of energies of nonbonded complexes. The library is furthermore equipped with geometry optimization, as well as modules for calculating molecular properties of general interest. Ideal gas thermodynamics is available and allows single structure as well as conformer (multistructure) averaged properties to be calculated. We offer the possibility to use several vibrational partition functions according to the nature of interactions being studied: for covalent systems, the traditional harmonic oscillator approximation is available; for nonbonded complexes, we systematically extended the partition function proposed by Grimme for all thermodynamic functions. The library is also capable of running Born-Oppenheimer molecular dynamics.

How to Catch the Ball: Fullerene Binding to the Corannulene Pincer

The corannulene pincer (also known in the literature as the buckycatcher) is a fascinating system that may encapsulate, among other molecules, the C60 and C70 fullerenes. These complexes are held together by strong π-stacking interactions. Although these are quantum mechanical effects, their description by quantum chemical methods has proved very hard. We used three semi-empirical methods, PM6-D3H4X, PM6-D3H+ and GFN2-xTB, to model the interactions. Binding to fullerenes was extended to all open conformations of the buckycatcher, and with the proper choice of solvation model and partition functions, we obtained Gibbs free energies of binding that deviated by 1.0–1.5 kcal/mol from the experimental data. Adding three-body dispersion to PM6-D3H+ led to even better agreement. These results agree better with the experimental data than calculations using higher-level methods at a significantly lower fraction of the computational cost. Furthermore, the formation of adducts with C60 was studied using dynamical simulations, which helped to build a more complete picture of the behavior of the corannulene pincer with fullerenes. We also investigated the use of exchange-binding models to recover more information on this system in solution. Though the final Gibbs free energies in solution were worsened, gas-phase enthalpies and entropies better mirrored the experimental data.

Anion Binding by Fluorescent Ureido-Hexahomotrioxacalix[3]arene Receptors: An NMR, Absorption and Emission Spectroscopic Study

Fluorescent receptors (4a-4c) based on (thio)ureido-functionalized hexahomotrioxacalix[3]arenes were synthesised and obtained in the partial cone conformation in solution. Naphthyl or pyrenyl fluorogenic units were introduced at the lower rim of the calixarene skeleton via a butyl spacer. The binding of biologically and environmentally relevant anions was studied with NMR, UV-vis absorption, and fluorescence titrations. Fluorescence of the pyrenyl receptor 4c displays both monomer and excimer fluorescence. The thermodynamics of complexation was determined in acetonitrile and was entropy-driven. Computational studies were also performed to bring further insight into the binding process. The data showed that association constants increase with the anion basicity, and AcO-, BzO- and F- were the best bound anions for all receptors. Pyrenylurea 4c is a slightly better receptor than naphthylurea 4a, and both are more efficient than naphthyl thiourea 4b. In addition, ureas 4a and 4c were also tested as ditopic receptors in the recognition of alkylammonium salts.

Human papillomavirus types 16 and 18 in carcinomas of the penis from Brazil

Human papillomavirus (HPV) type-16 DNA sequences were found in 26/53 (49%) and HPV type-18 in 5/53 (9%) of penile cancers. In only one specimen was HPV18 found on its own. HPV16 sequences were integrated into the host cell chromosomes although some monomeric and oligomeric free forms of DNA were detected in a few tissues. HPV18 was, as far as detectable, free and unintegrated in all tissues tested. HPV16 DNA was also detected in 40% of cases of carcinoma of the cervix in 19 women from the same social groups as the males. The viral DNA in the female cases was a mixture of integrated and free forms. The DNA binding protein (ICSP 11/12) of herpes simplex type 2 (HSV2) was not detected in 10 penile cancers tested with a monoclonal antibody.