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.

Selective scandium ion capture through coordination templating in a covalent organic framework

The use of coordination complexes within covalent organic frameworks can significantly diversify the structures and properties of this class of materials. Here we combined coordination chemistry and reticular chemistry by preparing frameworks that consist of a ditopic (p-phenylenediamine) and mixed tritopic moieties-an organic ligand and a scandium coordination complex of similar sizes and geometries, both bearing terminal phenylamine groups. Changing the ratio of organic ligand to scandium complex enabled the preparation of a series of crystalline covalent organic frameworks with tunable levels of scandium incorporation. Removal of scandium from the material with the highest metal content subsequently resulted in a 'metal-imprinted' covalent organic framework that exhibits a high affinity and capacity for Sc3+ ions in acidic environments and in the presence of competing metal ions. In particular, the selectivity of this framework for Sc3+ over common impurity ions such as La3+ and Fe3+ surpasses that of existing scandium adsorbents.

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.

Development of a Reactive Force Field for Simulating Photoinitiated Acrylate Polymerization

Light-driven and photocurable polymer-based additive manufacturing (AM) has enormous potential due to its excellent resolution and precision. Acrylated resins that undergo radical chain-growth polymerization are widely used in photopolymer AM due to their fast kinetics and often serve as a departure point for developing other resin materials for photopolymer-based AM technologies. For successful control of the photopolymer resins, the molecular basis of the acrylate free-radical polymerization has to be understood in detail. We present an optimized reactive force field (ReaxFF) for molecular dynamics (MD) simulations of acrylate polymer resins that captures radical polymerization thermodynamics and kinetics. The force field is trained against an extensive training set including density functional theory (DFT) calculations of reaction pathways along the radical polymerization from methyl acrylate to methyl butyrate, bond dissociation energies, and structures and partial charges of several molecules and radicals. We also found that it was critical to train the force field against an incorrect, nonphysical reaction pathway observed in simulations that used parameters not optimized for acrylate polymerization. The parameterization process utilizes a parallelized search algorithm, and the resulting model can describe polymer resin formation, crosslinking density, conversion rate, and residual monomers of the complex acrylate mixtures.

Coarse-grained modeling of polystyrene-modified CNTs and their interactions with lipid bilayers

In the present work, we describe Martini3 coarse-grained models of polystyrene and carboxyl-terminated polystyrene functionalized carbon nanotubes (CNTs) and investigate their interactions with lipid bilayers with and without cholesterol (CHOL) using molecular dynamics simulations. By changing the polystyrene chain length and grafting density at the end ring of the CNTs at two different nanotube concentrations, we observe the translocation of nanoparticles as well as changes in the lipid bilayer properties. Our results show that all developed models passively diffuse into the membranes without causing any damage to the membrane integrity, although high concentrations of CNTs induce structural and elastic changes in lipid bilayers. In the presence of CHOL, increasing CNT concentration results in decreased rates of CHOL transmembrane motions. On the other hand, CNTs are prone to lipid and polystyrene blockage, which affects their equilibrated configurations, and tilting behavior within the membranes. Hence, we demonstrate that polystyrene-functionalized CNTs are promising drug-carrier agents. However, polystyrene chain length and grafting density are important factors to consider to enhance the efficiency of drug delivery.

Production of Lipid Constructs by Design via Three-Dimensional Nanoprinting

Atomic force microscopy (AFM) in conjunction with microfluidic delivery was utilized to produce three-dimensional (3D) lipid structures following a custom design. While AFM is well-known for its spatial precision in imaging and 2D nanolithography, the development of AFM-based nanotechnology into 3D nanoprinting requires overcoming the technical challenges of controlling material delivery and interlayer registry. This work demonstrates the concept of 3D nanoprinting of amphiphilic molecules such as 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). Various formulations of POPC solutions were tested to achieve point, line, and layer-by-layer material delivery. The produced structures include nanometer-thick disks, long linear spherical caps, stacking grids, and organizational chiral architectures. The POPC molecules formed stacking bilayers in these constructions, as revealed by high-resolution structural characterizations. The 3D printing reached nanometer spatial precision over a range of 0.5 mm. The outcomes reveal the promising potential of our designed technology and methodology in the production of 3D structures from nanometer to continuum, opening opportunities in biomaterial sciences and engineering, such as in the production of 3D nanodevices, chiral nanosensors, and scaffolds for tissue engineering and regeneration.

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.

