Effect of ion to ligand ratio on the aqueous to organic relative solubility of a lanthanide-ligand complex.

In the solvent extraction of rare earth elements, mechanistic aspects remain unclear regarding where and how extractant molecules coordinate metal ions and transport them from the aqueous phase into the organic phase. Molecular dynamics simulations were used to examine how unprotonated di(2-ethylhexyl)phosphoric acid (DEHP-) ligands that coordinate the Gd3+ ion can transfer the ion across the water-organic interface. Using the umbrella sampling technique, potential of mean force profiles were constructed to quantify the relative solubility of the Gd3+ ion coordinated to 0-3 DEHP- ligands in either water, 1-octanol, or hexane solvents and at the water-organic interfaces. The simulations show the Gd-DEHP- complexes, at varying Ln-ligand ratios, preferentially solvate on water-organic interfaces. While the Gd(DEHP-)3 complex will diffuse past the aqueous-organic interface into the octanol solvent, it is thermodynamically preferred for the Gd(DEHP-)3 complex to remain in the water-hexane interface when there is no amphiphilic layer of excess ligand.

The influence of model building schemes and molecular dynamics sampling on QM-cluster models: the chorismate mutase case study.

Most QM-cluster models of enzymes are constructed based on X-ray crystal structures, which limits comparison to in vivo structure and mechanism. The active site of chorismate mutase from Bacillus subtilis and the enzymatic transformation of chorismate to prephenate is used as a case study to guide construction of QM-cluster models built first from the X-ray crystal structure, then from molecular dynamics (MD) simulation snapshots. The Residue Interaction Network ResidUe Selector (RINRUS) software toolkit, developed by our group to simplify and automate the construction of QM-cluster models, is expanded to handle MD to QM-cluster model workflows. Several options, some employing novel topological clustering from residue interaction network (RIN) information, are evaluated for generating conformational clustering from MD simulation. RINRUS then generates a statistical thermodynamic framework for QM-cluster modeling of the chorismate mutase mechanism via refining 250 MD frames with density functional theory (DFT). The 250 QM-cluster models sampled provide a mean ΔG‡ of 10.3 ± 2.6 kcal mol-1 compared to the experimental value of 15.4 kcal mol-1 at 25 °C. While the difference between theory and experiment is consequential, the level of theory used is modest and therefore "chemical" accuracy is unexpected. More important are the comparisons made between QM-cluster models designed from the X-ray crystal structure versus those from MD frames. The large variations in kinetic and thermodynamic properties arise from geometric changes in the ensemble of QM-cluster models, rather from the composition of the QM-cluster models or from the active site-solvent interface. The findings open the way for further quantitative and reproducible calibration in the field of computational enzymology using the model construction framework afforded with the RINRUS software toolkit.

Solution Structures of Europium Terpyridyl Complexes with Nitrate and Triflate Counterions in Acetonitrile.

Lanthanide-ligand complexes are key components of technological applications, and their properties depend on their structures in the solution phase, which are challenging to resolve experimentally or computationally. The coordination structure of the Eu3+ ion in different coordination environments in acetonitrile is examined using ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure (EXAFS) spectroscopy. AIMD simulations are conducted for the solvated Eu3+ ion in acetonitrile, both with or without a terpyridyl ligand, and in the presence of either triflate or nitrate counterions. EXAFS spectra are calculated directly from AIMD simulations and then compared to experimentally measured EXAFS spectra. In acetonitrile solution, both nitrate and triflate anions are shown to coordinate directly to the Eu3+ ion forming either ten- or eight-coordinate solvent complexes where the counterions are binding as bidentate or monodentate structures, respectively. Coordination of a terpyridyl ligand to the Eu3+ ion limits the available binding sites for the solvent and anions. In certain cases, the terpyridyl ligand excludes any solvent binding and limits the number of coordinated anions. The solution structure of the Eu-terpyridyl complex with nitrate counterions is shown to have a similar arrangement of Eu3+ coordinating molecules as the crystal structure. This study illustrates how a combination of AIMD and EXAFS can be used to determine how ligands, solvent, and counterions coordinate with the lanthanide ions in solution.

