Low effective mechanical advantage of giraffes' limbs during walking reveals trade-off between limb length and locomotor performance.

Giraffes (Giraffa camelopardalis) possess specialized locomotor morphology, namely elongate and gracile distal limbs. While this contributes to their overall height and enhances feeding behavior, we propose that the combination of long limb segments and modest muscle lever arms results in low effective mechanical advantage (EMA, the ratio of in-lever to out-lever moment arms), when compared with other cursorial mammals. To test this, we used a combination of experimentally measured kinematics and ground reaction forces (GRFs), musculoskeletal modeling, and inverse dynamics to calculate giraffe forelimb EMA during walking. Giraffes walk with an EMA of 0.34 (±0.05 SD), with no evident association with speed within their walking gait. Giraffe EMA was about four times lower than expectations extrapolated from other mammals, ranging from 0.03 to 297 kg, and this provides further evidence that EMA plateaus or even diminishes in mammals exceeding horse size. We further tested the idea that limb segment length is a factor which determines EMA, by modeling the GRF and muscle moment arms in the extinct giraffid Sivatherium giganteum and the other extant giraffid, Okapia johnstoni. Giraffa and Okapia shared similar EMA, despite a four to sixfold difference in body mass (Okapia EMA = 0.38). In contrast, Sivatherium, sharing a similar body mass with Giraffa, had greater EMA (0.59), which we propose reflects behavioral differences, such as a somewhat increased capability for athletic performance. Our modeling approach suggests that limb length is a determinant of GRF moment arm magnitude and that unless muscle moment arms scale isometrically with limb length, tall mammals are prone to low EMA.

Computational simulation of JAK/STAT signaling in somatic versus germline stem cells.

BACKGROUND: The Janus kinase/signal transducer and activator of transcription (JAK/STAT) signaling pathway regulates a variety of cellular processes. A major activation event in this pathway involves the phosphorylation of a tyrosine of STAT, converting unphosphorylated STAT (uSTAT) to phosphorylated STAT (pSTAT), an active transcription factor. In a noncanonical role, uSTAT contributes to the maintenance of heterochromatin stability. As such, an increase in pSTAT concurrently reduces uSTAT, resulting in heterochromatin loss, as observed in Drosophila somatic tissues. Paradoxically, an opposing phenomenon occurs in Drosophila male germline stem cells (GSCs), where the JAK/STAT pathway remains persistently active due to a continuous supply of ligands. Here, computational simulations were employed to dissect JAK/STAT pathway activation under different cellular contexts, mimicking somatic and germline cells. In these simulations, ordinary differential equations were leveraged to replicate the chemical reactions governing JAK/STAT signaling under different conditions. RESULTS: The outcomes indicate that transient ligand stimulation, typical in somatic tissues, led to a momentary reduction in uSTAT levels. Conversely, sustained ligand stimulation, a characteristic feature of the GSC niche, resulted in elevated uSTAT levels at equilibrium. CONCLUSION: The simulation suggests that the duration of ligand exposure could explain the observed opposite effects of JAK/STAT activation on heterochromatin in somatic versus GSCs.

Translating myosin-binding protein C and titin abnormalities to whole-heart function using a novel calcium-contraction coupling model.

Mutations in cardiac myosin-binding protein C (cMyBP-C) or titin may respectively lead to hypertrophic (HCM) or dilated (DCM) cardiomyopathies. The mechanisms leading to these phenotypes remain unclear because of the challenge of translating cellular abnormalities to whole-heart and system function. We developed and validated a novel computer model of calcium-contraction coupling incorporating the role of cMyBP-C and titin based on the key assumptions: 1) tension in the thick filament promotes cross-bridge attachment mechanochemically, 2) with increasing titin tension, more myosin heads are unlocked for attachment, and 3) cMyBP-C suppresses cross-bridge attachment. Simulated stationary calcium-tension curves, isotonic and isometric contractions, and quick release agreed with experimental data. The model predicted that a loss of cMyBP-C function decreases the steepness of the calcium-tension curve, and that more compliant titin decreases the level of passive and active tension and its dependency on sarcomere length. Integrating this cellular model in the CircAdapt model of the human heart and circulation showed that a loss of cMyBP-C function resulted in HCM-like hemodynamics with higher left ventricular end-diastolic pressures and smaller volumes. More compliant titin led to higher diastolic pressures and ventricular dilation, suggesting DCM-like hemodynamics. The novel model of calcium-contraction coupling incorporates the role of cMyBP-C and titin. Its coupling to whole-heart mechanics translates changes in cellular calcium-contraction coupling to changes in cardiac pump and circulatory function and identifies potential mechanisms by which cMyBP-C and titin abnormalities may develop into HCM and DCM phenotypes. This modeling platform may help identify distinct mechanisms underlying clinical phenotypes in cardiac diseases.

