Structure of the amino-terminal domain of phage 434 repressor at 2.0 A resolution

Protein stability prediction by fine-tuning a protein language model on a mega-scale dataset

Protein stability plays a crucial role in a variety of applications, such as food processing, therapeutics, and the identification of pathogenic mutations. Engineering campaigns commonly seek to improve protein stability, and there is a strong interest in streamlining these processes to enable rapid optimization of highly stabilized proteins with fewer iterations. In this work, we explore utilizing a mega-scale dataset to develop a protein language model optimized for stability prediction. ESMtherm is trained on the folding stability of 528k natural and de novo sequences derived from 461 protein domains and can accommodate deletions, insertions, and multiple-point mutations. We show that a protein language model can be fine-tuned to predict folding stability. ESMtherm performs reasonably on small protein domains and generalizes to sequences distal from the training set. Lastly, we discuss our model's limitations compared to other state-of-the-art methods in generalizing to larger protein scaffolds. Our results highlight the need for large-scale stability measurements on a diverse dataset that mirrors the distribution of sequence lengths commonly observed in nature.

Hybrid Vigor: Importance of Hybrid λ Phages in Early Insights in Molecular Biology

Laboratory-generated hybrids between phage λ and related phages played a seminal role in establishment of the λ model system, which, in turn, served to develop many of the foundational concepts of molecular biology, including gene structure and control. Important λ hybrids with phages 21 and 434 were the earliest of such phages. To understand the biology of these hybrids in full detail, we determined the complete genome sequences of phages 21 and 434. Although both genomes are canonical members of the λ-like phage family, they both carry unsuspected bacterial virulence gene types not previously described in this group of phages. In addition, we determined the sequences of the hybrid phages λ imm21, λ imm434, and λ h434 imm21. These sequences show that the replacements of λ DNA by nonhomologous segments of 21 or 434 DNA occurred through homologous recombination in adjacent sequences that are nearly identical in the parental phages. These five genome sequences correct a number of errors in published sequence fragments of the 21 and 434 genomes, and they point out nine nucleotide differences from Sanger's original λ sequence that are likely present in most extant λ strains in laboratory use today. We discuss the historical importance of these hybrid phages in the development of fundamental tenets of molecular biology and in some of the earliest gene cloning vectors. The 434 and 21 genomes reinforce the conclusion that the genomes of essentially all natural λ-like phages are mosaics of sequence modules from a pool of exchangeable segments.

Protein structure determination using a Riemannian approach

Protein NMR structure determination is one of the most extensively studied problems. Here, we adopt a novel method based on a matrix completion technique - the Riemannian approach - to rebuild the protein structure from the nuclear Overhauser effect distance restraints and the dihedral angle restraints. In comparison with the cyana method, the results generated via the Riemannian approach are more similar to the standard X-ray crystallographic structures as a result of the simple but powerful internal calculation processing function. In addition, our results demonstrate that the Riemannian approach has a comparable or even better performance than the cyana method on other structural assessment metrics, including the stereochemical quality and restraint violations. The Riemannian approach software is available at: https://github.com/xubiaopeng/Protein_Recon_MCRiemman.

Crystal structure of an anti-CRISPR protein, AcrIIA1

Clustered regularly interspaced short palindromic repeats (CRISPRs) and CRISPR-associated (Cas) proteins provide bacteria with RNA-based adaptive immunity against phage infection. To counteract this defense mechanism, phages evolved anti-CRISPR (Acr) proteins that inactivate the CRISPR-Cas systems. AcrIIA1, encoded by Listeria monocytogenes prophages, is the most prevalent among the Acr proteins targeting type II-A CRISPR-Cas systems and has been used as a marker to identify other Acr proteins. Here, we report the crystal structure of AcrIIA1 and its RNA-binding affinity. AcrIIA1 forms a dimer with a novel two helical-domain architecture. The N-terminal domain of AcrIIA1 exhibits a helix-turn-helix motif similar to transcriptional factors. When overexpressed in Escherichia coli, AcrIIA1 associates with RNAs, suggesting that AcrIIA1 functions via nucleic acid recognition. Taken together, the unique structural and functional features of AcrIIA1 suggest its distinct mode of Acr activity, expanding the diversity of the inhibitory mechanisms employed by Acr proteins.

Functional characterization of BcrR: a one-component transmembrane signal transduction system for bacitracin resistance

Bacitracin is a cell wall targeting antimicrobial with clinical and agricultural applications. With the growing mismatch between antimicrobial resistance and development, it is essential we understand the molecular mechanisms of resistance in order to prioritize and generate new effective antimicrobials. BcrR is a unique membrane-bound one-component system that regulates high-level bacitracin resistance in Enterococcus faecalis. In the presence of bacitracin, BcrR activates transcription of the bcrABD operon conferring resistance through a putative ATP-binding cassette (ABC) transporter (BcrAB). BcrR has three putative functional domains, an N-terminal helix-turn-helix DNA-binding domain, an intermediate oligomerization domain and a C-terminal transmembrane domain. However, the molecular mechanisms of signal transduction remain unknown. Random mutagenesis of bcrR was performed to generate loss- and gain-of-function mutants using transcriptional reporters fused to the target promoter PbcrA. Fifteen unique mutants were isolated across all three proposed functional domains, comprising 14 loss-of-function and one gain-of-function mutant. The gain-of-function variant (G64D) mapped to the putative dimerization domain of BcrR, and functional analyses indicated that the G64D mutant constitutively expresses the PbcrA-luxABCDE reporter. DNA-binding and membrane insertion were not affected in the five mutants chosen for further characterization. Homology modelling revealed putative roles for two key residues (R11 and S33) in BcrR activation. Here we present a new model of BcrR activation and signal transduction, providing valuable insight into the functional characterization of membrane-bound one-component systems and how they can coordinate critical bacterial responses, such as antimicrobial resistance.

