Multiscale Molecular Dynamics Simulations of the Homodimer Accessory Protein ORF7b of SARS-CoV-2.

The SARS-CoV-2 ORF7b protein has drawn attention for its potential role in viral pathogenesis, but its structural details and lateral membrane associations remain elusive. In this study, we conducted multiscale molecular dynamics simulations to provide detailed molecular insights of the protein's dimerization, which is crucial for unraveling its structural model of protein-protein interface important to regulating cellular immune response. To gain a deeper understanding of homodimer configurations, we employed a machine learning algorithm for structural-based clustering. Clusters were categorized into three distinct groups for both parallel and antiparallel orientations, highlighting the influence of the initial monomer conformation on dimer configurations. Analysis of hydrogen bonding and π-π and π-cation stacking interactions within clusters revealed variations in interactions between clusters. In parallel dimers, weak stacking interactions in the transmembrane (TM) region were observed. In contrast, antiparallel dimers exhibited strong hydrogen bonding and stacking interactions contributing to tight dimeric packing, both within and outside the TM domain. Overall, our study provides a comprehensive view of the structural dynamics of ORF7b homodimerization in both parallel and antiparallel orientations. These findings shed light on the molecular interactions involved in ORF7b dimerization, which are crucial for understanding its potential roles in SARS-CoV-2 pathogenesis. This knowledge could inform future research and therapeutic strategies targeting this viral protein.

Understanding how transmembrane domains regulate interactions between human BST-2 and the SARS-CoV-2 accessory protein ORF7a.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID, replicates at intracellular membranes. Bone marrow stromal antigen 2 (BST-2; tetherin) is an antiviral response protein that inhibits transport of viral particles after budding within infected cells. RNA viruses such as SARS-CoV-2 use various strategies to disable BST-2, including use of transmembrane 'accessory' proteins that interfere with BST-2 oligomerization. ORF7a is a small, transmembrane protein present in SARS-CoV-2 shown previously to alter BST-2 glycosylation and function. In this study, we investigated the structural basis for BST-2 ORF7a interactions, with a particular focus on transmembrane and juxtamembrane interactions. Our results indicate that transmembrane domains play an important role in BST-2 ORF7a interactions and mutations to the transmembrane domain of BST-2 can alter these interactions, particularly single-nucleotide polymorphisms in BST-2 that result in mutations such as I28S. Using molecular dynamics simulations, we identified specific interfaces and interactions between BST-2 and ORF7a to develop a structural basis for the transmembrane interactions. Differences in glycosylation are observed for BST-2 transmembrane mutants interacting with ORF7a, consistent with the idea that transmembrane domains play a key role in their heterooligomerization. Overall, our results indicate that ORF7a transmembrane domain interactions play a key role along with extracellular and juxtamembrane domains in modulating BST-2 function.

Leaflet Asymmetry Modeling in the Lipid Composition of Escherichia coli Cytoplasmic Membranes.

Lipid composition asymmetry between leaflets is important to cell function and plays a key role in the "positive inside" rule in transmembrane proteins. In this work, Escherichia coli inner plasma membrane models reflecting this asymmetry have been investigated at the early-log and stationary stages during the bacterial lifecycle using all-atom molecular dynamics simulations. The CHARMM36 lipid force field is used, and selected membrane properties are tested for variations between two leaflets and whole membranes. Our models include bacterial lipids with a cyclopropane moiety on the sn-2 acyl chain in the stationary membrane model. The PE/PG ratio for two leaflets reflects the "positive inside" rule of membrane proteins, set to 6.8 and 2.8 for the inner and outer leaflets of the two models, respectively. We are the first to model leaflet asymmetry in the lipid composition of E. coli cytoplasmic membranes and observe the effect on membrane properties in leaflets and whole membranes. Specifically, our results show that for the stationary phase bilayer, the surface area per lipid (SA/lipid) is larger, the thickness (2DC and DB) is smaller, the tilt angle is larger, the tilt modulus is smaller, and the deuterium order parameters (SCD) of sn-1 and sn-2 tails are lower, compared to the early-log stage. Moreover, the stationary stage bilayer has a positive spontaneous curvature, while the early-log stage has a near flat spontaneous curvature. For leaflet asymmetry, the inner leaflet has a larger SA/lipid, a smaller thickness, a smaller elastic tilt modulus (a larger tilt angle), and lower SCD, compared to the outer leaflet in both stages. Moreover, an asymmetric membrane involves a lipid tilt and a lateral extension, varying from a reference state of a pre-equilibrium membrane. This work encourages a more profound exploration of leaflet asymmetry in various other membrane models and how this might affect the structure and function of membrane-associated peptides and proteins.

All-Atom Modeling of Complex Cellular Membranes.

Cell membranes are composed of a variety of lipids and proteins where they interact with each other to fulfill their roles. The first step in modeling these interactions in molecular simulations is to have reliable mimetics of the membrane's lipid environment. This Feature Article presents our recent efforts to model complex cellular membranes using all-atom force fields. A short review of the CHARMM36 (C36) lipid force field and its recent update to incorporate the long-range dispersion is presented. Key examples of model membranes mimicking various species and organelles are given. These include single-celled organisms such as bacteria (E. coli., chlamydia, and P. aeruginosa) and yeast (plasma membrane, endoplasmic reticulum, and trans-Golgi network) and more advanced ones such as plants (soybean and Arabidopsis thaliana) and mammals (ocular lens, stratum corneum, and peripheral nerve myelin). Leaflet asymmetry in composition has also been applied to some of these models. With the increased lipid diversity in the C36 lipid FF, these complex models can better reflect the structural, mechanical, and dynamic properties of realistic membranes and open an opportunity to study biological processes involving other molecules.

Early structural development in melt-quenched polymer PTT from atomistic molecular dynamic simulations.

Molecular dynamics simulations are performed to study the initial structural development in poly(trimethylene terephthalate) (PTT) when quenched below its melting point. The development of local ordering has been observed in our simulations. The thermal properties, such as the glass transition temperature (T(g)) and the melting temperature (T(m)), determined from our simulations are in reasonable agreement with experimental values. It is found that, between these two temperatures, the number of local structures quickly increases during the thermal relaxation period soon after the system is quenched and starts to fluctuate afterwards. The formation and development of local structures is found to be driven mainly by the torsional and van der Waals forces and follows the classical nucleation-growth mechanism. The variation of local structures' fraction with temperature exhibits a maximum between T(g) and T(m), resembling the temperature dependence of the crystallization rate for most polymers. In addition, the backbone torsion distribution for segments within the local structures preferentially reorganizes to the trans-gauche-gauche-trans (t-g-g-t) conformation, the same as that in the crystalline state. As a consequence, we believe that such local structural ordering could be the baby nuclei that have been suggested to form in the early stage of polymer crystallization.