The molecular basis of dapsone activation of CYP2C9-catalyzed nonsteroidal anti-inflammatory drug oxidation.

Positive heterotropic cooperativity, or "activation," results in an instantaneous increase in enzyme activity in the absence of an increase in protein expression. Thus, cytochrome P450 (CYP) enzyme activation presents as a potential drug-drug interaction mechanism. It has been demonstrated previously that dapsone activates the CYP2C9-catalyzed oxidation of a number of nonsteroidal anti-inflammatory drugs in vitro. Here, we conducted molecular dynamics simulations (MDS) together with enzyme kinetic investigations and site-directed mutagenesis to elucidate the molecular basis of the activation of CYP2C9-catalyzed S-flurbiprofen 4'-hydroxylation and S-naproxen O-demethylation by dapsone. Supplementation of incubations of recombinant CYP2C9 with dapsone increased the catalytic efficiency of flurbiprofen and naproxen oxidation by 2.3- and 16.5-fold, respectively. MDS demonstrated that activation arises predominantly from aromatic interactions between the substrate, dapsone, and the phenyl rings of Phe114 and Phe476 within a common binding domain of the CYP2C9 active site, rather than involvement of a distinct effector site. Mutagenesis of Phe114 and Phe476 abrogated flurbiprofen and naproxen oxidation, and MDS and kinetic studies with the CYP2C9 mutants further identified a pivotal role of Phe476 in dapsone activation. MDS additionally showed that aromatic stacking interactions between two molecules of naproxen are necessary for binding in a catalytically favorable orientation. In contrast to flurbiprofen and naproxen, dapsone did not activate the 4'-hydroxylation of diclofenac, suggesting that the CYP2C9 active site favors cooperative binding of nonsteroidal anti-inflammatory drugs with a planar or near-planar geometry. More generally, the work confirms the utility of MDS for investigating ligand binding in CYP enzymes.

Mathematical models disentangle the role of IL-10 feedbacks in human monocytes upon proinflammatory activation.

Inflammation is one of the vital mechanisms through which the immune system responds to harmful stimuli. During inflammation, proinflammatory and anti-inflammatory cytokines interplay to orchestrate fine-tuned and dynamic immune responses. The cytokine interplay governs switches in the inflammatory response and dictates the propagation and development of the inflammatory response. Molecular pathways underlying the interplay are complex, and time-resolved monitoring of mediators and cytokines is necessary as a basis to study them in detail. Our understanding can be advanced by mathematical models that enable to analyze the system of interactions and their dynamical interplay in detail. We, therefore, used a mathematical modeling approach to study the interplay between prominent proinflammatory and anti-inflammatory cytokines with a focus on tumor necrosis factor and interleukin 10 (IL-10) in lipopolysaccharide-primed primary human monocytes. Relevant time-resolved data were generated by experimentally adding or blocking IL-10 at different time points. The model was successfully trained and could predict independent validation data and was further used to perform simulations to disentangle the role of IL-10 feedbacks during an acute inflammatory event. We used the insight to obtain a reduced predictive model including only the necessary IL-10-mediated feedbacks. Finally, the validated reduced model was used to predict early IL-10-tumor necrosis factor switches in the inflammatory response. Overall, we gained detailed insights into fine-tuning of inflammatory responses in human monocytes and present a model for further use in studying the complex and dynamic process of cytokine-regulated acute inflammation.

Cell-to-cell heterogeneity in ion channel conductance impacts substrate vulnerability to arrhythmia.