Cooperative "folding transition" in the sequence space facilitates function-driven evolution of protein families

In the protein sequence space, natural proteins form clusters of families which are characterized by their unique native folds whereas the great majority of random polypeptides are neither clustered nor foldable to unique structures. Since a given polypeptide can be either foldable or unfoldable, a kind of "folding transition" is expected at the boundary of a protein family in the sequence space. By Monte Carlo simulations of a statistical mechanical model of protein sequence alignment that coherently incorporates both short-range and long-range interactions as well as variable-length insertions to reproduce the statistics of the multiple sequence alignment of a given protein family, we demonstrate the existence of such transition between natural-like sequences and random sequences in the sequence subspaces for 15 domain families of various folds. The transition was found to be highly cooperative and two-state-like. Furthermore, enforcing or suppressing consensus residues on a few of the well-conserved sites enhanced or diminished, respectively, the natural-like pattern formation over the entire sequence. In most families, the key sites included ligand binding sites. These results suggest some selective pressure on the key residues, such as ligand binding activity, may cooperatively facilitate the emergence of a protein family during evolution. From a more practical aspect, the present results highlight an essential role of long-range effects in precisely defining protein families, which are absent in conventional sequence models.

Coevolutionary information, protein folding landscapes, and the thermodynamics of natural selection

The energy landscape used by nature over evolutionary timescales to select protein sequences is essentially the same as the one that folds these sequences into functioning proteins, sometimes in microseconds. We show that genomic data, physical coarse-grained free energy functions, and family-specific information theoretic models can be combined to give consistent estimates of energy landscape characteristics of natural proteins. One such characteristic is the effective temperature T(sel) at which these foldable sequences have been selected in sequence space by evolution. T(sel) quantifies the importance of folded-state energetics and structural specificity for molecular evolution. Across all protein families studied, our estimates for T(sel) are well below the experimental folding temperatures, indicating that the energy landscapes of natural foldable proteins are strongly funneled toward the native state.

The N-terminal domain of the repressor of Staphylococcus aureus phage Φ11 possesses an unusual dimerization ability and DNA binding affinity

Bacteriophage Φ11 uses Staphylococcus aureus as its host and, like lambdoid phages, harbors three homologous operators in between its two divergently oriented repressor genes. None of the repressors of Φ11, however, showed binding to all three operators, even at high concentrations. To understand why the DNA binding mechanism of Φ11 repressors does not match that of lambdoid phage repressors, we studied the N-terminal domain of the Φ11 lysogenic repressor, as it harbors a putative helix-turn-helix motif. Our data revealed that the secondary and tertiary structures of the N-terminal domain were different from those of the full-length repressor. Nonetheless, the N-terminal domain was able to dimerize and bind to the operators similar to the intact repressor. In addition, the operator base specificity, binding stoichiometry, and binding mechanism of this domain were nearly identical to those of the whole repressor. The binding affinities of the repressor and its N-terminal domain were reduced to a similar extent when the temperature was increased to 42°C. Both proteins also adequately dislodged a RNA polymerase from a Φ11 DNA fragment carrying two operators and a promoter. Unlike the intact repressor, the binding of the N-terminal domain to two adjacent operator sites was not cooperative in nature. Taken together, we suggest that the dimerization and DNA binding abilities of the N-terminal domain of the Φ11 repressor are distinct from those of the DNA binding domains of other phage repressors.

Improvement of structure-based potentials for protein folding by native and nonnative hydrogen bonds

Pure Gō models (where every native interaction equally stabilizes the folded state) have widely proved their convenience in the computational investigation of protein folding. However, a chemistry-based description of the real interactions also provides a desirable tune in the analysis of the folding process, and thus some hybrid Gō potentials that combine both aspects have been proposed. Among all the noncovalent interactions that contribute to protein folding, hydrogen bonds are the only ones with a partial covalent character. This feature makes them directional and, thus, more difficult to model as part of the coarse-grained descriptions that are typically employed in Gō models. Thanks to a simplified but rigorous representation of backbone hydrogen bonds that we have recently proposed, we present in this article a combined potential (Gō + backbone hydrogen bond) to study the thermodynamics of protein folding in the frame of very simple simulation models. We show that the explicit inclusion of hydrogen bonds leads to a systematic improvement in the description of protein folding. We discuss a representative set of examples (from two-state folders to downhill proteins, with different types of native structures) that reveal a relevant agreement with experimental data.

Generic coarse-grained model for protein folding and aggregation

A generic coarse-grained (CG) protein model is presented. The intermediate level of resolution (four beads per amino acid, implicit solvent) allows for accurate sampling of local conformations. It relies on simple interactions that emphasize structure, such as hydrogen bonds and hydrophobicity. Realistic alpha/beta content is achieved by including an effective nearest-neighbor dipolar interaction. Parameters are tuned to reproduce both local conformations and tertiary structures. The thermodynamics and kinetics of a three-helix bundle are studied. We check that the CG model is able to fold proteins with tertiary structures and amino acid sequences different from the one used for parameter tuning. By studying both helical and extended conformations we make sure the force field is not biased toward any particular secondary structure. The accuracy involved in folding not only the test protein but also other ones show strong evidence for amino acid cooperativity embedded in the model. Without any further adjustments or bias a realistic oligopeptide aggregation scenario is observed.

Direct stimulus perception and transcription activation by a membrane-bound DNA binding protein

