Search

In Gram-negative bacteria, the biogenesis of β-barrel outer membrane proteins is mediated by the β-barrel assembly machinery (BAM). The mechanism employed by BAM is complex and so far- incompletely understood. Here, we report the structures of BAM in nanodiscs, prepared using polar lipids and native membranes, where we observe an outward-open state. Mutations in the barrel domain of BamA reveal that plasticity in BAM is essential, particularly along the lateral seam of the barrel domain, which is further supported by molecular dynamics simulations that show conformational dynamics in BAM are modulated by the accessory proteins. We also report the structure of BAM in complex with EspP, which reveals an early folding intermediate where EspP threads from the underside of BAM and incorporates into the barrel domain of BamA, supporting a hybrid-barrel budding mechanism in which the substrate is folded into the membrane sequentially rather than as a single unit., The β-barrel assembly machinery (BAM) assists the folding and membrane insertion of bacterial outer membrane proteins. Here, the authors report structural characterization of BAM in lipid environment and in complex with the client protein EspP integrated into the barrel of BamA, providing insight into BAM mechanism of function.

Intrinsically disordered regions (IDRs) in proteins are often targets of combinatorial posttranslational modifications, which serve to regulate protein structure and function. Emerging evidence suggests that the N-terminal tails of G protein γ subunits, which are essential components of heterotrimeric G proteins, are intrinsically disordered, phosphorylation-dependent determinants of G protein signaling. Here, we found that the yeast Gγ subunit Ste18 underwent combinatorial, multisite phosphorylation events within its N-terminal IDR. G protein–coupled receptor (GPCR) activation and osmotic stress induced phosphorylation at Ser(7), whereas glucose and acid stress induced phosphorylation at Ser(3), which was a quantitative indicator of intracellular pH. Each site was phosphorylated by a distinct set of kinases, and phosphorylation of one site affected phosphorylation of the other, as determined through exposure to serial stimuli and through phosphosite mutagenesis. Lastly, we showed that phosphorylation resulted in changes in IDR structure and that different combinations of phosphorylation events modulated the activation rate and amplitude of the downstream mitogen-activated protein kinase Fus3. These data place Gγ subunits among intrinsically disordered proteins that undergo combinatorial posttranslational modifications that govern signaling pathway output.

The main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of Mpro, a cysteine protease, have been determined, facilitating structure-based drug design. Mpro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, Mpro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nu-cleophile Cys145 have been debated in previous studies of SARS-CoV Mpro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 Mpro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of Mpro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored Nδ (HD) and Nϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 Mpro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

Background In Gram-negative bacteria, type Va and Vc autotransporters are proteins that contain both a secreted virulence factor (the “passenger” domain) and a β-barrel that aids its export. While it is known that the folding and insertion of the β-barrel domain utilize the β-barrel assembly machinery (BAM) complex, how the passenger domain is secreted and folded across the membrane remains to be determined. The hairpin model states that passenger domain secretion occurs independently through the fully-formed and membrane-inserted β-barrel domain via a hairpin folding intermediate. In contrast, the BamA-assisted model states that the passenger domain is secreted through a hybrid of BamA, the essential subunit of the BAM complex, and the β-barrel domain of the autotransporter. Methods To ascertain the models' plausibility, we have used molecular dynamics to simulate passenger domain secretion for two autotransporters, EspP and YadA. Results We observed that each protein's β-barrel is unable to accommodate the secreting passenger domain in a hairpin configuration without major structural distortions. Additionally, the force required for secretion through EspP's β-barrel is more than that through the BamA β-barrel. Conclusions Secretion of autotransporters most likely occurs through an incompletely formed β-barrel domain of the autotransporter in conjunction with BamA. General significance Secretion of virulence factors is a process used by practically all pathogenic Gram-negative bacteria. Understanding this process is a necessary step towards limiting their infectious capacity.

Gaussian accelerated molecular dynamics (GaMD) is a recently developed enhanced sampling technique that provides efficient free energy calculations of biomolecules. Like the previous accelerated molecular dynamics (aMD), GaMD allows for "unconstrained" enhanced sampling without the need to set predefined collective variables and so is useful for studying complex biomolecular conformational changes such as protein folding and ligand binding. Furthermore, because the boost potential is constructed using a harmonic function that follows Gaussian distribution in GaMD, cumulant expansion to the second order can be applied to recover the original free energy profiles of proteins and other large biomolecules, which solves a long-standing energetic reweighting problem of the previous aMD method. Taken together, GaMD offers major advantages for both unconstrained enhanced sampling and free energy calculations of large biomolecules. Here, we have implemented GaMD in the NAMD package on top of the existing aMD feature and validated it on three model systems: alanine dipeptide, the chignolin fast-folding protein, and the M3 muscarinic G protein-coupled receptor (GPCR). For alanine dipeptide, while conventional molecular dynamics (cMD) simulations performed for 30 ns are poorly converged, GaMD simulations of the same length yield free energy profiles that agree quantitatively with those of 1000 ns cMD simulation. Further GaMD simulations have captured folding of the chignolin and binding of the acetylcholine (ACh) endogenous agonist to the M3 muscarinic receptor. The reweighted free energy profiles are used to characterize the protein folding and ligand binding pathways quantitatively. GaMD implemented in the scalable NAMD is widely applicable to enhanced sampling and free energy calculations of large biomolecules.