The solute carrier SPNS2 recruits PI(4,5)P2 to synergistically regulate transport of sphingosine-1-phosphate.

Solute carrier spinster homolog 2 (SPNS2), one of only four known major facilitator superfamily (MFS) lysolipid transporters in humans, exports sphingosine-1-phosphate (S1P) across cell membranes. Here, we explore the synergistic effects of lipid binding and conformational dynamics on SPNS2's transport mechanism. Using mass spectrometry, we discovered that SPNS2 interacts preferentially with PI(4,5)P2. Together with functional studies and molecular dynamics (MD) simulations, we identified potential PI(4,5)P2 binding sites. Mutagenesis of proposed lipid binding sites and inhibition of PI(4,5)P2 synthesis reduce S1P transport, whereas the absence of the N terminus renders the transporter essentially inactive. Probing the conformational dynamics of SPNS2, we show how synergistic binding of PI(4,5)P2 and S1P facilitates transport, increases dynamics of the extracellular gate, and stabilizes the intracellular gate. Given that SPNS2 transports a key signaling lipid, our results have implications for therapeutic targeting and also illustrate a regulatory mechanism for MFS transporters.

Polymyxin B1 in the Escherichia coli inner membrane: A complex story of protein and lipopolysaccharide-mediated insertion.

The rise in multi-drug resistant Gram-negative bacterial infections has led to an increased need for "last-resort" antibiotics such as polymyxins. However, the emergence of polymyxin-resistant strains threatens to bring about a post-antibiotic era. Thus, there is a need to develop new polymyxin-based antibiotics, but a lack of knowledge of the mechanism of action of polymyxins hinders such efforts. It has recently been suggested that polymyxins induce cell lysis of the Gram-negative bacterial inner membrane (IM) by targeting trace amounts of lipopolysaccharide (LPS) localized there. We use multiscale molecular dynamics (MD), including long-timescale coarse-grained (CG) and all-atom (AA) simulations, to investigate the interactions of polymyxin B1 (PMB1) with bacterial IM models containing phospholipids (PLs), small quantities of LPS, and IM proteins. LPS was observed to (transiently) phase separate from PLs at multiple LPS concentrations, and associate with proteins in the IM. PMB1 spontaneously inserted into the IM and localized at the LPS-PL interface, where it cross-linked lipid headgroups via hydrogen bonds, sampling a wide range of interfacial environments. In the presence of membrane proteins, a small number of PMB1 molecules formed interactions with them, in a manner that was modulated by local LPS molecules. Electroporation-driven translocation of PMB1 via water-filled pores was favored at the protein-PL interface, supporting the 'destabilizing' role proteins may have within the IM. Overall, this in-depth characterization of PMB1 modes of interaction reveals how small amounts of mislocalized LPS may play a role in pre-lytic targeting and provides insights that may facilitate rational improvement of polymyxin-based antibiotics.

Structure determination needs to go viral.

Viral diseases are expected to cause new epidemics in the future, therefore, it is essential to assess how viral diversity is represented in terms of deposited protein structures. Here, data were collected from the Protein Data Bank to screen the available structures of viruses of interest to WHO. Excluding SARS-CoV-2 and HIV-1, less than 50 structures were found per year, indicating a lack of diversity. Efforts to determine viral structures are needed to increase preparedness for future public health challenges.

Evolution of an Amniote-Specific Mechanism for Modulating Ubiquitin Signaling via Phosphoregulation of the E2 Enzyme UBE2D3.