Few membrane proteins with a role in transcriptional regulation have been studied, and none are able to perceive their respective stimuli and activate transcription of their regulons without the aid of auxiliary proteins. The bacitracin resistance regulator, BcrR, of Enterococcus faecalis is a membrane-bound DNA binding protein and is required for bacitracin-dependent expression of the bacitracin resistance genes, bcrABD. Here, we show that BcrR interacts directly with Zn2+ bacitracin (Kd = 2-5 micropM), but not metal-free bacitracin. A solution-based DNA binding assay demonstrated that the affinity of BcrR for its target DNA is much higher (Kd = 40 nM) than previously found for transmembrane regulators and is comparable to that of soluble DNA binding proteins. A construct of BcrR that lacked the transmembrane domain was unable to bind to DNA, indicating that membrane localization was important for DNA binding. Bacitracin did not cause a change in the DNaseI footprint of BcrR on the bcrA promoter, but in vitro transcription assays with BcrR proteoliposomes showed bacitracin-dependent activation of transcription. These findings demonstrate that BcrR is a bona fide one-component transmembrane signal transduction system, which perceives an extracellular stimulus (presence of bacitracin) and relays it to an intracellular transcriptional response independent of any auxiliary proteins.

Macromolecular modeling with rosetta

Advances over the past few years have begun to enable prediction and design of macromolecular structures at near-atomic accuracy. Progress has stemmed from the development of reasonably accurate and efficiently computed all-atom potential functions as well as effective conformational sampling strategies appropriate for searching a highly rugged energy landscape, both driven by feedback from structure prediction and design tests. A unified energetic and kinematic framework in the Rosetta program allows a wide range of molecular modeling problems, from fibril structure prediction to RNA folding to the design of new protein interfaces, to be readily investigated and highlights areas for improvement. The methodology enables the creation of novel molecules with useful functions and holds promise for accelerating experimental structural inference. Emerging connections to crystallographic phasing, NMR modeling, and lower-resolution approaches are described and critically assessed.

N15 Cro and lambda Cro: orthologous DNA-binding domains with completely different but equally effective homodimer interfaces

Bacteriophage Cro proteins bind to target DNA as dimers but do not all dimerize with equal strength, and differ in fold in the region of the dimer interface. We report the structure of the Cro protein from Enterobacteria phage N15 at 1.05 A resolution. The subunit fold contains five alpha-helices and is closely similar to the structure of P22 Cro (1.3 A backbone room mean square difference over 52 residues), but quite different from that of lambda Cro, a structurally diverged member of this family with a mixed alpha-helix/beta-sheet fold. N15 Cro crystallizes as a biological dimer with an extensive interface (1303 A(2) change in accessible surface area per dimer) and also dimerizes in solution with a K(d) of 5.1 +/- 1.5 microM. Its dimerization is much stronger than that of its structural homolog P22 Cro, which does not self-associate detectably in solution. Instead, the level of self-association and interfacial area for N15 Cro is similar to that of lambda Cro, even though these two orthologs do not share the same fold and have dimer interfaces that are qualitatively different in structure. The common Cro ancestor is thought to be an all-helical monomer similar to P22 Cro. We propose that two Cro descendants independently developed stronger dimerization by entirely different mechanisms.

Crystal structure of Mycobacterium tuberculosis LrpA, a leucine-responsive global regulator associated with starvation response

The bacterial leucine-responsive regulatory protein (Lrp) is a global transcriptional regulator that controls the expression of many genes during starvation and the transition to stationary phase. The Mycobacterium tuberculosis gene Rv3291c encodes a 150-amino acid protein (designated here as Mtb LrpA) with homology with Escherichia coli Lrp. The crystal structure of the native form of Mtb LrpA was solved at 2.1 A. The Mtb LrpA structure shows an N-terminal DNA-binding domain with a helix-turn-helix (HTH) motif, and a C-terminal regulatory domain. In comparison to the complex of E. coli AsnC with asparagine, the effector-binding pocket (including loop 100-106) in LrpA appears to be largely preserved, with hydrophobic substitutions consistent with its specificity for leucine. The effector-binding pocket is formed at the interface between adjacent dimers, with an opening to the core of the octamer as in AsnC, and an additional substrate-access channel opening to the outer surface of the octamer. Using electrophoretic mobility shift assays, purified Mtb LrpA protein was shown to form a protein-DNA complex with the lat promoter, demonstrating that the lat operon is a direct target of LrpA. Using computational analysis, a putative motif is identified in this region that is also present upstream of other operons differentially regulated under starvation. This study provides insights into the potential role of LrpA as a global regulator in the transition of M. tuberculosis to a persistent state.

How sensitive are nanosecond molecular dynamics simulations of proteins to changes in the force field?

The sensitivity of molecular dynamics simulations to variations in the force field has been examined in relation to a set of 36 structures corresponding to 31 proteins simulated by using different versions of the GROMOS force field. The three parameter sets used (43a1, 53a5, and 53a6) differ significantly in regard to the nonbonded parameters for polar functional groups and their ability to reproduce the correct solvation and partitioning behavior of small molecular analogues of the amino acid side chains. Despite the differences in the force field parameters no major differences could be detected in a wide range of structural properties such as the root-mean-square deviation from the experimental structure, radii of gyration, solvent accessible surface, secondary structure, or hydrogen bond propensities on a 5 to 10 ns time scale. The small differences that were observed correlated primarily with the presence of charged residues as opposed to residues that differed most between the parameter sets. The work highlights the variation that can be observed in nanosecond simulations of protein systems and implications of this for force field validation, as well as for the analysis of protein simulations in general.

Structure and reactivity of hydroxypropylphosphonic acid epoxidase in fosfomycin biosynthesis by a cation- and flavin-dependent mechanism