Analysis of the First Ion Coordination Sphere: A Toolkit to Analyze the Coordination Sphere of Ions.

Rapid and accurate approaches to characterizing the coordination structure of an ion are important for designing ligands and quantifying structure-property trends. Here, we introduce AFICS (Analysis of the First Ion Coordination Sphere), a tool written in Python 3 for analyzing the structural and geometric features of the first coordination sphere of an ion over the course of molecular dynamics simulations. The principal feature of AFICS is its ability to quantify the distortion a coordination geometry undergoes compared to uniform polyhedra. This work applies the toolkit to analyze molecular dynamics simulations of the well-defined coordination structure of aqueous Cr3+ along with the more ambiguous structure of aqueous Eu3+ chelated to ethylenediaminetetraacetic acid. The tool is targeted for analyzing ions with fluxional or irregular coordination structures (e.g., solution structures of f-block elements) but is generalized such that it may be applied to other systems.

Correction: The solution structures and relative stability constants of lanthanide-EDTA complexes predicted from computation.

Correction for 'The solution structures and relative stability constants of lanthanide-EDTA complexes predicted from computation' by Ravi D. O'Brien et al., Phys. Chem. Chem. Phys., 2022, 24, 10263-10271, https://doi.org/10.1039/D2CP01081J.

Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase.

Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.

Signal detection shapes ornament allometry in functionally convergent Caribbean Anolis and Southeast Asian Draco lizards.

Visual ornaments have long been assumed to evolve hyper-allometry as an outcome of sexual selection. Yet growing evidence suggests many sexually selected morphologies can exhibit other scaling patterns with body size, including hypo-allometry. The large conspicuous throat fan, or dewlap, of arboreal Caribbean Anolis lizards was one ornament previously thought to conform to the classical expectation of hyper-allometry. We re-evaluated this classic example alongside a second arboreal group of lizards that has also independently evolved a functionally equivalent dewlap, the Southeast Asian Draco lizards. Across multiple closely related species in both genera, the Anolis and Draco dewlaps were either isometric or had hypo-allometric scaling patterns. In the case of the Anolis dewlap, variation in dewlap allometry was predicted by the distance of conspecifics and the light environment in which the dewlap was typically viewed. Signal efficacy, therefore, appears to have driven the evolution of hypo-allometry in what was originally thought to be a sexually selected ornament with hyper-allometry. Our findings suggest that other elaborate morphological structures used in social communication might similarly exhibit isometric or hypo-allometric scaling patterns because of environmental constraints on signal detection.

The solution structures and relative stability constants of lanthanide-EDTA complexes predicted from computation.

Ligand selectivity to specific lanthanide (Ln) ions is key to the separation of rare earth elements from each other. Ligand selectivity can be quantified with relative stability constants (measured experimentally) or relative binding energies (calculated computationally). The relative stability constants of EDTA (ethylenediaminetetraacetic acid) with La3+, Eu3+, Gd3+, and Lu3+ were predicted from relative binding energies, which were quantified using electronic structure calculations with relativistic effects and based on the molecular structures of Ln-EDTA complexes in solution from density functional theory molecular dynamics simulations. The protonation state of an EDTA amine group was varied to study pH ∼7 and ∼11 conditions. Further, simulations at 25 °C and 90 °C were performed to elucidate how structures of Ln-EDTA complexes varying with temperature are related to complex stabilities at different pH conditions. Relative stability trends are predicted from computation for varying Ln3+ ions (La, Eu, Gd, Lu) with a single ligand (EDTA at pH ∼11), as well as for a single Ln3+ ion (La) with varying ligands (EDTA at pH ∼7 and ∼11). Changing the protonation state of an EDTA amine site significantly changes the solution structure of the Ln-EDTA complex resulting in a reduction of the complex stability. Increased Ln-ligand complex stability is correlated to reduced structural variations in solution upon an increase in temperature.

Cheminformatic quantum mechanical enzyme model design: A catechol-O-methyltransferase case study.

