Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “binding affinities”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A substrate-induced gating mechanism is conserved among Gram-positive IgA1 metalloproteases

The mucosal adaptive immune response is dependent on the production of IgA antibodies and particularly IgA1, yet opportunistic bacteria have evolved mechanisms to specifically block this response by producing IgA1 proteases (IgA1Ps). Our lab was the first to describe the structures of a metal-dependent IgA1P (metallo-IgA1P) produced from Gram-positive Streptococcus pneumoniae both in the absence and presence of its IgA1 substrate through cryo-EM single particle reconstructions. This prior study revealed an active-site gating mechanism reliant on substrate-induced conformational changes to the enzyme that begged the question of whether such a mechanism is conserved among the wider Gram-positive metallo-IgA1P subfamily of virulence factors. Here, we used cryo-EM to characterize the metallo-IgA1P of a more distantly related family member from Gemella haemolysans, an emerging opportunistic pathogen implicated in meningitis, endocarditis, and more recently bacteremia in the elderly. While the substrate-free structures of these two metallo-IgA1Ps exhibit differences in the relative starting positions of the domain responsible for gating substrate, the enzymes have similar domain orientations when bound to IgA1. Together with biochemical studies that indicate these metallo-IgA1Ps have similar binding affinities and activities, these data indicate that metallo-IgA1P binding requires the specific IgA1 substrate to open the enzymes for access to their active site and thus, largely conform to an “induced fit” model.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Recognition of nuclear export signals by CRM1 carrying the oncogenic E571K mutation

The E571K mutation of CRM1 is highly prevalent in some cancers, but its mechanism of tumorigenesis is unclear. Glu571 of CRM1 is located in its nuclear export signal (NES)-binding groove, suggesting that binding of select NESs may be altered. We generated HEK 293 cells with either monoallelic CRM1WT/E571K or biallelic CRM1E571K/E571K using CRISPR/Cas9. We also combined analysis of binding affinities and structures of 27 diverse NESs for wild-type and E571K CRM1 with structure-based bioinformatics. While most NESs bind the two CRM1 similarly, NESs from Mek1, eIF4E-transporter, and RPS2 showed >10-fold affinity differences. These NESs have multiple charged side chains binding close to CRM1 position 571, but this feature alone was not sufficient to predict different binding to CRM1(E571K). Consistent with eIF4E-transporter NES binding weaker to CRM1(E571K), eIF4E-transporter was mislocalized in tumor cells carrying CRM1(E571K). This serves as proof of concept that understanding how CRM1(E571K) affects NES binding provides a platform for identifying cargoes that are mislocalized in cancer upon CRM1 mutation. Finally, we showed that large affinity changes seen with some NES peptides (of Mek1 and RPS2) do not always translate to the full-length cargoes, suggesting limitations with current NES prediction methods. Therefore, comprehensive studies like ours are imperative to identify CRM1 cargoes with real pathogenic potential.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Differentiating Ligand Tailoring and Cation Incorporation as Strategies for Tuning Heterobimetallic Cerium Complexes

Tuning of redox-active complexes featuring metals with high coordination numbers by incorporation of secondary redox-inactive cations has received far less attention than it deserves. Here, appending moderate steric bulk to a tripodal ligand framework has been tested for its influence on secondary-cation-driven structural and electrochemical tuning of cerium, a lanthanide that tends to adopt high coordination numbers. A quasi -C 3 -symmetric cerium(III) complex denoted [Ce] has been prepared that features pendant benzyloxy groups, and this work demonstrates that this species offers a site capable of binding single Na + or Ca 2+ ions. Electrochemical and UV-visible spectroscopic studies reveal equilibrium binding affinity of [Ce] for Na + in acetonitrile solvent, contrasting with tight binding of all cations in all other previously studied systems of this type. The modulated cation binding can be attributed to the bulky benzyloxy groups, which impact the thermodynamics of cation binding but do not impede the formation of cerium centers with coordination number 8 upon binding of either Na + or Ca2 + . Finally, the Ce(IV/III) reduction potential was found to be tunable under the equilibrium binding conditions, highlighting the potentially significant role that controlled structural changes can play in modulating the solution properties of heterobimetallic complexes.

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

TwoFold: Highly accurate structure and affinity prediction for protein-ligand complexes from sequences