Determining whether an ectopic depolarization will lead to a self-perpetuating arrhythmia is of critical importance in determining arrhythmia risk, so it is necessary to understand what factors impact substrate vulnerability. This study sought to explore the impact of cell-to-cell heterogeneity in ion channel conductance on substrate vulnerability to arrhythmia by measuring the duration of the vulnerable window in computational models of one-dimensional cables of ventricular cardiomyocytes. We began by using a population of uniform cable models to determine the mechanisms underlying the vulnerable window phenomenon. We found that in addition to the known importance of GNa, the conductances GCa,L and GKr also play a minor role in determining the vulnerable window duration. We also found that a steeper slope of the repolarizing action potential during the vulnerable window correlated with a shorter vulnerable window duration in uniform cables. We applied our understanding from these initial simulations to an investigation of the vulnerable window in heterogeneous cable models. The heterogeneous cables displayed a great deal of intra-cable variation in vulnerable window duration, highly sensitive to the cardiomyocytes in the local environment of the ectopic stimulus. Coupling strength modulated not only the magnitude of the vulnerable window duration but also the extent of intra-tissue variability in vulnerable window duration.NEW & NOTEWORTHY We investigate the impact of cell-to-cell heterogeneity in ion channel conductance on substrate vulnerability to arrhythmia by measuring the vulnerable window duration in computational cardiomyocyte cable models. We demonstrate a wide range of intra-cable variability in vulnerable window duration (VWD) and show how this is changed by ion channel block and coupling strength perturbations.

Multifunctional Chiral Chemically-Powered Micropropellers for Cargo Transport and Manipulation.

Practical applications of synthetic self-propelled nano and microparticles for microrobotics, targeted drug delivery, and manipulation at the nanoscale are rapidly expanding. However, fabrication limitations often hinder progress, resulting in relatively simple shapes and limited functionality. Here, taking advantage of 3D nanoscale printing, chiral micropropellers powered by the hydrogen peroxide reduction reaction are fabricated. Due to their chirality, the propellers exhibit multifunctional behavior controlled by an applied magnetic field: spinning in place (loitering), directed migration in the prescribed direction, capture, and transport of polymer cargo particles. Design parameters of the propellers are optimized by computation modeling based on mesoscale molecular dynamics. It is predicted by computer simulations, and confirmed experimentally, that clockwise rotating propellers attract each other and counterclockwise repel. These results shed light on how chirality and shape optimization enhance the functionality of synthetic autonomous micromachines.

Electrophysiological patterns and structural substrates of Brugada syndrome: Critical appraisal and computational analyses.

Brugada syndrome (BrS) is a cardiac electrophysiological disease with unknown etiology, associated with sudden cardiac death. Symptomatic patients are treated with implanted cardiac defibrillator, but no risk stratification strategy is effective in patients that are at low to medium arrhythmic risk. Cardiac computational modeling is an emerging tool that can be used to verify the hypotheses of pathogenesis and inspire new risk stratification strategies. However, to obtain reliable results computational models must be validated with consistent experimental data. We reviewed the main electrophysiological and structural variables from BrS clinical studies to assess which data could be used to validate a computational approach. Activation delay in the epicardial right ventricular outflow tract is a consistent finding, as well as increased fibrosis and subclinical alterations of right ventricular functional and morphological parameters. The comparison between other electrophysiological variables is hindered by methodological differences between studies, which we commented. We conclude by presenting a recent theory unifying electrophysiological and structural substrate in BrS and illustrate how computational modeling could help translation to risk stratification.

Computer model of IL-6-dependent rheumatoid arthritis in F759 mice.

The interleukin-6 (IL-6) amplifier, which describes the simultaneous activation of signal transducer and activator of transcription 3 (STAT3) and NF-κb nuclear factor kappa B (NF-κB), in synovial fibroblasts causes the infiltration of immune cells into the joints of F759 mice. The result is a disease that resembles human rheumatoid arthritis. However, the kinetics and regulatory mechanisms of how augmented transcriptional activation by STAT3 and NF-κB leads to F759 arthritis is unknown. We here show that the STAT3-NF-κB complex is present in the cytoplasm and nucleus and accumulates around NF-κB binding sites of the IL-6 promoter region and established a computer model that shows IL-6 and IL-17 (interleukin 17) signaling promotes the formation of the STAT3-NF-κB complex followed by its binding on promoter regions of NF-κB target genes to accelerate inflammatory responses, including the production of IL-6, epiregulin, and C-C motif chemokine ligand 2 (CCL2), phenotypes consistent with in vitro experiments. The binding also promoted cell growth in the synovium and the recruitment of T helper 17 (Th17) cells and macrophages in the joints. Anti-IL-6 blocking antibody treatment inhibited inflammatory responses even at the late phase, but anti-IL-17 and anti-TNFα antibodies did not. However, anti-IL-17 antibody at the early phase showed inhibitory effects, suggesting that the IL-6 amplifier is dependent on IL-6 and IL-17 stimulation at the early phase, but only on IL-6 at the late phase. These findings demonstrate the molecular mechanism of F759 arthritis can be recapitulated in silico and identify a possible therapeutic strategy for IL-6 amplifier-dependent chronic inflammatory diseases.

