Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “binding”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Adsorption of Ag, Au, Cu, and Ni on MoS 2 : theory and experiment

Abstract Here, we present results of a computational and experimental study of adsorption of various metals on MoS2. In particular, we analyzed the binding mechanism of four metallic elements (Ag, Au, Cu, Ni) on MoS2. Among these elements, Ni exhibits the strongest binding and lowest mobility on the surface of MoS2. On the other hand, Au and Ag bond very weakly to the surface and have very high mobilities. Our calculations for Cu show that its bonding and surface mobility are between these two groups. Experimentally, Ni films exhibit a composition characterized by randomly oriented nanoscale clusters. This is consistent with the larger cohesive energy of Ni atoms as compared with their binding energy with MoS2, which is expected to result in 3D clusters. In contrast, Au and Ag tend to form atomically flat plateaued structures on MoS2, which is contrary to their larger cohesive energy as compared to their weak binding with MoS2. Cu displays a surface morphology somewhat similar to Ni, featuring larger nanoscale clusters. However, unlike Ni, in many cases Cu exhibits small plateaued surfaces on these clusters. This suggests that Cu likely has two competing mechanisms that cause it to span the behaviors seen in the Ni and Au/Ag film morphologies. These results indicate that calculations of the initial binding conditions could be useful for predicting film morphologies. In addition, out calculations show that the adsorption of adatoms with odd electron number like Ag, Au, and Cu results in 100% spin-polarization and integer magnetic moment of the system. Adsorption of Ni adatoms, with even electron number, does not induce a magnetic transition.

Harms, Haley↗
DOE OSTI2020

TBDB: a database of structurally annotated T-box riboswitch:tRNA pairs