We describe our development of ab initio protein-ligand binding pose prediction models based on transformers and binding affinity prediction models based on the neural tangent kernel (NTK). Folding both protein and ligand, the TwoFold models achieve efficient and quality predictions matching state-of-the-art implementations while additionally reconstructing protein structures. In conclusion, solving NTK models points to a new use case for highly optimized linear solver benchmarking codes on HPC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Structurally silent peptide anchor modifications allosterically modulate T cell recognition in a receptor-dependent manner

Presentation of peptides by class I MHC proteins underlies T cell immune responses to pathogens and cancer. The association between peptide binding affinity and immunogenicity has led to the engineering of modified peptides with improved MHC binding, with the hope that these peptides would be useful for eliciting cross-reactive immune responses directed toward their weak binding, unmodified counterparts. Increasing evidence, however, indicates that T cell receptors (TCRs) can perceive such anchor-modified peptides differently than wild-type (WT) peptides, although the scope of discrimination is unclear. We show here that even modifications at primary anchors that have no discernible structural impact can lead to substantially stronger or weaker T cell recognition depending on the TCR. Surprisingly, the effect of peptide anchor modification can be sensed by a TCR at regions distant from the site of modification, indicating a through-protein mechanism in which the anchor residue serves as an allosteric modulator for TCR binding. Our findings emphasize caution in the use and interpretation of results from anchor-modified peptides and have implications for how anchor modifications are accounted for in other circumstances, such as predicting the immunogenicity of tumor neoantigens. Our data also highlight an important need to better understand the highly tunable dynamic nature of class I MHC proteins and the impact this has on various forms of immune recognition.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Decoupling of catalysis and transition state analog binding from mutations throughout a phosphatase revealed by high-throughput enzymology

Using high-throughput microfluidic enzyme kinetics (HT-MEK), we measured over 9,000 inhibition curves detailing impacts of 1,004 single-site mutations throughout the alkaline phosphatase PafA on binding affinity for two transition state analogs (TSAs), vanadate and tungstate. As predicted by catalytic models invoking transition state complementary, mutations to active site and active-site-contacting residues had highly similar impacts on catalysis and TSA binding. Unexpectedly, most mutations to more distal residues that reduced catalysis had little or no impact on TSA binding and many even increased tungstate affinity. These disparate effects can be accounted for by a model in which distal mutations alter the enzyme’s conformational landscape, increasing the occupancy of microstates that are catalytically less effective but better able to accommodate larger transition state analogs. In support of this ensemble model, glycine substitutions (rather than valine) were more likely to increase tungstate affinity (but not more likely to impact catalysis), presumably due to increased conformational flexibility that allows previously disfavored microstates to increase in occupancy. These results indicate that residues throughout an enzyme provide specificity for the transition state and discriminate against analogs that are larger only by tenths of an Ångström. Thus, engineering enzymes that rival the most powerful natural enzymes will likely require consideration of distal residues that shape the enzyme’s conformational landscape and fine-tune active-site residues. Biologically, the evolution of extensive communication between the active site and remote residues to aid catalysis may have provided the foundation for allostery to make it a highly evolvable trait.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Selective acetate recognition and sensing using SWCNTs functionalized with croconamides

Croconamides are a new class of dual-hydrogen bond donors that have attracted interest in anion recognition and sensing applications due to their strong H-bond donation capability. In this work, we report chemiresistive sensors for the detection of anions based on poly(4-vinylpyridine) (P4VP)-wrapped single-walled carbon nanotubes (SWCNTs) functionalized with croconamide-based selectors. Model structures of asymmetric N,N-substituted croconamide selectors with N,N-substitutional groups of alkyl cationic pyridinium and electron-withdrawing moieties such as aniline (6) and 3,5-bis(trifluoromethyl)phenyl (7) groups were synthesized and their binding affinities for various anions were characterized using UV–vis and 1 H NMR titrations. Switchable binding of oxoanions was observed with selector 7, wherein HSO 4 – formed hydrogen bonding interactions and acetate (AcO – ) induced deprotonation of the NH- protons of the croconamide. Large sensitivity transitions were obtained toward AcO – for both P4VP-6-SWCNT (S = 95.68 ± 2.11 %) and P4VP-7-SWCNT (S = 140.91 ± 3.68 %) at 83.33 mM, which was attributed to deprotonation of the croconamide selectors upon addition of AcO – . In conclusion, three serial sensors were assembled in a sensor array that rapid and selective screening of AcO – was demonstrated against other oxoanions such as HSO 4 – and H 2 PO 4 – in real-time.

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

Efficacy and Delivery of Novel FAST Agents for Coronaviruses