To accurately simulate the inner workings of an enzyme active site with quantum mechanics (QM), not only must the reactive species be included in the model but also important surrounding residues, solvent, or coenzymes involved in crafting the microenvironment. Our lab has been developing the Residue Interaction Network Residue Selector (RINRUS) toolkit to utilize interatomic contact network information for automated, rational residue selection and QM-cluster model generation. Starting from an x-ray crystal structure of catechol-O-methyltransferase, RINRUS was used to construct a series of QM-cluster models. The reactant, product, and transition state of the methyl transfer reaction were computed for a total of 550 models, and the resulting free energies of activation and reaction were used to evaluate model convergence. RINRUS-designed models with only 200-300 atoms are shown to converge. RINRUS will serve as a cornerstone for improved and automated cheminformatics-based enzyme model design.

Deep-time convergent evolution in animal communication presented by shared adaptations for coping with noise in lizards and other animals.

Convergence in communication appears rare compared with other forms of adaptation. This is puzzling, given communication is acutely dependent on the environment and expected to converge in form when animals communicate in similar habitats. We uncover deep-time convergence in territorial communication between two groups of tropical lizards separated by over 140 million years of evolution: the Southeast Asian Draco and Caribbean Anolis. These groups have repeatedly converged in multiple aspects of display along common environmental gradients. Robot playbacks to free-ranging lizards confirmed that the most prominent convergence in display is adaptive, as it improves signal detection. We then provide evidence from a sample of the literature to further show that convergent adaptation among highly divergent animal groups is almost certainly widespread in nature. Signal evolution is therefore curbed towards the same set of adaptive solutions, especially when animals are challenged with the problem of communicating effectively in noisy environments.

Quantifying Inter-Residue Contacts through Interaction Energies.

The validity and accuracy of protein modeling is dependent on constructing models that account for the inter-residue interactions crucial for protein structure and function. Residue interaction networks derived from interatomic van der Waals contacts have previously demonstrated usefulness toward designing protein models, but there has not yet been evidence of a connection between network-predicted interaction strength and quantitative interaction energies. This work evaluates the intraprotein contact networks of five proteins against ab initio interaction energies computed using symmetry-adapted perturbation theory. To more appropriately capture the local chemistry of the protein, we deviate from traditional protein network analysis to redefine the interacting nodes in terms of main chain and side chain functional groups rather than complete amino acids. While there is no simple correspondence between the features of the contact network and actual interaction strength, random forest models constructed from minimal structural, network, and chemical descriptors are capable of accurately predicting interaction energy. The results of this work serve as a foundation for the development and improvement of functional group-based contact networks.

Captopril mitigates splenomegaly and myelofibrosis in the Gata1low murine model of myelofibrosis.

Allogeneic stem cell transplantation is currently the only curative therapy for primary myelofibrosis (MF), while the JAK2 inhibitor, ruxolitinib. Has been approved only for palliation. Other therapies are desperately needed to reverse life-threatening MF. However, the cell(s) and cytokine(s) that promote MF remain unclear. Several reports have demonstrated that captopril, an inhibitor of angiotensin-converting enzyme that blocks the production of angiotensin II (Ang II), mitigates fibrosis in heart, lung, skin and kidney. Here, we show that captopril can mitigate the development of MF in the Gata1low mouse model of primary MF. Gata1low mice were treated with 79 mg/kg/d captopril in the drinking water from 10 to 12 months of age. At 13 months of age, bone marrows were examined for fibrosis, megakaryocytosis and collagen expression; spleens were examined for megakaryocytosis, splenomegaly and collagen expression. Treatment of Gata1low mice with captopril in the drinking water was associated with normalization of the bone marrow cellularity; reduced reticulin fibres, splenomegaly and megakaryocytosis; and decreased collagen expression. Our findings suggest that treating with the ACE inhibitors captopril has a significant benefit in overcoming pathological changes associated with MF.

A transition state "trapped"? QM-cluster models of engineered threonyl-tRNA synthetase.