The biosynthesis of fosfomycin, an oxirane antibiotic in clinical use, involves a unique epoxidation catalyzed by (S)-2-hydroxypropylphosphonic acid epoxidase (HPPE). The reaction is essentially dehydrogenation of a secondary alcohol. A high-resolution crystallographic analysis reveals that the HPPE subunit displays a two-domain combination. The C-terminal or catalytic domain has the cupin fold that binds a divalent cation, whereas the N-terminal domain carries a helix-turn-helix motif with putative DNA-binding helices positioned 34 A apart. The structure of HPPE serves as a model for numerous proteins, of ill-defined function, predicted to be transcription factors but carrying a cupin domain at the C terminus. Structure-reactivity analyses reveal conformational changes near the catalytic center driven by the presence or absence of ligand, that HPPE is a Zn(2+)/Fe(2+)-dependent epoxidase, proof that flavin mononucleotide is required for catalysis, and allow us to propose a simple mechanism that is compatible with previous experimental data. The participation of the redox inert Zn(2+) in the mechanism is surprising and indicates that Lewis acid properties of the metal ions are sufficient to polarize the substrate and, aided by flavin mononucleotide reduction, facilitate the epoxidation.

Refinement of protein cores and protein–peptide interfaces using a potential scaling approach

Refinement of side chain conformations in protein model structures and at the interface of predicted protein-protein or protein-peptide complexes is an important step during protein structural modelling and docking. A common approach for side chain prediction is to assume a rigid protein main chain for both docking partners and search for an optimal set of side chain rotamers to optimize the steric fit. However, depending on the target-template similarity in the case of comparative protein modelling and on the accuracy of an initially docked complex, the main chain template structure is only an approximation of a realistic target main chain. An inaccurate rigid main chain conformation can in turn interfere with the prediction of side chain conformations. In the present study, a potential scaling approach (PS-MD) during a molecular dynamics (MD) simulation that also allows the inclusion of explicit solvent has been used to predict side chain conformations on semi-flexible protein main chains. The PS-MD method converges much faster to realistic protein-peptide interface structures or protein core structures than standard MD simulations. Depending on the accuracy of the protein main chain, it also gives significantly better results compared with the standard rotamer search method.

High-resolution crystal structure of the restriction-modification controller protein C.AhdI from Aeromonas hydrophila

Restriction-modification (R-M) systems serve to protect the host bacterium from invading bacteriophage. The multi-component system includes a methyltransferase, which recognizes and methylates a specific DNA sequence, and an endonuclease which recognises the same sequence and cleaves within or close to this site. The endonuclease will only cleave DNA that is unmethylated at the specific site, thus host DNA is protected while non-host DNA is cleaved. However, following DNA replication, expression of the endonuclease must be delayed until the host DNA is appropriately methylated. In many R-M systems, this regulation is achieved at the transcriptional level via the controller protein, or C-protein. We have solved the first X-ray structure of an R-M controller protein, C.AhdI, to 1.69 A resolution using selenomethionine MAD. C.AhdI is part of a Type IIH R-M system from the pathogen Aeromonas hydrophila. The structure reveals an all-alpha protein that contains a classical helix-turn-helix (HTH) domain and can be assigned to the Xre family of transcriptional regulators. Unlike its monomeric structural homologues, an extended helix generates an interface that results in dimerisation of the free protein. The dimer is electrostatically polarised and a positively charged surface corresponds to the position of the DNA recognition helices of the HTH domain. Comparison with the structure of the lambda cI ternary complex suggests that C.AhdI activates transcription through direct contact with the sigma70 subunit of RNA polymerase.

Structure and DNA-binding properties of the cytolysin regulator CylR2 from Enterococcus faecalis

Enterococcus faecalis is one of the major causes for hospital-acquired antibiotic-resistant infections. It produces an exotoxin, called cytolysin, which is lethal for a wide range of Gram-positive bacteria and is toxic to higher organisms. Recently, the regulation of the cytolysin operon was connected to autoinduction by a quorum-sensing mechanism involving the CylR1/CylR2 two-component regulatory system. We report here the crystal structure of CylR2 and its properties in solution as determined by heteronuclear NMR spectroscopy. The structure reveals a rigid dimer containing a helix-turn-helix DNA-binding motif as part of a five-helix bundle that is extended by an antiparallel beta-sheet. We show that CylR2 is a DNA-binding protein that binds specifically to a 22 bp fragment of the cytolysin promoter region. NMR chemical shift perturbation experiments identify surfaces involved in DNA binding and are in agreement with a model for the CylR2/DNA complex that attributes binding specificity to a complex network of CylR2/DNA interactions. Our results propose a mechanism where repression is achieved by CylR2 obstruction of the promoter preventing biosynthesis of the cytolysin operon transcript.

Monovalent Cations Regulate DNA Sequence Recognition by 434 Repressor

The bacteriophage 434 repressor distinguishes between its six naturally occurring binding sites using indirect readout. In indirect readout, sequence-dependent differences in the structure and flexibility of non-contacted bases in a protein's DNA-binding site modulate the affinity of DNA for protein. The conformation and flexibility of a DNA sequence can be influenced by the interaction of the DNA bases or backbone with solution components. We examined the effect of changing the cation-type present in solution on the stability and structure of 434 repressor complexes with wild-type and mutant OR1 and OR3, binding sites that differ in their contacted and non-contacted base sequences. We find that the affinity of repressor for OR1, but not for OR3, depends remarkably on the type and concentration of monovalent cation. Moreover, the formation of a stable, specific repressor-OR1 complex requires the presence of monovalent cations; however, repressor-OR3 complex formation has no such requirement. Changing monovalent cation type alters the ability of repressor to protect OR1, but not OR3, from *OH radical cleavage. Altering the relative length of the poly(dA) x poly(dT) tract in the non-contacted regions of the OR1 and OR3 can reverse the cation sensitivity of repressor's affinities for these two sites. Taken together these findings show that cation-dependent alterations in DNA structure underlies indirect readout of DNA sequence by 434 repressor and perhaps other proteins.

Secondary Structure Switching in Cro Protein Evolution