Polystyrene-modified carbon nanotubes: Promising carriers in targeted drug delivery

To design drug-delivery agents for therapeutic and diagnostic applications, understanding the mechanisms by which covalently functionalized carbon nanotubes penetrate and interact with cell membranes is of great importance. Here, we report all-atom molecular dynamics results from polystyrene and carboxyl-terminated polystyrene-modified carbon nanotubes and show their translocation behavior across a model lipid bilayer together with their potential to deliver a molecule of the drug ibuprofen into the cell. Our results indicate that functionalized carbon nanotubes are internalized by the membrane in hundreds of nanoseconds and that drug loading increases the internalization speed further. Both loaded and unloaded tubes cross the closest leaflet of the bilayer by nonendocytic pathways, and for the times studied, the drug molecule remains trapped inside the pristine tube while remaining attached at the end of polystyrene-modified tube. On the other hand, carboxyl-terminated polystyrene functionalization allows the drug to be completely released into the lower leaflet of the bilayer without imposing damage to the membrane. This study shows that polystyrene functionalization is a promising alternative and facilitates drug delivery as a benchmark case.

Production of novel SARS-CoV-2 Spike truncations in Chinese hamster ovary cells leads to high expression and binding to antibodies

SARS-CoV-2 Spike is a key protein that mediates viral entry into cells and elicits antibody responses. Its importance in infection, diagnostics, and vaccinations has created a large demand for purified Spike for clinical and research applications. Spike is difficult to express, prompting modifications to the protein and expression platforms to improve yields. Alternatively, the Spike receptor-binding domain (RBD) is commonly expressed with higher titers, though it has lower sensitivity in serological assays. Here, we improve transient Spike expression in Chinese hamster ovary (CHO) cells. We demonstrate that Spike titers increase significantly over the expression period, maximizing at 14 mg L-1 on day 7. In comparison, RBD titers peak at 54 mg L-1 on day 3. Next, we develop eight Spike truncations (T1-T8) in pursuit of truncation with high expression and antibody binding. The truncations T1 and T4 express at 130 and 73 mg L-1 , respectively, which are higher than our RBD titers. Purified proteins were evaluated for binding to antibodies raised against full-length Spike. T1 has similar sensitivity as Spike against a monoclonal antibody and even outperforms Spike for a polyclonal antibody. These results suggest that T1 is a promising Spike alternative for use in various applications.

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.

Impact of Surface Polarity on Lipid Assembly under Spatial Confinement

Molecular dynamics (MD) simulations in the MARTINI model are used to study the assembly of 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) molecules under spatial confinement, such as during solvent evaporation from ultrasmall (femtoliter quantity) droplets. The impact of surface polarity on molecular assembly is discussed in detail. To the best of our knowledge, this work represents the first of its kind. Our results reveal that solvent evaporation gives rise to the formation of well-defined stacks of lipid bilayers in a smectic alignment. These smectic mesophases form on both polar and nonpolar surfaces but with a notable distinction. On polar surfaces, the director of the stack is oriented perpendicular to the support surface. By contrast, the stacks orient at an angle on the nonpolar surfaces. The packing of head groups on surfaces and lipid molecular mobility exhibits significant differences as surface polarity changes. The role of glycerol in the assembly and stability is also revealed. The insights revealed from the simulation have a significant impact on additive manufacturing, biomaterials, model membranes, and engineering protocells. For example, POPC assemblies via evaporation of ultrasmall droplets were produced and characterized. The trends compare well with the bilayer stack models. The surface polarity influences the local morphology and structures at the interfaces, which could be rationalized via the molecule-surface interactions observed from simulations.

A Computational Algorithm to Assess the Physiochemical Determinants of T Cell Receptor Dissociation Kinetics

The rational design of T Cell Receptors (TCRs) for immunotherapy has stagnated due to a limited understanding of the dynamic physiochemical features of the TCR that elicit an immunogenic response. The physiochemical features of the TCR-peptide major histocompatibility complex (pMHC) bond dictate bond lifetime which, in turn, correlates with immunogenicity. Here, we: i) characterize the force-dependent dissociation kinetics of the bond between a TCR and a set of pMHC ligands using Steered Molecular Dynamics (SMD); and ii) implement a machine learning algorithm to identify which physiochemical features of the TCR govern dissociation kinetics. Our results demonstrate that the total number of hydrogen bonds between the CDR2β-MHC⍺(β), CDR1α-Peptide, and CDR3β-Peptide are critical features that determine bond lifetime.