Impact of surrounding tissue-type and peri-electrode gap in stereoelectroencephalography guided (SEEG) radiofrequency thermocoagulation (RF-TC): a computational study.

PURPOSE: To use computational modeling to provide a complete and logical description of the electrical and thermal behavior during stereoelectroencephalography-guided (SEEG) radiofrequency thermo-coagulation (RF-TC). METHODS: A coupled electrical-thermal model was used to obtain the temperature distributions in the tissue during RF-TC. The computer model was first validated by an ex vivo model based on liver fragments and later used to study the impact of three different factors on the coagulation zone size: 1) the difference in the tissue surrounding the electrode (gray/white matter), 2) the presence of a peri-electrode gap occupied by cerebrospinal fluid (CSF), and 3) the energy setting used (power-duration). RESULTS: The model built for the experimental validation was able to predict both the evolution of impedance and the short diameter of the coagulation zone (error < 0.01 mm) reasonably well but overestimated the long diameter by 2 - 3 mm. After adapting the model to clinical conditions, the simulation showed that: 1) Impedance roll-off limited the coagulation size but involved overheating (around 100 °C); 2) The type of tissue around the contacts (gray vs. white matter) had a moderate impact on the coagulation size (maximum difference 0.84 mm), and 3) the peri-electrode gap considerably altered the temperature distributions, avoided overheating, although the diameter of the coagulation zone was not very different from the no-gap case (<0.2 mm). CONCLUSIONS: This study showed that computer modeling, especially subject- and scenario-specific modeling, can be used to estimate in advance the electrical and thermal performance of the RF-TC in brain tissue.

Formation of Amyloid-Like Conformational States of ß-Structured Membrane Proteins on the Example of OMPF Porin from the Yersinia pseudotuberculosis Outer Membrane.

The work presents results of the in vitro and in silico study of formation of amyloid-like structures under harsh denaturing conditions by non-specific OmpF porin of Yersinia pseudotuberculosis (YpOmpF), a membrane protein with ß-barrel conformation. It has been shown that in order to obtain amyloid-like porin aggregates, preliminary destabilization of its structure in a buffer solution with acidic pH at elevated temperature followed by long-term incubation at room temperature is necessary. After heating at 95°C in a solution with pH 4.5, significant conformational rearrangements are observed in the porin molecule at the level of tertiary and secondary structure of the protein, which are accompanied by the increase in the content of total ß-structure and sharp decrease in the value of characteristic viscosity of the protein solution. Subsequent long-term exposure of the resulting unstable intermediate YpOmpF at room temperature leads to formation of porin aggregates of various shapes and sizes that bind thioflavin T, a specific fluorescent dye for the detection of amyloid-like protein structures. Compared to the initial protein, early intermediates of the amyloidogenic porin pathway, oligomers, have been shown to have increased toxicity to the Neuro-2aCCL-131™ mouse neuroblastoma cells. The results of computer modeling and analysis of the changes in intrinsic fluorescence during protein aggregation suggest that during formation of amyloid-like aggregates, changes in the structure of YpOmpF affect not only the areas with an internally disordered structure corresponding to the external loops of the porin, but also main framework of the molecule, which has a rigid spatial structure inherent to ß-barrel.

Simulating left atrial arrhythmias with an interactive N-body model.