Genetic variation in the enzymes that catalyze posttranslational modification of proteins is a potentially important source of phenotypic variation during evolution. Ubiquitination is one such modification that affects turnover of virtually all of the proteins in the cell in addition to roles in signaling and epigenetic regulation. UBE2D3 is a promiscuous E2 enzyme, which acts as an ubiquitin donor for E3 ligases that catalyze ubiquitination of developmentally important proteins. We have used protein sequence comparison of UBE2D3 orthologs to identify a position in the C-terminal α-helical region of UBE2D3 that is occupied by a conserved serine in amniotes and by alanine in anamniote vertebrate and invertebrate lineages. Acquisition of the serine (S138) in the common ancestor to modern amniotes created a phosphorylation site for Aurora B. Phosphorylation of S138 disrupts the structure of UBE2D3 and reduces the level of the protein in mouse embryonic stem cells (ESCs). Substitution of S138 with the anamniote alanine (S138A) increases the level of UBE2D3 in ESCs as well as being a gain of function early embryonic lethal mutation in mice. When mutant S138A ESCs were differentiated into extraembryonic primitive endoderm, levels of the PDGFRα and FGFR1 receptor tyrosine kinases were reduced and primitive endoderm differentiation was compromised. Proximity ligation analysis showed increased interaction between UBE2D3 and the E3 ligase CBL and between CBL and the receptor tyrosine kinases. Our results identify a sequence change that altered the ubiquitination landscape at the base of the amniote lineage with potential effects on amniote biology and evolution.

Modifying the catalytic preference of alpha-amylase toward n-alkanes for bioremediation purposes using in silico strategies.

Since the beginning of oil exploration, whole ecosystems have been affected by accidents and bad practices involving petroleum compounds. In this sense, bioremediation stands out as the cheapest and most eco-friendly alternatives to reverse the damage done in oil-impacted areas. However, more efforts must be made to engineer enzymes that could be used in the bioremediation process. Interestingly, a recent work described that α-amylase, one of the most evolutionary conserved enzymes, was able to promiscuously degrade n-alkanes, a class of molecules abundant in the petroleum admixture. Considering that α-amylase is expressed in almost all known organisms, and employed in numerous biotechnological processes, using it can be a great leap toward more efficient applications of enzyme or microorganism-consortia bioremediation approaches. In this work, we employed a strict computational approach to design new α-amylase mutants with potentially enhanced catalytic efficiency toward n-alkanes. Using in silico techniques, such as molecular docking, molecular dynamics, metadynamics, and residue-residue interaction networks, we generated mutants potentially more efficient for degrading n-alkanes, L183Y, and N314A. Our results indicate that the new mutants have an increased binding rate for tetradecane, the longest n-alkane previously tested, which can reside in the catalytic center for more extended periods. Additionally, molecular dynamics and network analysis showed that the new mutations have no negative impact on protein structure than the WT. Our results aid in solidifying this enzyme as one more tool in the petroleum bioremediation toolbox.

Making it Rain: Cloud-Based Molecular Simulations for Everyone.

We present a user-friendly front-end for running molecular dynamics (MD) simulations using the OpenMM toolkit on the Google Colab framework. Our goals are (1) to highlight the usage of a cloud-computing scheme for educational purposes for a hands-on approach when learning MD simulations and (2) to exemplify how low-income research groups can perform MD simulations in the microsecond time scale. We hope this work facilitates teaching and learning of molecular simulation throughout the community.

The Lazy Life of Lipid-Linked Oligosaccharides in All Life Domains.

Lipid-linked oligosaccharides (LLOs) play an important role in the N-glycosylation pathway as the donor substrate of oligosaccharyltransferases (OSTs), which are responsible for the en bloc transfer of glycan chains onto a nascent polypeptide. The lipid component of LLO in both eukarya and archaea consists of a dolichol, and an undecaprenol in prokarya, whereas the number of isoprene units may change between species. Given the potential relevance of LLOs and their related enzymes to diverse biotechnological applications, obtaining reliable LLO models from distinct domains of life could support further studies on complex formation and their processing by OSTs, as well as protein engineering on such systems. In this work, molecular modeling techniques, such as quantum mechanics calculations, molecular dynamics simulations, and metadynamics were employed to study eukaryotic (Glc3-Man9-GlcNAc2-PP-Dolichol), bacterial (Glc1-GalNAc5-Bac1-PP-Undecaprenol), and archaeal (Glc1-Man1-Gal1-Man1-Glc1-Gal1-Glc1-P-Dolichol) LLOs in membrane bilayers. Microsecond molecular dynamics simulations and metadynamics calculations of LLOs revealed that glycan chains are more prone to interact with the membrane lipid head groups, while the PP linkages are positioned at the lipid phosphate head groups level. The dynamics of isoprenoid chains embedded within the bilayer are described, and membrane dynamics and related properties are also investigated. Overall, there are similarities regarding the structure and dynamics of the eukaryotic, the bacterial, and the archaeal LLOs in bilayers, which can support the comprehension of their association with OSTs. These data may support future studies on the transferring mechanism of the oligosaccharide chain to an acceptor protein.

