Generative artificial intelligence performs rudimentary structural biology modeling.

Natural language-based generative artificial intelligence (AI) has become increasingly prevalent in scientific research. Intriguingly, capabilities of generative pre-trained transformer (GPT) language models beyond the scope of natural language tasks have recently been identified. Here we explored how GPT-4 might be able to perform rudimentary structural biology modeling. We prompted GPT-4 to model 3D structures for the 20 standard amino acids and an α-helical polypeptide chain, with the latter incorporating Wolfram mathematical computation. We also used GPT-4 to perform structural interaction analysis between nirmatrelvir and its target, the SARS-CoV-2 main protease. Geometric parameters of the generated structures typically approximated close to experimental references. However, modeling was sporadically error-prone and molecular complexity was not well tolerated. Interaction analysis further revealed the ability of GPT-4 to identify specific amino acid residues involved in ligand binding along with corresponding bond distances. Despite current limitations, we show the capacity of natural language generative AI to perform basic structural biology modeling and interaction analysis with atomic-scale accuracy.

Research involving the recently deceased: ethics questions that must be answered.

Research involving recently deceased humans that are physiologically maintained following declaration of death by neurologic criteria-or 'research involving the recently deceased'-can fill a translational research gap while reducing harm to animals and living human subjects. It also creates new challenges for honouring the donor's legacy, respecting the rights of donor loved ones, resource allocation and public health. As this research model gains traction, new empirical ethics questions must be answered to preserve public trust in all forms of tissue donation and in the practice of medicine while respecting the legacy of the deceased and the rights of donor loved ones. This article suggests several topics for immediate investigation to understand the attitudes and experiences of researchers, clinical collaborators, donor loved ones and the public to ensure research involving the recently deceased advances ethically.

Ceramide as an endothelial cell surface receptor and a lung-specific lipid vascular target for circulating ligands.

The vascular endothelium from individual organs is functionally specialized, and it displays a unique set of accessible molecular targets. These serve as endothelial cell receptors to affinity ligands. To date, all identified vascular receptors have been proteins. Here, we show that an endothelial lung-homing peptide (CGSPGWVRC) interacts with C16-ceramide, a bioactive sphingolipid that mediates several biological functions. Upon binding to cell surfaces, CGSPGWVRC triggers ceramide-rich platform formation, activates acid sphingomyelinase and ceramide production, without the associated downstream apoptotic signaling. We also show that the lung selectivity of CGSPGWVRC homing peptide is dependent on ceramide production in vivo. Finally, we demonstrate two potential applications for this lipid vascular targeting system: i) as a bioinorganic hydrogel for pulmonary imaging and ii) as a ligand-directed lung immunization tool against COVID-19. Thus, C16-ceramide is a unique example of a lipid-based receptor system in the lung vascular endothelium targeted in vivo by circulating ligands such as CGSPGWVRC.

Integrated metabolic and genetic analysis reveals distinct features of human differentiated thyroid cancer.

