Using DelPhi to compute electrostatic potentials and assess their contribution to interactions

The Vibrio cholerae minor pilin TcpB mediates uptake of the cholera toxin phage CTXφ

Virulent strains of the bacterial pathogen Vibrio cholerae cause the diarrheal disease cholera by releasing cholera toxin into the small intestine. V. cholerae acquired its cholera toxin genes by lysogenic infection with the filamentous bacteriophage CTXφ. CTXφ uses its minor coat protein pIII, located in multiple copies at the phage tip, to bind to the V. cholerae toxin-coregulated pilus (TCP). However, the molecular details of this interaction and the mechanism of phage internalization are not well-understood. The TCP filament is a polymer of major pilins, TcpA, and one or more minor pilin, TcpB. TCP are retractile, with both retraction and assembly initiated by TcpB. Consistent with these roles in pilus dynamics, we hypothesized that TcpB controls both binding and internalization of CTXφ. To test this hypothesis, we determined the crystal structure of the C-terminal half of TcpB and characterized its interactions with CTXφ pIII. We show that TcpB is a homotrimer in its crystallographic form as well as in solution and is present in multiple copies at the pilus tip, which likely facilitates polyvalent binding to pIII proteins at the phage tip. We further show that recombinant forms of TcpB and pIII interact in vitro, and both TcpB and anti-TcpB antibodies block CTXφ infection of V. cholerae Finally, we show that CTXφ uptake requires TcpB-mediated retraction. Our data support a model whereby CTXφ and TCP bind in a tip-to-tip orientation, allowing the phage to be drawn into the V. cholerae periplasm as an extension of the pilus filament.

Comparative Analysis of Electrostatic Models for Ligand Docking

The precise modeling of molecular interactions remains an important goal among molecular modeling techniques. Some of the challenges in the field include the precise definition of a Hamiltonian for biomolecular systems, together with precise parameters derived from Molecular Mechanics Force Fields, for example. The problem is even more challenging when interaction energies from different species are computed, such as the interaction energy involving a ligand and a protein, given that small differences must be computed from large energies. Here we evaluated the effects of the electrostatic model for ligand binding energy evaluation in the context of ligand docking. For this purpose, a classical Coulomb potential with distance-dependent dielectrics was compared with a Poisson-Boltzmann (PB) model for electrostatic potential computation, based on DelPhi calculations. We found that, although the electrostatic energies were highly correlated for the Coulomb and PB models, the ligand pose and the enrichment of actual ligands against decoy compounds, were improved when binding energies were computed using PB as compared to the Coulomb model. We observed that the electrostatic energies computed with the Coulomb model were, on average, ten times larger than the energies computed with the PB model, suggesting a strong overestimation of the polar interactions in the Coulomb model. We also found that a slightly smoothed Lennard-Jones potential combined with the PB model resulted in a good compromise between ligand sampling and energetic scoring.

Functional interplay between the RK motif and linker segment dictates Oct4-DNA recognition

The POU family transcription factor Oct4 plays pivotal roles in regulating pluripotency and somatic cell reprogramming. Previous studies have indicated an important role for major groove contacts in Oct4–DNA recognition; however, the contributions of the RK motif in the POUh domain and the linker segment joining the two DNA-binding domains remain poorly understood. Here, by combining molecular modelling and functional assays, we find that the RK motif is essential for Oct4–DNA association by recognizing the narrowed DNA minor groove. Intriguingly, computational simulations reveal that the function of the RK motif may be finely tuned by H-bond interactions with the partially disordered linker segment and that breaking these interactions significantly enhances the DNA binding and reprogramming activities of Oct4. These findings uncover a self-regulatory mechanism for specific Oct4–DNA recognition and provide insights into the functional crosstalk at the molecular level that may illuminate mechanistic studies of the Oct protein family and possibly transcription factors in the POU family. Our gain-of-function Oct4 mutants might also be useful tools for use in reprogramming and regenerative medicine.

Structural insights of the ssDNA binding site in the multifunctional endonuclease AtBFN2 from Arabidopsis thaliana