BACKGROUND: Heart disease and strokes are leading global killers. While atrial arrhythmias are not deadly by themselves, they can disrupt blood flow in the heart, causing blood clots. These clots can travel to the brain, causing strokes, or to the coronary arteries, causing heart attacks. Additionally, prolonged periods of elevated heart rates can lead to structural and functional changes in the heart, ultimately leading to heart failure if untreated. The left atrium, with its more complex topology, is the primary site for complex arrhythmias. Much remains unknown about the causes of these arrhythmias, and computer modeling is employed to study them. METHODS: We use N-body modeling techniques and parallel computing to build an interactive model of the left atrium. Through user input, individual muscle attributes can be adjusted, and ectopic events can be placed to induce arrhythmias in the model. Users can test ablation scenarios to determine the most effective way to eliminate these arrhythmias. RESULTS: We set up muscle conditions that either spontaneously generate common arrhythmias or, with a properly timed and located ectopic event, induce an arrhythmia. These arrhythmias were successfully eliminated with simulated ablation. CONCLUSIONS: We believe the model could be useful to doctors, researchers, and medical students studying left atrial arrhythmias.

Fully angularly resolved 3D microrheology with optical tweezers.

Microrheology with optical tweezers (MOT) is an all-optical technique that allows the user to investigate a materials' viscoelastic properties at microscopic scales, and is particularly useful for those materials that feature complex microstructures, such as biological samples. MOT is increasingly being employed alongside 3D imaging systems and particle tracking methods to generate maps showing not only how properties may vary between different points in a sample but also how at a single point the viscoelastic properties may vary with direction. However, due to the diffraction limited shape of focussed beams, optical traps are inherently anisotropic in 3D. This can result in a significant overestimation of the fluids' viscosity in certain directions. As such, the rheological properties can only be accurately probed along directions parallel or perpendicular to the axis of trap beam propagation. In this work, a new analytical method is demonstrated to overcome this potential artefact. This is achieved by performing principal component analysis on 3D MOT data to characterise the trap, and then identify the frequency range over which trap anisotropy influences the data. This approach is initially applied to simulated data for a Newtonian fluid where the trap anisotropy induced maximum error in viscosity is reduced from ~ 150% to less than 6%. The effectiveness of the method is corroborated by experimental MOT measurements performed with water and gelatine solutions, thus confirming that the microrheology of a fluid can be extracted reliably across a wide frequency range and in any arbitrary direction. This work opens the door to fully spatially and angularly resolved 3D mapping of the rheological properties of soft materials over a broad frequency range.

Multiscale Computer Model of the Spinal Dorsal Horn Reveals Changes in Network Processing Associated with Chronic Pain.

Pain-related sensory input is processed in the spinal dorsal horn (SDH) before being relayed to the brain. That processing profoundly influences whether stimuli are correctly or incorrectly perceived as painful. Significant advances have been made in identifying the types of excitatory and inhibitory neurons that comprise the SDH, and there is some information about how neuron types are connected, but it remains unclear how the overall circuit processes sensory input or how that processing is disrupted under chronic pain conditions. To explore SDH function, we developed a computational model of the circuit that is tightly constrained by experimental data. Our model comprises conductance-based neuron models that reproduce the characteristic firing patterns of spinal neurons. Excitatory and inhibitory neuron populations, defined by their expression of genetic markers, spiking pattern, or morphology, were synaptically connected according to available qualitative data. Using a genetic algorithm, synaptic weights were tuned to reproduce projection neuron firing rates (model output) based on primary afferent firing rates (model input) across a range of mechanical stimulus intensities. Disparate synaptic weight combinations could produce equivalent circuit function, revealing degeneracy that may underlie heterogeneous responses of different circuits to perturbations or pathologic insults. To validate our model, we verified that it responded to the reduction of inhibition (i.e., disinhibition) and ablation of specific neuron types in a manner consistent with experiments. Thus validated, our model offers a valuable resource for interpreting experimental results and testing hypotheses in silico to plan experiments for examining normal and pathologic SDH circuit function.SIGNIFICANCE STATEMENT We developed a multiscale computer model of the posterior part of spinal cord gray matter (spinal dorsal horn), which is involved in perceiving touch and pain. The model reproduces several experimental observations and makes predictions about how specific types of spinal neurons and synapses influence projection neurons that send information to the brain. Misfiring of these projection neurons can produce anomalous sensations associated with chronic pain. Our computer model will not only assist in planning future experiments, but will also be useful for developing new pharmacotherapy for chronic pain disorders, connecting the effect of drugs acting at the molecular scale with emergent properties of neurons and circuits that shape the pain experience.

Phosphatidylinositol phosphates modulate interactions between the StarD4 sterol trafficking protein and lipid membranes.

