SARS-CoV-2 omicron spike simulations: broad antibody escape, weakened ACE2 binding, and modest furin cleavage.

The recent emergence of the omicron variant of the SARS-CoV-2 virus with large numbers of mutations has raised concern about a potential new surge in infections. Here we use molecular dynamics to study the biophysics of the interface of the BA1 and BA2 omicron spike protein binding to (i) the ACE2 receptor protein, (ii) antibodies from all known binding regions, and (iii) the furin binding domain. Our simulations suggest that while there is a significant reduction of antibody (Ab) binding strength corresponding to escape, the omicron spikes pay a cost in terms of weaker receptor binding as measured by interfacial hydrogen bonds (H-bond). The furin cleavage domain (FCD) is the same or weaker binding than the delta variant, suggesting lower fusogenicity resulting in less viral load and disease intensity than the delta variant. IMPORTANCE The BA1 and BA2 and closely related BA2.12.2 and BA.5 omicron variants of SARS-CoV-2 dominate the current global infection landscape. Given the high number of mutations, particularly those which will lead to antibody escape, it is important to establish accurate methods that can guide developing health policy responses that identify at a fundamental level whether omicron and its variants are more threatening than its predecessors, especially delta. The importance of our work is to demonstrate that simple in silico simulations can predict biochemical binding details of the omicron spike protein that have epidemiological consequences, especially for binding to the cells and for fusing the viral membrane with the cells. In each case, we predicted weaker binding of the omicron spike, which agreed with subsequent experimental results. Future virology experiments will be needed to test these predictions further.

Using machine learning to detect coronaviruses potentially infectious to humans.

Establishing the host range for novel viruses remains a challenge. Here, we address the challenge of identifying non-human animal coronaviruses that may infect humans by creating an artificial neural network model that learns from spike protein sequences of alpha and beta coronaviruses and their binding annotation to their host receptor. The proposed method produces a human-Binding Potential (h-BiP) score that distinguishes, with high accuracy, the binding potential among coronaviruses. Three viruses, previously unknown to bind human receptors, were identified: Bat coronavirus BtCoV/133/2005 and Pipistrellus abramus bat coronavirus HKU5-related (both MERS related viruses), and Rhinolophus affinis coronavirus isolate LYRa3 (a SARS related virus). We further analyze the binding properties of BtCoV/133/2005 and LYRa3 using molecular dynamics. To test whether this model can be used for surveillance of novel coronaviruses, we re-trained the model on a set that excludes SARS-CoV-2 and all viral sequences released after the SARS-CoV-2 was published. The results predict the binding of SARS-CoV-2 with a human receptor, indicating that machine learning methods are an excellent tool for the prediction of host expansion events.

A biophysical model for plant cell plate maturation based on the contribution of a spreading force.

Plant cytokinesis, a fundamental process of plant life, involves de novo formation of a "cell plate" partitioning the cytoplasm of dividing cells. Cell plate formation is directed by orchestrated delivery, fusion of cytokinetic vesicles, and membrane maturation to form a nascent cell wall by timely deposition of polysaccharides. During cell plate maturation, the fragile membrane network transitions to a fenestrated sheet and finally a young cell wall. Here, we approximated cell plate sub-structures with testable shapes and adopted the Helfrich-free energy model for membranes, including a stabilizing and spreading force, to understand the transition from a vesicular network to a fenestrated sheet and mature cell plate. Regular cell plate development in the model was possible, with suitable bending modulus, for a two-dimensional late stage spreading force of 2-6 pN/nm, an osmotic pressure difference of 2-10 kPa, and spontaneous curvature between 0 and 0.04 nm-1. With these conditions, stable membrane conformation sizes and morphologies emerged in concordance with stages of cell plate development. To reach a mature cell plate, our model required the late-stage onset of a spreading/stabilizing force coupled with a concurrent loss of spontaneous curvature. Absence of a spreading/stabilizing force predicts failure of maturation. The proposed model provides a framework to interrogate different players in late cytokinesis and potentially other membrane networks that undergo such transitions. Callose, is a polysaccharide that accumulates transiently during cell plate maturation. Callose-related observations were consistent with the proposed model's concept, suggesting that it is one of the factors involved in establishing the spreading force.

Bottom-up synthesis of protein-based nanomaterials from engineered ß-solenoid proteins.

