Proton-pumping rhodopsins (PPRs) are photoactive retinal-binding proteins that transport ions across biological membranes in response to light. These proteins are interesting for lightharvesting applications in bioenergy production, in optogenetics applications in neuroscience, and as fluorescent sensors of membrane potential. Little is known, however, about how the protein sequence determines the considerable variation in spectral properties of PPRs from different biological niches or how to engineer these properties in a given PPR. Here we report a comprehensive study of amino acid substitutions in the retinal binding pocket of Gloeobacter violacaeus rhodopsin (GR) that tune its spectral properties. Directed evolution generated 70 GR variants with absorption maxima shifted by up to +/-80 nm, extending the protein's light absorption significantly beyond the range of known natural PPRs. While proton pumping activity was disrupted in many of the spectrally shifted variants, we identified single tuning mutations that incurrred blue and red shifts of 42 nm and 22 nm, respectively, that did not disrupt proton pumping. Blue-shifting mutations were distributed evenly along the retinal molecule while redshifting mutations were clustered near the residue K257, which forms a covalent bond with retinal through a Schiff base linkage. Thirty-four of the identified tuning mutations are not found in known microbial rhodopsins. We discovered a subset of red-shifted GRs that exhibit high levels of fluorescence relative to the wild-type protein.

Folding and cellular localization of many proteins of Gram-negative bacteria rely on a network of chaperones and secretion systems. Among them is the lipase-specific foldase Lif, a membrane-bound steric chaperone that tightly binds (K D = 29 nM) and mediates folding of the lipase LipA, a virulence factor of the pathogenic bacterium P. aeruginosa. Lif consists of five-domains, including a mini domain MD1 essential for LipA folding. However, the molecular mechanism of Lif-assisted LipA folding remains elusive. Here, we show in in vitro experiments using a soluble form of Lif (sLif) that isolated MD1 inhibits sLif-assisted LipA activation. Furthermore, the ability to activate LipA is lost in the variant sLif Y99A , in which the evolutionary conserved amino acid Y99 from helix α1 of MD1 is mutated to alanine. This coincides with an approximately threefold reduced affinity of the variant to LipA together with increased flexibility of sLif Y99A in the complex as determined by polarization-resolved fluorescence spectroscopy. We have solved the NMR solution structures of P. aeruginosa MD1 and variant MD1 Y99A revealing a similar fold indicating that a structural modification is likely not the reason for the impaired activity of variant sLif Y99A. Molecular dynamics simulations of the sLif:LipA complex in connection with rigidity analyses suggest a long-range network of interactions spanning from Y99 of sLif to the active site of LipA, which might be essential for LipA activation. These findings provide important details about the putative mechanism for LipA activation and point to a general mechanism of protein folding by multi-domain steric chaperones. The Gram-negative human pathogen Pseudomonas aeruginosa produces a wide range of extracellular enzymes 1,2 , among them the lipase LipA, a secreted putative virulence factor 3-5. For its conversion into an enzymatically active conformation, LipA requires the assistance of an inner membrane-bound chaperone named lipase-specific foldase (Lif) 6. On the folding pathway, LipA can adopt several structurally different intermediates: an inactive and unfolded molten globule-like conformation 7 , a near-natively folded pre-active conformation 8 and two folded conformations that differ in the structure of the α-helical lid covering the active site, with the folded closed

Lipases are essential and widely used biocatalysts. Hence, the production of lipases requires a detailed understanding of the molecular mechanism of its folding and secretion. Lipase A from Pseudomonas aeruginosa, PaLipA, constitutes a prominent example that has additional relevance because of its role as a virulence factor in many diseases. PaLipA requires the assistance of a membrane-integrated steric chaperone, the lipasespecific foldase Lif, to achieve its enzymatically active state. However, the molecular mechanism of how Lif activates its cognate lipase has remained elusive. Here, we show by molecular dynamics simulations at the atomistic level and potential of mean force computations that Lif catalyzes the activation process of PaLipA by structurally stabilizing an intermediate PaLipA conformation, particularly a β-sheet in the region of residues 17-30, such that the opening of PaLipA's lid domain is facilitated. This opening allows substrate access to PaLipA's catalytic site. A surprising and so far not fully understood aspect of our study is that the open state of PaLipA is unstable compared to the closed one according to our computational and in vitro biochemical results. We thus speculate that further interactions of PaLipA with the Xcp secretion machinery and/or components of the extracellular matrix contribute to the remaining activity of secreted PaLipA.

Maximum entropy methods (MEMs) determine posterior distributions by combining experimental data with prior information. MEMs are frequently used to reconstruct conformational ensembles of molecular systems for experimental information and initial molecular ensembles. We performed time-resolved Forster resonance energy transfer (FRET) experiments to probe the interdye distance distributions of the lipase-specific foldase Lif in the apo state, which likely has highly flexible, disordered, and/or ordered structural elements. Distance distributions estimated from ensembles of molecular dynamics (MD) simulations serve as prior information, and FRET experiments, analyzed within a Bayesian framework to recover distance distributions, are used for optimization. We tested priors obtained by MD with different force fields (FFs) tailored to ordered (FF99SB, FF14SB, and FF19SB) and disordered proteins (IDPSFF and FF99SBdisp). We obtained five substantially different posterior ensembles. As in our FRET experiments the noise is characterized by photon counting statistics, for a validated dye model, MEM can quantify consistencies between experiment and prior or posterior ensembles. However, posterior populations of conformations are uncorrelated to structural similarities for individual structures selected from different prior ensembles. Therefore, we assessed MEM simulating varying priors in synthetic experiments with known target ensembles. We found that (i) the prior and experimental information must be carefully balanced for optimal posterior ensembles to minimize perturbations of populations by overfitting and (ii) only ensemble-integrated quantities like inter-residue distance distributions or density maps can be reliably obtained but not ensembles of atomistic structures. This is because MEM optimizes ensembles but not individual structures. This result for a highly flexible system suggests that structurally varying priors calculated from varying prior ensembles, e.g., generated with different FFs, may serve as an ad hoc estimate for MEM reconstruction robustness.