The multi S1/P1 nuclease AtBFN2 (EC 3.1.30.1) encoded by the Arabidopsis thaliana At1g68290 gene is a glycoprotein that digests RNA, ssDNA, and dsDNA. AtBFN2 depends on three zinc ions for cleaving DNA and RNA at 3′-OH to yield 5′-nucleotides. In addition, AtBFN2′s enzymatic activity is strongly glycan dependent. Plant Zn²⁺-dependent endonucleases present a unique fold, and belong to the Phospholipase C (PLC)/P1 nuclease superfamily. In this work, we present the first complete, ligand-free, AtBFN2 crystal structure, along with sulfate, phosphate and ssDNA co-crystal structures. With these, we were able to provide better insight into the glycan structure and possible enzymatic mechanism. In comparison with other nucleases, the AtBFN2/ligand-free and AtBFN2/PO₄ models suggest a similar, previously proposed, catalytic mechanism. Our data also confirm that the phosphate and vanadate can inhibit the enzyme activity by occupying the active site. More importantly, the AtBFN2/A₅T structure reveals a novel and conserved secondary binding site, which seems to be important for plant Zn²⁺-dependent endonucleases. Based on these findings, we propose a rational ssDNA binding model, in which the ssDNA wraps itself around the protein and the attached surface glycan, in turn, reinforces the binding complex.

Structural characterization of CFA/III and Longus type IVb pili from enterotoxigenic Escherichia coli

The type IV pili are helical filaments found on many Gram-negative pathogenic bacteria, with multiple diverse roles in pathogenesis, including microcolony formation, adhesion, and twitching motility. Many pathogenic enterotoxigenic Escherichia coli (ETEC) isolates express one of two type IV pili belonging to the type IVb subclass: CFA/III or Longus. Here we show a direct correlation between CFA/III expression and ETEC aggregation, suggesting that these pili, like the Vibrio cholerae toxin-coregulated pili (TCP), mediate microcolony formation. We report a 1.26-Å resolution crystal structure of CofA, the major pilin subunit from CFA/III. CofA is very similar in structure to V. cholerae TcpA but possesses a 10-amino-acid insertion that replaces part of the α2-helix with an irregular loop containing a 3(10)-helix. Homology modeling suggests a very similar structure for the Longus LngA pilin. A model for the CFA/III pilus filament was generated using the TCP electron microscopy reconstruction as a template. The unique 3(10)-helix insert fits perfectly within the gap between CofA globular domains. This insert, together with differences in surface-exposed residues, produces a filament that is smoother and more negatively charged than TCP. To explore the specificity of the type IV pilus assembly apparatus, CofA was expressed heterologously in V. cholerae by replacing the tcpA gene with that of cofA within the tcp operon. Although CofA was synthesized and processed by V. cholerae, no CFA/III filaments were detected, suggesting that the components of the type IVb pilus assembly system are highly specific to their pilin substrates.

Stability of domain structures in multi-domain proteins

Multi-domain proteins have many advantages with respect to stability and folding inside cells. Here we attempt to understand the intricate relationship between the domain-domain interactions and the stability of domains in isolation. We provide quantitative treatment and proof for prevailing intuitive ideas on the strategies employed by nature to stabilize otherwise unstable domains. We find that domains incapable of independent stability are stabilized by favourable interactions with tethered domains in the multi-domain context. Stability of such folds to exist independently is optimized by evolution. Specific residue mutations in the sites equivalent to inter-domain interface enhance the overall solvation, thereby stabilizing these domain folds independently. A few naturally occurring variants at these sites alter communication between domains and affect stability leading to disease manifestation. Our analysis provides safe guidelines for mutagenesis which have attractive applications in obtaining stable fragments and domain constructs essential for structural studies by crystallography and NMR.

On the role of electrostatics in protein-protein interactions