Biomolecular self-assembly is an emerging bottom-up approach for the synthesis of novel nanomaterials. DNA and viruses have both been used to create scaffolds but the former lacks chemical diversity and the latter lack spatial control. To date, the use of protein scaffolds to template materials on the nanoscale has focused on amyloidogenic proteins that are known to form fibrils or two-protein systems where a second protein acts as a cross-linker. We previously developed a unique approach for self-assembly of nanomaterials based on engineering ß-solenoid proteins (BSPs) to polymerize into micrometer-length fibrils. BSPs have highly regular geometries, tunable lengths, and flat surfaces that are amenable to engineering and functionalization. Here, we present a newly engineered BSP based on the antifreeze protein of the beetle Rhagium inquisitor (RiAFP-m9), which polymerizes into stable fibrils under benign conditions. Gold nanoparticles were used to functionalize the RiAFP-m9 fibrils as well as those assembled from the previously described SBAFP-m1 protein. Cysteines incorporated into the sequences provide site-specific gold attachment. Additionally, silver was deposited on the gold-labelled fibrils by electroless plating to create nanowires. These results bolster prospects for programable self-assembly of BSPs to create scaffolds for functional nanomaterials.

In silico stress-strain measurements on self-assembled protein lattices.

Due to their large mechanical strength and potential for functionalization, beta-solenoid proteins show promise as building blocks in biomaterials applications such as two- and three-dimensional scaffolds. We have designed simulation models of two-dimensional square and honeycomb protein lattices by covalently linking a beta-solenoid protein, the spruce budworm antifreeze protein (SBAFP), to symmetric protein multimers. Periodic boundary conditions applied to the simulation cell allow for the simulation of an infinite lattice. We use molecular dynamics to strain the lattice by deforming the simulation cell and measuring the resulting stress tensor. We evaluate the linear portion of stress-strain curves to extract the corresponding bulk and shear elastic moduli. When strained at a rate of 0.3 nm ps-1, the lattices yield a bulk modulus of approximately 3 GPa. This large elastic modulus demonstrates that 2-dimensional structures designed from beta-solenoid proteins can be expected to retain the exceptional material strength of their building blocks.

High Tensile Strength of Engineered ß-Solenoid Fibrils via Sonication and Pulling.

We present estimates of ultimate tensile strength (UTS) for two engineered ß-solenoid protein mutant fibril structures (spruce budworm and Rhagium inquisitor antifreeze proteins) derived from sonication-based measurements and from force pulling molecular dynamics simulations, both in water. Sonication experiments generate limiting scissioned fibrils with a well-defined length-to-width correlation for the mutant spruce budworm protein and the resultant UTS estimate is 0.66 ± 0.08 GPa. For fibrils formed from engineered R. inquisitor antifreeze protein, depending upon geometry, we estimate UTSs of 3.5 ± 3.2-5.5 ± 5.1 GPa for proteins with interfacial disulfide bonds, and 1.6 ± 1.5-2.5 ± 2.3 GPa for the reduced form. The large error bars for the R. inquisitor structures are intrinsic to the broad distribution of limiting scission lengths. Simulations provide pulling velocity-dependent UTSs increasing from 0.2 to 1 GPa in the available speed range, and 1.5 GPa extrapolated to the speeds expected in the sonication experiments. Simulations yield low-velocity values for the Young's modulus of 6.0 GPa. Without protein optimization, these mechanical parameters are similar to those of spider silk and Kevlar, but in contrast to spider silk, these proteins have a precisely known sequence-structure relationship.

Structures of Bacteroides fragilis uridine 5'-diphosphate-N-acetylglucosamine (UDP-GlcNAc) acyltransferase (BfLpxA).

Uridine 5'-diphosphate-N-acetylglucosamine (UDP-GlcNAc) acyltransferase (LpxA) catalyzes a reversible reaction for adding an O-acyl group to the GlcNAc in UDP-GlcNAc in the first step of lipid A biosynthesis. Lipid A constitutes a major component of lipopolysaccharides, also referred to as endotoxins, which form the outer monolayer of the outer membrane of Gram-negative bacteria. Ligand-free and UDP-GlcNAc-bound crystal structures of LpxA from Bacteroides fragilis NCTC 9343, the most common pathogenic bacteria found in abdominal abscesses, have been determined and are presented here. The enzyme crystallizes in a cubic space group, with the crystallographic threefold axis generating the biological functional homotrimer and with each monomer forming a nine-rung left-handed ß-helical (LßH) fold in the N-terminus followed by an α-helical motif in the C-terminus. The structure is highly similar to LpxA from other organisms. Yet, despite sharing a similar LßH structure with LpxAs from Escherichia coli and others, previously unseen calcium ions are observed on the threefold axis in B. fragilis LpxA to help stabilize the trimeric assembly.