In a recent study [Science, 2015, 347, 6224], protein engineering was used to design a core within the enzyme threonyl-tRNA synthetase (ThrRS) capable of stabilizing the coplanar transition state conformation of an inserted noncanonical p-biphenylalanine (BiPhe) residue. Using the X-ray crystal structures of the preliminary (Protein Data Bank entries 4S02, 4S0J, 4S0L, 4S0I, and 4S0K) and final (PDB entry 4S03) ThrRS proteins, fully quantum mechanical (QM) cluster models were constructed and analyzed. Density functional theory and molecular dynamics computations were performed to investigate the energetic profiles of BiPhe dihedral rotation within the ThrRS models. For the 4S03 model, results indicate that steric and hydrophobic forces of the residues surrounding BiPhe eliminate the coplanar transition state entirely. Molecular dynamics simulations were carried out that confirmed the extent of BiPhe rotational flexibility, and provided additional information on barrier heights of full BiPhe rotation. Transition states of near-coplanar biphenyl rings of BiPhe were found for the 4S0I and 4S0K models, but are not likely persistent on any observable timescale. The dihedral angle of the biphenyl moiety is thermally allowed to fluctuate within the ThrRS protein core models by a range of 17°-26°. BiPhe-residue interaction counts (RICs) were used to compare the interaction differences among the different ThrRS cores. The RICs demonstrate how BiPhe is compacted within the 4S03 core, resulting in the experimentally observed "trapped" coplanar transition state analogue. This work presents a unique application of QM-cluster models towards studying the inner workings of proteins, and suggests avenues that computational chemistry can be used to further guide bioengineering.

Ecological Release from Aquatic Predation Is Associated with the Emergence of Marine Blenny Fishes onto Land.

An ecological release from competition or predation is a frequent adaptive explanation for the colonization of novel environments, but empirical data are limited. On the island of Rarotonga, several blenny fish species appear to be in the process of colonizing land. Anecdotal observations have implied that aquatic predation is an important factor in prompting such amphibious fish behavior. We provide evidence supporting this hypothesis by demonstrating that amphibious blennies shift their abundance up and down the shoreline to remain above predatory fishes that periodically move into intertidal areas during high tide. A predation experiment using blenny mimics confirmed a high risk of aquatic predation for blennies, significantly higher than predation experienced on land. These data suggest that predation has played an active role in promoting terrestrial activity in amphibious blennies and provide a rare example of how ecological release from predation could drive the colonization of a novel environment.

Circulating Tumor DNA to Monitor Therapy for Aggressive B-Cell Lymphomas.

OPINION STATEMENT: The goal of therapy for aggressive B-cell lymphomas such as diffuse large B-cell lymphoma (DLBCL) and Burkitt lymphoma (BL) is to achieve cure. Combination chemotherapy with rituximab cures most patients, but those with recurrent disease have a poor prognosis. Medical imaging scans such as computed tomography (CT) and positron emission tomography (PET) are the principal methods to assess response and monitor for disease relapse after therapy but are fundamentally limited by risks of radiation, cost, and a lack of tumor specificity. Novel sequencing-based DNA monitoring methods are capable of quantifying small amounts of circulating tumor DNA (ctDNA) before, during, and after therapy for mature B-cell lymphomas. Detection of ctDNA encoding clonal rearranged variable-diversity-joining (VDJ) receptor gene sequences has demonstrated improved analytical sensitivity and enhanced tumor specificity compared to imaging scans in DLBCL, offering broad clinical applicability across a range of aggressive B-cell lymphomas. Molecular monitoring of ctDNA has vaulted into the spotlight as a promising non-invasive tool with immediate clinical impact on monitoring for recurrence after therapy prior to clinical symptoms. As these clinical observations are validated, ctDNA monitoring needs to be investigated as a tool for response-adapted therapy and as a marker of minimal residual disease upon completion of therapy in aggressive B-cell lymphomas. Molecular monitoring of ctDNA holds tremendous promise that may ultimately transform our ability to monitor disease in aggressive B-cell lymphomas.

Repeated evolution and the impact of evolutionary history on adaptation.