The role of Zn2+, dimerization and N-glycosylation in the interaction of Auxin-Binding Protein 1 (ABP1) with different auxins.

Auxin is critical for plant growth and development. The main natural auxin is indole-3-acetic acid (IAA), whereas 1-naphthalene acetic acid (NAA) is a synthetic form. Auxin-Binding Protein 1 (ABP1) specifically binds auxins, presumably playing roles as receptor in nontranscriptional cell responses. ABP1 structure was previously established from maize at 1.9 Å resolution. To gain further insight on ABP1 structural biology, this study was carried out employing molecular dynamics simulations of the complete models of the oligomeric glycosylated proteins from maize and Arabidopsis thaliana with or without auxins. In maize, both Zn2+ coordination and glycosylation promoted conformational stability and most of such stabilization effect was located on the N-terminal region. The α-helix of C-terminal regions in ABP1 of both species unfolded during simulations, assuming a more extended structure in maize. In Arabidopsis, the helix appeared more stable, being preserved in most of the monomeric simulations and unfolding when the protein was in the dimeric form. In Arabidopsis ABP1 bound to IAA or NAA, glycosylation structures arranged around the protein, covering the putative site of entrance or egress of auxin. NAA bound protein folding was more similar to the crystal structure showing higher stability compared to that of IAA bound. The molecular structural differences of ABP1 found between the species and auxin types indicate that this auxin-binding protein shows functional specificities in dicots and monocots, as well as in auxin type binding.

In silico Investigation of the PglB Active Site Reveals Transient Catalytic States and Octahedral Metal Ion Coordination.

The last step of the bacterial N-glycosylation pathway involves PglB, an oligosaccharyltransferase, which is responsible for the en bloc transfer of a fully assembled oligosaccharide chain to a protein possessing the extended motif D/E-X-N-X-S/T. Recently, this molecule had its full structure elucidated, enabling the description of its domains and the proposition of a catalytic mechanism. By employing molecular dynamics simulations, we were able to evaluate structural aspects of PglB, suggesting prevalent motions that may bring insights into the mechanism of the glycosylated peptide detachment. Additionally, we identified transient states at the catalytic site, in which the previously described carboxamide twisting mechanism was observed. Aided by quantum mechanics calculations for each different conformational states of the catalytic site, we determined the presence of an octahedral metal coordination, along with the presence of one water molecule at the catalytic site.

Atomic model and micelle dynamics of QS-21 saponin.

QS-21 is a saponin extracted from Quillaja saponaria, widely investigated as a vaccine immunoadjuvant. However, QS-21 use is mainly limited by its chemical instability, significant variety in molecular composition and low tolerance dose in mammals. Also, this compound tends to form micelles in a concentration-dependent manner. Here, we aimed to characterize its conformation and the process of micelle formation, both experimentally and computationally. Therefore, molecular dynamics (MD) simulations were performed in systems containing different numbers of QS-21 molecules in aqueous solution, in order to evaluate the spontaneous micelle formation. The applied methodology allowed the generation of micelles whose sizes were shown to be in high agreement with small-angle X-ray scattering (SAXS). Furthermore, the ester linkage between fucose and acyl chain was less solvated in the micellar form, suggesting a reduction in hydrolysis. This is the first atomistic interpretation of previous experimental data, the first micellar characterization of saponin micelles by SAXS and first tridimensional model of a micelle constituted of saponins, contributing to the understanding of the molecular basis of these compounds.

Molecular Crowding Alters the Interactions of Polymyxin Lipopeptides within the Periplasm of E. coli: Insights from Molecular Dynamics.