There is substantial evidence for extensive nonvesicular sterol transport in cells. For example, lipid transfer by the steroidogenic acute regulator-related proteins (StarD) containing a StarT domain has been shown to involve several pathways of nonvesicular trafficking. Among the soluble StarT domain-containing proteins, StarD4 is expressed in most tissues and has been shown to be an effective sterol transfer protein. However, it was unclear whether the lipid composition of donor or acceptor membranes played a role in modulating StarD4-mediated transport. Here, we used fluorescence-based assays to demonstrate a phosphatidylinositol phosphate (PIP)-selective mechanism by which StarD4 can preferentially extract sterol from liposome membranes containing certain PIPs (especially, PI(4,5)P2 and to a lesser degree PI(3,5)P2). Monophosphorylated PIPs and other anionic lipids had a smaller effect on sterol transport. This enhancement of transport was less effective when the same PIPs were present in the acceptor membranes. Furthermore, using molecular dynamics (MD) simulations, we mapped the key interaction sites of StarD4 with PIP-containing membranes and identified residues that are important for this interaction and for accelerated sterol transport activity. We show that StarD4 recognizes membrane-specific PIPs through specific interaction with the geometry of the PIP headgroup as well as the surrounding membrane environment. Finally, we also observed that StarD4 can deform membranes upon longer incubations. Taken together, these results suggest a mechanism by which PIPs modulate cholesterol transfer activity via StarD4.

Structural insights into the role of the WW2 domain on tandem WW-PPxY motif interactions of oxidoreductase WWOX.

Class I WW domains are present in many proteins of various functions and mediate protein interactions by binding to short linear PPxY motifs. Tandem WW domains often bind peptides with multiple PPxY motifs, but the interplay of WW-peptide interactions is not always intuitive. The WW domain-containing oxidoreductase (WWOX) harbors two WW domains: an unstable WW1 capable of PPxY binding and stable WW2 that cannot bind PPxY. The WW2 domain has been suggested to act as a WW1 domain chaperone, but the underlying mechanism of its chaperone activity remains to be revealed. Here, we combined NMR, isothermal calorimetry, and structural modeling to elucidate the roles of both WW domains in WWOX binding to its PPxY-containing substrate ErbB4. Using NMR, we identified an interaction surface between these two domains that supports a WWOX conformation compatible with peptide substrate binding. Isothermal calorimetry and NMR measurements also indicated that while binding affinity to a single PPxY motif is marginally increased in the presence of WW2, affinity to a dual-motif peptide increases 10-fold. Furthermore, we found WW2 can directly bind double-motif peptides using its canonical binding site. Finally, differential binding of peptides in mutagenesis experiments was consistent with a parallel N- to C-terminal PPxY tandem motif orientation in binding to the WW1-WW2 tandem domain, validating structural models of the interaction. Taken together, our results reveal the complex nature of tandem WW-domain organization and substrate binding, highlighting the contribution of WWOX WW2 to both protein stability and target binding.

AI protein structure prediction-based modeling and mutagenesis of a protostome receptor and peptide ligands reveal key residues for their interaction.

The protostome leucokinin (LK) signaling system, including LK peptides and their G protein-coupled receptors, has been characterized in several species. Despite the progress, molecular mechanisms governing LK peptide-receptor interactions remain to be elucidated. Previously, we identified a precursor protein for Aplysia leucokinin-like peptides (ALKs) that contains the greatest number of amidated peptides among LK precursors in all species identified so far. Here, we identified the first ALK receptor from Aplysia, ALKR. We used cell-based IP1 activation assays to demonstrate that two ALK peptides with the most copies, ALK1 and ALK2, activated ALKR with high potencies. Other endogenous ALK-derived peptides bearing the FXXWX-amide motif also activated ALKR to various degrees. Our examination of cross-species activity of ALKs with the Anopheles LK receptor was consistent with a critical role for the FXXWX-amide motif in receptor activity. Furthermore, we showed, through alanine substitution of ALK1, the highly conserved phenylalanine (F), tryptophan (W), and C-terminal amidation were each essential for receptor activation. Finally, we used an artificial intelligence-based protein structure prediction server (Robetta) and Autodock Vina to predict the ligand-bound conformation of ALKR. Our model predicted several interactions (i.e., hydrophobic interactions, hydrogen bonds, and amide-pi stacking) between ALK peptides and ALKR, and several of our substitution and mutagenesis experiments were consistent with the predicted model. In conclusion, our results provide important information defining possible interactions between ALK peptides and their receptors. The workflow utilized here may be useful for studying other ligand-receptor interactions for a neuropeptide signaling system, particularly in protostomes.