We report the solution structure of the Cro protein from bacteriophage P22. Comparisons of its sequence and structure to those of lambda Cro strongly suggest an alpha-to-beta secondary structure switching event during Cro evolution. The folds of P22 Cro and lambda Cro share a three alpha helix fragment comprising the N-terminal half of the domain. However, P22 Cro's C terminus folds as two helices, while lambda Cro's folds as a beta hairpin. The all-alpha fold found for P22 Cro appears to be ancestral, since it also occurs in cI proteins, which are anciently duplicated paralogues of Cro. PSI-BLAST and transitive homology analyses strongly suggest that the sequences of P22 Cro and lambda Cro are globally homologous despite encoding different folds. The alpha+beta fold of lambda Cro therefore likely evolved from its all-alpha ancestor by homologous secondary structure switching, rather than by nonhomologous replacement of both sequence and structure.

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.

Retroevolution of lambda Cro toward a stable monomer

The Cro protein from bacteriophage lambda has a dimeric alpha+beta fold that evolved from an ancestral all-alpha monomer. The sequence mutations responsible for this dramatic structural evolution are unknown. Here we use analysis of sequence alignments to show that Ala-33, a small side chain in the hydrophobic "ball-and-socket" dimer interface of lambda Cro, was a much larger tryptophan side chain at a previous point in evolution. The retroevolutionary lambda Cro-A33W mutant shows a 10-fold reduction in dimerization affinity relative to the wild type as well as a large increase in monomer thermal stability (Delta T(m) > 10 degrees C), apparently due to partial filling of the hydrophobic socket from within the same monomer. An additional mutation in the dimer interface, F58D, almost completely abolishes detectable dimerization while maintaining the high monomer stability. The secondary structure content of the monomerized versions of lambda Cro is similar to that of the wild-type protein, and the tertiary structure of the monomer appears relatively well defined. These results (i) support a model in which the ball-and-socket dimer interface of lambda Cro was created by altered volume mutations within a limited branch of the Cro lineage and (ii) suggest the possibility that the evolution of the alpha+beta dimer from an all-alpha monomer proceeded through an alpha+beta monomer intermediate.

Structural characterisation and functional significance of transient protein-protein interactions

Protein-protein complexes that dissociate and associate readily, often depending on the physiological condition or environment, play an important role in many biological processes. In order to characterise these "transient" protein-protein interactions, two sets of complexes were collected and analysed. The first set consists of 16 experimentally validated "weak" transient homodimers, which are known to exist as monomers and dimers at physiological concentration, with dissociation constants in the micromolar range. A set of 23 functionally validated transient (i.e. intracellular signalling) heterodimers comprise the second set. This set includes complexes that are more stable, with nanomolar binding affinities, and require a molecular trigger to form and break the interaction. In comparison to more stable homodimeric complexes, the weak homodimers demonstrate smaller contact areas between protomers and the interfaces are more planar and polar on average. The physicochemical and geometrical properties of these weak homodimers more closely resemble those of non-obligate hetero-oligomeric complexes, whose components can exist either as monomers or as complexes in vivo. In contrast to the weak transient dimers, "strong" transient dimers often undergo large conformational changes upon association/dissociation and are characterised with larger, less planar and sometimes more hydrophobic interfaces. From sequence alignments we find that the interface residues of the weak transient homodimers are generally more conserved than surface residues, consistent with being constrained to maintain the protein-protein interaction during evolution. Protein families that include members with different oligomeric states or structures are identified, and found to exhibit a lower sequence conservation at the interface. The results are discussed in terms of the physiological function and evolution of protein-protein interactions.

Ab initio prediction of helical segments in polypeptides

An ab initio method has been developed to predict helix formation for polypeptides. The approach relies on the systematic analysis of overlapping oligopeptides to determine the helical propensity for individual residues. Detailed atomistic level modeling, including entropic contributions, and solvation/ionization energies calculated through the solution of the Poisson-Boltzmann equation, is utilized. The calculation of probabilities for helix formation is based on the generation of ensembles of low energy conformers. The approach, which is easily amenable to parallelization, is shown to perform very well for several benchmark polypeptide systems, including the bovine pancreatic trypsin inhibitor, the immunoglobulin binding domain of protein G, the chymotrypsin inhibitor 2, the R69 N-terminal domain of phage 434 repressor, and the wheat germ agglutinin.

Prolonged stationary-phase incubation selects for lrp mutations in Escherichia coli K-12

Evolution by natural selection occurs in cultures of Escherichia coli maintained under carbon starvation stress. Mutants of increased fitness express a growth advantage in stationary phase (GASP) phenotype, enabling them to grow and displace the parent as the majority population. The first GASP mutation was identified as a loss-of-function allele of rpoS, encoding the stationary-phase global regulator, sigma(S) (M. M. Zambrano, D. A. Siegele, M. A. Almirón, A. Tormo, and R. Kolter, Science 259:1757-1760, 1993). We now report that a second global regulator, Lrp, can also play a role in stationary-phase competition. We found that a mutant that took over an aged culture of an rpoS strain had acquired a GASP mutation in lrp. This GASP allele, lrp-1141, encodes a mutant protein lacking the critical glycine in the turn of the helix-turn-helix DNA-binding domain. The lrp-1141 allele behaves as a null mutation when in single copy and is dominant negative when overexpressed. Hence, the mutant protein appears to retain stability and the ability to dimerize but lacks DNA-binding activity. We also demonstrated that a lrp null allele generated by a transposon insertion has a fitness gain identical to that of the lrp-1141 allele, verifying that cells lacking Lrp activity have a competitive advantage during prolonged starvation. Finally, we tested by genetic analysis the hypothesis that the lrp-1141 GASP mutation confers a fitness gain by enhancing amino acid catabolism during carbon starvation. We found that while amino acid catabolism may play a role, it is not necessary for the lrp GASP phenotype, and hence the lrp GASP phenotype is due to more global physiological changes.

Constructing side chains on near-native main chains for ab initio protein structure prediction