BACKGROUND: Differentiated thyroid cancer (DTC) affects thousands of lives worldwide each year. Typically, DTC is a treatable disease with a good prognosis. Yet, some patients are subjected to partial or total thyroidectomy and radioiodine therapy to prevent local disease recurrence and metastasis. Unfortunately, thyroidectomy and/or radioiodine therapy often worsen(s) quality of life and might be unnecessary in indolent DTC cases. On the other hand, the lack of biomarkers indicating a potential metastatic thyroid cancer imposes an additional challenge to managing and treating patients with this disease. AIM: The presented clinical setting highlights the unmet need for a precise molecular diagnosis of DTC and potential metastatic disease, which should dictate appropriate therapy. MATERIALS AND METHODS: In this article, we present a differential multi-omics model approach, including metabolomics, genomics, and bioinformatic models, to distinguish normal glands from thyroid tumours. Additionally, we are proposing biomarkers that could indicate potential metastatic diseases in papillary thyroid cancer (PTC), a sub-class of DTC. RESULTS: Normal and tumour thyroid tissue from DTC patients had a distinct yet well-defined metabolic profile with high levels of anabolic metabolites and/or other metabolites associated with the energy maintenance of tumour cells. The consistency of the DTC metabolic profile allowed us to build a bioinformatic classification model capable of clearly distinguishing normal from tumor thyroid tissues, which might help diagnose thyroid cancer. Moreover, based on PTC patient samples, our data suggest that elevated nuclear and mitochondrial DNA mutational burden, intra-tumour heterogeneity, shortened telomere length, and altered metabolic profile reflect the potential for metastatic disease. DISCUSSION: Altogether, this work indicates that a differential and integrated multi-omics approach might improve DTC management, perhaps preventing unnecessary thyroid gland removal and/or radioiodine therapy. CONCLUSIONS: Well-designed, prospective translational clinical trials will ultimately show the value of this integrated multi-omics approach and early diagnosis of DTC and potential metastatic PTC.

Preclinical optimization of a diode laser-based clamp-free partial nephrectomy in a large animal model.

Kidney cancer is a common urologic malignancy with either laparoscopic (LPN) or robotic partial nephrectomy as therapeutic options of choice for localized tumors. However, renal resection and suturing are challenging steps of the procedure that can lead to complications such as prolonged warm ischemia, bleeding, and urinary fistulas. LPN with a diode laser is an efficient technique due to its cutting and/or coagulation attributes. Surprisingly, key laser features such as wavelength and power remain undefined. Using a large porcine model, we evaluated the laser range of wavelength and power in a clamp-free LPN and compared it to the established gold-standard LPN technique (i.e., cold-cutting and suturing). By analyzing surgery duration, bleeding, presence of urine leak, tissue damage related to the resected renal fragment and the remaining organ, hemoglobin levels, and renal function, we show that an optimized experimental diode laser clamp-free LPN (wavelength, 980 nm; power, 15 W) had shorter surgery time with less bleeding, and better postoperative renal function recovery when compared to the well-established technique. Together, our data indicate that partial nephrectomy with a diode laser clamp-free LPN technique is an improved alternative to the gold-standard technique. Therefore, translational clinical trials towards human patient applications are readily feasible.

Modeling of ACE2 and antibodies bound to SARS-CoV-2 provides insights into infectivity and immune evasion.

Given the COVID-19 pandemic, there is interest in understanding ligand-receptor features and targeted antibody-binding attributes against emerging SARS-CoV-2 variants. Here, we developed a large-scale structure-based pipeline for analysis of protein-protein interactions regulating SARS-CoV-2 immune evasion. First, we generated computed structural models of the Spike protein of 3 SARS-CoV-2 variants (B.1.1.529, BA.2.12.1, and BA.5) bound either to a native receptor (ACE2) or to a large panel of targeted ligands (n = 282), which included neutralizing or therapeutic monoclonal antibodies. Moreover, by using the Barnes classification, we noted an overall loss of interfacial interactions (with gain of new interactions in certain cases) at the receptor-binding domain (RBD) mediated by substituted residues for neutralizing complexes in classes 1 and 2, whereas less destabilization was observed for classes 3 and 4. Finally, an experimental validation of predicted weakened therapeutic antibody binding was performed in a cell-based assay. Compared with the original Omicron variant (B.1.1.529), derivative variants featured progressive destabilization of antibody-RBD interfaces mediated by a larger set of substituted residues, thereby providing a molecular basis for immune evasion. This approach and findings provide a framework for rapidly and efficiently generating structural models for SARS-CoV-2 variants bound to ligands of mechanistic and therapeutic value.

A modeling-based approach to optimize COVID-19 vaccine dosing schedules for improved protection.

