Enhancement of SARS-CoV-2 Infection via Crosslinking of Adjacent Spike Proteins by N-Terminal Domain-Targeting Antibodies.

The entry of SARS-CoV-2 into host cells is mediated by the interaction between the spike receptor-binding domain (RBD) and host angiotensin-converting enzyme 2 (ACE2). Certain human antibodies, which target the spike N-terminal domain (NTD) at a distant epitope from the host cell binding surface, have been found to augment ACE2 binding and enhance SARS-CoV-2 infection. Notably, these antibodies exert their effect independently of the antibody fragment crystallizable (Fc) region, distinguishing their mode of action from previously described antibody-dependent infection-enhancing (ADE) mechanisms. Building upon previous hypotheses and experimental evidence, we propose that these NTD-targeting infection-enhancing antibodies (NIEAs) achieve their effect through the crosslinking of neighboring spike proteins. In this study, we present refined structural models of NIEA fragment antigen-binding region (Fab)-NTD complexes, supported by molecular dynamics simulations and hydrogen-deuterium exchange mass spectrometry (HDX-MS). Furthermore, we provide direct evidence confirming the crosslinking of spike NTDs by NIEAs. Collectively, our findings advance our understanding of the molecular mechanisms underlying NIEAs and their impact on SARS-CoV-2 infection.

Celastrol suppresses humoral immune responses and autoimmunity by targeting the COMMD3/8 complex.

Celastrol, a bioactive molecule extracted from the Tripterygium wilfordii plant, has been shown to exhibit anti-inflammatory properties. However, its mechanism of action has not been fully elucidated. Here, we show that celastrol suppresses humoral immune responses and autoimmunity by disabling a protein complex consisting of copper metabolism MURR1 domain-containing (COMMD) 3 and COMMD8 (COMMD3/8 complex), a signaling adaptor for chemoattractant receptors. Having demonstrated the involvement of the COMMD3/8 complex in a mouse model of rheumatoid arthritis, we identified celastrol as a compound that covalently bound to and dissociated the COMMD3/8 complex. Celastrol inhibited B cell migration, reduced antibody responses, and blocked arthritis progression, recapitulating deficiency of the COMMD3/8 complex. These effects of celastrol were abolished in mice expressing a celastrol-resistant mutant of the COMMD3/8 complex. These findings establish that celastrol exerts immunosuppressive activity by targeting the COMMD3/8 complex. Our study suggests that the COMMD3/8 complex is a potentially druggable target in autoimmune diseases and points to celastrol as a lead pharmacologic candidate in this capacity.

Flexible, Functional, and Familiar: Characteristics of SARS-CoV-2 Spike Protein Evolution.

The SARS-CoV-2 S protein is a major point of interaction between the virus and the human immune system. As a consequence, the S protein is not a static target but undergoes rapid molecular evolution. In order to more fully understand the selection pressure during evolution, we examined residue positions in the S protein that vary greatly across closely related viruses but are conserved in the subset of viruses that infect humans. These "evolutionarily important" residues were not distributed evenly across the S protein but were concentrated in two domains: the N-terminal domain and the receptor-binding domain, both of which play a role in host cell binding in a number of related viruses. In addition to being localized in these two domains, evolutionary importance correlated with structural flexibility and inversely correlated with distance from known or predicted host receptor-binding residues. Finally, we observed a bias in the composition of the amino acids that make up such residues toward more human-like, rather than virus-like, sequence motifs.

Methods for sequence and structural analysis of B and T cell receptor repertoires.

B cell receptors (BCRs) and T cell receptors (TCRs) make up an essential network of defense molecules that, collectively, can distinguish self from non-self and facilitate destruction of antigen-bearing cells such as pathogens or tumors. The analysis of BCR and TCR repertoires plays an important role in both basic immunology as well as in biotechnology. Because the repertoires are highly diverse, specialized software methods are needed to extract meaningful information from BCR and TCR sequence data. Here, we review recent developments in bioinformatics tools for analysis of BCR and TCR repertoires, with an emphasis on those that incorporate structural features. After describing the recent sequencing technologies for immune receptor repertoires, we survey structural modeling methods for BCR and TCRs, along with methods for clustering such models. We review downstream analyses, including BCR and TCR epitope prediction, antibody-antigen docking and TCR-peptide-MHC Modeling. We also briefly discuss molecular dynamics in this context.