An unprecedented fully reduced {MoV 60} polyoxometalate: from an all-inorganic molecular light-absorber model to improved photoelectronic performance

Fully reduced polyoxometalates are predicted to give rise to a broad and strong absorption spectrum, suitable energy levels, and unparalleled electronic and optical properties. However, they are not available to date. Here, an unprecedented fully reduced polyoxomolybdate cluster, namely Na8[MoV 60O140(OH)28]·19H2O {MoV 60}, was successfully designed and obtained under hydrothermal conditions, which is rare and is the largest fully reduced polyoxometalate reported so far. The MoV 60 molecule describes one Keggin {ε-Mo12} encapsulated in an unprecedented {Mo24} cage, giving rise to a double truncated tetrahedron quasi-nesting architecture, which is further face-capped by another four {Mo6} tripods. Its crystalline stability in air, solvent tolerance, and photosensitivity were all shown. As a cheap and robust molecular light-absorber model possessing wide light absorption, MoV 60 was applied to build a co-sensitized solar cell photoelectronic device along with N719 dyes and the optimal power conversion efficiency was 28% higher than that of single-dye sensitization. These results show that MoV 60 polyoxometalate could serve as an ideal model for the design and synthesis of all-inorganic molecular light-absorbers for other light-driven processes in the future.

Using Molecular Dynamics Simulations to Interrogate T Cell Receptor Non-Equilibrium Kinetics

•

Insights into the atomic-scale interaction of the T Cell Receptor with the peptide Major Histocompatibility Complex.

•

Investigation of the physiochemical features that correspond with T Cell Receptor recognition during dynamic dissociation.

•

Implications of force-dependent non-equilibrium kinetics on T Cell Receptor mechanosensing.

An atomic-scale mechanism of T Cell Receptor (TCR) mechanosensing of peptides in the binding groove of the peptide-major histocompatibility complex (pMHC) may inform the design of novel TCRs for immunotherapies. Using steered molecular dynamics simulations, our study demonstrates that mutations to peptides in the binding groove of the pMHC – which are known to discretely alter the T cell response to an antigen – alter the MHC conformation at equilibrium. This subsequently impacts the overall strength (duration and length) of the TCR-pMHC bond under constant load. Moreover, physiochemical features of the TCR-pMHC dynamic bond strength, such as hydrogen bonds and Lennard-Jones contacts, correlate with the immunogenic response elicited by the specific peptide in the MHC groove. Thus, formation of transient TCR-pMHC bonds is characteristic of immunogenic peptides, and steered molecular dynamics simulations can be used in the overall design strategy of TCRs for immunotherapies.

Comparing the Expense and Accuracy of Methods to Simulate Atomic Vibrations in Rubrene

Atomic vibrations can inform about materials properties from hole transport in organic semiconductors to correlated disorder in metal-organic frameworks. Currently, there are several methods for predicting these vibrations using simulations, but the accuracy-efficiency tradeoffs have not been examined in depth. In this study, rubrene is used as a model system to predict atomic vibrational properties using six different simulation methods: density functional theory, density functional tight binding, density functional tight binding with a Chebyshev polynomial-based correction, a trained machine learning model, a pretrained machine learning model called ANI-1, and a classical forcefield model. The accuracy of each method is evaluated by comparison to the experimental inelastic neutron scattering spectrum. All methods discussed here show some accuracy across a wide energy region, though the Chebyshev-corrected tight-binding method showed the optimal combination of high accuracy with low expense. We then offer broad simulation guidelines to yield efficient, accurate results for inelastic neutron scattering spectrum prediction.

SARS-CoV-2 spike binding to ACE2 is stronger and longer ranged due to glycan interaction

Highly detailed steered molecular dynamics simulations are performed on differently glycosylated receptor binding domains of the severe acute respiratory syndrome coronavirus-2 spike protein. The binding strength and the binding range increase with glycosylation. The interaction energy rises very quickly when pulling the proteins apart and only slowly drops at larger distances. We see a catch-slip-type behavior whereby interactions during pulling break and are taken over by new interactions forming. The dominant interaction mode is hydrogen bonds, but Lennard-Jones and electrostatic interactions are relevant as well.