The cell envelope of Gram-negative bacteria is a crowded tripartite architecture that separates the cell interior from the external environment. Two membranes encapsulate the aqueous periplasm, which contains the cell wall. Little is known about the mechanisms via which antimicrobial peptides move through the periplasm from the outer membrane to their site of action, the inner membrane. We utilize all-atom molecular dynamics to study two antimicrobial peptides, polymyxins B1 and E, within models of the E. coli periplasm crowded to different extents. In a simple chemical environment, both PMB1 and PME bind irreversibly to the cell wall. The presence of specific macromolecules leads to competition with the polymyxins for cell wall interaction sites, resulting in polymyxin dissociation from the cell wall. Chemical complexity also impacts interactions between polymyxins and Braun's lipoprotein; thus, the interaction modes of lipoprotein antibiotics within the periplasm are dependent upon the nature of the other species present.

Computational microbiology of bacteria: Advancements in molecular dynamics simulations.

Microbiology is traditionally considered within the context of wet laboratory methodologies. Computational techniques have a great potential to contribute to microbiology. Here, we describe our loose definition of "computational microbiology" and provide a short survey focused on molecular dynamics simulations of bacterial systems that fall within this definition. It is our contention that increased compositional complexity and realistic levels of molecular crowding within simulated systems are key for bridging the divide between experimental and computational microbiology.

Unrestrained conformational characterization of Stenocereus eruca saponins in aqueous and nonaqueous solvents.

Saponins are secondary metabolites that have a plethora of biological activities. However, the absence of knowledge of their 3D structures is a major drawback for structural-based strategies in medicinal chemistry. To address this problem, the current work presents structural models of Stenocereus eruca saponins, named erucasaponin A and stellatoside B. These compounds were constructed on the basis of a combination of unrestrained molecular dynamics (MD) simulations and NOESY data, in both pyridine and water. The models obtained in this way offer a robust description of the saponin dynamics in solution and support the use of submicrosecond MD simulations in describing and predicting glycoconjugate conformations.

Donor-strand exchange drives assembly of the TasA scaffold in Bacillus subtilis biofilms.

Many bacteria in nature exist in multicellular communities termed biofilms, where cells are embedded in an extracellular matrix that provides rigidity to the biofilm and protects cells from chemical and mechanical stresses. In the Gram-positive model bacterium Bacillus subtilis, TasA is the major protein component of the biofilm matrix, where it has been reported to form functional amyloid fibres contributing to biofilm structure and stability. Here, we present electron cryomicroscopy structures of TasA fibres, which show that, rather than forming amyloid fibrils, TasA monomers assemble into fibres through donor-strand exchange, with each subunit donating a ß-strand to complete the fold of the next subunit along the fibre. Combining electron cryotomography, atomic force microscopy, and mutational studies, we show how TasA fibres congregate in three dimensions to form abundant fibre bundles that are essential for B. subtilis biofilm formation. Our study explains the previously observed biochemical properties of TasA and shows how a bacterial extracellular globular protein can assemble from monomers into ß-sheet-rich fibres, and how such fibres assemble into bundles in biofilms.

Lipids mediate supramolecular outer membrane protein assembly in bacteria.

ß Barrel outer membrane proteins (OMPs) cluster into supramolecular assemblies that give function to the outer membrane (OM) of Gram-negative bacteria. How such assemblies form is unknown. Here, through photoactivatable cross-linking into the Escherichia coli OM, coupled with simulations, and biochemical and biophysical analysis, we uncover the basis for OMP clustering in vivo. OMPs are typically surrounded by an annular shell of asymmetric lipids that mediate higher-order complexes with neighboring OMPs. OMP assemblies center on the abundant porins OmpF and OmpC, against which low-abundance monomeric ß barrels, such as TonB-dependent transporters, are packed. Our study reveals OMP-lipid-OMP complexes to be the basic unit of supramolecular OMP assembly that, by extending across the entire cell surface, couples the requisite multifunctionality of the OM to its stability and impermeability.

Uncovering cryptic pockets in the SARS-CoV-2 spike glycoprotein.

