Characterizing the regulatory Fas (CD95) epitope critical for agonist antibody targeting and CAR-T bystander function in ovarian cancer.

Receptor clustering is the most critical step to activate extrinsic apoptosis by death receptors belonging to the TNF superfamily. Although clinically unsuccessful, using agonist antibodies, the death receptors-5 remains extensively studied from a cancer therapeutics perspective. However, despite its regulatory role and elevated function in ovarian and other solid tumors, another tumor-enriched death receptor called Fas (CD95) remained undervalued in cancer immunotherapy until recently, when its role in off-target tumor killing by CAR-T therapies was imperative. By comprehensively analyzing structure studies in the context of the binding epitope of FasL and various preclinical Fas agonist antibodies, we characterize a highly significant patch of positively charged residue epitope (PPCR) in its cysteine-rich domain 2 of Fas. PPCR engagement is indispensable for superior Fas agonist signaling and CAR-T bystander function in ovarian tumor models. A single-point mutation in FasL or Fas that interferes with the PPCR engagement inhibited apoptotic signaling in tumor cells and T cells. Furthermore, considering that clinical and immunological features of the autoimmune lymphoproliferative syndrome (ALPS) are directly attributed to homozygous mutations in FasL, we reveal differential mechanistic details of FasL/Fas clustering at the PPCR interface compared to described ALPS mutations. As Fas-mediated bystander killing remains vital to the success of CAR-T therapies in tumors, our findings highlight the therapeutic analytical design for potentially effective Fas-targeting strategies using death agonism to improve cancer immunotherapy in ovarian and other solid tumors.

A dynamic biomimetic model of the membrane-bound CD4-CD3-TCR complex during pMHC disengagement.

The coordinated (dis)engagement of the membrane-bound T cell receptor (TCR)-CD3-CD4 complex from the peptide-major histocompatibility complex (pMHC) is fundamental to TCR signal transduction and T cell effector function. As such, an atomic-scale understanding would not only enhance our basic understanding of the adaptive immune response but would also accelerate the rational design of TCRs for immunotherapy. In this study, we explore the impact of the CD4 coreceptor on the TCR-pMHC (dis)engagement by constructing a molecular-level biomimetic model of the CD3-TCR-pMHC and CD4-CD3-TCR-pMHC complexes within a lipid bilayer. After allowing the system complexes to equilibrate (engage), we use steered molecular dynamics to dissociate (disengage) the pMHC. We find that 1) the CD4 confines the pMHC closer to the T cell by 1.8 nm at equilibrium; 2) CD4 confinement shifts the TCR along the MHC binding groove engaging a different set of amino acids and enhancing the TCR-pMHC bond lifetime; 3) the CD4 translocates under load increasing the interaction strength between the CD4-pMHC, CD4-TCR, and CD4-CD3; and 4) upon dissociation, the CD3-TCR complex undergoes structural oscillation and increased energetic fluctuation between the CD3-TCR and CD3-lipids. These atomic-level simulations provide mechanistic insight on how the CD4 coreceptor impacts TCR-pMHC (dis)engagement. More specifically, our results provide further support (enhanced bond lifetime) for a force-dependent kinetic proofreading model and identify an alternate set of amino acids in the TCR that dominate the TCR-pMHC interaction and could thus impact the design of TCRs for immunotherapy.

A molecular dynamics investigation of N-glycosylation effects on T-cell receptor kinetics.

The binding interaction between the T-cell receptor (TCR) and peptide-major histocompatibility complex (pMHC) is modulated by several factors (known and unknown), however, investigations into effects of glycosylation are limited. A fully glycosylated computational model of the TCR bound to the pMHC is developed to investigate the effects of glycosylation on dissociation kinetics from the pMHC. Here, we examine the effects of N-glycosylation on TCR-pMHC bond strength using steered molecular dynamic simulations. N-glycosylation is a post-translational modification that adds sugar moieties to molecules and can modulate the activity of several immune molecules. Using a TCR-pMHC pair found in melanoma as a case study, our study demonstrates that N-glycosylation of the TCR-pMHC alters the proteins' conformation; increases the bond lifetime; and increases the number of hydrogen bonds and Lennard-Jones Contacts involved in the TCR-pMHC bond. We find that weak glycan-protein or glycan-glycan interactions impact the equilibrated structure of the TCR and pMHC leading to an increase in the overall bond strength of the TCR-pMHC complex including the duration and energetic strength under constant load. These results indicate that N-glycosylation plays an important role in the TCR-pMHC bond and should be considered in future computational and experimental studies.Communicated by Ramaswamy H. Sarma.