Structural Modeling of Lymphocyte Receptors and Their Antigens.

Structural modeling plays a key role in protein function prediction on a genome-wide scale. For B and T lymphocyte receptors, the critical functional question is: which antigens and epitopes are targeted? With emerging B cell receptor (BCR) and T cell receptor (TCR) sequencing methods improving in both breadth and depth, there is a growing need for methods that can help answer this question. Since lymphocyte-antigen recognition depends on complementarity, structural modeling is likely to play an important role in understanding antigen specificity and affinity. In the case of BCRs, such modeling methods have a long history in the study and design of antibodies. However, for TCRs there are relatively few publicly available modeling tools, and, to our knowledge, none that incorporate interaction between TCRs and peptide-MHC (pMHC) complexes. Here, we provide a web-based tool, ImmuneScape ( https://sysimm.org/immune-scape/ ), to carry out TCR-pMHC modeling as a first step toward structure-based function prediction.

Understanding Conformational Dynamics of Complex Lipid Mixtures Relevant to Biology.

This is a perspective article entitled "Frontiers in computational biophysics: understanding conformational dynamics of complex lipid mixtures relevant to biology" which is following a CECAM meeting with the same name.

Prediction of Thylakoid Lipid Binding Sites on Photosystem II.

The thylakoid membrane has a unique lipid composition, consisting mostly of galactolipids. These thylakoid lipids have important roles in photosynthesis. Here, we investigate to what extent these lipids bind specifically to the Photosystem II complex. To this end, we performed coarse-grain MD simulations of the Photosystem II complex embedded in a thylakoid membrane with realistic composition. Based on >85 µs simulation time, we find that monogalactosyldiacylglycerol and sulfoquinovosyldiacylglycerol lipids are enriched in the annular shell around the protein, and form distinct binding sites. From the analysis of residue contacts, we conclude that electrostatic interactions play an important role in stabilizing these binding sites. Furthermore, we find that chlorophyll a has a prevalent role in the coordination of the lipids. In addition, we observe lipids to diffuse in and out of the plastoquinone exchange cavities, allowing exchange of cocrystallized lipids with the bulk membrane and suggesting a more open nature of the plastoquinone exchange cavity. Together, our data provide a wealth of information on protein-lipid interactions for a key protein in photosynthesis.

Design and Properties of Genetically Encoded Probes for Sensing Macromolecular Crowding.

Cells are highly crowded with proteins and polynucleotides. Any reaction that depends on the available volume can be affected by macromolecular crowding, but the effects of crowding in cells are complex and difficult to track. Here, we present a set of Förster resonance energy transfer (FRET)-based crowding-sensitive probes and investigate the role of the linker design. We investigate the sensors in vitro and in vivo and by molecular dynamics simulations. We find that in vitro all the probes can be compressed by crowding, with a magnitude that increases with the probe size, the crowder concentration, and the crowder size. We capture the role of the linker in a heuristic scaling model, and we find that compression is a function of size of the probe and volume fraction of the crowder. The FRET changes observed in Escherichia coli are more complicated, where FRET-increases and scaling behavior are observed solely with probes that contain the helices in the linker. The probe with the highest sensitivity to crowding in vivo yields the same macromolecular volume fractions as previously obtained from cell dry weight. The collection of new probes provides more detailed readouts on the macromolecular crowding than a single sensor.

Exchange pathways of plastoquinone and plastoquinol in the photosystem II complex.

Plastoquinone (PLQ) acts as an electron carrier between photosystem II (PSII) and the cytochrome b6f complex. To understand how PLQ enters and leaves PSII, here we show results of coarse grained molecular dynamics simulations of PSII embedded in the thylakoid membrane, covering a total simulation time of more than 0.5 ms. The long time scale allows the observation of many spontaneous entries of PLQ into PSII, and the unbinding of plastoquinol (PLQol) from the complex. In addition to the two known channels, we observe a third channel for PLQ/PLQol diffusion between the thylakoid membrane and the PLQ binding sites. Our simulations point to a promiscuous diffusion mechanism in which all three channels function as entry and exit channels. The exchange cavity serves as a PLQ reservoir. Our simulations provide a direct view on the exchange of electron carriers, a key step of the photosynthesis machinery.

Molecular Dynamics of Photosystem II Embedded in the Thylakoid Membrane.

Photosystem II (PSII) is one of the key protein complexes in photosynthesis. We introduce a coarse grained model of PSII and present the analysis of 60 µs molecular dynamics simulations of PSII in both monomeric and dimeric form, embedded in a thylakoid membrane model that reflects its native lipid composition. We describe in detail the setup of the protein complex and the many natural cofactors and characterize their mobility. Overall we find that the protein subunits and cofactors are more flexible toward the periphery of the complex as well as near the PLQ exchange cavity and at the dimer interface. Of all cofactors, ß-carotenes show the highest mobility. Some of the ß-carotenes diffuse in and out of the protein complex via the thylakoid membrane. In contrast with the PSII dimer, the monomeric form adopts a tilted conformation in the membrane, with strong interactions between the soluble PsbO subunit and the glycolipid headgroups. Interestingly, the tilted conformation causes buckling of the membrane. Together, our results provide an unprecedented view of PSII dynamics on a microsecond time scale. Our data may be used as basis for the interpretation of experimental data as well as for theoretical models describing exciton energy transfer.

Characterization of thylakoid lipid membranes from cyanobacteria and higher plants by molecular dynamics simulations.

The thylakoid membrane is mainly composed of non-common lipids, so called galactolipids. Despite the importance of these lipids for the function of the photosynthetic reaction centers, the molecular organization of these membranes is largely unexplored. Here we use multiscale molecular dynamics simulations to characterize the thylakoid membrane of both cyanobacteria and higher plants. We consider mixtures of up to five different galactolipids plus phosphatidylglycerol to represent these complex membranes. We find that the different lipids generally mix well, although nanoscale heterogeneities are observed especially in case of the plant membrane. The fluidity of the cyanobacterial membrane is markedly reduced compared to the plant membrane, even considering elevated temperatures at which thermophilic cyanobacteria are found. We also find that the plant membrane more readily undergoes a phase transformation to an inverted hexagonal phase. We furthermore characterized the conformation and dynamics of the cofactors plastoquinone and plastoquinol, revealing of the fast flip-flop rates for the non-reduced form. Together, our results provide a molecular view on the dynamical organization of the thylakoid membrane.

Lipid organization of the plasma membrane.

The detailed organization of cellular membranes remains rather elusive. Based on large-scale molecular dynamics simulations, we provide a high-resolution view of the lipid organization of a plasma membrane at an unprecedented level of complexity. Our plasma membrane model consists of 63 different lipid species, combining 14 types of headgroups and 11 types of tails asymmetrically distributed across the two leaflets, closely mimicking an idealized mammalian plasma membrane. We observe an enrichment of cholesterol in the outer leaflet and a general non-ideal lateral mixing of the different lipid species. Transient domains with liquid-ordered character form and disappear on the microsecond time scale. These domains are coupled across the two membrane leaflets. In the outer leaflet, distinct nanodomains consisting of gangliosides are observed. Phosphoinositides show preferential clustering in the inner leaflet. Our data provide a key view on the lateral organization of lipids in one of life's fundamental structures, the cell membrane.

Martini Force Field Parameters for Glycolipids.

We present an extension of the Martini coarse-grained force field to glycolipids. The glycolipids considered here are the glycoglycerolipids monogalactosyldiacylglycerol (MGDG), sulfoquinovosyldiacylglycerol (SQDG), digalactosyldiacylglycerol (DGDG), and phosphatidylinositol (PI) and its phosphorylated forms (PIP, PIP2), as well as the glycosphingolipids galactosylceramide (GCER) and monosialotetrahexosylganglioside (GM1). The parametrization follows the same philosophy as was used previously for lipids, proteins, and carbohydrates focusing on the reproduction of partitioning free energies of small compounds between polar and nonpolar solvents. Bonded parameters are optimized by comparison to lipid conformations sampled with an atomistic force field, in particular with respect to the representation of the most populated states around the glycosidic linkage. Simulations of coarse-grained glycolipid model membranes show good agreement with atomistic simulations as well as experimental data available, especially concerning structural properties such as electron densities, area per lipid, and membrane thickness. Our coarse-grained model opens the way to large scale simulations of biological processes in which glycolipids are important, including recognition, sorting, and clustering of both external and membrane bound proteins.