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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

CoCo-MD: A Simple and Effective Method for the Enhanced Sampling of Conformational Space.

CoCo ("complementary coordinates") is a method for ensemble enrichment based on principal component analysis (PCA) that was developed originally for the investigation of NMR data. Here we investigate the potential of the CoCo method, in combination with molecular dynamics simulations (CoCo-MD), to be used more generally for the enhanced sampling of conformational space. Using the alanine penta-peptide as a model system, we find that an iterative workflow, interleaving short multiple-walker MD simulations with long-range jumps through conformational space informed by CoCo analysis, can increase the rate of sampling of conformational space up to 10 times for the same computational effort (total number of MD timesteps). Combined with the reservoir-REMD method, free energies can be readily calculated. An alternative, approximate but fast and practically useful, alternative approach to unbiasing CoCo-MD generated data is also described. Applied to cyclosporine A, we can achieve far greater conformational sampling than has been reported previously, using a fraction of the computational resource. Simulations of the maltose binding protein, begun from the "open" state, effectively sample the "closed" conformation associated with ligand binding. The PCA-based approach means that optimal collective variables to enhance sampling need not be defined in advance by the user but are identified automatically and are adaptive, responding to the characteristics of the developing ensemble. In addition, the approach does not require any adaptations to the associated MD code and is compatible with any conventional MD package.

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.

Dynamics of DDB2-DDB1 complex under different naturally-occurring mutants in Xeroderma Pigmentosum disease.

BACKGROUND: Xeroderma Pigmentosum (XP) is a disease caused by mutations in the nucleotide excision repair (NER) pathway. Patients with XP exhibit a high propensity to skin cancers and some subtypes of XP can even present neurological impairments. During NER, DDB2 (XPE), in complex with DDB1 (DDB-Complex), performs the DNA lesion recognition. However, not much is known about how mutations found in XP patients affect the DDB2 structure and complex assembly. Thus, we searched for structural evidence associated with the role of three naturally occurring mutations found in XPE patients: R273H, K244E, and L350P. METHODS: Each mutant was individually constructed and submitted to multiple molecular dynamics simulations, done in triplicate for each designed system. Additionally, Dynamic Residue Interaction Networks were designed for each system and analyzed parallel with the simulations. RESULTS: DDB2 mutations promoted loss of flexibility in the overall protein structure, producing a different conformational behavior in comparison to the WT, especially in the region comprising residues 354 to 371. Furthermore, the DDB-complex containing the mutated forms of DDB2 showed distinct behaviors for each mutant: R273H displayed higher structural instability when complexed; L350P affected DDB1 protein-protein binding with DDB2; and K244E, altered the complex binding trough different ways than L350P. CONCLUSIONS: The data gathered throughout the analyses helps to enlighten the structural basis for how naturally occurring mutations found in XPE patients impact on DDB2 and DDB1 function. GENERAL SIGNIFICANCE: Our data influence not only on the knowledge of XP but on the DNA repair mechanisms of NER itself.

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.