Determining structure and action mechanism of LBF14 by molecular simulation.

The recently discovered, membrane-active peptide LBF14 contains several non-proteinogenic amino acids and is able to transform vesicles into tubule networks. The exact membrane interaction mechanism and detailed secondary structure are yet to be determined. We performed molecular dynamics simulations of LBF14 and let it fold de novo into its ensemble of native secondary structures. Histidine protonation state effects on secondary structure were investigated. An MD simulation of the peptide with a lipid bilayer was performed. Simulation results were compared to circular dichroism and electron paramagnetic resonance data of previous studies. LBF14 contains a conserved helical section in an otherwise random structure. Helical stability is influenced by histidine protonation. The peptide localized to the polar layer of the membrane, consistent with experimental results. While the overall secondary structure is unaffected by membrane interaction, Ramachandran plot analysis yielded two distinct peptide conformations during membrane interaction. This conformational change was accompanied by residue repositioning within the membrane. LBF14 only affected the local order in the membrane, and had no measurable effect on pressure. The simulation results are consistent with the previously proposed membrane interaction mechanism of LBF14 and can additionally explain the local interaction mechanism. Communicated by Ramaswamy H. Sarma.

Effects of N-Glycosylation on the Structure, Function, and Stability of a Plant-Made Fc-Fusion Anthrax Decoy Protein.

Protein N-glycosylation is an important post-translational modification and has influences on a variety of biological processes at the cellular and molecular level, making glycosylation a major study aspect for glycoprotein-based therapeutics. To achieve a comprehensive understanding on how N-glycosylation impacts protein properties, an Fc-fusion anthrax decoy protein, viz rCMG2-Fc, was expressed in Nicotiana benthamiana plant with three types of N-glycosylation profiles. Three variants were produced by targeting protein to plant apoplast (APO), endoplasmic reticulum (ER) or removing the N-glycosylation site by a point mutation (Agly). Both the APO and ER variants had a complex-type N-glycan (GnGnXF) as their predominant glycans. In addition, ER variant had a higher concentration of mannose-type N-glycans (50%). The decoy protein binds to the protective antigen (PA) of anthrax through its CMG2 domain and inhibits toxin endocytosis. The protein expression, sequence, N-glycosylation profile, binding kinetics to PA, toxin neutralization efficiency, and thermostability were determined experimentally. In parallel, we performed molecular dynamics (MD) simulations of the predominant full-length rCMG2-Fc glycoform for each of the three N-glycosylation profiles to understand the effects of glycosylation at the molecular level. The MAN8 glycoform from the ER variant was additionally simulated to resolve differences between the APO and ER variants. Glycosylation showed strong stabilizing effects on rCMG2-Fc during in planta accumulation, evidenced by the over 2-fold higher expression and less protein degradation observed for glycosylated variants compared to the Agly variant. Protein function was confirmed by toxin neutralization assay (TNA), with effective concentration (EC50) rankings from low to high of 67.6 ng/ml (APO), 83.15 ng/ml (Agly), and 128.9 ng/ml (ER). The binding kinetics between rCMG2-Fc and PA were measured with bio-layer interferometry (BLI), giving sub-nanomolar affinities regardless of protein glycosylation and temperatures (25 and 37°C). The protein thermostability was examined utilizing the PA binding ELISA to provide information on EC50 differences. The fraction of functional ER variant decayed after overnight incubation at 37°C, and no significant change was observed for APO or Agly variants. In MD simulations, the MAN8 glycoform exhibits quantitatively higher distance between the CMG2 and Fc domains, as well as higher hydrophobic solvent accessible surface areas (SASA), indicating a possibly higher aggregation tendency of the ER variant. This study highlights the impacts of N-glycosylation on protein properties and provides insight into the effects of glycosylation on protein molecular dynamics.