BACKGROUND: Whether natural selection can erase the imprint of past evolutionary history from phenotypes has been a topic of much debate. A key source of evidence that present-day selection can override historically contingent effects comes from the repeated evolution of similar adaptations in different taxa. Yet classic examples of repeated evolution are often among closely related taxa, suggesting the likelihood that similar adaptations evolve is contingent on the length of time separating taxa. To resolve this, we performed a meta-analysis of published reports of repeated evolution. RESULTS: Overall, repeated evolution was far more likely to be documented among closely related than distantly related taxa. However, not all forms of adaptation seemed to exhibit the same pattern. The evolution of similar behavior and physiology seemed frequent in distantly related and closely related taxa, while the repeated evolution of morphology was heavily skewed towards closely related taxa. Functionally redundant characteristics-alternative phenotypes that achieve the same functional outcome-also appeared less contingent. CONCLUSIONS: If the literature provides a reasonable reflection of the incidence of repeated evolution in nature, our findings suggest that natural selection can overcome contingent effects to an extent, but it depends heavily on the aspect of the phenotype targeted by selection.

Heterogeneous atypical cell populations are present in blood of metastatic breast cancer patients.

INTRODUCTION: Circulating tumor cells (CTCs) are commonly isolated from the blood by targeting the epithelial cell adhesion molecule (EpCAM) through positive selection. However, EpCAM can be downregulated during metastatic progression, or it can be initially not present. We designed the present prospective trial to characterize CTCs as well as other circulating cell populations in blood samples from women with metastatic breast cancer without EpCAM-dependent enrichment and/or isolation technology. METHODS: A total of 32 patients with metastatic breast cancer were enrolled, and blood samples were processed using a previously described negative depletion immunomagnetic methodology. Samples from healthy volunteers were run as controls (n = 5). Multistep sequential labeling was performed to label and fix cell-surface markers followed by permeabilization for cytokeratins (CK) 8, 18 and 19. Multiparametric flow cytometry (FCM) analysis was conducted using a BD LSR II flow cytometer or a BD FACSAria II or FACSAria III cell sorter. Immunocytochemical staining on postenrichment specimens for DAPI, EpCAM, CD45, CK, epidermal growth factor receptor and vimentin was performed. Expression of these markers was visualized using confocal microscopy (CM). RESULTS: CD45-negative/CK-positive (CD45- CK+) populations with EpCAM + and EpCAM - expression were identified with both FCM and CM from the negatively enriched patient samples. In addition, EpCAM + and EpCAM - populations that were CK + and coexpressing the pan-hematopoietic marker CD45 were also noted. There were more CK + EpCAM - events/ml than CK + EpCAM + events/ml in both the CD45- and CD45+ fractions (both statistically significant at P ≤ 0.0005). The number of CK + CD45- and CK + CD45+ events per milliliter in blood samples (regardless of EpCAM status) was higher in patient samples than in normal control samples (P ≤ 0.0005 and P ≤ 0.026, respectively). Further, a significant fraction of the CK + CD45+ events also expressed CD68, a marker associated with tumor-associated macrophages. Higher levels of CD45-CK + EpCAM - were associated with worse overall survival (P = 0.0292). CONCLUSIONS: Metastatic breast cancer patients have atypical cells that are CK + EpCAM - circulating in their blood. Because a substantial number of these patients do not have EpCAM + CTCs, additional studies are needed to evaluate the role of EpCAM - circulating cells as a prognostic and predictive marker.

Tissue Eosinophilia in B-cell Lymphoma: An Underrecognized Phenomenon.