We proposed to test and develop advanced delivery for novel agents from our collaborators Facile Accelerated Specific Therapeutics (FAST) platform to reduce coronavirus replication. Sachi Bioworks Inc., Prof. Anushree Chatterjee, and Prof. Prashant Nagpal at the University of Colorado Boulder have developed a bioinformatics and synthesis pipeline to produce sequence specific theranostic agents (agents that can be therapies and/or diagnostics) that are inherently transported into the cytoplasm of mammalian host cells and sequence-specifically interfere in nucleic acid replication. The agent comprises a small nanoparticle (2-5 nm) chosen for ideal cellular transport and/or imaging conjugated to a short, synthetic DNA analog oligomer designed for binding to one or more target viral sequences. The sequence specific binding of the FAST agent to its target prevents nucleic acid replication due to its high affinity binding. While the small nanoparticle facilitates delivery in vitro, we plan to package the FAST agents into a larger nanoparticle (80-300 nm) for future in vivo delivery applications. Our team at Sandia has expertise encapsulating biomolecules including protein, DNA, and RNA into solid lipid nanoparticles (LNP) and lipid coated mesoporous silica nanoparticles (LC-MSN) and shown successful delivery in mouse models to multiple tissues. Our team focused on formulation parameters for FAST agents into lipid nanoparticles (LNP) and lipid coated mesoporous silica nanoparticles (LC-MSN) for enhanced delivery and/or efficacy and in vivo translation. We used lipid formulas that have been shown in literature to facility in vitro and more importantly, in vivo delivery. In our work discussed below, we successfully demonstrate loading and release of FAST agents on silica core and stable LC-MSN in a reasonable size range for in vivo testing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

From small changes to big gains: pyridinium-based tetralactam macrocycle for enhanced sugar recognition in water

The complex distribution of functional groups in carbohydrates, coupled with their strong solvation in water, makes them challenging targets for synthetic receptors. Despite extensive research into various molecular frameworks, most synthetic carbohydrate receptors have exhibited low affinities, and their interactions with sugars in aqueous environments remain poorly understood. In this work, we present a simple pyridinium-based hydrogen-bonding receptor derived from a subtle structural modification of a well-known tetralactam macrocycle. This small structural change resulted in a dramatic enhancement of glucose binding affinity, increasing from 56 M −1 to 3001 M −1 . Remarkably, the performance of our synthetic lectin surpasses that of the natural lectin, concanavalin A, by over fivefold. X-ray crystallography of the macrocycle–glucose complex reveals a distinctive hydrogen bonding pattern, which allows for a larger surface overlap between the receptor and glucose, contributing to the enhanced affinity. Furthermore, this receptor possesses allosteric binding sites, which involve chloride binding and trigger receptor aggregation. This unique allosteric process reveals the critical role of structural flexibility in this hydrogen-bonding receptor for the effective recognition of sugars. We also demonstrate the potential of this synthetic lectin as a highly sensitive glucose sensor in aqueous solutions.

Zhai, Canjia↗
DOE OSTI2022

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. In this study, we aim to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

A ligand discovery toolbox for the WWE domain family of human E3 ligases

The WWE domain is a relatively under-researched domain found in twelve human proteins and characterized by a conserved tryptophan-tryptophan-glutamate (WWE) sequence motif. Six of these WWE domain-containing proteins also contain domains with E3 ubiquitin ligase activity. The general recognition of poly-ADP-ribosylated substrates by WWE domains suggests a potential avenue for development of Proteolysis-Targeting Chimeras (PROTACs). Here, we present novel crystal structures of the HUWE1, TRIP12, and DTX1 WWE domains in complex with PAR building blocks and their analogs, thus enabling a comprehensive analysis of the PAR binding site structural diversity. Furthermore, we introduce a versatile toolbox of biophysical and biochemical assays for the discovery and characterization of novel WWE domain binders, including fluorescence polarization-based PAR binding and displacement assays, 15 N-NMR-based binding affinity assays and 19 F-NMR-based competition assays. Through these assays, we have characterized the binding of monomeric iso -ADP-ribose ( iso -ADPr) and its nucleotide analogs with the aforementioned WWE proteins. Finally, we have utilized the assay toolbox to screen a small molecule fragment library leading to the successful discovery of novel ligands targeting the HUWE1 WWE domain.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Local Chain Dynamics in Sequence-Controlled Polymers as a Tunable Handle for Rare Earth Sequestration