While the development of different vaccines slowed the dissemination of SARS-CoV-2, the occurrence of breakthrough infections has continued to fuel the COVID-19 pandemic. To secure at least partial protection in the majority of the population through 1 dose of a COVID-19 vaccine, delayed administration of boosters has been implemented in many countries. However, waning immunity and emergence of new variants of SARS-CoV-2 suggest that such measures may induce breakthrough infections due to intermittent lapses in protection. Optimizing vaccine dosing schedules to ensure prolonged continuity in protection could thus help control the pandemic. We developed a mechanistic model of immune response to vaccines as an in silico tool for dosing schedule optimization. The model was calibrated with clinical data sets of acquired immunity to COVID-19 mRNA vaccines in healthy and immunocompromised participants and showed robust validation by accurately predicting neutralizing antibody kinetics in response to multiple doses of COVID-19 mRNA vaccines. Importantly, by estimating population vulnerability to breakthrough infections, we predicted tailored vaccination dosing schedules to minimize breakthrough infections, especially for immunocompromised individuals. We identified that the optimal vaccination schedules vary from CDC-recommended dosing, suggesting that the model is a valuable tool to optimize vaccine efficacy outcomes during future outbreaks.

Integrated metabolic and genetic analysis reveals distinct features of primary differentiated thyroid cancer and its metastatic potential in humans.

Differentiated thyroid cancer (DTC) affects thousands of lives worldwide every year. Typically, DTC is a treatable disease with a good prognosis. Yet, some patients are subjected to partial or total thyroidectomy and radioiodine therapy to prevent local disease recurrence and metastasis. Unfortunately, thyroidectomy and/or radioiodine therapy often worsen(s) the quality of life and might be unnecessary in indolent DTC cases. This clinical setting highlights the unmet need for a precise molecular diagnosis of DTC, which should dictate appropriate therapy. Here we propose a differential multi-omics model approach to distinguish normal gland from thyroid tumor and to indicate potential metastatic diseases in papillary thyroid cancer (PTC), a sub-class of DTC. Based on PTC patient samples, our data suggest that elevated nuclear and mitochondrial DNA mutational burden, intratumor heterogeneity, shortened telomere length, and altered metabolic profile reflect the potential for metastatic disease. Specifically, normal and tumor thyroid tissues from these patients had a distinct yet well-defined metabolic profile with high levels of anabolic metabolites and/or other metabolites associated with the energy maintenance of tumor cells. Altogether, this work indicates that a differential and integrated multi-omics approach might improve DTC management, perhaps preventing unnecessary thyroid gland removal and/or radioiodine therapy. Well-designed, prospective translational clinical trials will ultimately show the value of this targeted molecular approach. TRANSLATIONAL RELEVANCE: In this article, we propose a new integrated metabolic, genomic, and cytopathologic methods to diagnose Differentiated Thyroid Cancer when the conventional methods failed. Moreover, we suggest metabolic and genomic markers to help predict high-risk Papillary Thyroid Cancer. Both might be important tools to avoid unnecessary surgery and/or radioiodine therapy that can worsen the quality of life of the patients more than living with an indolent Thyroid nodule.

Dysregulation of the PRUNE2/PCA3 genetic axis in human prostate cancer: from experimental discovery to validation in two independent patient cohorts.

