CG2AT2 : an enhanced fragment-based approach for serial multi-scale molecular dynamics simulations

Coarse-grained molecular dynamics provides a means for simulating the assembly and interactions of macromolecular complexes at a reduced level of representation, thereby allowing both longer timescale and larger sized simulations. Here, we describe an enhanced fragment-based protocol for converting macromolecular complexes from coarse-grained to atomistic resolution, for further refinement and analysis. While the focus is upon systems that comprise an integral membrane protein embedded in a phospholipid bilayer, the technique is also suitable for membrane-anchored and soluble protein/nucleotide complexes. Overall, this provides a method for generating an accurate and well-equilibrated atomic-level description of a macromolecular complex. The approach is evaluated using a diverse test set of 11 system configurations of varying size and complexity. Simulations are assessed in terms of protein stereochemistry, conformational drift, lipid/protein interactions, and lipid dynamics

A Pathway Model to Understand the Evolution of Spike Protein Binding to ACE2 in SARS-CoV-2 Variants

After the SARS-CoV-2 Wuhan variant that gave rise to the pandemic, other variants named Delta, Omicron, and Omicron-2 sequentially became prevalent, with mutations spread around the viral genome, including on the spike (S) protein; in order to understand the resultant in gains in infectivity, we interrogated in silico both the equilibrium binding and the binding pathway of the virus’ receptor-binding domain (RBD) to the angiotensin-converting enzyme 2 (ACE2) receptor. We interrogated the molecular recognition between the RBD of different variants and ACE2 through supervised molecular dynamics (SuMD) and classic molecular dynamics (MD) simulations to address the effect of mutations on the possible S protein binding pathways. Our results indicate that compensation between binding pathway efficiency and stability of the complex exists for the Omicron BA.1 receptor binding domain, while Omicron BA.2′s mutations putatively improved the dynamic recognition of the ACE2 receptor, suggesting an evolutionary advantage over the previous strains

Insights into membrane protein–lipid interactions from free energy calculations

Integral membrane proteins are regulated by specific interactions with lipids from the surrounding bilayer. The structures of protein–lipid complexes can be determined through a combination of experimental and computational approaches, but the energetic basis of these interactions is difficult to resolve. Molecular dynamics simulations provide the primary computational technique to estimate the free energies of these interactions. We demonstrate that the energetics of protein–lipid interactions may be reliably and reproducibly calculated using three simulation-based approaches: potential of mean force calculations, alchemical free energy perturbation, and well-tempered metadynamics. We employ these techniques within the framework of a coarse-grained force field and apply them to both bacterial and mammalian membrane protein–lipid systems. We demonstrate good agreement between the different techniques, providing a robust framework for their automated implementation within a pipeline for annotation of newly determined membrane protein structures

Garden varieties: how attractive are recommended garden plants to butterflies?

One way the public can engage in insect conservation is through wildlife gardening, including the growing of insect-friendly flowers as sources of nectar. However, plant varieties differ in the types of insects they attract. To determine which garden plants attracted which butterflies, we counted butterflies nectaring on 11 varieties of summer-flowering garden plants in a rural garden in East Sussex, UK. These plants were all from a list of 100 varieties considered attractive to British butterflies, and included the five varieties specifically listed by the UK charity Butterfly Conservation as best for summer nectar. A total of 2659 flower visits from 14 butterfly and one moth species were observed. We performed a principal components analysis which showed contrasting patterns between the species attracted to Origanum vulgare and Buddleia davidii. The “butterfly bush” Buddleia attracted many nymphalines, such as the peacock, Inachis io, but very few satyrines such as the gatekeeper, Pyronia tithonus, which mostly visited Origanum. Eupatorium cannibinum had the highest Simpson’s Diversity score of 0.75, while Buddleia and Origanum were lower, scoring 0.66 and 0.50 respectively. No one plant was good at attracting all observed butterfly species, as each attracted only a subset of the butterfly community. We conclude that to create a butterfly-friendly garden, a variety of plant species are required as nectar sources for butterflies. Furthermore, garden plant recommendations can probably benefit from being more precise as to the species of butterfly they attract

Can sacrificial feeding areas protect aquatic plants from herbivore grazing? Using behavioural ecology to inform wildlife management