Chain dynamics govern the intricate behaviors of proteins, underpinning functions such as catalysis, recognition, and stimulus response, and are an increasingly appreciated aspect of structure–function relationships. Analogously, manipulating chain dynamics and structure in abiotic polymers via sequence control is an exciting, yet underexplored, strategy for improving material functions. In this work, we report a systematic study relating the sequence of polymeric sequestrants to their structure and dynamics, as well as to their binding affinity and selectivity for model substrates, rare earth elements (REEs). A series of sequence-controlled polymers with metal chelating, solubilizing, and structure forming monomers was synthesized via multiblock polymerization, yielding compositionally identical polymers with spectroscopically resolved domains and distinct morphologies. Using a combination of small-angle X-ray scattering and 19 F NMR relaxometry measurements, we connected differences in polymer structure and dynamics to polymer sequence variables such as the patchiness (density) of the structure forming monomer and the location of the chelating monomer. Furthermore, we found that, relative to calcium, all polymers in the series collapse more and have slower dynamics when binding REEs (lanthanum and lutetium) , though the extent of these effects were sequence-dependent and localized to specific domains within the polymer. Notably, sequence-controlled polymers that exhibited the largest conformational and dynamic changes upon binding REEs also bound REEs with the greatest affinity and modest selectivity. Collectively, these results correlate monomer patterning with dynamics, morphology, and REE binding performance en route to the development of efficient and selective macromolecular chelators.

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

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MPro WT ) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373), and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ~400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decrease in free energy (ΔG) and thermal stability (T m ), respectively, relative to MPro WT . The H41A mutation results in a 20-fold increase in the dissociation constant (K d ), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and T m , respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MPro H41A , whereas no significant change is observed in binding to MPro WT . Structures of the four inhibitor complexes with MPro 1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MPro WT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MPro WT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MPro C145A also results in a significant, but smaller, increase in K d and decrease in ΔG and T m , relative to NMV.

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

Influence of Steric and Electronic Properties of P2 Groups on Covalent Inhibitor Binding to SARS-CoV-2 Main Protease

The main protease (MPro) of SARS-CoV-2 is a critical enzyme required for viral replication, making it a prime target for antiviral drug development. Covalent inhibitors, which form a stable interaction with the catalytic C145, have demonstrated strong inhibition of MPro, but the influence of steric and electronic properties of P2 substituents, designed to engage the S2 substrate-binding subsite within the MPro active site, on inhibitor binding affinity remains underexplored. Here, in this study, we design and characterize two hybrid covalent inhibitors, BBH-3 and BBH-4, and present their X-ray crystallographic structures in complex with MPro, providing molecular insights into how their distinct P2 groups, a dichlorobenzyl moiety in BBH-3 and an adamantyl substituent in BBH-4, affect binding conformation and active site adaptability. Comparative structural analyses with previously characterized inhibitors, including BBH-2 and Mcule-5948770040, reveal how the P2 bulkiness and electronic properties influence active site dynamics, particularly through interactions with the S2 and S5 subsites. The P2 group of BBH-3 induces conformational shifts in the S2 helix and the S5 loop, while BBH-4 displaces M49, stabilizing its binding through hydrophobic interactions. Isothermal titration calorimetry further elucidates the impact of P2 modifications on inhibitor affinity, revealing a delicate balance between enthalpic and entropic contributions. The data demonstrate that BBH-3 exhibits less favorable binding, affirming that dichlorobenzyl substitution at the P2 position has a more negative impact on the affinity for MPro than bulky saturated cyclic groups. This underscores the feature that MPro active site malleability may be accompanied by a conformational strain.

SARS-CoV-2↗
DOE OSTI2025

Structural insights into isoform-specific RAS-PI3Kα interactions and the role of RAS in PI3Kα activation

Abstract Mutations in RAS and PI3Kα are major drivers of human cancer. Their interaction plays a crucial role in activating PI3Kα and amplifying the PI3K-AKT-mTOR pathway. Disrupting RAS-PI3Kα interaction enhances survival in lung and skin cancer models and reduces tumor growth and angiogenesis, although the structural details of this interaction remain unclear. Here, we present structures of KRAS, RRAS2, and MRAS bound to the catalytic subunit (p110α) of PI3Kα, elucidating the interaction interfaces and local conformational changes upon complex formation. Structural and mutational analyses highlighted key residues in RAS and PI3Kα impacting binding affinity and revealed isoform-specific differences at the interaction interface in RAS and PI3K isoforms, providing a rationale for their differential affinities. Notably, in the RAS-p110α complex structures, RAS interaction with p110α is limited to the RAS-binding domain and does not involve the kinase domain. This study underscores the pivotal role of the RAS-PI3Kα interaction in PI3Kα activation and provides a blueprint for designing PI3Kα isoform-specific inhibitors to disrupt this interaction.

