The Pattern of GH Action in the Mouse Brain.

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 IP 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 because 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 (C60H28) 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.

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.

Acid Rain and Flue Gas: Quantum Chemical Hydrolysis of NO2.

Despite decades of efforts, much is still unknown about the hydrolysis of nitrogen dioxide (NO2 ), 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 NO2 . 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.

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.

Local complete active space second-order perturbation theory using pair natural orbitals (PNO-CASPT2).

We present a CASPT2 method which exploits local approximations to achieve linear scaling of the computational effort with the molecular size, provided the active space is small and local. The inactive orbitals are localized, and the virtual space for each electron pair is spanned by a domain of pair-natural orbitals (PNOs). The configuration space is internally contracted, and the PNOs are defined for uniquely defined orthogonal pairs. Distant pair energies are obtained by multipole approximations, so that the number of configurations that are explicitly treated in the CASPT2 scales linearly with molecular size (assuming a constant active space). The PNOs are generated using approximate amplitudes obtained in a pair-specific semi-canonical basis of projected atomic orbitals (PAOs). The evaluation and transformation of the two-electron integrals use the same parallel local density fitting techniques as recently described for linear-scaling PNO-LMP2 (local second-order Møller-Plesset perturbation theory). The implementation of the amplitude equations, which are solved iteratively, employs the local integrated tensor framework. The efficiency and accuracy of the method are tested for excitation energies and correlation energies. It is demonstrated that the errors introduced by the local approximations are very small. They can be well controlled by few parameters for the distant pair approximation, initial PAO domains, and the PNO domains.

Analysis of association of clinical aspects and IL1B tagSNPs with severe preeclampsia.

OBJECTIVE: This study investigates the association between IL1B genotypes using a tag SNP (single polymorphism) approach, maternal and environmental factors in Brazilian women with severe preeclampsia. METHODS: A case-control study with a total of 456 patients (169 preeclamptic women and 287 controls) was conducted in the two reference maternity hospitals of Sergipe state, Northeast Brazil. A questionnaire was administered and DNA was extracted to genotype the population for four tag SNPs of the IL1Beta: rs 1143643, rs 1143633, rs 1143634 and rs 1143630. Haplotype association analysis and p-values were calculated using the THESIAS test. Odds ratio (OR) estimation, confidence interval (CI) and multivariate logistic regression were performed. RESULTS: High pregestational body mass index (pre-BMI), first gestation, cesarean section, more than six medical visits, low level of consciousness on admission and TC and TT genotype in rs1143630 of IL1Beta showed association with the preeclamptic group in univariate analysis. After multivariate logistic regression pre-BMI, first gestation and low level of consciousness on admission remained associated. CONCLUSION: We identified an association between clinical variables and preeclampsia. Univariate analysis suggested that inflammatory process-related genes, such as IL1B, may be involved and should be targeted in further studies. The identification of the genetic background involved in preeclampsia host response modulation is mandatory in order to understand the preeclampsia process.

Similarities and differences between WHO criteria and two other approaches for maternal near miss diagnosis.

OBJECTIVES: To evaluate the similarities, differences and diagnostic aspects between World Health Organization (WHO) criteria and two other maternal near miss (MNM) diagnostic tools. METHODS: A cross-sectional study was conducted from June 2011 to May 2012 in two reference maternity hospitals in Aracaju, Brazil. Prospective case identification and data collection were performed and patients were classified as an MNM case according to WHO, Waterstone and literature-based criteria. The diagnostic properties and concordance of literature-based and Waterstone criteria were calculated using WHO criteria as standard. RESULTS: Of a total of 20 435 patients, 19 239 women did not have potentially life-threatening conditions, there were 17 maternal deaths, and 77 MNM cases based on the WHO criteria. Waterstone and literature-based criteria identified 404 and 959 MNM cases, respectively, most of them related to hypertensive disorders and haemorrhage. The sensitivity, specificity and accuracy in diagnosing MNM cases using Waterstone and literature-based criteria were above 90%, but Waterstone sensitivity was 48.1%. The similarities between the Waterstone and literature-based criteria were very weak compared to WHO criteria, with a positive percentage concordance below 9%. CONCLUSIONS: Although using WHO guidelines to detect MNM cases can be difficult when implemented in low-resource settings, the results from this study reinforce the importance of this tool in detecting the truly severe cases. Waterstone and literature-based criteria are not suitable for identifying indubitable MNM. However, they remain useful as a preliminary step to select potentially severe cases, mainly those related to hypertension and haemorrhage.

The prevalence of severe maternal morbidity and near miss and associated factors in Sergipe, Northeast Brazil.

BACKGROUND: The investigation of severe maternal morbidity (SAMM) and maternal near miss (NM) and associated risk factors is important for the global reduction of maternal mortality. This study investigated the prevalence of SAMM and NM cases and the associated risk factors in two reference maternity hospitals in a capital city in Northeast-Brazil. METHODS: A cross-sectional study with a nested case-control component was conducted from June-2011 to May-2012. Case identification was prospective and data collection was performed according to WHO criteria and definitions. Odds ratio with confidence intervals and multivariate analysis were used whenever possible. RESULTS: There were 16,243 deliveries, 1,102 SAMM cases, 77 NM cases and 17 maternal deaths. The maternal NM outcome ratio was 5.8 cases/1,000 live births (LB); the total prevalence of SAMM + NM was 72.6 cases/1,000 LB, the maternal near miss: mortality ratio was 4.5cases/1 maternal death (18% of mortality index). Management-based criteria were the most common events for NM (87.1%) and hypertensive disorders for SAMM (67.5%). Higher age, previous abortion and caesarean delivery, the non-adhesion to antenatal care, current caesarean delivery and bad perinatal results were associated with SAMM/NM. In the multivariate analysis, patient's status, previous caesarian and abortion and level of consciousness were significant when analyzed together. CONCLUSIONS: SAMM and NM situations were prevalent in the studied population and some risk factors seem to be associated with the event, particularly previous gestational antecedents. Protocols based on SAMM/NM situations can save lives and decrease maternal mortality.

Methods for the analysis of complex fluorescence decays: sum of Becquerel functions versus sum of exponentials.

Ensemble fluorescence decays are usually analyzed with a sum of exponentials. However, broad continuous distributions of lifetimes, either unimodal or multimodal, occur in many situations. A simple and flexible fitting function for these cases that encompasses the exponential is the Becquerel function. In this work, the applicability of the Becquerel function for the analysis of complex decays of several kinds is tested. For this purpose, decays of mixtures of four different fluorescence standards (binary, ternary and quaternary mixtures) are measured and analyzed. For binary and ternary mixtures, the expected sum of narrow distributions is well recovered from the Becquerel functions analysis, if the correct number of components is used. For ternary mixtures, however, satisfactory fits are also obtained with a number of Becquerel functions smaller than the true number of fluorophores in the mixture, at the expense of broadening the lifetime distributions of the fictitious components. The quaternary mixture studied is well fitted with both a sum of three exponentials and a sum of two Becquerel functions, showing the inevitable loss of information when the number of components is large. Decays of a fluorophore in a heterogeneous environment, known to be represented by unimodal and broad continuous distributions (as previously obtained by the maximum entropy method), are also measured and analyzed. It is concluded that these distributions can be recovered by the Becquerel function method with an accuracy similar to that of the much more complex maximum entropy method. It is also shown that the polar (or phasor) plot is not always helpful for ascertaining the degree (and kind) of complexity of a fluorescence decay.