Background: We have previously shown that the long non-coding (lnc)RNA prostate cancer associated 3 (PCA3; formerly prostate cancer antigen 3) functions as a trans-dominant negative oncogene by targeting the previously unrecognized prostate cancer suppressor gene PRUNE2 (a homolog of the Drosophila prune gene), thereby forming a functional unit within a unique allelic locus in human cells. Here, we investigated the PCA3/PRUNE2 regulatory axis from early (tumorigenic) to late (biochemical recurrence) genetic events during human prostate cancer progression. Methods: The reciprocal PCA3 and PRUNE2 gene expression relationship in paired prostate cancer and adjacent normal prostate was analyzed in two independent retrospective cohorts of clinically annotated cases post-radical prostatectomy: a single-institutional discovery cohort (n=107) and a multi-institutional validation cohort (n=497). We compared the tumor gene expression of PCA3 and PRUNE2 to their corresponding expression in the normal prostate. We also serially examined clinical/pathological variables including time to disease recurrence. Results: We consistently observed increased expression of PCA3 and decreased expression of PRUNE2 in prostate cancer compared with the adjacent normal prostate across all tumor grades and stages. However, there was no association between the relative gene expression levels of PCA3 or PRUNE2 and time to disease recurrence, independent of tumor grades and stages. Conclusions: We concluded that upregulation of the lncRNA PCA3 and targeted downregulation of the protein-coding PRUNE2 gene in prostate cancer could be early (rather than late) molecular events in the progression of human prostate tumorigenesis but are not associated with biochemical recurrence. Further studies of PCA3/PRUNE2 dysregulation are warranted. Funding: We received support from the Human Tissue Repository and Tissue Analysis Shared Resource from the Department of Pathology of the University of New Mexico School of Medicine and a pilot award from the University of New Mexico Comprehensive Cancer Center. RP and WA were supported by awards from the Levy-Longenbaugh Donor-Advised Fund and the Prostate Cancer Foundation. EDN reports research fellowship support from the Brazilian National Council for Scientific and Technological Development (CNPq), Brazil, and the Associação Beneficente Alzira Denise Hertzog Silva (ABADHS), Brazil. This work has been funded in part by the NCI Cancer Center Support Grants (CCSG; P30) to the University of New Mexico Comprehensive Cancer Center (CA118100) and the Rutgers Cancer Institute of New Jersey (CA072720).

A modeling-based approach to optimize COVID-19 vaccine dosing schedules for improved protection.

While the development of different vaccines has slowed the dissemination of SARS-CoV-2, the occurrence of breakthrough infections continues to fuel the pandemic. As a strategy to secure at least partial protection, with a single dose of a given COVID-19 vaccine to maximum possible fraction of the population, delayed administration of subsequent doses (or boosters) has been implemented in many countries. However, waning immunity and emergence of new variants of SARS-CoV-2 suggest that such measures may jeopardize the attainment of herd immunity due to intermittent lapses in protection. Optimizing vaccine dosing schedules could thus make the difference between periodic occurrence of breakthrough infections or effective control of the pandemic. To this end, we have developed a mechanistic mathematical model of adaptive immune response to vaccines and demonstrated its applicability to COVID-19 mRNA vaccines as a proof-of-concept for future outbreaks. The model was thoroughly calibrated against multiple clinical datasets involving immune response to SARS-CoV-2 infection and mRNA vaccines in healthy and immunocompromised subjects (cancer patients undergoing therapy); the model showed robust clinical validation by accurately predicting neutralizing antibody kinetics, a correlate of vaccine-induced protection, in response to multiple doses of mRNA vaccines. Importantly, we estimated population vulnerability to breakthrough infections and predicted tailored vaccination dosing schedules to maximize protection and thus minimize breakthrough infections, based on the immune status of a sub-population. We have identified a critical waiting window for cancer patients (or, immunocompromised subjects) to allow recovery of the immune system (particularly CD4+ T-cells) for effective differentiation of B-cells to produce neutralizing antibodies and thus achieve optimal vaccine efficacy against variants of concern, especially between the first and second doses. Also, we have obtained optimized dosing schedules for subsequent doses in healthy and immunocompromised subjects, which vary from the CDC-recommended schedules, to minimize breakthrough infections. The developed modeling tool is based on generalized adaptive immune response to antigens and can thus be leveraged to guide vaccine dosing schedules during future outbreaks.

Leveraging the scientific findings to develop therapeutic strategies for dormant breast cancer cells.