Effective wildlife management is needed for conservation, economic and human well-being objectives. However, traditional population control methods are frequently ineffective, unpopular with stakeholders, may affect non-target species, and can be both expensive and impractical to implement. New methods which address these issues and offer effective wildlife management are required. We used an individual-based model to predict the efficacy of a sacrificial feeding area in preventing grazing damage by mute swans (Cygnus olor) to adjacent river vegetation of high conservation and economic value. The accuracy of model predictions was assessed by a comparison with observed field data, whilst prediction robustness was evaluated using a sensitivity analysis. We used repeated simulations to evaluate how the efficacy of the sacrificial feeding area was regulated by (i) food quantity, (ii) food quality, and (iii) the functional response of the forager. Our model gave accurate predictions of aquatic plant biomass, carrying capacity, swan mortality, swan foraging effort, and river use. Our model predicted that increased sacrificial feeding area food quantity and quality would prevent the depletion of aquatic plant biomass by swans. When the functional response for vegetation in the sacrificial feeding area was increased, the food quantity and quality in the sacrificial feeding area required to protect adjacent aquatic plants were reduced. Our study demonstrates how the insights of behavioural ecology can be used to inform wildlife management. The principles that underpin our model predictions are likely to be valid across a range of different resource-consumer interactions, emphasising the generality of our approach to the evaluation of strategies for resolving wildlife management problems

Incubation of ovine scrapie with environmental matrix results in biological and biochemical changes of PrPSc over time

Ovine scrapie can be transmitted via environmental reservoirs. A pool of ovine scrapie isolates were incubated on soil for one day or thirteen months and eluted prion was used to challenge tg338 mice transgenic for ovine PrP. After one-day incubation on soil, two PrPSc phenotypes were present: G338 or Apl338ii. Thirteen months later some divergent PrPSc phenotypes were seen: a mixture of Apl338ii with either G338 or P338, and a completely novel PrPSc deposition, designated Cag338. The data show that prolonged ageing of scrapie prions within an environmental matrix may result in changes in the dominant PrPSc biological/biochemical properties

“Hummingbird” floral traits interact synergistically to discourage visitation by bumble bee foragers

Pollination syndromes are suites of floral traits presumed to reflect adaptations to attract and utilize a “primary” type of animal pollinator. However, syndrome traits may also function to deter “secondary” flower visitors that reduce plant fitness through their foraging activities. Here we use the hummingbird-pollinated plant species Mimulus cardinalis as a model to investigate the potential deterrent effects of classic bird syndrome traits on bumble bee foragers. To establish that M. cardinalis flowers elicit an avoidance response in bees, we assessed the choice behavior of individual foragers on a mixed experimental array of M. cardinalis and its bee-pollinated sister species M. lewisii. As expected, bees showed a strong preference against M. cardinalis flowers (only 22% of total bee visits were to M. cardinalis), but surprisingly also showed a high degree of individual specialization (95.2% of total plant transitions were between conspecifics). To determine M. cardinalis floral traits that discourage bee visitation, we then assessed foraging responses of individuals to M. cardinalis-like and M. lewisii-like floral models differing in color, orientation, reward, and combinations thereof. Across experiments, M. cardinalis-like trait combinations consistently produced a higher degree of flower avoidance behavior and individual specialization than expected based on bee responses to each trait in isolation. We then conducted a series of flower discrimination experiments to assess the ability of bees to utilize traits and trait combinations associated with each species. Relative to M. lewisii-like alternatives, M. cardinalis-like traits alone had a minimal effect on bee foraging proficiency but together increased the time bees spent searching for rewarding flowers from 1.49 to 2.65 s per visit. Collectively, our results show that M. cardinalis flowers impose foraging costs on bumble bees sufficient to discourage visitation and remarkably, generate such costs through synergistic color-orientation and color-reward trait interactions. Floral syndromes therefore represent complex adaptations to multiple pollinator groups, rather than simply the primary pollinator

In silico and structural analyses demonstrate that intrinsic protein motions guide T cell receptor complementarity determining region loop flexibility

T-cell immunity is controlled by T cell receptor (TCR) binding to peptide major histocompatibility complexes (pMHCs). The nature of the interaction between these two proteins has been the subject of many investigations because of its central role in immunity against pathogens, cancer, in autoimmunity, and during organ transplant rejection. Crystal structures comparing unbound and pMHC-bound TCRs have revealed flexibility at the interaction interface, particularly from the perspective of the TCR. However, crystal structures represent only a snapshot of protein conformation that could be influenced through biologically irrelevant crystal lattice contacts and other factors. Here, we solved the structures of three unbound TCRs from multiple crystals. Superposition of identical TCR structures from different crystals revealed some conformation differences of up to 5 Å in individual complementarity determining region (CDR) loops that are similar to those that have previously been attributed to antigen engagement. We then used a combination of rigidity analysis and simulations of protein motion to reveal the theoretical potential of TCR CDR loop flexibility in unbound state. These simulations of protein motion support the notion that crystal structures may only offer an artifactual indication of TCR flexibility, influenced by crystallization conditions and crystal packing that is inconsistent with the theoretical potential of intrinsic TCR motions