Contributions of the international plant science community to the fight against infectious diseases in humans-part 2: Affordable drugs in edible plants for endemic and re-emerging diseases

The fight against infectious diseases often focuses on epidemics and pandemics, which demand urgent resources and command attention from the health authorities and media. However, the vast majority of deaths caused by infectious diseases occur in endemic zones, particularly in developing countries, placing a disproportionate burden on underfunded health systems and often requiring international interventions. The provision of vaccines and other biologics is hampered not only by the high cost and limited scalability of traditional manufacturing platforms based on microbial and animal cells, but also by challenges caused by distribution and storage, particularly in regions without a complete cold chain. In this review article, we consider the potential of molecular farming to address the challenges of endemic and re-emerging diseases, focusing on edible plants for the development of oral drugs. Key recent developments in this field include successful clinical trials based on orally delivered dried leaves of Artemisia annua against malarial parasite strains resistant to artemisinin combination therapy, the ability to produce clinical-grade protein drugs in leaves to treat infectious diseases and the long-term storage of protein drugs in dried leaves at ambient temperatures. Recent FDA approval of the first orally delivered protein drug encapsulated in plant cells to treat peanut allergy has opened the door for the development of affordable oral drugs that can be manufactured and distributed in remote areas without cold storage infrastructure and that eliminate the need for expensive purification steps and sterile delivery by injection.

Contributions of the international plant science community to the fight against human infectious diseases - part 1: epidemic and pandemic diseases

Infectious diseases, also known as transmissible or communicable diseases, are caused by pathogens or parasites that spread in communities by direct contact with infected individuals or contaminated materials, through droplets and aerosols, or via vectors such as insects. Such diseases cause ˜17% of all human deaths and their management and control places an immense burden on healthcare systems worldwide. Traditional approaches for the prevention and control of infectious diseases include vaccination programmes, hygiene measures and drugs that suppress the pathogen, treat the disease symptoms or attenuate aggressive reactions of the host immune system. The provision of vaccines and biologic drugs such as antibodies is hampered by the high cost and limited scalability of traditional manufacturing platforms based on microbial and animal cells, particularly in developing countries where infectious diseases are prevalent and poorly controlled. Molecular farming, which uses plants for protein expression, is a promising strategy to address the drawbacks of current manufacturing platforms. In this review article, we consider the potential of molecular farming to address healthcare demands for the most prevalent and important epidemic and pandemic diseases, focussing on recent outbreaks of high-mortality coronavirus infections and diseases that disproportionately affect the developing world.

Porous Aromatic Framework Nanosheets Anchored with Lewis Pairs for Efficient and Recyclable Heterogeneous Catalysis

Lewis pairs (LPs) with outstanding performance for nonmetal-mediated catalysis reactions have high fundamental interest and remarkable application prospects. However, their solubility characteristics lead to instability and deactivation upon recycling. Here, the layered porous aromatic framework (PAF-6), featuring two kinds of Lewis base sites (NPiperazine and NTriazine), is exfoliated into few-layer nanosheets to form the LP entity with the Lewis acid. After comparison with various porous networks and verification by density functional theory (DFT) calculations, the NTriazine atom in the specific spatial environment is determined to preferably coordinate with the electron-deficient boron compound in a sterically hindered pattern. LP-bare porous product displays high catalytic activity for the hydrogenation of both olefin and imine compounds, and demonstrates ≈100% activity after 10 successful cycles in hydrogenation reactions. Considering the natural advantage of porous organic frameworks to construct LP groups opens up novel prospects for preparing other nonmetallic heterogeneous catalysts for efficient and recyclable catalysis.

Development and simulation of fully glycosylated molecular models of ACE2-Fc fusion proteins and their interaction with the SARS-CoV-2 spike protein binding domain

We develop fully glycosylated computational models of ACE2-Fc fusion proteins which are promising targets for a COVID-19 therapeutic. These models are tested in their interaction with a fragment of the receptor-binding domain (RBD) of the Spike Protein S of the SARS-CoV-2 virus, via atomistic molecular dynamics simulations. We see that some ACE2 glycans interact with the S fragments, and glycans are influencing the conformation of the ACE2 receptor. Additionally, we optimize algorithms for protein glycosylation modelling in order to expedite future model development. All models and algorithms are openly available.