Response to Extreme Temperatures of Mesoporous Silica MCM-41: Porous Structure Transformation Simulation and Modification of Gas Adsorption Properties.

Molecular dynamics (MD) and Monte Carlo (MC) simulations were applied together for the first time to reveal the porous structure transformation mechanisms of mesoporous silica MCM-41 subjected to temperatures up to 2885 K. Silica was experimentally characterized to inform the models and enable prediction of changes in gas adsorption/separation properties. MD simulations suggest that the pore closure process is activated by a collective diffusion of matrix atoms into the porous region, accompanied by bond reformation at the surface. Degradation is kinetically limited, such that complete pore closure is postponed at high heating rates. We experimentally observe decreased gas adsorption with increasing temperature in mesoporous silica heated at fixed rates, due to pore closure and structural degradation consistent with simulation predictions. Applying the Kissinger equation, we find a strong correlation between the simulated pore collapse temperatures and the experimental values which implies an activation energy of 416 ± 17 kJ/mol for pore closure. MC simulations give the adsorption and selectivity for thermally treated MCM-41, for N2, Ar, Kr, and Xe at room temperature within the 1-10 000 kPa pressure range. Relative to pristine MCM-41, we observe that increased surface roughness due to decreasing pore size amplifies the difference of the absolute adsorption amount differently for different adsorbate molecules. In particular, we find that adsorption of strongly interacting molecules can be enhanced in the low-pressure region while adsorption of weakly interacting molecules is inhibited. This then results in higher selectivity in binary mixture adsorption in mesoporous silica.

A quantum chemistry study of curvature effects on boron nitride nanotubes/nanosheets for gas adsorption.

Quantum chemistry calculations were performed to investigate the effect of the surface curvature of a Boron Nitride (BN) nanotube/nanosheet on gas adsorption. Curved boron nitride layers with different curvatures interacting with a number of different gases including noble gases, oxygen, and water on both their convex and concave sides of the surface were studied using density functional theory (DFT) with a high level dispersion corrected functional. Potential energy surfaces of the gas molecules interacting with the selected BN surfaces were investigated. In addition, the charge distribution and electrostatic potential contour of the selected BN surfaces are discussed. The results reveal how the curvature of the BN surfaces affects gas adsorption. In particular, small curvatures lead to a slight difference in the physisorption energy, while large curvatures present distinct potential energy surfaces, especially for the short-range repulsion.

Molecular dynamics of different polymer blends containing poly(2,6-dimethyl-1,4-phenylene ether).

Detailed atomistic molecular dynamics simulations were performed to investigate the behavior of two different binary blends, a miscible system poly(2,6-dimethyl-1,4-phenylene ether)-polystyrene (PPE-PS) and an immiscible system poly(2,6-dimethyl-1,4-phenylene ether)-poly(methyl methacrylate) (PPE-PMMA). We compared these two blends to study how PPE behaves when blended with different polymers. In both cases, the structure and phase behavior of polymer melts were studied by means of radial distribution functions (RDFs). Radii of gyration illustrate the static properties. Packing features of the benzene rings were observed in PPE and PS, both PS and PPE were well dispersed over the whole time scale of simulation. Furthermore, there was a tendency for aggregation of PMMA chains in PPE-PMMA systems. The mean squared displacements of monomers and whole chains describe the mobility of polymers in various systems.

Simulation of the effects of chain architecture on the sorption of ethylene in polyethylene.

An osmotic ensemble hyperparallel tempering technique has been developed to study the solubility of ethylene in amorphous linear low-density polyethylene of different chain architectures. The NERD united-atom force field (Nath, Escobedo, and de Pablo revised united-atom force field) is used in all simulations. We have investigated the effect of polyethylene chain length and branching on ethylene solubility. In this study, we have considered short-chain branching of amorphous linear low-density ethylene-1-hexene copolymers under typical polymerization reactor conditions. It is observed that, in the polymer, ethylene prefers to reside in the vicinity of polymer chain ends. This clustering causes a decrease in ethylene solubility with polymer chain length. When short-chain branches are introduced to a linear polymer chain, however, the chain-end clustering effect is counteracted by a higher density, thereby leading to an ethylene solubility almost identical to that in the linear polymer.