Tissue eosinophilia is seldom reported in B-cell lymphoma. It poses diagnostic challenges and frequently leads to the consideration of other diagnoses, particularly T-cell lymphomas. The scarce literature underscores the need for in-depth studies to enhance awareness and understanding of this phenomenon. We investigated 54 cases of B-cell lymphoma with notable tissue eosinophils, analyzing clinical information, hematoxylin and eosin staining, immunohistochemistry, and PCR-based clonality analysis. Nodal marginal zone lymphoma (NMZL) emerged as the most prevalent type (n=26), followed by B-cell lymphoma, not otherwise specified (n=13), diffuse large B-cell lymphoma (n=10), follicular lymphoma (n=2), chronic lymphocytic leukemia/small lymphocytic lymphoma (n=1), extranodal marginal zone lymphoma (n=1), and primary cutaneous marginal zone lymphoma (n=1). Shared features across different lymphoma types, best exemplified by NMZL, included plasmacytic differentiation (57.7%), increased vascularity (84.6%) with a tendency for perivascular distribution of neoplastic cells, and a tumor microenvironment abundant in T cells and histiocytes; some cases showed increased PD-1-positive cells. These features often raise consideration of angioimmunoblastic T-cell lymphoma. Along with clonality analysis, features supporting the diagnosis of B-cell lymphoma included cytological atypia in B cells rather than T cells, and the lack of follicular dendritic cell meshwork expansion. In addition, diffuse large B-cell lymphoma frequently exhibited interfollicular distribution and monocytoid appearance, indicating the possibility of transformed NMZL. Collectively, tissue eosinophilia can occur in diverse B-cell lymphomas but is most prevalent in tumors with a postgerminal stage of differentiation.

Splenic manifestations of chronic autoimmune disorder: a report of five cases with histiocytic necrotizing change in four cases.

AIMS: Autoimmune diseases (AD) are associated with lymphadenopathy and splenomegaly. Changes in the spleen have not been characterized completely in AD; we describe splenectomy specimens from five patients with chronic AD, highlighting the presence of necrotizing histiocytosis. METHODS AND RESULTS: Of the patients (three males and two females; mean 40 years), four had systemic lupus erythematosus; one had rheumatoid arthritis. All had moderate splenomegaly (213-803 g, mean 421 g). Four cases exhibited necrosis with apoptosis and karyorrhectic debris occurring in the white pulp and minimal acute inflammation; one showed florid follicular hyperplasia. Splenic involvement ranged from focal to extensive. Plasma cells were negative for IgG4. Haematoxylin bodies were not identified. Stains for infectious organisms were negative. Immunohistochemical studies showed that lymphocytes surrounding the necrosis were a mixture of CD4(+) and CD8(+) T cells; CD123-positive plasmacytoid dendritic cells were not present, and staining for kappa and lambda light chains showed no clonality. 16S rDNA PCR was performed; no amplification was seen in three of four cases tested for bacteria specific rDNA. Epstein-Barr virus-encoded RNA (EBER) in situ hybridization studies highlighted rare positive cells in four cases. CONCLUSIONS: Splenomegaly in AD is thought to be hyperplasic, but we present four cases showing histiocytic necrosis, a finding which should be considered part of the spectrum of AD in the spleen.

Impact of presurgical breast magnetic resonance imaging (MRI) on surgical planning - a retrospective analysis from a private radiology group.

Presurgical bilateral breast MRI is being increasingly utilized in newly diagnosed breast cancer patients to identify the presence of other potentially malignant lesions and thereby guide surgical planning. The authors have performed a retrospective analysis of presurgical MRIs in 445 breast cancer patients, consecutively seen in a community practice setting, to determine what risk factors may be associated with bilateral synchronous breast cancer. Clinical, imaging, surgical, and pathological data were collected on these newly diagnosed patients, who underwent presurgical breast MRIs from November 1, 2008 through October 31, 2009. The MRIs detected 84 incremental (additional) malignancies in 66 (14.8%) patients, including 22 patients in the contralateral breast (4.9%) and 48 patients who had ipsilateral incremental malignancies (10.8%). The presurgical breast MRIs had a sensitivity of 99% and a specificity of 86%. Based on the imaging study and the subsequent biopsy(s), 105 patients (23.6%) had changes in their surgical planning. For those with incremental breast malignancies, the synchronous bilateral breast cancer patients were older (67.5 ± 13.5 years, p = 0.01). Patients with ipsilateral incremental malignancies were younger (55.6 ± 10.4 years of age, p = 0.03). A strong family history and ductal carcinoma in situ (DCIS) as the index lesion also appeared to be associated with increased bilateral malignancy detection. The authors conclude that although presurgical MRI remains controversial, our data support its role in surgical planning for newly diagnosed breast cancer patients.