Breast cancer (BC) metastasis can occur decades before clinical diagnosis. During this time, the cancer cells (BCCs) can remain dormant for decades. This type of dormancy also occurs during remission where the dormant BCCs adapt cycling quiescence within the tissue microenvironment. BC shows preference for the bone marrow (BM), resulting in poor prognosis. The BM provides a challenge due to the complex niche between the peripheral interface and endosteum. The process of dormancy begins upon entry into the marrow with the changes facilitated through crosstalk between the cancer cells and tissue niche. More importantly, dormancy can occur at any time during the disease process, including the time during treatment. This perspective discusses the challenges posed by the marrow microenvironment to develop treatment. The article discusses the complex mechanisms at each compartment within the marrow niche and the added negative issue of toxicity to the endogenous stem cells.

Genomic landscape of lymphatic malformations: a case series and response to the PI3Kα inhibitor alpelisib in an N-of-1 clinical trial.

Background: Lymphatic malformations (LMs) often pose treatment challenges due to a large size or a critical location that could lead to disfigurement, and there are no standardized treatment approaches for either refractory or unresectable cases. Methods: We examined the genomic landscape of a patient cohort of LMs (n = 30 cases) that underwent comprehensive genomic profiling using a large-panel next-generation sequencing assay. Immunohistochemical analyses were completed in parallel. Results: These LMs had low mutational burden with hotspot PIK3CA mutations (n = 20) and NRAS (n = 5) mutations being most frequent, and mutually exclusive. All LM cases with Kaposi sarcoma-like (kaposiform) histology had NRAS mutations. One index patient presented with subacute abdominal pain and was diagnosed with a large retroperitoneal LM harboring a somatic PIK3CA gain-of-function mutation (H1047R). The patient achieved a rapid and durable radiologic complete response, as defined in RECIST1.1, to the PI3Kα inhibitor alpelisib within the context of a personalized N-of-1 clinical trial (NCT03941782). In translational correlative studies, canonical PI3Kα pathway activation was confirmed by immunohistochemistry and human LM-derived lymphatic endothelial cells carrying an allele with an activating mutation at the same locus were sensitive to alpelisib treatment in vitro, which was demonstrated by a concentration-dependent drop in measurable impedance, an assessment of cell status. Conclusions: Our findings establish that LM patients with conventional or kaposiform histology have distinct, yet targetable, driver mutations. Funding: R.P. and W.A. are supported by awards from the Levy-Longenbaugh Fund. S.G. is supported by awards from the Hugs for Brady Foundation. This work has been funded in part by the NCI Cancer Center Support Grants (CCSG; P30) to the University of Arizona Cancer Center (CA023074), the University of New Mexico Comprehensive Cancer Center (CA118100), and the Rutgers Cancer Institute of New Jersey (CA072720). B.K.M. was supported by National Science Foundation via Graduate Research Fellowship DGE-1143953. Clinical trial number: NCT03941782.

The epigenetic reader PHF21B modulates murine social memory and synaptic plasticity-related genes.

Synaptic dysfunction is a manifestation of several neurobehavioral and neurological disorders. A major therapeutic challenge lies in uncovering the upstream regulatory factors controlling synaptic processes. Plant homeodomain (PHD) finger proteins are epigenetic readers whose dysfunctions are implicated in neurological disorders. However, the molecular mechanisms linking PHD protein deficits to disease remain unclear. Here, we generated a PHD finger protein 21B-depleted (Phf21b-depleted) mutant CRISPR mouse model (hereafter called Phf21bΔ4/Δ4) to examine Phf21b's roles in the brain. Phf21bΔ4/Δ4 animals exhibited impaired social memory. In addition, reduced expression of synaptic proteins and impaired long-term potentiation were observed in the Phf21bΔ4/Δ4 hippocampi. Transcriptome profiling revealed differential expression of genes involved in synaptic plasticity processes. Furthermore, we characterized a potentially novel interaction of PHF21B with histone H3 trimethylated lysine 36 (H3K36me3), a histone modification associated with transcriptional activation, and the transcriptional factor CREB. These results establish PHF21B as an important upstream regulator of synaptic plasticity-related genes and a candidate therapeutic target for neurobehavioral dysfunction in mice, with potential applications in human neurological and psychiatric disorders.