NFX1-123: A potential therapeutic target in cervical cancer.

NFX1-123 is a splice variant isoform of the NFX1 gene. It is highly expressed in cervical cancers caused by HPV, and NFX1-123 is a protein partner with the HPV oncoprotein E6. Together, NFX1-123 and E6 affect cellular growth, longevity, and differentiation. The expression status of NFX1-123 in cancers beyond cervical and head and neck cancers, and its potential as therapeutic target, have not been investigated. TSVdb of TCGA was used to quantify NFX1-123 expression in 24 cancers compared with normal tissues. The NFX1-123 protein structure was predicted and then submitted to retrieve suitable drug molecules. The top four compounds, found to bind in silico to NFX1-123, were tested experimentally to determine their effects on NFX1-123-related cellular growth, survival, and migration. 46% of cancers (11 of 24 had significant differences in NFX1-123 expression, with nine having had greater NFX1-123 expression, when compared with adjacent normal tissues. Bioinformatics and proteomic predictive analysis modeled the three-dimensional structure of NFX1-123, and drug libraries were screened for high-binding affinity compounds using this modeled structure. Seventeen drugs with binding energies ranging from -1.3 to -10 Kcal/mol were identified. The top four compounds were used to treat HPV- and HPV+ cervical cancer cell lines, three of which (Ropitoin, R428 and Ketoconazole) reduced NFX1-123 protein levels, inhibited cellular growth, survival, and migration, and enhanced the cytotoxicity of Cisplatin. These findings highlight cancers expressing high levels of NFX1-123, and drugs that target it, may reduce cellular growth, survival, and migration, making NFX1-123 a potential novel therapeutic target.

A mathematical model of SARS-CoV-2 immunity predicts paxlovid rebound.

Nirmatrelvir/ritonavir (Paxlovid), an oral antiviral medication targeting SARS-CoV-2, remains an important treatment for COVID-19. Initial studies of nirmatrelvir/ritonavir were performed in SARS-CoV-2 unvaccinated patients without prior confirmed SARS-CoV-2 infection; however, most individuals have now either been vaccinated and/or have experienced SARS-CoV-2 infection. After nirmatrelvir/ritonavir became widely available, reports surfaced of "Paxlovid rebound," a phenomenon in which symptoms (and SARS-CoV-2 test positivity) would initially resolve, but after finishing treatment, symptoms and test positivity would return. We used a previously described parsimonious mathematical model of immunity to SARS-CoV-2 infection to model the effect of nirmatrelvir/ritonavir treatment in unvaccinated and vaccinated patients. Model simulations show that viral rebound after treatment occurs only in vaccinated patients, while unvaccinated (SARS-COV-2 naïve) patients treated with nirmatrelvir/ritonavir do not experience any rebound in viral load. This work suggests that an approach combining parsimonious models of the immune system could be used to gain important insights in the context of emerging pathogens.

Investigating the Orientation of an Interfacially Adsorbed Monoclonal Antibody and Its Fragments Using Neutron Reflection.

Interfacial adsorption is a molecular process occurring during the production, purification, transport, and storage of antibodies, with a direct impact on their structural stability and subsequent implications on their bioactivities. While the average conformational orientation of an adsorbed protein can be readily determined, its associated structures are more complex to characterize. Neutron reflection has been used in this work to investigate the conformational orientations of the monoclonal antibody COE-3 and its Fab and Fc fragments at the oil/water and air/water interfaces. Rigid body rotation modeling was found to be suitable for globular and relatively rigid proteins such as the Fab and Fc fragments but less so for relatively flexible proteins such as full COE-3. Fab and Fc fragments adopted a 'flat-on' orientation at the air/water interface, minimizing the thickness of the protein layer, but they adopted a substantially tilted orientation at the oil/water interface with increased layer thickness. In contrast, COE-3 was found to adsorb in tilted orientations at both interfaces, with one fragment protruding into the solution. This work demonstrates that rigid-body modeling can provide additional insights into protein layers at various interfaces relevant to bioprocess engineering.