Science & Technology - Other Topics↗
DOE OSTI2025

Structure–Function Relationships in Sequence-Controlled Copolymers for Rare Earth Element Chelation

The ability to tune material function through primary sequence is a defining feature of biological macromolecules, allowing precise control over structure and target interactions in complex aqueous environments. However, translating sequence–structure–function relationships to synthetic macromolecules is challenging due to their dispersity in sequence, conformation, and composition. Here, we report systematic studies of amphiphilic polymer chelators designed to probe how composition and patterning influence binding affinity and selectivity for rare earth elements (REEs), a series of technologically relevant metals with challenging separation profiles. A library of copolymers varying hydrophobic monomer composition and patterning was synthesized via reversible addition–fragmentation chain transfer (RAFT) polymerization, spanning statistical, gradient, and block architectures. REE binding was quantified using a high-throughput colorimetric assay, and reconstruction of polymer ensembles using kinetic stochastic simulations enabled quantitative comparisons of sequence heterogeneity, linking local monomer colocalization to emergent REE binding. Further, we investigated the role of different hydrophobic comonomers in tuning metal coordination, with binding trends linked to structural features that influence binding site desolvation. Complementary dynamic light scattering (DLS) and small-angle X-ray scattering (SAXS) measurements showed that both polymer and monomer architecture modulate metal-induced conformational changes, and that multichain assembly behavior emerges beyond critical hydrophobic thresholds. Sequence control also altered REE selectivity, with nonmonotonic differences observed across compositionally identical polymers with different sequence architectures. Together, these findings establish design principles that connect polymer sequence and structure to binding performance, guiding the design of macromolecular chelators with enhanced affinity and selectivity for applications in separations, sensing, and catalysis.

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

The effect of protein mutations on drug binding suggests ensuing personalised drug selection

The advent of personalised medicine promises a deeper understanding of mechanisms and therefore therapies. However, the connection between genomic sequences and clinical treatments is often unclear. We studied 50 breast cancer patients belonging to a population-cohort in the state of Qatar. From Sanger sequencing, we identified several new deleterious mutations in the estrogen receptor 1 gene (ESR1). The effect of these mutations on drug treatment in the protein target encoded by ESR1, namely the estrogen receptor, was achieved via rapid and accurate protein–ligand binding affinity interaction studies which were performed for the selected drugs and the natural ligand estrogen. Four nonsynonymous mutations in the ligand-binding domain were subjected to molecular dynamics simulation using absolute and relative binding free energy methods, leading to the ranking of the efficacy of six selected drugs for patients with the mutations. Our study shows that a personalised clinical decision system can be created by integrating an individual patient’s genomic data at the molecular level within a computational pipeline which ranks the efficacy of binding of particular drugs to variant proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Computing the Relative Affinity of Chlorophylls a and b to Light-Harvesting Complex II

In plants and algae, the primary antenna protein bound to photosystem II is light-harvesting complex II (LHCII), a pigment–protein complex that binds eight chlorophyll (Chl) a molecules and six Chl b molecules. Chl a and Chl b differ only in that Chl a has a methyl group (–CH 3 ) on one of its pyrrole rings, while Chl b has a formyl group (–CHO) at that position. This blue-shifts the Chl b absorbance relative to Chl a . It is not known how the protein selectively binds the right Chl type at each site. Knowing the selection criteria would allow the design of light-harvesting complexes that bind different Chl types, modifying an organism to utilize the light of different wavelengths. The difference in the binding affinity of Chl a and Chl b in pea and spinach LHCII was calculated using multiconformation continuum electrostatics and free energy perturbation. Both methods have identified some Chl sites where the bound Chl type ( a or b ) has a significantly higher affinity, especially when the protein provides a hydrogen bond for the Chl b formyl group. However, the Chl a sites often have little calculated preference for one Chl type, so they are predicted to bind a mixture of Chl a and b . The electron density of the spinach LHCII was reanalyzed, which, however, confirmed that there is negligible Chl b in the Chl a -binding sites. Finally, it is suggested that the protein chooses the correct Chl type during folding, segregating the preferred Chl to the correct binding site.

chemical calculations↗