Genetic and Structural Analysis of SARS-CoV-2 Spike Protein for Universal Epitope Selection.

Evaluation of immunogenic epitopes for universal vaccine development in the face of ongoing SARS-CoV-2 evolution remains a challenge. Herein, we investigate the genetic and structural conservation of an immunogenically relevant epitope (C662-C671) of spike (S) protein across SARS-CoV-2 variants to determine its potential utility as a broad-spectrum vaccine candidate against coronavirus diseases. Comparative sequence analysis, structural assessment, and molecular dynamics simulations of C662-C671 epitope were performed. Mathematical tools were employed to determine its mutational cost. We found that the amino acid sequence of C662-C671 epitope is entirely conserved across the observed major variants of SARS-CoV-2 in addition to SARS-CoV. Its conformation and accessibility are predicted to be conserved, even in the highly mutated Omicron variant. Costly mutational rate in the context of energy expenditure in genome replication and translation can explain this strict conservation. These observations may herald an approach to developing vaccine candidates for universal protection against emergent variants of coronavirus.

Translational Modeling Identifies Synergy between Nanoparticle-Delivered miRNA-22 and Standard-of-Care Drugs in Triple-Negative Breast Cancer.

PURPOSE: Downregulation of miRNA-22 in triple-negative breast cancer (TNBC) is associated with upregulation of eukaryotic elongation 2 factor kinase (eEF2K) protein, which regulates tumor growth, chemoresistance, and tumor immunosurveillance. Moreover, exogenous administration of miRNA-22, loaded in nanoparticles to prevent degradation and improve tumor delivery (termed miRNA-22 nanotherapy), to suppress eEF2K production has shown potential as an investigational therapeutic agent in vivo. METHODS: To evaluate the translational potential of miRNA-22 nanotherapy, we developed a multiscale mechanistic model, calibrated to published in vivo data and extrapolated to the human scale, to describe and quantify the pharmacokinetics and pharmacodynamics of miRNA-22 in virtual patient populations. RESULTS: Our analysis revealed the dose-response relationship, suggested optimal treatment frequency for miRNA-22 nanotherapy, and highlighted key determinants of therapy response, from which combination with immune checkpoint inhibitors was identified as a candidate strategy for improving treatment outcomes. More importantly, drug synergy was identified between miRNA-22 and standard-of-care drugs against TNBC, providing a basis for rational therapeutic combinations for improved response CONCLUSIONS: The present study highlights the translational potential of miRNA-22 nanotherapy for TNBC in combination with standard-of-care drugs.

Spatial and temporal dynamics of HDACs class IIa following mild traumatic brain injury in adult rats.

The fundamental role of epigenetic regulatory mechanisms involved in neuroplasticity and adaptive responses to traumatic brain injury (TBI) is gaining increased recognition. TBI-induced neurodegeneration is associated with several changes in the expression-activity of various epigenetic regulatory enzymes, including histone deacetylases (HDACs). In this study, PET/CT with 6-([18F]trifluoroacetamido)-1- hexanoicanilide ([18F]TFAHA) to image spatial and temporal dynamics of HDACs class IIa expression-activity in brains of adult rats subjected to a weight drop model of diffuse, non-penetrating, mild traumatic brain injury (mTBI). The mTBI model was validated by histopathological and immunohistochemical analyses of brain tissue sections for localization and magnitude of expression of heat-shock protein-70 kDa (HSP70), amyloid precursor protein (APP), cannabinoid receptor-2 (CB2), ionized calcium-binding adapter protein-1 (IBA1), histone deacetylase-4 and -5 (HDAC4 and HDAC5). In comparison to baseline, the expression-activities of HDAC4 and HDAC5 were downregulated in the hippocampus, nucleus accumbens, peri-3rd ventricular part of the thalamus, and substantia nigra at 1-3 days post mTBI, and remained low at 7-8 days post mTBI. Reduced levels of HDAC4 and HDAC5 expression observed in neurons of these brain regions post mTBI were associated with the reduced nuclear and neuropil levels of HDAC4 and HDAC5 with the shift to perinuclear localization of these enzymes. These results support the rationale for the development of therapeutic strategies to upregulate expression-activity of HDACs class IIa post-TBI. PET/CT (MRI) with [18F]TFAHA can facilitate the development and clinical translation of unique therapeutic approaches to upregulate the expression and activity of HDACs class IIa enzymes in the brain after TBI.