Is there value in constructing side chains while searching protein conformational space during an ab initio simulation? If so, what is the most computationally efficient method for constructing these side chains? To answer these questions, four published approaches were used to construct side chain conformations on a range of near-native main chains generated by ab initio protein structure prediction methods. The accuracy of these approaches was compared with a naive approach that selects the most frequently observed rotamer for a given amino acid to construct side chains. An all-atom conditional probability discriminatory function is useful at selecting conformations with overall low all-atom root mean square deviation (r.m.s.d.) and the discrimination improves on sets that are closer to the native conformation. In addition, the naive approach performs as well as more sophisticated methods in terms of the percentage of chi(1) angles built accurately and the all-atom r. m.s.d., between the native and near-native conformations. The results suggest that the naive method would be extremely useful for fast and efficient side chain construction on vast numbers of conformations for ab initio prediction of protein structure.

DNA-induced conformational changes in bacteriophage 434 repressor

Although bacteriophage 434 repressor binds to its specific DNA sites only as a dimer, formation of the dimers in solution occurs at concentrations three orders of magnitude higher than those needed to bind the 434 operator DNA. Our results suggest that both specific and non-specific DNA induce conformational changes in repressor that lead to formation of repressor dimers. The repressor conformational changes induced by DNA occur at concentrations much lower than those needed for binding of repressor, suggesting that the alternative conformations of repressor persist even if the protein is not in direct contact with DNA. Hence, DNA acts in a "catalytic" fashion to induce a steady-state amount of an alternative repressor conformation that has an enhanced affinity for its specific binding site. These findings suggest that the repressor conformer induced by non-specific DNA is the form of the repressor that is optimized for searching for DNA binding sites along non-specific DNA. Upon finding a binding site, the repressor protein undergoes an additional conformational change that allows it to "lock-on" to its specific site.

Folding propensities of synthetic peptide fragments covering the entire sequence of phage 434 Cro protein

The phage 434 Cro protein, the N-terminal domain of its repressor (R1-69) and that of phage lambda (lambda6-85) constitute a group of small, monomeric, single-domain folding units consisting of five helices with striking structural similarity. The intrinsic helix stabilities in lambda6-85 have been correlated to its rapid folding behavior, and a residual hydrophobic cluster found in R1-69 in 7 M urea has been proposed as a folding initiation site. To understand the early events in the folding of 434 Cro, and for comparison with R1-69 and lambda6-85, we examined the conformational behavior of five peptides covering the entire 434 Cro sequence in water, 40% (by volume) TFE/water, and 7 M urea solutions using CD and NMR. Each peptide corresponds to a helix and adjacent residues as identified in the native 434 Cro NMR and crystal structures. All are soluble and monomeric in the solution conditions examined except for the peptide corresponding to the 434 Cro helix 4, which has low water solubility. Helix formation is observed for the 434 Cro helix 1 and helix 2 peptides in water, for all the peptides in 40% TFE and for none in 7 M urea. NMR data indicate that the helix limits in the peptides are similar to those in the native protein helices. The number of side-chain NOEs in water and TFE correlates with the helix content, and essentially none are observed in 7 M urea for any peptide, except that for helix 5, where a hydrophobic cluster may be present. The low intrinsic folding propensities of the five helices could account for the observed stability and folding behavior of 434 Cro and is, at least qualitatively, in accord with the results of the recently described diffusion-collision model incorporating intrinsic helix propensities.

Molecular dynamics simulation in solvent of the bacteriophage 434 cI repressor protein DNA binding domain amino acids (R1-69) in complex with its cognate operator (OR1) DNA sequence

We investigated protein/DNA interactions, using molecular dynamics simulations computed between a 10 Angstom water layer model of the 434 cI Repressor protein DNA binding domain (DBD) amino acids (R1-69) and DNA of operator (OR1) and its flanks consisting of 28 nucleotide base pairs. Hydrogen bonding interactions were monitored. In addition, van der Waals and electrostatic interaction energies were calculated. Amino acids of the 434 cI repressor DNA recognition helix 3 formed both direct and water mediated hydrogen bonds at cognate codon-anticodon nucleotide base and backbone sites within the OR1 DNA major groove halfsites and flanking regions. In addition, hydrophilic amino acids within the loop between helix 3 and helix 4 have strong electrostatic attraction to codon-anticodon nucleotides located within the central nucleotides of the minor groove between the OR1 major groove halfsites. These interactions together induced significant structural changes in the operator DNA manifested by overtwisting of the central nucleotide base pairs and narrowing of the minor groove between the DNA major groove halfsites. Finally, these findings offer a code for site specific DNA recognition by the 434 cI repressor protein.

Side-chain and backbone flexibility in protein core design

We have developed a computational approach for the design and prediction of hydrophobic cores that includes explicit backbone flexibility. The program consists of a two-stage combination of a genetic algorithm and monte carlo sampling using a torsional model of the protein. Backbone structures are evaluated either by a canonical force-field or a constraining potential that emphasizes the preservation of local geometry. The utility of the method for protein design and engineering is explored by designing three novel hydrophobic core variants of the protein 434 cro. We use the new method to evaluate these and previously designed 434 cro variants, as well as a series of phage T4 lysozyme variants. In order to properly evaluate the influence of backbone flexibility, we have also analyzed the effects of varying amounts of side-chain flexibility on the performance of fixed backbone methods. Comparison of results using a fixed versus flexible backbone reveals that, surprisingly, the two methods are almost equivalent in their abilities to predict relative experimental stabilities, but only when full side-chain flexibility is allowed. The prediction of core side-chain structure can vary dramatically between methods. In some, but not all, cases the flexible backbone method is a better predictor of structure. The development of a flexible backbone approach to core design is particularly important for attempts at de novo protein design, where there is no prior knowledge of a precise backbone structure.

A thermodynamic study of the 434-repressor N-terminal domain and of its covalently linked dimers