The role of electrostatics in protein-protein interactions and binding is reviewed in this paper. A brief outline of the computational modeling, in the framework of continuum electrostatics, is presented and the basic electrostatic effects occurring upon the formation of the complex are discussed. The effect of the salt concentration and pH of the water phase on protein-protein binding free energy is demonstrated which indicates that the increase of the salt concentration tends to weaken the binding, an observation that is attributed to the optimization of the charge-charge interactions across the interface. It is pointed out that the pH-optimum (pH of optimal binding affinity) varies among the protein-protein complexes, and perhaps is a result of their adaptation to particular subcellular compartments. The similarities and differences between hetero- and homo-complexes are outlined and discussed with respect to the binding mode and charge complementarity.

Crystal structure of the Vibrio cholerae colonization factor TcpF and identification of a functional immunogenic site

Vibrio cholerae relies on two main virulence factors--toxin-coregulated pilus (TCP) and cholera toxin--to cause the gastrointestinal disease cholera. TCP is a type IV pilus that mediates bacterial autoagglutination and colonization of the intestine. TCP is encoded by the tcp operon, which also encodes TcpF, a protein of unknown function that is secreted by V. cholerae in a TCP-dependent manner. Although TcpF is not required for TCP biogenesis, a tcpF mutant has a colonization defect in the infant mouse cholera model that is as severe as a pilus mutant. Furthermore, TcpF antisera protect against V. cholerae infection. TcpF has no apparent sequence homology to any known protein. Here, we report the de novo X-ray crystal structure of TcpF and the identification of an epitope that is critical for its function as a colonization factor. A monoclonal antibody recognizing this epitope is protective against V. cholerae challenge and adds to the protection provided by an anti-TcpA antibody. These data suggest that TcpF has a novel function in V. cholerae colonization and define a region crucial for this function.

Structure of a trimeric bacterial microcompartment shell protein, EtuB, associated with ethanol utilization in Clostridium kluyveri

It has been suggested that ethanol metabolism in the strict anaerobe Clostridium kluyveri occurs within a metabolosome, a subcellular proteinaceous bacterial microcompartment. Two bacterial microcompartment shell proteins [EtuA (ethanol utilization shell protein A) and EtuB] are found encoded on the genome clustered with the genes for ethanol utilization. The function of the bacterial microcompartment is to facilitate fermentation by sequestering the enzymes, substrates and intermediates. Recent structural studies of bacterial microcompartment proteins have revealed both hexamers and pentamers that assemble to generate the pseudo-icosahedral bacterial microcompartment shell. Some of these shell proteins have pores on their symmetry axes. Here we report the structure of the trimeric bacterial microcompartment protein EtuB, which has a tandem structural repeat within the subunit and pseudo-hexagonal symmetry. The pores in the EtuB trimer are within the subunits rather than between symmetry related subunits. We suggest that the evolutionary advantage of this is that it releases the pore from the rotational symmetry constraint allowing more precise control of the fluxes of asymmetric molecules, such as ethanol, across the pore. We also model EtuA and demonstrate that the two proteins have the potential to interact to generate the casing for a metabolosome.

Comparing 2-nt 3' overhangs against blunt-ended siRNAs: a systems biology based study

In this study, we formulate a computational reaction model following a chemical kinetic theory approach to predict the binding rate constant for the siRNA-RISC complex formation reaction. The model allowed us to study the potency difference between 2-nt 3' overhangs against blunt-ended siRNA molecules in an RNA interference (RNAi) system. The rate constant predicted by this model was fed into a stochastic simulation of the RNAi system (using the Gillespie stochastic simulator) to study the overall potency effect. We observed that the stochasticity in the transcription/translation machinery has no observable effects in the RNAi pathway. Sustained gene silencing using siRNAs can be achieved only if there is a way to replenish the dsRNA molecules in the cell. Initial findings show about 1.5 times more blunt-ended molecules will be required to keep the mRNA at the same reduced level compared to the 2-nt overhang siRNAs. However, the mRNA levels jump back to saturation after a longer time when blunt-ended siRNAs are used. We found that the siRNA-RISC complex formation reaction rate was 2 times slower when blunt-ended molecules were used pointing to the fact that the presence of the 2-nt overhangs has a greater effect on the reaction in which the bound RISC complex cleaves the mRNA.