In Silico Measurements of Twist and Bend Moduli for ß-Solenoid Protein Self-Assembly Units.

We compute potentials of mean force for bend and twist deformations via force pulling and umbrella sampling experiments for four ß-solenoid proteins (BSPs) that show promise in nanotechnology applications. In all cases, we find quasi-Hooke's law behavior until the point of rupture. Bending moduli show modest anisotropy for two-sided and three-sided BSPs, and little anisotropy for a four-sided BSP. There is a slight clockwise/counterclockwise asymmetry in the twist potential of mean force, showing greater stiffness when the applied twist follows the intrinsic twist. When we extrapolate to beam theory appropriate for amyloid fibrils of the BSPs, we find bend/twist moduli which are somewhat smaller than those in the literature for other amyloid fibrils. Twist persistence lengths are on the order of a micron, and bend persistence lengths are several microns. Provided the intrinsic twist can be reversed, these results support the usage of BSPs in biomaterials applications.

Simulated cytoskeletal collapse via tau degradation.

We present a coarse-grained two dimensional mechanical model for the microtubule-tau bundles in neuronal axons in which we remove taus, as can happen in various neurodegenerative conditions such as Alzheimers disease, tauopathies, and chronic traumatic encephalopathy. Our simplified model includes (i) taus modeled as entropic springs between microtubules, (ii) removal of taus from the bundles due to phosphorylation, and (iii) a possible depletion force between microtubules due to these dissociated phosphorylated taus. We equilibrate upon tau removal using steepest descent relaxation. In the absence of the depletion force, the transverse rigidity to radial compression of the bundles falls to zero at about 60% tau occupancy, in agreement with standard percolation theory results. However, with the attractive depletion force, spring removal leads to a first order collapse of the bundles over a wide range of tau occupancies for physiologically realizable conditions. While our simplest calculations assume a constant concentration of microtubule intercalants to mediate the depletion force, including a dependence that is linear in the detached taus yields the same collapse. Applying percolation theory to removal of taus at microtubule tips, which are likely to be the protective sites against dynamic instability, we argue that the microtubule instability can only obtain at low tau occupancy, from 0.06-0.30 depending upon the tau coordination at the microtubule tips. Hence, the collapse we discover is likely to be more robust over a wide range of tau occupancies than the dynamic instability. We suggest in vitro tests of our predicted collapse.

Sequence-dependent stability test of a left-handed ß-helix motif.

The left-handed ß-helix (LHBH) is an intriguing, rare structural pattern in polypeptides that has been implicated in the formation of amyloid aggregates. We used accurate all-atom replica-exchange molecular dynamics (REMD) simulations to study the relative stability of diverse sequences in the LHBH conformation. Ensemble-average coordinates from REMD served as a scoring criterion to identify sequences and threadings optimally suited to the LHBH, as in a fold recognition paradigm. We examined the repeatability of our REMD simulations, finding that single simulations can be reliable to a quantifiable extent. We find expected behavior for the positive and negative control cases of a native LHBH and intrinsically disordered sequences, respectively. Polyglutamine and a designed hexapeptide repeat show remarkable affinity for the LHBH motif. A structural model for misfolded murine prion protein was also considered, and showed intermediate stability under the given conditions. Our technique is found to be an effective probe of LHBH stability, and promises to be scalable to broader studies of this and potentially other novel or rare motifs. The superstable character of the designed hexapeptide repeat suggests theoretical and experimental follow-ups.

Massive conformation change in the prion protein: Using dual-basin structure-based models to find misfolding pathways.

We employ all-atom structure-based models with a force field with multiple energetic basins for the C-terminal (residues 166-226) of the mammalian prion protein. One basin represents the known alpha-helical (αH) structure while the other represents the same residues in a left-handed beta-helical (LHBH) conformation. The LHBH structure has been proposed to help describe one class of in vitro grown fibrils, as well as possibly self-templating the conversion of normal cellular prion protein to the infectious form. Yet, it is unclear how the protein may make this global rearrangement. Our results demonstrate that the conformation changes are not strongly limited by large-scale geometry modification and that there may exist an overall preference for the LHBH conformation. Furthermore, our model presents novel intermediate trapping conformations with twisted LHBH structure.