Porous Organic Frameworks Featured by Distinct Confining Fields for the Selective Hydrogenation of Biomass-Derived Ketones

The asymmetric hydrogenation of biomass-derived molecules for the preparation of single enantiomer compounds is an effective method to reduce the rapid consumption of fossil resources. Porous organic frameworks (POFs) with pure organic surfaces may provide unusual confinement effects for organic substrates in chiral catalysis. Here, a series of POF catalysts are designed with chiral active centers decorated into sharply defined one-dimensional channels with diameters in the range of 1.2-2.9 nm. Due to the synergistic effect originating from the conjugated inner wall, the POF material (aperture size 2.4 nm) concentrates over 90% of aromatic species into the porous architecture, and its affinity is one or two orders of magnitude higher than those of classical porous solids. As determined by PBE+D3 calculation, the phenyl fragment reveals strong π-π interaction for steric hindrance around the metal active site to achieve stronger asymmetric induction. Therefore, this POF catalyst achieves high conversion (>99% yield) and enantioselectivity (>99% ee) for various substrates. The advantages of using the POF platform as a chiral catalyst can provide new perspectives on POF-based solid-state host-guest chemistry and asymmetric heterogeneous catalysis.

Confining Liquids inside Carbon Nanotubes: Accelerated Molecular Dynamics with Spliced, Soft-Core Potentials and Simulated Annealing

Understanding emergent phenomena of fluids under physical confinement requires the development of advanced tools for rapid and accurate simulation of their physiochemical properties. Simulating liquid molecules commensurate in size with the nanoscale enclosures that confine them is a key challenge. We demonstrate an accelerated molecular dynamics simulation technique that combines soft-core potentials (SCP) and simulated annealing (SA) to analyze confined liquids. This integrated SCP/SA method relies on a new spliced soft-core potential (SSCP), which enables tunable accuracy with respect to the target hard-core potential (HCP). SCP/SA enables the packing of enclosures with bulk material in a controlled, thermodynamically consistent manner. The enhanced SSCP accuracy is a critical feature of SCP/SA, enabling a smooth transition between the SCP and the HCP at a desired SCP hardness. We applied SCP/SA to the problem of filling a carbon nanotube (CNT) in periodic boundary conditions with a popular ionic liquid (IL), 1-butyl-3-methylimidazolium hexafluorophosphate [BMIM⁺][PF₆^-]. We performed a series of triplicate simulations on systems with varying CNT diameter and charge to demonstrate SCP/SA's versatility. Beyond this IL/CNT system, the SCP/SA simulation framework has a broad range of potential applications, not limited to nanoscale enclosures and interfaces, including both solid-state and biological systems.

New Means to Control Molecular Assembly

While self-assembly of molecules is relatively well-known and frequently utilized in chemical synthesis and material science, controlled assembly of molecules represents a new concept and approach. The present work demonstrates the concept of controlled molecular assembly using a non-spherical biomolecule, heparosan tetrasaccharide (MW = 1.099 kD). The key to controlled assembly is the fact that ultra-small solution droplets exhibit different evaporation dynamics from those of larger ones. Using an independently controlled microfluidic probe in an atomic force microscope, sub-femtoliter aqueous droplets containing designed molecules produce well-defined features with dimensions as small as tens of nanometers. The initial shape of the droplet and the concentration of solute within the droplet dictate the final assembly of molecules due to the ultrafast evaporation rate and dynamic spatial confinement of the droplets. The level of control demonstrated in this work brings us closer to programmable synthesis for chemistry and materials science which can be used to develop vehicles for drug delivery three-dimensional nanoprinting in additive manufacturing.

Harnessed Dopant Block Copolymers Assist Decorating Membrane Pores: A Dissipative Particle Dynamics Study