The isolated N-terminal 1-69 domain of the 434-phage repressor, R69, and its covalently linked (head-to-tail and tail-to-tail) dimers have been studied by differential scanning microcalorimetry (DSC) and CD. At neutral solvent conditions the R69 domain maintains its native structure, both in isolated form and within the dimers. The stability of the domain depends highly upon pH within the acidic range, thus at pH 2 and low ionic strength R69 is already partially unfolded at room temperature. The thermodynamic parameters of unfolding calculated from the DSC data are typical for small globular proteins. At neutral pH and moderate ionic strength, the domains of the dimers behave as two independent units with unfolding parameters similar to those of the isolated domain, which means that linking two R69 domains, either by a long peptide linker or by a designed C-terminal disulfide bridge, does not induce any cooperation between them.

Validation of NMR side-chain conformations by packing calculations

The precision and accuracy of protein structures determined by nuclear magnetic resonance (NMR) spectroscopy depend on the completeness of input experimental data set. Typically, rather than a single structure, an ensemble of up to 20 equally representative conformers is generated and routinely deposited in the Protein Database. There are substantially more experimentally derived restraints available to define the main-chain coordinates than those of the side chains. Consequently, the side-chain conformations among the conformers are more variable and less well defined than those of the backbone. Even when a side chain is determined with high precision and is found to adopt very similar orientations among all the conformers in the ensemble, it is possible that its orientation might still be incorrect. Thus, it would be helpful if there were a method to assess independently the side-chain orientations determined by NMR. Recently, homology modeling by side-chain packing algorithms has been shown to be successful in predicting the side-chain conformations of the buried residues for a protein when the main-chain coordinates and sequence information are given. Since the main-chain coordinates determined by NMR are consistently more reliable than those of the side-chains, we have applied the side-chain packing algorithms to predict side-chain conformations that are compatible with the NMR-derived backbone. Using four test cases where the NMR solution structures and the X-ray crystal structure of the same protein are available, we demonstrate that the side-chain packing method can provide independent validation for the side-chain conformations of NMR structures. Comparison of the side-chain conformations derived by side-chain packing prediction and by NMR spectroscopy demonstrates that when there is agreement between the NMR model and the predicted model, on average 78% of the time the X-ray structure also concurs. While the side-chain packing method can confirm the reliable residue conformations in NMR models, more importantly, it can also identify the questionable residue conformations with an accuracy of 60%. This validation method can serve to increase the confidence level for potential users of structural models determined by NMR.

An evolutionary link between sporulation and prophage induction in the structure of a repressor:anti-repressor complex

Spore formation is an extreme response of some bacteria to adversity. In Bacillus subtilis the proteins of the sin, sporulation inhibition, region form a component of an elaborate molecular circuitry that regulates the commitment to sporulation. SinR is a tetrameric repressor protein that binds to the promoters of genes essential for entry into sporulation and prevents their transcription. This repression is overcome through the activity of SinI, which disrupts the SinR tetramer through the formation of a SinI-SinR heterodimer. The interactions governing this curious quaternary transition are revealed in the crystal structure of the SinI-SinR complex. The most striking, and unexpected, finding is that the tertiary structure of the DNA-binding domain of SinR is identical with that of the corresponding domains of the repressor proteins, CI and Cro, of bacteriophage 434 that regulate lysis/lysogeny. This structural similarity greatly exceeds that between SinR and any bacterial protein or between the 434 repressor proteins and their homologues in the closely related bacteriophage lambda. The close evolutionary relationship implied by the structures of SinR and the 434 repressors provokes both comparison of their functions and a speculative consideration of the intriguing possibility of an evolutionary link between the two adaptive responses, sporulation and prophage induction.

Carboxyl-terminal domain dimer interface mutant 434 repressors have altered dimerization and DNA binding specificities

Strong dimerization of the repressor, mediated by the carboxyl (C)-terminal domain, is a prerequisite for forming a specific complex with DNA and cooperative DNA binding to form tetramers. We have generated a computer model of the C-terminal domain of the 434 repressor based on the crystal structure of the homologous UmuD' protein. This model predicts that residues in the primary sequence between 93 and 168 contribute to the dimer interface. We changed several amino acid residues located in this region. Gel filtration and crosslinking assays were used to characterize the strength and specificity of dimerization of the purified repressor C-terminal domain dimer interface mutants. These results indicate that amino acid residues K121, H139, D161 and N163 contribute to the strength and/or specificity of dimerization. The relative affinity of the bacteriophage 434 repressor for 434 operators is determined, in part, by the repressor's ability to detect sequence-dependent structural alterations in the non-contacted region at the center of an operator site. We find that the relative ability of C-terminal domain dimer interface mutant repressors to dimerize does not necessarily predict their relative abilities to bind DNA, and that these proteins are deficient in detecting non-contacted base-dependent differences in operator strength. Our results show that the structure of the DNA in complex with these mutant proteins differs from that found in wild-type repressor-operator complexes, even though the sites of these mutations lie in a separate domain from that which contacts the DNA. These observations demonstrate that the structural integrity of the C-terminal domain dimer interface is required to appropriately orient the DNA binding information contained within the DNA-contacting N-terminal domain.

Sequence specificity of bacteriophage 434 repressor-operator complexation

The binding affinity of the bacteriophage 434 repressor for its DNA operator depends strongly on the nature of two central base-pairs that are not in contact with the dimeric protein. In order to investigate the origin of this sequence specificity, we carried out molecular modelling of five model operators with central TA, AT, CG, GC and IC sequences. The five oligomers were studied both before and after complexation with the N-terminal binding domain of the 434 repressor. The relative importance of nucleic acid flexibility on operator-repressor binding was studied via a low frequency normal mode analysis using an internal coordinate method that we developed recently. The results suggest a higher twisting flexibility for TA and AT central steps than for CG, GC or IC steps, but the differences appear to be too small to account for the strength of repressor binding. An energetic analysis of the model operator-repressor complexes reveals rather that the preference for A.T and T.A base-pairs is electrostatic in origin and is linked to the presence of cationic Arg43 side-chains of repressor. This conclusion is supported by comparison with an R43A mutant and correlates with the sequence dependence of the electrostatic potential in the central minor groove of the operators.