Left handed beta helix models for mammalian prion fibrils.

We propose models for in vitro grown mammalian prion protein fibrils based upon left handed beta helices formed both from the N-terminal and C-terminal regions of the proteinase resistant infectious prion core. The C-terminal threading onto a beta-helical structure is almost uniquely determined by fixing the cysteine disulfide bond on a helix corner. In comparison to known left handed helical peptides, the resulting model structures have similar stability attributes including relatively low root mean square deviations in all atom molecular dynamics, substantial side-chain-to-side-chain hydrogen bonding, good volume packing fraction, and low hydrophilic/hydrophobic frustration. For the N-terminus, we propose a new threading of slightly more than two turns, which improves upon the above characteristics relative to existing three turn beta-helical models. The N-terminal and C-terminal beta helices can be assembled into eight candidate models for the fibril repeat units, held together by large hinge (order 30 residues) domain swapping, with three amenable to fibril promoting domain swapping via a small (five residue) hinge on the N-terminal side. Small concentrations of the metastable C-terminal beta helix in vivo might play a significant role in templating the infectious conformation and in enhancing conversion kinetics for inherited forms of the disease and explain resistance (for canines) involving hypothesized coupling to the methionine 129 sulfur known to play a role in human disease.

Theoretical study of DNA damage recognition via electron transfer from the [4Fe-4S] complex of MutY.

The mechanism of site-specific recognition of DNA by proteins has been a long-standing issue. The DNA glycosylase MutY, for instance, must find the rare 8-oxoguanine-adenine mismatches among the large number of basepairs in the DNA. This protein has a [4Fe-4S] cluster, which is highly conserved in species as diverse as Escherichia Coli and Homo sapiens. The mixed-valent nature of this cluster suggests that charge transfer may play a role in MutY's function. We have studied the energetics of the charge transfer in Bacillus stearothermophilus MutY-DNA complex using multiscale calculation including density functional theory and molecular dynamics. The [4Fe-4S] cluster in MutY is found to undergo 2+ to 3+ oxidation when coupling to DNA through hole transfer, especially when MutY is near an oxoguanine modified base (oxoG). Employing the Marcus theory for electron transfer, we find near optimal Frank-Condon factors for electron transfer from MutY to oxoguanine modified base. MutY has modest selectivity for oxoguanine over guanine due to the difference in oxidation potential. The tunneling matrix element is significantly reduced with the mutation R149W, whereas the mutation L154F reduces the tunneling matrix element as well as the Frank-Condon factor. Both L154F and R149W mutations are known to dramatically reduce or eliminate repair efficiency. We suggest a scenario where the charge transfer leads to a stabilization of the specific binding conformation, which is likely the recognition mode, thus enabling it to find the damaged site efficiently.

A mechanism for copper inhibition of infectious prion conversion.

We employ ab initio electronic structure calculations to obtain two structural models for copper bound in the strongest binding site of the noninfectious form of the prion protein. The models are compatible with available experimental constraints from electron spin resonance data. The bending of the peptide backbone attendant with the copper binding is not compatible with the requisite straight beta-strand backbone structure for the same sequence contained in two recently proposed models of the prion protein structure in its infectious form. We hypothesize that copper binding at this site is protective against conversion to the infectious form, discuss experimental data that appear to support and conflict with our hypothesis, and propose tests using recombinant prion protein, genetically modified cultured neurons, and transgenic mice.

Prion disease: exponential growth requires membrane binding.

A hallmark feature of prions, whether in mammals or yeast and fungi, is exponential growth associated with fission or autocatalysis of protein aggregates. We have employed a rigorous kinetic analysis to recent data from transgenic mice lacking a glycosylphosphatidylinositol membrane anchor to the normal cellular PrP(C) protein, which show that toxicity requires the membrane binding. We find as well that the membrane is necessary for exponential growth of prion aggregates; without it, the kinetics is simply the quadratic-in-time growth characteristic of linear elongation as observed frequently in in vitro amyloid growth experiments with other proteins. This requires both: i), a substantial intercellular concentration of anchorless PrP(C), and ii), a concentration of small scrapies seeding aggregates from the inoculum, which remains relatively constant with time and exceeds the concentration of large polymeric aggregates. We also can explain via this analysis why mice heterozygous for the anchor-full/anchor-free PrP(C) proteins have more rapid incubation than mice heterozygous for anchor-full/null PrP(C), and contrast the mammalian membrane associated fission or autocatalysis with the membrane free fission of yeast and fungal prions.