Self-assembly of asymmetric block copolymers (BCPs) around active pore edges has emerged as an important strategy to produce smart membranes with tunable pathways for solute transport. However, thus far, it is still challenging to manipulate pore shape and functionality for directional transformation under external stimuli. Here, a versatile strategy by mesoscale simulations to design stimuli-responsive pores with various edge decorations in hybrid membranes is reported. Dopant BCPs are used as decorators to stabilize pore edges and extend their function in reconfiguring pores in response to repeated membrane stretching/shrinking caused by external stimuli. The decoration morphologies are predictable since the assemblies of dopant BCPs around pore edges are closely related to their self-assemblies in solution. The coassembly between different BCPs in the hybrid membrane for the control of pore morphology is featured, and the parameter settings, including block incompatibility and molecular architecture for the construction of a specific pore, are determined. Results show that harnessed dopant BCPs in the hybrid membrane can enhance pore formation and induce directional pore shape and functionality transformation. Diversified pore decorations exhibit potential that can be further explored in selective solute transport and the design of stimuli-responsive smart nanodevices.

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 (EC₅₀) 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 EC₅₀ 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.

Computational modelling of atomic layer etching of chlorinated germanium surfaces by argon

The atomic layer etching of chlorinated germanium surfaces under argon bombardment was simulated using molecular dynamics with a newly fitted Tersoff potential. The chlorination energy determines the threshold energy for etching and the number of etched atoms in the bombardment phase. Etch rate is determined by bombardment energy.

Atomistic modeling of La3+ doping segregation effect on nanocrystalline yttria-stabilized zirconia

The effect of La3+ doping on the structure and ionic conductivity change in nanocrystalline yttria-stabilized zirconia (YSZ) was studied using a combination of Monte Carlo and molecular dynamics simulations. The simulation revealed the segregation of La3+ at eight tilt grain boundary (GB) structures and predicted an average grain boundary (GB) energy decrease of 0.25 J m-2, which is close to the experimental values reported in the literature. Cation stabilization was found to be the main reason for the GB energy decrease, and energy fluctuations near the grain boundary are smoothed out with La3+ segregation. Both dynamic and energetic analysis on the Σ13(510)/[001] GB structure revealed La3+ doping hinders O2- diffusion in the GB region, where the diffusion coefficient monotonically decreases with increasing La3+ doping concentration. The effect was attributed to the increase in the site-dependent migration barriers for O2- hopping caused by segregated La3+, which also leads to anisotropic diffusion at the GB.

Molecular Dynamics Modeling of Methylene Blue-DOPC Lipid Bilayer Interactions

We present a coarse-grained MARTINI model for methylene blue (MB) and investigate the interactions of MB with dioleylphosphatidylcholine (DOPC) lipid bilayers by molecular dynamics simulations. Our results show that the charge state of MB and the oxidation degree of the DOPC bilayer play critical roles on membrane properties. Oxidation of the DOPC bilayer significantly increases permeability of water and MB molecules, irrespective of the charge state of MB. The most significant changes in membrane properties are obtained for peroxidized lipid bilayers in the presence of cationic MB, with ∼11% increase in the membrane area per lipid head group and ∼7 and 44% reduction in membrane thickness and lateral diffusivity, respectively.

Controllable Multigeometry Nanoparticles via Cooperative Assembly of Amphiphilic Diblock Copolymer Blends with Asymmetric Architectures

Multigeometry nanoparticles with high complexity in composition and structure have attracted significant attention for enhanced functionality. We assess a simple but versatile strategy to construct hybrid nanoparticles with subdivided geometries through the cooperative assembly of diblock copolymer blends with asymmetric architectures. We report the formation of multicompartmental, vesicular, cylindrical, and spherical structures from pure AB systems. Then, we explore the assemblies of binary AB/AC blends, where the two incompatible, hydrophobic diblock copolymers subdivide into self-assembled local geometries, and the complexity of the obtained morphologies increases. We expand the strategy to ternary AB/AC/AD systems by tuning the effect of phase separation of different hydrophobic domains on the surface or internal region of the nanoparticle. The kinetic control of the coassembly in the initial stage is crucial for controlling the final morphology. The interactions of copolymers with different block lengths and chemistries enable the stabilization of interfaces, rims and ends of subdomains in the hybrid multigeometry nanoparticles. With further exploration of size and shape, the dependence of local geometry on the volume fraction is discussed. We show an efficient approach for controllable multigeometry nanoparticle construction that will be useful for multifunctional and hierarchical nanomaterials.

Structural analysis of human glycoprotein butyrylcholinesterase using atomistic molecular dynamics: The importance of glycosylation site ASN241