Molecular ecology and evolution of Streptococcus thermophilus bacteriophages--a review

Bacteriophages attacking Streptococcus thermophilus, a lactic acid bacterium used in milk fermentation, are a threat to the dairy industry. These small isometric-headed phages possess double-stranded DNA genomes of 31 to 45 kb. Yoghurt-derived phages exhibit a limited degree of variability, as defined by restriction pattern and host range, while a large diversity of phage types have been isolated from cheese factories. Despite this diversity all S. thermophilus phages, virulent and temperate, belong to a single DNA homology group. Several mechanisms appear to create genetic variability in this phage group. Site-specific deletions, one type possibly mediated by a viral recombinase/integrase, which transformed a temperate into a virulent phage, were observed. Recombination as a result of superinfection of a lysogenic host has been reported. Comparative DNA sequencing identified up to 10% sequence diversity due to point mutations. Genome sequencing of the prototype temperate phage phi Sfi21 revealed many predicted proteins which showed homology with phages from Lactococcus lactis suggesting horizontal gene transfer. Homology with phages from evolutionary unrelated bacteria like E. coli (e.g. lambdoid phage 434 and P1) and Mycobacterium phi L5 was also found. Due to their industrial importance, the existence of large phage collections, and the whole phage genome sequencing projects which are currently underway, the S. thermophilus phages may present an interesting experimental system to study bacteriophage evolution.

Alignment algorithm for homology modeling and threading

A DNA/protein sequence comparison is a popular computational tool for molecular biologists. Finding a good alignment implies an evolutionary and/or functional relationship between proteins or genomic loci. Sequential similarity between two proteins indicates their structural resemblance, providing a practical approach for structural modeling, when structure of one of these proteins is known. The first step in the homology modeling is a construction of an accurate sequence alignment. The commonly used alignment algorithms do not provide an adequate treatment of the structurally mismatched residues in locally dissimilar regions. We propose a simple modification of the existing alignment algorithm which treats these regions properly and demonstrate how this modification improves sequence alignments in real proteins.

Thermodynamic model of secondary structure for alpha-helical peptides and proteins

A thermodynamic model describing formation of alpha-helices by peptides and proteins in the absence of specific tertiary interactions has been developed. The model combines free energy terms defining alpha-helix stability in aqueous solution and terms describing immersion of every helix or fragment of coil into a micelle or a nonpolar droplet created by the rest of protein to calculate averaged or lowest energy partitioning of the peptide chain into helical and coil fragments. The alpha-helix energy in water was calculated with parameters derived from peptide substitution and protein engineering data and using estimates of nonpolar contact areas between side chains. The energy of nonspecific hydrophobic interactions was estimated considering each alpha-helix or fragment of coil as freely floating in the spherical micelle or droplet, and using water/cyclohexane (for micelles) or adjustable (for proteins) side-chain transfer energies. The model was verified for 96 and 36 peptides studied by 1H-nmr spectroscopy in aqueous solution and in the presence of micelles, respectively ([set 1] and [set 2]) and for 30 mostly alpha-helical globular proteins ([set 3]). For peptides, the experimental helix locations were identified from the published medium-range nuclear Overhauser effects detected by 1H-nmr spectroscopy. For sets 1, 2, and 3, respectively, 93, 100, and 97% of helices were identified with average errors in calculation of helix boundaries of 1.3, 2.0, and 4.1 residues per helix and an average percentage of correctly calculated helix-coil states of 93, 89, and 81%, respectively. Analysis of adjustable parameters of the model (the entropy and enthalpy of the helix-coil transition, the transfer energy of the helix backbone, and parameters of the bound coil), determined by minimization of the average helix boundary deviation for each set of peptides or proteins, demonstrates that, unlike micelles, the interior of the effective protein droplet has solubility characteristics different from that for cyclohexane, does not bind fragments of coil, and lacks interfacial area.

Characterization of the lysogeny DNA module from the temperate Streptococcus thermophilus bacteriophage phi Sfi21

Phage phi Sfi21, the only temperate Streptococcus thermophilus phage from our phage collection, showed extensive DNA homology with virulent phages from lytic group I. Southern blot hybridizations demonstrated that the phi Sfi21-specific DNA was clustered in an approximately 6.6-kb-long region, the putative lysogeny module. Sequence analysis and database research identified an integrase within this module; orf 203 with homology to an anonymous orf 258 from the temperate lactococcal phage BK5-T; orf 127 and orf 122 with weak homology to the N- and C-terminal parts, respectively, of the cl-like repressor from lactococcal phages Tuc2009 and BK5-T; orf 75 with homology to a repressor protein from lambdoid phage 434 and an anti-repressor ant with homology to phage P1. The molecular arrangement of the predicted orfs in phage phi Sfi21 was very similar to that of the lactococcal phage BK5-T. The transition from phi Sfi21-specific DNA into DNA shared with virulent phages was abrupt and flanked at one side by notable DNA repeats. Sequence analysis identified a holin protein to the left of the lysogeny module. A site-specific deletion of 2.4 kb, which reproducibly transformed phi Sfi21 into a lytic phage, was localized in the lysogeny module. It was flanked at both sides by conspicuous DNA repeats. One repeat region reflected the DNA around the attP site, while the other reflected the putative genetic switch region between repressor and anti-repressor genes. S. thermophilus host Sfi1 transformed with a plasmid containing int and orf 203 showed resistance to superinfection by heterologous phages, but not by the homologous phi Sfi21. Part of the int gene could be deleted without loss of this activity, while a deletion in orf 203 resulted in loss of the phage resistance. We speculate on the possibility of a bipartite immunity system for the control of lysogeny in phi Sfi21.