Abstract T-box riboswitches constitute a large family of tRNA-binding leader sequences that play a central role in gene regulation in many gram-positive bacteria. Accurate inference of the tRNA binding to T-box riboswitches is critical to predict their cis-regulatory activity. However, there is no central repository of information on the tRNA binding specificities of T-box riboswitches, and de novo prediction of binding specificities requires advanced knowledge of computational tools to annotate riboswitch secondary structure features. Here, we present the T-box Riboswitch Annotation Database (TBDB, https://tbdb.io), an open-access database with a collection of 23,535 T-box riboswitch sequences, spanning the major phyla of 3,632 bacterial species. Among structural predictions, the TBDB also identifies specifier sequences, cognate tRNA binding partners, and downstream regulatory targets. To our knowledge, the TBDB presents the largest collection of feature, sequence, and structural annotations carried out on this important family of regulatory RNA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structures of CTCF–DNA complexes including all 11 zinc fingers

The CCCTC-binding factor (CTCF) binds tens of thousands of enhancers and promoters on mammalian chromosomes by means of its 11 tandem zinc finger (ZF) DNA-binding domain. In addition to the 12–15-bp CORE sequence, some of the CTCF binding sites contain 5' upstream and/or 3' downstream motifs. Here, we describe two structures for overlapping portions of human CTCF, respectively, including ZF1–ZF7 and ZF3–ZF11 in complex with DNA that incorporates the CORE sequence together with either 3' downstream or 5' upstream motifs. Like conventional tandem ZF array proteins, ZF1–ZF7 follow the right-handed twist of the DNA, with each finger occupying and recognizing one triplet of three base pairs in the DNA major groove. ZF8 plays a unique role, acting as a spacer across the DNA minor groove and positioning ZF9–ZF11 to make cross-strand contacts with DNA. We ascribe the difference between the two subgroups of ZF1–ZF7 and ZF8–ZF11 to residues at the two positions -6 and -5 within each finger, with small residues for ZF1–ZF7 and bulkier and polar/charged residues for ZF8–ZF11. ZF8 is also uniquely rich in basic amino acids, which allows salt bridges to DNA phosphates in the minor groove. Highly specific arginine–guanine and glutamine–adenine interactions, used to recognize G:C or A:T base pairs at conventional base-interacting positions of ZFs, also apply to the cross-strand interactions adopted by ZF9–ZF11. The differences between ZF1–ZF7 and ZF8–ZF11 can be rationalized structurally and may contribute to recognition of high-affinity CTCF binding sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural basis for cloverleaf RNA-initiated viral genome replication

The genomes of positive-strand RNA viruses serve as a template for both protein translation and genome replication. In enteroviruses, a cloverleaf RNA structure at the 5' end of the genome functions as a switch to transition from viral translation to replication by interacting with host poly(C)-binding protein 2 (PCBP2) and the viral 3CD pro protein. We determined the structures of cloverleaf RNA from coxsackievirus and poliovirus. Cloverleaf RNA folds into an H-type four-way junction and is stabilized by a unique adenosine-cytidine-uridine (A•C-U) base triple involving the conserved pyrimidine mismatch region. The two PCBP2 binding sites are spatially proximal and are located on the opposite end from the 3CD pro binding site on cloverleaf. We determined that the A•C-U base triple restricts the flexibility of the cloverleaf stem–loops resulting in partial occlusion of the PCBP2 binding site, and elimination of the A•C-U base triple increases the binding affinity of PCBP2 to the cloverleaf RNA. Based on the cloverleaf structures and biophysical assays, we propose a new mechanistic model by which enteroviruses use the cloverleaf structure as a molecular switch to transition from viral protein translation to genome replication.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Structural and functional insights into the interaction between the bacteriophage T4 DNA processing proteins gp32 and Dda

Abstract Bacteriophage T4 is a classic model system for studying the mechanisms of DNA processing. A key protein in T4 DNA processing is the gp32 single-stranded DNA-binding protein. gp32 has two key functions: it binds cooperatively to single-stranded DNA (ssDNA) to protect it from nucleases and remove regions of secondary structure, and it recruits proteins to initiate DNA processes including replication and repair. Dda is a T4 helicase recruited by gp32, and we purified and crystallized a gp32–Dda–ssDNA complex. The low-resolution structure revealed how the C-terminus of gp32 engages Dda. Analytical ultracentrifugation analyses were consistent with the crystal structure. An optimal Dda binding peptide from the gp32 C-terminus was identified using surface plasmon resonance. The crystal structure of the Dda–peptide complex was consistent with the corresponding interaction in the gp32–Dda–ssDNA structure. A Dda-dependent DNA unwinding assay supported the structural conclusions and confirmed that the bound gp32 sequesters the ssDNA generated by Dda. The structure of the gp32–Dda–ssDNA complex, together with the known structure of the gp32 body, reveals the entire ssDNA binding surface of gp32. gp32–Dda–ssDNA complexes in the crystal are connected by the N-terminal region of one gp32 binding to an adjacent gp32, and this provides key insights into this interaction.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Mechanistic studies of small molecule ligands selective to RNA single G bulges

Abstract Small-molecule RNA binders have emerged as an important pharmacological modality. A profound understanding of the ligand selectivity, binding mode, and influential factors governing ligand engagement with RNA targets is the foundation for rational ligand design. Here, we report a novel class of coumarin derivatives exhibiting selective binding affinity towards single G RNA bulges. Harnessing the computational power of all-atom Gaussian accelerated molecular dynamics simulations, we unveiled a rare minor groove binding mode of the ligand with a key interaction between the coumarin moiety and the G bulge. This predicted binding mode is consistent with results obtained from structure-activity relationship studies and transverse relaxation measurements by nuclear magnetic resonance spectroscopy. We further generated 444 molecular descriptors from 69 coumarin derivatives and identified key contributors to the binding events, such as charge state and planarity, by lasso (least absolute shrinkage and selection operator) regression. Our work deepened the understanding of RNA-small molecule interactions and integrated a new framework for the rational design of selective small-molecule RNA binders.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Stereoselective recognition of morphine enantiomers by μ -opioid receptor

Stereospecific recognition of chiral molecules plays a crucial role in biological systems. The μ-opioid receptor (MOR) exhibits binding affinity towards (-)-morphine, a well-established gold standard in pain management, while it shows minimal binding affinity for the (+)-morphine enantiomer, resulting in a lack of analgesic activity. Understanding how MOR stereoselectively recognizes morphine enantiomers has remained a puzzle in neuroscience and pharmacology for over half-a-century due to the lack of direct observation techniques. To unravel this mystery, we constructed the binding and unbinding processes of morphine enantiomers with MOR via molecular dynamics simulations to investigate the thermodynamics and kinetics governing MOR's stereoselective recognition of morphine enantiomers. Our findings reveal that the binding of (-)-morphine stabilizes MOR in its activated state, exhibiting a deep energy well and a prolonged residence time. In contrast, (+)-morphine fails to sustain the activation state of MOR. Furthermore, the results suggest that specific residues, namely D114 2.50 and D147 3.32 , are deprotonated in the active state of MOR bound to (-)-morphine. This work highlights that the selectivity in molecular recognition goes beyond binding affinities, extending into the realm of residence time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Density functionals combined with van der Waals corrections for graphene adsorbed on layered materials

Standard density functionals like Perdew-Burke-Ernzerhof (PBE) or Strongly Constrained and Appropriately Normed (SCAN) need a correction to account for long-range van der Waals (vdW) interaction. The damped Zaremba-Kohn model (dZK) starts from a formula for the vdW interaction of a distant atom with a solid surface, both with known dielectric properties, damps this formula at short range, and then treats an adsorbed molecule or atomic layer as a collection of renormalized atoms. Here, we extend this model to graphene adsorbed on semiconducting layered materials [bulk graphite and hBN, and multilayer transition metal dichalcogenides (TMD)] by including the C 4 asymptotic term and multiple electrostatic image effects due to the two surfaces of the substrate slabs in the vdW calculations. The resulting SCAN-vdW-dZK and PBE-vdW-dZK give approximately the same results for the systems considered here, in agreement with available reference values. The predicted binding energies are roughly 25% lower than those from SCAN+rVV10 (revised Vydrov and Van Voorhis 2010), and ~15% lower than those from PBE+rVV10. Since SCAN+rVV10 usually overbinds, the predicted binding energies by SCAN-vdW-dZK and PBE-vdW-dZK are expected to be closer to the true values. The predicted equilibrium binding distances from SCAN-vdW-dZK and PBE-vdW-dZK are slightly larger (~ 0.1 Å) than those from SCAN+rVV10, and close to those from SCAN. The binding energy depends upon the number of substrate layers more strongly in vdW-dZK than in rVV10. The C 4 -term contributions can be 40% of the total vdW interactions, and the C 5 term contributes about 10%. The effects of images and the back surfaces of slabs can contribute about 4–10%. The vdW interaction energy power laws from the vdW-dZK model for graphene adsorbed on multilayer MoS 2 show slowly varying decay to the pairwise exponent –4 with increasing separation $\textit{D}$, very similar to those obtained from the random-phase approximation and renormalization group approaches by Ambrosetti et al. Both the PBE-vdW-dZK and SCAN-vdW-dZK give a greater increase in interlayer binding energy when the TMD substrate changes from monolayer to four-layer for the graphene/TMD adsorption systems. This is consistent with the relevant electron-energy-loss spectroscopy experimental results showing increased dielectric response from the substrate with increasing substrate layer number.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Biexcitons are bound in CsPbBr 3 perovskite nanocrystals

Here, we study the energetics of quasiparticle excitations in CsPbBr 3 perovskite nanocrystals using path integral molecular dynamics simulations. Employing detailed molecular models, we elucidate the interplay of anharmonic lattice degrees of freedom, dielectric confinement, and electronic correlation on exciton and biexciton binding energies of a range of nanocrystal sizes. We find generally good agreement with some experimental observations of binding energies and additionally explain the observed size-dependent Stokes shift. The explicit model calculations are compared with simplified approximations to rationalize the lattice contributions to binding. We find that polaron formation significantly reduces exciton binding energies, whereas these effects are negligible for biexciton interactions. While experimentally the binding energy of biexcitons is uncertain, based on our study, we conclude that biexcitons are bound in CsPbBr 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Designing polymeric membranes with coordination chemistry for high-precision ion separations

State-of-the-art polymeric membranes are unable to perform the high-precision ion separations needed for technologies essential to a circular economy and clean energy future. Coordinative interactions are a mechanism to increase sorption of a target species into a membrane, but the effects of these interactions on membrane permeability and selectivity are poorly understood. We use a multilayered polymer membrane to assess how ion-membrane binding energies affect membrane permeability of similarly sized cations: Cu 2+ , Ni 2+ , Zn 2+ , Co 2+ , and Mg 2+ . We report that metals with higher binding energy to iminodiacetate groups of the polymer more selectively permeate through the membrane in multisalt solutions than single-salt solutions. In contrast, weaker binding species are precluded from diffusing into the polymer membrane, which leads to passage proportional to binding energy and independent of membrane thickness. Our findings demonstrate that selectivity of polymeric membranes can markedly increase by tailoring ion-membrane binding energy and minimizing membrane thickness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Clostridioides difficile utilizes siderophores as an iron source and FhuDBGC contributes to ferrichrome uptake

ABSTRACT Clostridioides difficileremains a public health threat commonly observed following antibiotic use. Due to the importance of iron in many cell processes, most bacteria, includingC. difficile, have multiple mechanisms of acquiring iron. Previous studies have examined ferrous iron uptake inC. difficile; here we focus on the role of siderophores. In a growth assay, we show thatC. difficilecan use a variety of siderophores as the sole iron source. InC. difficile,two ABC transporters induced under low iron conditions are predicted siderophore importers: FhuDBGC and YclNOPQ. We hypothesized that these transporters are responsible for the uptake of the siderophores we tested. To investigate the specificity of these transporters, we purified the substrate binding proteins and examined siderophore binding using thermal shift. We demonstrate increased stability between one siderophore binding protein, FhuD, and the siderophore ferrichrome, suggesting a binding interaction. This specificity correlates with the inability of anΔfhuDBGCmutant to grow efficiently under iron limiting conditions in the presence of ferrichrome. WhileC. difficileused additional siderophores in our growth experiments, we did not observe increased thermal stability between the receptor proteins and any of the other siderophores tested, suggesting these siderophores do not bind these receptors and other siderophore import mechanisms remain to be elucidated. Redundancy in iron acquisition is a microbial survival adaptation to cope with the constant battle for iron within a host. Greater knowledge about howC. difficileacquires iron will provide insight about howC. difficilecolonizes and persists in the colon. IMPORTANCE This study is the first example ofC. difficilegrowing with siderophores as the sole iron source and describes the characterization of the ferric hydroxamate uptake ABC transporter (FhuDBGC). This transporter shows specificity to the siderophore ferrichrome. While not required for pathogenesis, this transporter highlights the redundancy in iron acquisition mechanisms thatC. difficileuses to compete for iron during an infection.

Microbiology↗
DOE OSTI2021

Peptidoglycan Contribution to the B Cell Superantigen Activity of Staphylococcal Protein A

Staphylococcus aureus causes reiterative and chronic persistent infections. This can be explained by the formidable ability of this pathogen to escape immune surveillance mechanisms. Cells of S. aureus display the abundant staphylococcal protein A (SpA). SpA binds to immunoglobulin (Ig) molecules and coats the bacterial surface to prevent phagocytic uptake. SpA also binds and cross-links variable heavy 3 (V H 3) idiotype (IgM) B cell receptors, promoting B cell expansion and the secretion of nonspecific V H 3-IgM via a mechanism requiring CD4 + T cell help. SpA binding to antibodies is mediated by the N-terminal Ig-binding domains (IgBDs). The so-called region X, uncharacterized LysM domain, and C-terminal LPXTG sorting signal for peptidoglycan attachment complete the linear structure of the protein. Here, we report that both the LysM domain and the LPXTG motif sorting signal are required for the B cell superantigen activity of SpA in a mouse model of infection. SpA molecules purified from staphylococcal cultures are sufficient to exert B cell superantigen activity and promote immunoglobulin secretion as long as they carry intact LysM and LPXTG motif domains with bound peptidoglycan fragments. The LysM domain binds the glycan chains of peptidoglycan fragments, whereas the LPXTG motif is covalently linked to wall peptides lacking glycan. These findings emphasize the complexity of SpA interactions with B cell receptors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Bindee

Bindee is a clang tool that outputs a simple pybind11 template given a C++ file for efficient generation of C++-Python bindings. Bindee is intended to be a helper tool for minimizing initial user effort and safeguarding against common runtime errors. Bindee relies on two open-source software to produce bindings. Clang's LibTooling enables bindee to traverse a C++ file's AST to pick out bindable variables and functions, or "bindees." PyBind11 is templated, header-only library for generating C++-Python bindings for variables and functions. Operating purely in C++, bindee does not require learning any new API for accomplishing its task. Additionally, picking the correct pybind11 API for a given bindable element is handled without interaction from the user. Any user input is denoted by '@TEXT@' string substitution. Bindee is capable of generating modular bindings for public class methods, public class variables, enumerations, and free functions and variables. Bindings for templates are also supported. For those familiar with pybind11, bindee does not handle trampolines, C++ extensions through lambdas, or custom type casters.

Kimkno, JasonS.↗
DOE OSTI2023

Structure of coagulation factor VIII bound to a patient-derived anti-C1 domain antibody inhibitor

The development of pathogenic antibody inhibitors against coagulation factor VIII (FVIII) occurs in ~30% of congenital hemophilia A patients receiving FVIII replacement therapy as well as in all cases of acquired hemophilia A. KM33 is an anti-C1 domain antibody inhibitor previously isolated from a severe hemophilia A patient. In addition to potently blocking FVIII binding to von Willebrand factor and phospholipid surfaces, KM33 disrupts FVIII binding to lipoprotein receptor-related protein 1 (LRP1), which drives FVIII hepatic clearance and antigen presentation in dendritic cells. Here, we report on the structure of FVIII bound to NB33, a recombinant derivative of KM33, by single-particle cryo-electron microscopy. Furthermore, structural analysis reveals the NB33 epitope localizes to FVIII residues R2090-S2094 and I2158-R2159 which constitute membrane-binding loops in the C1 domain. Further analysis reveals multiple FVIII lysine and arginine residues, previously shown to mediate binding to LRP1, dock onto an acidic cleft at the NB33 variable domain interface, thus blocking a putative LRP1 binding site. Together, these results demonstrate a novel mechanism of FVIII inhibition by a patient-derived antibody inhibitor and provide structural evidence toward engineering FVIII with reduced LRP1-mediated clearance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A generalized Flory-Stockmayer kinetic theory of connectivity percolation and rigidity percolation of cytoskeletal networks

Actin networks are essential for living cells to move, reproduce, and sense their environments. The dynamic and rheological behavior of actin networks is modulated by actin-binding proteins such as α-actinin, Arp2/3, and myosin. There is experimental evidence that actin-binding proteins modulate the cooperation of myosin motors by connecting the actin network. In this work, we present an analytical mean field model, using the Flory-Stockmayer theory of gelation, to understand how different actin-binding proteins change the connectivity of the actin filaments as the networks are formed. We follow the kinetics of the networks and estimate the concentrations of actin-binding proteins that are needed to reach connectivity percolation as well as to reach rigidity percolation. We find that Arp2/3 increases the actomyosin connectivity in the network in a non-monotonic way. We also describe how changing the connectivity of actomyosin networks modulates the ability of motors to exert forces, leading to three possible phases of the networks with distinctive dynamical characteristics: a sol phase, a gel phase, and an active phase. Thus, changes in the concentration and activity of actin-binding proteins in cells lead to a phase transition of the actin network, allowing the cells to perform active contraction and change their rheological properties.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Characterizing inhibitors of human AP endonuclease 1

AP endonuclease 1 (APE1) processes DNA lesions including apurinic/apyrimidinic sites and 3´-blocking groups, mediating base excision repair and single strand break repair. Much effort has focused on developing specific inhibitors of APE1, which could have important applications in basic research and potentially lead to clinical anticancer agents. We used structural, biophysical, and biochemical methods to characterize several reported inhibitors, including 7-nitroindole-2-carboxylic acid (CRT0044876), given its small size, reported potency, and widespread use for studying APE1. Intriguingly, NMR chemical shift perturbation (CSP) experiments show that CRT0044876 and three similar indole-2-carboxylic acids bind a pocket distal from the APE1 active site. A crystal structure confirms these findings and defines the pose for 5-nitroindole-2-carboxylic acid. However, dynamic light scattering experiments show the indole compounds form colloidal aggregates that could bind (sequester) APE1, causing nonspecific inhibition. Endonuclease assays show the compounds lack significant APE1 inhibition under conditions (detergent) that disrupt aggregation. Thus, binding of the indole-2-carboxylic acids at the remote pocket does not inhibit APE1 repair activity. Myricetin also forms aggregates and lacks APE1 inhibition under aggregate-disrupting conditions. Two other reported compounds (MLS000552981, MLS000419194) inhibit APE1 in vitro with low micromolar IC50 and do not appear to aggregate in this concentration range. However, NMR CSP experiments indicate the compounds do not bind specifically to apo- or Mg 2+ -bound APE1, pointing to a non-specific mode of inhibition, possibly DNA binding. Our results highlight methods for rigorous interrogation of putative APE1 inhibitors and should facilitate future efforts to discover compounds that specifically inhibit this important repair enzyme.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Computationally Guided and Experimentally Validated Design of Custom Chelators for Critical Mineral Recovery

Selective, high throughput separation of target critical metals from complex environments such as fly ash leachates and mining process streams presents a significant challenge for economical production. Custom chelators and sorbents are an attractive technology for selective metal extraction, however it can be difficult to predict their performance, and significant experimental efforts are often required to develop chelating technologies. Here, we present a computational strategy focused on modelling chelator-metal binding interactions and benchmark these results versus experimental data. A computational pipeline combining forcefield, semiempirical, and meta-GGA methods with a thermodynamic framework optimized for error cancellation has been developed to predict binding energies of chelator complexes towards critical mineral recovery applications. This approach, originally validated on [2.2.2] cryptates binding mono- and divalent cations, demonstrated robust predictive capabilities with an R2 of 0.850 against experimental aqueous binding energies. The workflow includes metadynamics for exploring high-dimensional potential energy surfaces and a cluster-continuum model for accurate yet computationally efficient solvation modeling. Error cancellation between solvation energies of free and chelator-coordinated ions enables faster convergence, even with finite cluster sizes. Initial studies on the cryptates revealed consistent metal-ligand coordination patterns, with systematic variations influenced by ion size and charge, highlighting key structural features linked to binding selectivity. Further studies of a proprietary chelator have resulted in identification of previously unreported selectivity towards economically significant metals, which in-house experiments have confirmed, demonstrating the feasibility of this approach. By applying this methodology to new chelators targeting critical minerals such as lithium, cobalt, nickel and other strategic metals, we aim to accelerate the discovery of next-generation chelators for efficient recovery, recycling, and separation processes. This computational framework serves as the backbone of a high-throughput design pipeline tailored for sustainable resource utilization and may be applied to a wide range of systems to meet experimental needs.

computational materials↗
DOE OSTI2022

A Coarse-Grained Model of Affinity Maturation Indicates the Importance of B-Cell Receptor Avidity in Epitope Subdominance

The elicitation of broadly neutralizing antibodies (bnAbs) is a major goal in the design of vaccines against rapidly-mutating viruses. In the case of influenza, many bnAbs that target conserved epitopes on the stem of the hemagglutinin protein (HA) have been discovered. However, these antibodies are rare, are not boosted well upon reinfection, and often have low neutralization potency, compared to strain-specific antibodies directed to the HA head. Different hypotheses have been proposed to explain this phenomenon. We use a coarse-grained computational model of the germinal center reaction to investigate how B-cell receptor binding valency affects the growth and affinity maturation of competing B-cells. We find that receptors that are unable to bind antigen bivalently, and also those that do not bind antigen cooperatively, have significantly slower rates of growth, memory B-cell production, and, under certain conditions, rates of affinity maturation. The corresponding B-cells are predicted to be outcompeted by B-cells that bind bivalently and cooperatively. We use the model to explore strategies for a universal influenza vaccine, e.g., how to boost the concentrations of the slower growing cross-reactive antibodies directed to the stem. The results suggest that, upon natural reinfections subsequent to vaccination, the protectiveness of such vaccines would erode, possibly requiring regular boosts. Collectively, our results strongly support the importance of bivalent antibody binding in immunodominance, and suggest guidelines for developing a universal influenza vaccine.

60 APPLIED LIFE SCIENCES↗