The COVID-19 pandemic has prompted a rapid response in vaccine and drug development. Herein, we modeled a complete membrane-embedded SARS-CoV-2 spike glycoprotein and used molecular dynamics simulations with benzene probes designed to enhance discovery of cryptic pockets. This approach recapitulated lipid and host metabolite binding sites previously characterized by cryo-electron microscopy, revealing likely ligand entry routes, and uncovered a novel cryptic pocket with promising druggable properties located underneath the 617-628 loop. A full representation of glycan moieties was essential to accurately describe pocket dynamics. A multi-conformational behavior of the 617-628 loop in simulations was validated using hydrogen-deuterium exchange mass spectrometry experiments, supportive of opening and closing dynamics. The pocket is the site of multiple mutations associated with increased transmissibility found in SARS-CoV-2 variants of concern including Omicron. Collectively, this work highlights the utility of the benzene mapping approach in uncovering potential druggable sites on the surface of SARS-CoV-2 targets.

Polymyxin B1 within the E. coli cell envelope: insights from molecular dynamics simulations.

Polymyxins are used as last-resort antibiotics, where other treatments have been ineffectual due to antibiotic resistance. However, resistance to polymyxins has also been now reported, therefore it is instructive to characterise at the molecular level, the mechanisms of action of polymyxins. Here we review insights into these mechanisms from molecular dynamics simulations and discuss the utility of simulations as a complementary technique to  experimental methodologies.

The hitchhiker's guide to the periplasm: Unexpected molecular interactions of polymyxin B1 in E. coli.

The periplasm of Gram-negative bacteria is a complex, highly crowded molecular environment. Little is known about how antibiotics move across the periplasm and the interactions they experience. Here, atomistic molecular dynamics simulations are used to study the antibiotic polymyxin B1 within models of the periplasm, which are crowded to different extents. We show that PMB1 is likely to be able to "hitchhike" within the periplasm by binding to lipoprotein carriers-a previously unreported passive transport route. The simulations reveal that PMB1 forms both transient and long-lived interactions with proteins, osmolytes, lipids of the outer membrane, and the cell wall, and is rarely uncomplexed when in the periplasm. Furthermore, it can interfere in the conformational dynamics of native proteins. These are important considerations for interpreting its mechanism of action and are likely to also hold for other antibiotics that rely on diffusion to cross the periplasm.

Structural Basis for Silicic Acid Uptake by Higher Plants.

Many of the world's most important food crops such as rice, barley and maize accumulate silicon (Si) to high levels, resulting in better plant growth and crop yields. The first step in Si accumulation is the uptake of silicic acid by the roots, a process mediated by the structurally uncharacterised NIP subfamily of aquaporins, also named metalloid porins. Here, we present the X-ray crystal structure of the archetypal NIP family member from Oryza sativa (OsNIP2;1). The OsNIP2;1 channel is closed in the crystal structure by the cytoplasmic loop D, which is known to regulate channel opening in classical plant aquaporins. The structure further reveals a novel, five-residue extracellular selectivity filter with a large diameter. Unbiased molecular dynamics simulations show a rapid opening of the channel and visualise how silicic acid interacts with the selectivity filter prior to transmembrane diffusion. Our results will enable detailed structure-function studies of metalloid porins, including the basis of their substrate selectivity.

Development of GROMOS-Compatible Parameter Set for Simulations of Chalcones and Flavonoids.

Chalcones and flavonoids constitute a large family of plant secondary metabolites that have been explored as a potential source of novel pharmaceutical products. While the simulation of these compounds by molecular dynamics (MD) can be a valuable strategy to assess their conformational properties and so further develop their role in drug discovery, there are no set of force field parameters specifically designed and experimentally validated for their conformational description in condensed phase. So the current work developed a new parameter set for MD simulations of these compounds' main scaffolds under GROMOS force field. We employed a protocol adjusting the atomic charges and torsional parameters to the respective quantum mechanical derived dipole moments and dihedrals rotational profiles, respectively. Experimental properties of organic liquids were used as references to the calculated values to validate the parameters. Additionally, metadynamics simulations were performed to evaluate the conformational space of complex chalcones and flavonoids, while NOE contacts during simulations were measured and compared to experimental data. Accordingly, the employed protocol allowed us to obtain force field parameters that reproduce well the target data and may be expected to contribute in more accurate computational studies on the biological/therapeutical role of such molecules.