Human butyrylcholinesterase (BChE) is a glycoprotein capable of bioscavenging toxic compounds such as organophosphorus (OP) nerve agents. For commercial production of BChE, it is practical to synthesize BChE in non–human expression systems, such as plants or animals. However, the glycosylation profile in these systems is significantly different from the human glycosylation profile, which could result in changes in BChE’s structure and function. From our investigation, we found that the glycan attached to ASN₂₄₁ is both structurally and functionally important due to its close proximity to the BChE tetramerization domain and the active site gorge. To investigate the effects of populating glycosylation site ASN₂₄₁, monomeric human BChE glycoforms were simulated with and without site ASN₂₄₁ glycosylated. Our simulations indicate that the structure and function of human BChE are significantly affected by the absence of glycan 241.

Development and Application of a Coarse-Grained Model for PNIPAM by Iterative Boltzmann Inversion and Its Combination with Lattice Boltzmann Hydrodynamics

The polymer poly(N-isopropylacrylamide) (PNIPAM) is studied using a novel combination of multiscale modeling methodologies. We develop an iterative Boltzmann inversion potential of concentrated PNIPAM solutions and combine it with lattice Boltzmann as a Navier-Stokes equation solver for the solvent. We study in detail the influence of the methodology on statics and dynamics of the system. The combination is successful and significantly simpler and faster than other mapping techniques for polymer solution while keeping the correct hydrodynamics. The model can semiquantitatively describe the correct phase behavior and polymer dynamics.

Controllable multicompartment morphologies from cooperative self-assembly of copolymer-copolymer blends

Multicompartment nanostructures, such as microcapsules with clearly separated shell and core, are not easily accessible by conventional block copolymer self-assembly. We assess a versatile computational strategy through cooperative assembly of diblock copolymer blends to generate spherical and cylindrical compartmentalized micelles with intricate structures and morphologies. The co-assembly strategy combines the advantages of polymer blending and incompatibility-induced phase separation. Following this strategy, various nanoassemblies of pure AB, binary AB/AC and ternary AB/AC/AD systems such as compartmentalized micelles with sponge-like, Janus, capsule-like and onion-like morphologies can be obtained. The formation and structural adjustment of microcapsule micelles, in which the shell or core can be occupied by either pure or mixed diblock copolymers, were explored. The mechanism involving the separation of shell and core copolymers is attributed to the stretching force differences of copolymers which drive the arrangement of different copolymers in a pathway to minimize the total interfacial energy. Moreover, by adjusting block interactions, an efficient approach is exhibited for regulating the shell or core composition and morphology in microcapsule micelles, such as the transition from the "pure shell/mixed core" morphology to the "mixed shell/pure core" morphology in the AB/AC/AD micelle. This mesoscale simulation study identifies the key factors governing co-assembly of diblock copolymer blends and provides bottom-up insights towards the design and optimization of new multicompartment micelles.

Tunable Permeability of Cross-Linked Microcapsules from pH-Responsive Amphiphilic Diblock Copolymers: A Dissipative Particle Dynamics Study

Using dissipative particle dynamics simulation, we probe the tunable permeability of cross-linked microcapsules made from pH-sensitive diblock copolymers poly(ethylene oxide)-b-poly(N,N-diethylamino-2-ethyl methacrylate) (PEO-b-PDEAEMA). We first examine the self-assembly of non-cross-linked microcapsules and their pH-responsive collapse and then explore the effects of cross-linking and block interaction on the swelling or deswelling of cross-linked microcapsules. Our results reveal a preferential loading of hydrophobic dicyclopentadiene (DCPD) molecules in PEO-b-PDEAEMA copolymers. Upon reduction of pH, non-cross-linked microcapsules fully decompose into small wormlike clusters as a result of large self-repulsions of protonated copolymers. With increasing degree of cross-linking, the morphology of the microcapsule becomes more stable to pH change. The highly cross-linked microcapsule shell undergoes significant local polymer rearrangement in acidic solution, which eliminates the amphiphilicility and therefore enlarges the permeability of the shell. The responsive cross-linked shell experiences a disperse-to-buckle configurational transition upon reduction of pH, which is effective for the steady or pulsatile regulation of shell permeability. The swelling rate of the cross-linked shell is dependent on both electrostatic and nonelectrostatic interactions between the pH-sensitive groups as well as the other groups. Our study highlights the combination of cross-linking structure and block interactions in stabilizing microcapsules and tuning their selective permeability.

