Simulating Gram-Negative Bacterial Outer Membrane: A Coarse Grain Model

PyCGTOOL: Automated Generation of Coarse-Grained Molecular Dynamics Models from Atomistic Trajectories

Development of coarse-grained (CG) molecular dynamics models is often a laborious process which commonly relies upon approximations to similar models, rather than systematic parametrization. PyCGTOOL automates much of the construction of CG models via calculation of both equilibrium values and force constants of internal coordinates directly from atomistic molecular dynamics simulation trajectories. The derivation of bespoke parameters from atomistic simulations improves the quality of the CG model compared to the use of generic parameters derived from other molecules, while automation greatly reduces the time required. The ease of configuration of PyCGTOOL enables the rapid investigation of multiple atom-to-bead mappings and topologies. Although we present PyCGTOOL used in combination with the GROMACS molecular dynamics engine its use of standard trajectory input libraries means that it is in principle compatible with other software. The software is available from the URL https://github.com/jag1g13/pycgtool as the following doi: 10.5281/zenodo.259330 .

Two Outer Membrane Proteins Contribute to Caulobacter crescentus Cellular Fitness by Preventing Intracellular S-Layer Protein Accumulation

Surface layers, or S-layers, are two-dimensional protein arrays that form the outermost layer of many bacteria and archaea. They serve several functions, including physical protection of the cell from environmental threats. The high abundance of S-layer proteins necessitates a highly efficient export mechanism to transport the S-layer protein from the cytoplasm to the cell exterior. Caulobacter crescentus is unique in that it has two homologous, seemingly redundant outer membrane proteins, RsaF_a and RsaF_b, which together with other components form a type I protein translocation pathway for S-layer export. These proteins have homology to Escherichia coli TolC, the outer membrane channel of multidrug efflux pumps. Here we provide evidence that, unlike TolC, RsaF_a and RsaF_b are not involved in either the maintenance of membrane stability or the active export of antimicrobial compounds. Rather, RsaF_a and RsaF_b are required to prevent intracellular accumulation and aggregation of the S-layer protein RsaA; deletion of RsaF_a and RsaF_b led to a general growth defect and lowered cellular fitness. Using Western blotting, transmission electron microscopy, and transcriptome sequencing (RNA-seq), we show that loss of both RsaF_a and RsaF_b led to accumulation of insoluble RsaA in the cytoplasm, which in turn caused upregulation of a number of genes involved in protein misfolding and degradation pathways. These findings provide new insight into the requirement for RsaF_a and RsaF_b in cellular fitness and tolerance to antimicrobial agents and further our understanding of the S-layer export mechanism on both the transcriptional and translational levels in C. crescentusIMPORTANCE Decreased growth rate and reduced cell fitness are common side effects of protein production in overexpression systems. Inclusion bodies typically form inside the cell, largely due to a lack of sufficient export machinery to transport the overexpressed proteins to the extracellular environment. This phenomenon can conceivably also occur in natural systems. As one example of a system evolved to prevent intracellular protein accumulation, our study demonstrates that Caulobacter crescentus has two homologous outer membrane transporter proteins that are involved in S-layer export. This is an interesting case study that demonstrates how bacteria can evolve redundancy to ensure adequate protein export functionality and maintain high cellular fitness. Moreover, we provide evidence that these two outer membrane proteins, although being the closest C. crescentus homologs to TolC in E. coli, do not process TolC functionality in C. crescentus.

Methodology for identification of pore forming antimicrobial peptides from soy protein subunits β-conglycinin and glycinin

Antimicrobial peptides (AMPs) inactivate microbial cells through pore formation in cell membrane. Because of their different mode of action compared to antibiotics, AMPs can be effectively used to combat drug resistant bacteria in human health. AMPs can also be used to replace antibiotics in animal feed and immobilized on food packaging films. In this research, we developed a methodology based on mechanistic evaluation of peptide-lipid bilayer interaction to identify AMPs from soy protein. Production of AMPs from soy protein is an attractive, cost-saving alternative for commercial consideration, because soy protein is an abundant and common protein resource. This methodology is also applicable for identification of AMPs from any protein. Initial screening of peptide segments from soy glycinin (11S) and soy β-conglycinin (7S) subunits was based on their hydrophobicity, hydrophobic moment and net charge. Delicate balance between hydrophilic and hydrophobic interactions is necessary for pore formation. High hydrophobicity decreases the peptide solubility in aqueous phase whereas high hydrophilicity limits binding of the peptide to the bilayer. Out of several candidates chosen from the initial screening, two peptides satisfied the criteria for antimicrobial activity, viz. (i) lipid-peptide binding in surface state and (ii) pore formation in transmembrane state of the aggregate. This method of identification of antimicrobial activity via molecular dynamics simulation was shown to be robust in that it is insensitive to the number of peptides employed in the simulation, initial peptide structure and force field. Their antimicrobial activity against Listeria monocytogenes and Escherichia coli was further confirmed by spot-on-lawn test.

Molecular dynamics simulations of membrane proteins and their interactions: from nanoscale to mesoscale

Molecular dynamics simulations provide a computational tool to probe membrane proteins and systems at length scales ranging from nanometers to close to a micrometer, and on microsecond timescales. All atom and coarse-grained simulations may be used to explore in detail the interactions of membrane proteins and specific lipids, yielding predictions of lipid binding sites in good agreement with available structural data. Building on the success of protein-lipid interaction simulations, larger scale simulations reveal crowding and clustering of proteins, resulting in slow and anomalous diffusional dynamics, within realistic models of cell membranes. Current methods allow near atomic resolution simulations of small membrane organelles, and of enveloped viruses to be performed, revealing key aspects of their structure and functionally important dynamics.

Lateral opening in the intact β-barrel assembly machinery captured by cryo-EM

The β-barrel assembly machinery (BAM) is a ∼203 kDa complex of five proteins (BamA-E), which is essential for viability in E. coli. BAM promotes the folding and insertion of β-barrel proteins into the outer membrane via a poorly understood mechanism. Several current models suggest that BAM functions through a 'lateral gating' motion of the β-barrel of BamA. Here we present a cryo-EM structure of the BamABCDE complex, at 4.9 Å resolution. The structure is in a laterally open conformation showing that gating is independent of BamB binding. We describe conformational changes throughout the complex and interactions between BamA, B, D and E, and the detergent micelle that suggest communication between BAM and the lipid bilayer. Finally, using an enhanced reconstitution protocol and functional assays, we show that for the outer membrane protein OmpT, efficient folding in vitro requires lateral gating in BAM.