Mathematical modeling of cancer immunotherapy for personalized clinical translation.

Encouraging advances are being made in cancer immunotherapy modeling, especially in the key areas of developing personalized treatment strategies based on individual patient parameters, predicting treatment outcomes and optimizing immunotherapy synergy when used in combination with other treatment approaches. Here we present a focused review of the most recent mathematical modeling work on cancer immunotherapy with a focus on clinical translatability. It can be seen that this field is transitioning from pure basic science to applications that can make impactful differences in patients' lives. We discuss how researchers are integrating experimental and clinical data to fully inform models so that they can be applied for clinical predictions, and present the challenges that remain to be overcome if widespread clinical adaptation is to be realized. Lastly, we discuss the most promising future applications and areas that are expected to be the focus of extensive upcoming modeling studies.

Structural models of SARS-CoV-2 Omicron variant in complex with ACE2 receptor or antibodies suggest altered binding interfaces.

There is enormous ongoing interest in characterizing the binding properties of the SARS-CoV-2 Omicron Variant of Concern (VOC) (B.1.1.529), which continues to spread towards potential dominance worldwide. To aid these studies, based on the wealth of available structural information about several SARS-CoV-2 variants in the Protein Data Bank (PDB) and a modeling pipeline we have previously developed for tracking the ongoing global evolution of SARS-CoV-2 proteins, we provide a set of computed structural models (henceforth models) of the Omicron VOC receptor-binding domain (omRBD) bound to its corresponding receptor Angiotensin-Converting Enzyme (ACE2) and a variety of therapeutic entities, including neutralizing and therapeutic antibodies targeting previously-detected viral strains. We generated bound omRBD models using both experimentally-determined structures in the PDB as well as machine learningbased structure predictions as starting points. Examination of ACE2-bound omRBD models reveals an interdigitated multi-residue interaction network formed by omRBD-specific substituted residues (R493, S496, Y501, R498) and ACE2 residues at the interface, which was not present in the original Wuhan-Hu-1 RBD-ACE2 complex. Emergence of this interaction network suggests optimization of a key region of the binding interface, and positive cooperativity among various sites of residue substitutions in omRBD mediating ACE2 binding. Examination of neutralizing antibody complexes for Barnes Class 1 and Class 2 antibodies modeled with omRBD highlights an overall loss of interfacial interactions (with gain of new interactions in rare cases) mediated by substituted residues. Many of these substitutions have previously been found to independently dampen or even ablate antibody binding, and perhaps mediate antibody-mediated neutralization escape ( e.g ., K417N). We observe little compensation of corresponding interaction loss at interfaces when potential escape substitutions occur in combination. A few selected antibodies ( e.g ., Barnes Class 3 S309), however, feature largely unaltered or modestly affected protein-protein interfaces. While we stress that only qualitative insights can be obtained directly from our models at this time, we anticipate that they can provide starting points for more detailed and quantitative computational characterization, and, if needed, redesign of monoclonal antibodies for targeting the Omicron VOC Spike protein. In the broader context, the computational pipeline we developed provides a framework for rapidly and efficiently generating retrospective and prospective models for other novel variants of SARS-CoV-2 bound to entities of virological and therapeutic interest, in the setting of a global pandemic.