Molecular dynamics simulations of ternary lipid bilayers containing plant sterol and glucosylceramide

An atomic-level molecular dynamics simulation was carried out to study the effects of a plant sterol (sitosterol) and glucosylceramide (GlcCer) on a 1,2-dilinoleoylposphocholine (DLiPC) membrane. Initially, a membrane containing 50mol% sitosterol was compared with that containing the same ratio of cholesterol. These simulations showed differential condensing and ordering effects of sitosterol and cholesterol, with cholesterol being slightly more efficient than sitosterol in packing the membrane more tightly to a liquid ordered phase. By incorporation of 9.3% GlcCer on DLiPC/sterol membrane no notable change was observed in terms of area per lipid, bilayer thickness, order parameter and lateral diffusion. Some clusters of GlcCer/sterol were observed at higher ratio of GlcCer (15.5%), supporting the existence of GlcCer/sitosterol-enriched L_o-domains on the nanometer scale in the plant lipid mixture.

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 N₂, 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.

Understanding the Interaction of Pluronics L61 and L64 with a DOPC Lipid Bilayer: An Atomistic Molecular Dynamics Study

We investigate the interactions of Pluronics L61 and L64 with a dioleylphosphatidylcholine (DOPC) lipid bilayer by atomistic molecular dynamics simulations using the all-atom OPLS force field. Our results show that the initial configuration of the polymer with respect to the bilayer determines its final conformation within the bilayer. When the polymer is initially placed at the lipid/water interface, we observe partial insertion of the polymer in a U-shaped conformation. On the other hand, when the polymer is centered at the bilayer, it stabilizes to a transmembrane state, which facilitates water transport across the bilayer. We show that membrane thickness decreases while its fluidity increases in the presence of Pluronics. When the polymer concentration inside the bilayer is high, pore formation is initiated with L64. Our results show good agreement with existing experimental data and reveal that the hydrophilic/lipophilic balance of the polymer plays a critical role in the interaction mechanisms as well as in the dynamics of Pluronics with and within the bilayer.

Reactive Molecular Dynamics Simulations of the Silanization of Silica Substrates by Methoxysilanes and Hydroxysilanes

We perform reactive molecular dynamics simulations of monolayer formation by silanes on hydroxylated silica substrates. Solutions composed of alkylmethoxysilanes or alkylhydroxysilanes in hexane are placed in contact with a hydroxylated silica surface and simulated using a reactive force field (ReaxFF). In particular, we have modeled the deposition of butyl-, octyl-, and dodecyltrimethoxysilane to observe the dependence of alkylsilyl chain length on monolayer formation. We additionally modeled silanization using dodecyltrihydroxysilane, which allows for the comparison of two grafting mechanisms of alkoxysilanes: (1) direct condensation of alkoxysilane with surface-bound silanols and (2) a two-step hydrolysis-condensation mechanism. To emulate an infinite reservoir of reactive solution far away from the substrate, we have developed a method in which new precursor molecules are periodically added to a region of the simulation box located away from the surface. It is determined that the contact angle of alkyl tails bound to the surface is dependent on their grafting density. During the early stages of grafting alkoxy- and hydroxysilanes to the substrate, a preference is shown for silanes to condense with silanols further from the substrate surface and also close to neighboring surface-bound silanols. The kinetics of silica silanization by hydroxysilanes was observed to be much faster than for methoxysilanes. However, the as-deposited hydroxysilane monolayers show similar morphological characteristics to those formed by methoxysilanes.

Molecular modeling of lipid probes and their influence on the membrane

In this review a number of Molecular Dynamics simulation studies are discussed which focus on the understanding of the behavior of lipid probes in biomembranes. Experiments often use specialized probe moieties or molecules to report on the behavior of a membrane and try to gain information on the membrane as a whole from the probe lipids as these probes are the only things an experiment sees. Probes can be used to make NMR, EPR and fluorescence accessible to the membrane and use fluorescent or spin-active moieties for this purpose. Clearly membranes with and without probes are not identical which makes it worthwhile to elucidate the differences between them with detailed atomistic simulations. In almost all cases these differences are confined to the local neighborhood of the probe molecules which are sparsely used and generally present as single molecules. In general, the behavior of the bulk membrane lipids can be qualitatively understood from the probes but in most cases their properties cannot be directly quantitatively deduced from the probe behavior. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.