Structure of infectious prions: stabilization by domain swapping.

A candidate structure for the minimal prion infectious unit is a recently discovered protein oligomer modeled as a beta-helical prion trimer (BPT); BPTs can stack to form cross-beta fibrils and may provide insight into protein aggregates of other amyloid diseases. However, the BPT lacks a clear intermonomer binding mechanism. Here we propose an alternative domain-swapped trimeric prion (DSTP) model and show with molecular dynamics (MD) that the DSTP has more favorable intermonomer hydrogen bonding and proline dihedral strain energy than the BPT. This new structural proposal may be tested by lysine and N terminus fluorescent resonance energy transfer (FRET) either directly on recombinant prion protein amyloid aggregates or on synthetic constructs that contain the proline/lysine-rich hinge region critical for domains to swap. In addition, the domain swapping may provide 1) intrinsic entanglement, which can contribute to the remarkable temperature stability of the infectious prion structure and help explain the absence of PrP(Sc) monomers, 2) insight into why specific prolines are potentially relevant to three inherited forms of prion disease, and 3) a simple explanation of prion strains assuming the strain is encoded in the monomer number of the oligomers.

Lateral organization of cholesterol molecules in lipid-cholesterol assemblies.

We present results of an off-lattice simulation of a two-component planar system, as a model for lateral organization of cholesterol molecules in lipid-cholesterol assemblies. We explore the existence of "superlattice" structures even in fluid systems, in the absence of an underlying translational long-range order, and study their coupling to hexatic or bond-orientational order. We discuss our results in context of geometric superlattice theories and "condensation complexes" in understanding a variety of experiments in artificial lipid-cholesterol assemblies.

Modeling amyloid beta-peptide insertion into lipid bilayers.

Inspired by recent suggestions that the Alzheimer's amyloid beta peptide (Abeta) can insert into cell membranes and form harmful ion channels, we model insertion of the 40- and 42-residue forms of the peptide into cell membranes using a Monte Carlo code which is specific at the amino acid level. We examine insertion of the regular Abeta peptide as well as mutants causing familial Alzheimer's disease, and find that all but one of the mutants change the insertion behavior by causing the peptide to spend more simulation steps in only one leaflet of the bilayer. We also find that Abeta42, because of the extra hydrophobic residues relative to Abeta40, is more likely to adopt this conformation than Abeta40 in both wild-type and mutant forms. We argue qualitatively why these effects happen. Here, we present our results and develop the hypothesis that this partial insertion increases the probability of harmful channel formation. This hypothesis can partly explain why these mutations are neurotoxic simply due to peptide insertion behavior. We further apply this model to various artificial Abeta mutants which have been examined experimentally, and offer testable experimental predictions contrasting the roles of aggregation and insertion with regard to toxicity of Abeta mutants. These can be used through further experiments to test our hypothesis.

Simulations of oligomeric intermediates in prion diseases.

We extend our previous stochastic cellular automata-based model for two-dimensional (areal) aggregation of prion proteins on neuronal surfaces. The new anisotropic model allows us to simulate both strong beta-sheet and weaker attachment bonds between proteins. Constraining binding directions allows us to generate aggregate structures with the hexagonal lattice symmetry found in recently observed in vitro experiments. We argue that these constraints on rules may correspond to underlying steric constraints on the aggregation process. We find that monomer-dominated growth of the areal aggregate is too slow to account for some observed doubling-time-to-incubation-time ratios inferred from data, and so consider aggregation dominated by relatively stable but noninfectious oligomeric intermediates. We compare a kinetic theory analysis of oligomeric aggregation to spatially explicit simulations of the process. We find that with suitable rules for misfolding of oligomers, possibly due to water exclusion by the surrounding aggregate, the resulting oligomeric aggregation model maps onto our previous monomer aggregation model. Therefore it can produce some of the same attractive features for the description of prion incubation time data. We propose experiments to test the oligomeric aggregation model.