Protein-protein Docking Simulations Research Articles

POP1 might be recruiting its type-Ia interface for NLRP3-mediated PYD-PYD interaction: Insights from MD simulation.

Inflammasomes are multiprotein caspase-activating complexes that enhance the maturation and release of proinflammatory cytokines (IL-1β and IL-18) in response to the invading pathogen and/or host-derived cellular stress. These are assembled by the sensory proteins (viz NLRC4, NLRP1, NLRP3, and AIM-2), adaptor protein (ASC), and effector molecule procaspase-1. In NLRP3-mediated inflammasome activation, ASC acts as a mediator between NLRP3 and procaspase-1 for the transmission of signals. A series of homotypic protein-protein interactions (NLRP3PYD :ASCPYD and ASCCARD :CASP1CARD ) propagates the downstream signaling for the production of proinflammatory cytokines. Pyrin-only protein 1 (POP1) is known to act as the regulator of inflammasome. It modulates the ASC-mediated inflammasome assembly by interacting with pyrin domain (PYD) of ASC. However, despite similar electrostatic surface potential, the interaction of POP1 with NLRP3PYD is obscured till date. Herein, to explore the possible PYD-PYD interactions between NLRP3PYD and POP1, a combined approach of protein-protein docking and molecular dynamics simulation was adapted. The current study revealed that POP1's type-Ia interface and type-Ib interface of NLRP3PYD might be crucial for 1:1 PYD-PYD interaction. In addition to type-I mode of interaction, we also observed type-II and type-III interaction modes in two different dynamically stable heterotrimeric complexes (POP1-NLRP3-NLRP3 and POP1-NLRP3-POP1). The inter-residual/atomic distance calculation exposed several critical residues that possibly govern the said interaction, which need further investigation. Overall, the findings of this study will shed new light on hitherto concealed molecular mechanisms underlying NLRP3-mediated inflammasome, which will have strong future therapeutic implications.

Molecular and structural characterization of novel cystatins from the taiga tick Ixodes persulcatus

Cystatins are cysteine peptidase inhibitors that in ticks mediate processes such as blood feeding and digestion. The ixodid tick Ixodes persulcatus is endemic to the Eurasia, where it is the principal vector of Lyme borreliosis. To date, no I. persulcatus cystatin has been characterized. In the present work, we describe three novel cystatins from I. persulcatus, named JpIpcys2a, JpIpcys2b and JpIpcys2c. In addition, the potential of tick cystatins as cross-protective antigens was evaluated by vaccination of hamsters using BrBmcys2c, a cystatin from Rhipicephalus microplus, against I. persulcatus infestation. Sequence analysis showed that motifs that are characteristic of cystatins type 2 are fully conserved in JpIpcys2b, while mutations are present in both JpIpcys2a and JpIpcys2c. Protein-protein docking simulations further revealed that JpIpcys2a, JpIpcys2b and JpIpcys2c showed conserved binding sites to human cathepsins L, all of them covering the active site cleft. Cystatin transcripts were detected in different I. persulcatus tissues and instars, showing their ubiquitous expression during I. persulcatus development. Serological analysis showed that although hamsters immunized with BrBmcys2c developed a humoral immune response, this response was not adequate to protect against a heterologous challenge with I. persulcatus adult ticks. The lack of cross-protection provided by BrBmcys2c immunization is perhaps linked to the fact that cystatins cluster into multigene protein families that are expressed differentially and exhibit functional redundancy. How to target such small proteins that are secreted in low quantities remains a challenge in the development of suitable anti-tick vaccine antigens.

Computational and biochemical characterization of two partially overlapping interfaces and multiple weak-affinity K-Ras dimers

Recent studies found that membrane-bound K-Ras dimers are important for biological function. However, the structure and thermodynamic stability of these complexes remained unknown because they are hard to probe by conventional approaches. Combining data from a wide range of computational and experimental approaches, here we describe the structure, dynamics, energetics and mechanism of assembly of multiple K-Ras dimers. Utilizing a range of techniques for the detection of reactive surfaces, protein-protein docking and molecular simulations, we found that two largely polar and partially overlapping surfaces underlie the formation of multiple K-Ras dimers. For validation we used mutagenesis, electron microscopy and biochemical assays under non-denaturing conditions. We show that partial disruption of a predicted interface through charge reversal mutation of apposed residues reduces oligomerization while introduction of cysteines at these positions enhanced dimerization likely through the formation of an intermolecular disulfide bond. Free energy calculations indicated that K-Ras dimerization involves direct but weak protein-protein interactions in solution, consistent with the notion that dimerization is facilitated by membrane binding. Taken together, our atomically detailed analyses provide unique mechanistic insights into K-Ras dimer formation and membrane organization as well as the conformational fluctuations and equilibrium thermodynamics underlying these processes.

Subunit orientation in the Escherichia coli enterobactin biosynthetic EntA-EntE complex revealed by a two-hybrid approach.

The siderophore enterobactin is synthesized by the enzymes EntA-F and EntH in the Escherichia coli cytoplasm. We previously reported invitro evidence of an interaction between tetrameric EntA and monomeric EntE. Here we used bacterial adenylate cyclase two-hybrid (BACTH) assays to demonstrate that the E.coli EntA-EntE interaction occurs intracellularly. Furthermore, to obtain information on subunit orientation in the EntA-EntE complex, we fused BACTH reporter fragments T18 and T25 to EntA and EntE in both N-terminal and C-terminal orientations. To validate functionality of our fusion proteins, we performed Chrome Azurol S (CAS) assays using E.coli entE(-) and entA(-) knockout strains transformed with our BACTH constructs. We found that transformants expressing N-terminal and C-terminal T18/T25 fusions to EntE exhibited CAS signals, indicating that these constructs could rescue the entE(-) phenotype. While expression of EntA with N-terminal T18/T25 fusions exhibited CAS signals, C-terminal fusions did not, presumably due to disruption of the EntA tetramer invivo. Bacterial growth assays supported our CAS findings. Co-transformation of functional T18/T25 fusions into cya(-)E.coli BTH101cells resulted in positive BACTH signals only when T18/T25 fragments were fused to the N-termini of both EntA and EntE. Co-expression of N-terminally fused EntA with C-terminally fused EntE resulted in no detectable BACTH signal. Analysis of protein expression by Western blotting confirmed that the loss of BACTH signal was not due to impaired expression of fusion proteins. Based on our results, we propose that the N-termini of EntA and EntE are proximal in the intracellular complex, while the EntA N-terminus and EntE C-terminus are distal. A protein-protein docking simulation using SwarmDock was in agreement with our experimental observations.

Identification of TAX2 peptide as a new unpredicted anti-cancer agent.

The multi-modular glycoprotein thrombospondin-1 (TSP-1) is considered as a key actor within the tumor microenvironment. Besides, TSP-1 binding to CD47 is widely reported to regulate cardiovascular function as it promotes vasoconstriction and angiogenesis limitation. Therefore, many studies focused on targeting TSP-1:CD47 interaction, aiming for up-regulation of physiological angiogenesis to enhance post-ischemia recovery or to facilitate engraftment. Thus, we sought to identify an innovative selective antagonist for TSP-1:CD47 interaction. Protein-protein docking and molecular dynamics simulations were conducted to design a novel CD47-derived peptide, called TAX2. TAX2 binds TSP-1 to prevent TSP-1:CD47 interaction, as revealed by ELISA and co-immunoprecipitation experiments. Unexpectedly, TAX2 inhibits in vitro and ex vivo angiogenesis features in a TSP-1-dependent manner. Consistently, our data highlighted that TAX2 promotes TSP-1 binding to CD36-containing complexes, leading to disruption of VEGFR2 activation and downstream NO signaling. Such unpredicted results prompted us to investigate TAX2 potential in tumor pathology. A multimodal imaging approach was conducted combining histopathological staining, MVD, MRI analysis and μCT monitoring for tumor angiography longitudinal follow-up and 3D quantification. TAX2 in vivo administrations highly disturb syngeneic melanoma tumor vascularization inducing extensive tumor necrosis and strongly inhibit growth rate and vascularization of human pancreatic carcinoma xenografts in nude mice.

Protein-Protein Docking and Molecular Dynamics Simulations Elucidated Binding Modes of FUBI-p62 UBA Complex

Protein-Protein Docking and Molecular Dynamics Simulations Elucidated Binding Modes of FUBI-p62 UBA Complex

Molecular Dynamics and Regulation of Butyrylcholinesterase Cholinergic Activity by RNA Binding Proteins

RNA recognition motif (RRM) is the most abundant RNA-binding domain, mainly known to involve in regulating post-transcriptional processes like intron splicing, stability and mRNA intracellular translocation. As RNA metabolism in neurons attributes a significant proportionality towards neurodegenerative diseases, potential role of RRM domains has become more conspicuous. Here, through flexible protein-protein docking and molecular dynamics simulation approaches, we postulate biochemically auspicious interactions occurring between Butyrylcholinesterase (BuChE) and RRM domains of transactive response DNA binding protein, embryonic lethal abnormal vision like protein 4 and heterogeneous nuclear ribonucleoprotein A1 protein. Through a comprehensive analysis of these interactions, we observed an exclusive binding behavior for RRM domains. Evidently, upon binding to RRMs, omega loop of BuChE attains a closed conformation and masks the access of substrate to the catalytic triad and oxyanion hole. Moreover, prominent adjustments were detected in RRM-bound BuChE at the peripheral anionic, choline binding and proton transfer sites. We propose that interaction of RNA-binding proteins to BuChE may decrease its ability to hydrolyze multiple choline esters, which may contribute in delayed progression of poor cognition during neurodegenerative disorders.

Identification of Ser/Thr kinase and forkhead associated domains in Mycobacterium ulcerans: characterization of novel association between protein kinase Q and MupFHA.

Background Mycobacterium ulcerans, the causative agent of Buruli ulcer in humans, is unique among the members of Mycobacterium genus due to the presence of the virulence determinant megaplasmid pMUM001. This plasmid encodes multiple virulence-associated genes, including mup011, which is an uncharacterized Ser/Thr protein kinase (STPK) PknQ.Methodology/Principal FindingsIn this study, we have characterized PknQ and explored its interaction with MupFHA (Mup018c), a FHA domain containing protein also encoded by pMUM001. MupFHA was found to interact with PknQ and suppress its autophosphorylation. Subsequent protein-protein docking and molecular dynamic simulation analyses showed that this interaction involves the FHA domain of MupFHA and PknQ activation loop residues Ser170 and Thr174. FHA domains are known to recognize phosphothreonine residues, and therefore, MupFHA may be acting as one of the few unusual FHA-domain having overlapping specificity. Additionally, we elucidated the PknQ-dependent regulation of MupDivIVA (Mup012c), which is a DivIVA domain containing protein encoded by pMUM001. MupDivIVA interacts with MupFHA and this interaction may also involve phospho-threonine/serine residues of MupDivIVA.Conclusions/SignificanceTogether, these results describe novel signaling mechanisms in M. ulcerans and show a three-way regulation of PknQ, MupFHA, and MupDivIVA. FHA domains have been considered to be only pThr specific and our results indicate a novel mechanism of pSer as well as pThr interaction exhibited by MupFHA. These results signify the need of further re-evaluating the FHA domain –pThr/pSer interaction model. MupFHA may serve as the ideal candidate for structural studies on this unique class of modular enzymes.

Molecular Interaction Between PPTI and Shakv1.1 Potassium Channel Explored by Docking and Molecular Dynamics Simulation

PPTI, proposed as a potential dendrotoxin, which isolated from the venom of Pseudocerastes persicus has 58 amino acid residues and cross-linked by three disulfide bonds. We explored the interaction between PPTI and ShaKv1.1 potassium channel by successive application of protein-protein docking and molecular dynamic simulation. The analysis of the complexes revealed that the Lys5 residue of PPTI plugged its side chain into the channel selectivity filter. Molecular dynamic simulation in membrane environment of the complex structure resulted in docking stage confirmed the stability of the complex.

Structural Analysis of Respirasomes in Electron Transfer Pathway of Acidithiobacillus ferrooxidans: A Computer-Aided Molecular Designing Study

Acidithiobacillus ferrooxidans obtains its metabolic energy by reducing extracellular ferrous iron with either downhill or uphill electron transfer pathway. The downhill electron transfer pathway has been substantially explored in recent years to underpin the mechanism of iron respiration but, there exists a wide gap in our present understanding on how these proteins are organized as a supercomplex and what sort of atomic level interactions governs their stability in the iron respiratory chain. In the present study, we aimed at unraveling the structural basis of supermolecular association of respirasomes using protein threading, protein-protein docking, and molecular dynamics (MD) simulation protocols. Our results revealed that Phe312 of outer membrane cytochrome c plays a crucial role in diffusing electrons from heme C group to Asp73 of rusticyanin. In line with the previous experimental results, His143 of rusticyanin was found to have a stable interaction with Glu121 of periplasmic cytochrome c4. Cytochrome c4 interacts with subunit B of cytochrome c oxidase through Lys146 and Thr148 of the conserved hydrophobic/aromatic motif 145-WKWTFSY-151 to attain stability during simulation. Phe468 of cytochrome c oxidase was found indispensable for stabilizing heme aa3 during MD simulation. Taken together, we conclude that the molecular interactions of charged and hydrophobic amino acids present on the surface of each respirasome form a hypothetical electron wire in the iron respiratory supercomplex of A. ferrooxidans.

Structural, phylogenetic and docking studies of D-amino acid oxidase activator (DAOA), a candidate schizophrenia gene.

BackgroundSchizophrenia is a neurodegenerative disorder that occurs worldwide and can be difficult to diagnose. It is the foremost neurological disorder leading to suicide among patients in both developed and underdeveloped countries. D-amino acid oxidase activator (DAOA), also known as G72, is directly implicated in the glutamateric hypothesis of schizophrenia. It activates D-amino acid oxidase, which oxidizes D-serine, leading to modulation of the N-methyl-D-aspartate receptor.MethodsMODELLER (9v10) was utilized to generate three dimensional structures of the DAOA candidate gene. The HOPE server was used for mutational analysis. The Molecular Evolutionary Genetics Analysis (MEGA5) tool was utilized to reconstruct the evolutionary history of the candidate gene DAOA. AutoDock was used for protein-ligand docking and Gramm-X and PatchDock for protein-protein docking.ResultsA suitable template (1ZCA) was selected by employing BLASTp on the basis of 33% query coverage, 27% identity and E-value 4.9. The Rampage evaluation tool showed 91.1% favored region, 4.9% allowed region and 4.1% outlier region in DAOA. ERRAT demonstrated that the predicted model had a 50.909% quality factor. Mutational analysis of DAOA revealed significant effects on hydrogen bonding and correct folding of the DAOA protein, which in turn affect protein conformation. Ciona was inferred as the outgroup. Tetrapods were in their appropriate clusters with bifurcations. Human amino acid sequences are conserved, with chimpanzee and gorilla showing more than 80% homology and bootstrap value based on 1000 replications. Molecular docking analysis was employed to elucidate the binding mode of the reported ligand complex for DAOA. The docking experiment demonstrated that DAOA is involved in major amino acid interactions: the residues that interact most strongly with the ligand C28H28N3O5PS2 are polar but uncharged (Gln36, Asn38, Thr 122) and non-polar hydrophobic (Ile119, Ser171, Ser21, Ala31). Protein-protein docking simulation demonstrated two ionic bonds and one hydrogen bond involving DAOA. Lys-7 of the receptor protein interacted with Lys-163 and Asp-2037. Tyr-03 interacted with Arg-286 of the ligand protein and formed a hydrogen bond.ConclusionThe predicted interactions might serve to inhibit the disease-related allele. It is assumed that current bioinformatics methods will contribute significantly to identifying, analyzing and curing schizophrenia. There is an urgent need to develop effective drugs for schizophrenia, and tools for examining candidate genes more accurately and efficiently are required.

Towards Fast and Accurate Calculation of Protein pKa Values Exploiting Various Degrees of Conformational Flexibility

Virtually every biological process is pH dependent. Many fundamental biological phenomena including protein folding, enzyme catalysis, protein-protein interactions and pathological conditions are profoundly influenced by the pH of their environment. Accurate prediction of the pKa values of residues that can adopt variable protonation states would be a significant step towards probing the effects of pH on proteins. Despite the progress in pKa prediction algorithms, developing a method that can incorporate extensive protein conformational flexibility while retaining a relatively small computational resource footprint remains a significant challenge. We developed a fast and accurate method to predict pKas of residues that commonly exhibit variable protonation states in proteins. The algorithm utilizes Metropolis Monte Carlo strategy to incorporate conformational flexibility through extensive amino-acid side-chain sampling and is based on Rosetta score function with an explicit term to account for protonation probabilities of individual amino-acids. We tested the technique by predicting pKas for an assembled dataset comprising 306 residues from 44 proteins resulting in a root-mean-square deviation (RMSD) of 0.81 from experimental values. We analyzed the effects of employing increasing levels of conformational flexibility by (1) sampling the side-chains of neighboring residues and (2) using a generated ensemble of 50 diverse backbone conformers resulting in 77% and 79% predictions < 1 pH unit from experimental values respectively. Using additional degrees of freedom allowed capture of vital hydrogen-bonding and charge-charge interactions, but resulted in structural rearrangements negating pKa shifts in some cases. The method yielded good results when used to predict large pKa shifts in point mutants from Staphylococcal nuclease. Finally, we employed our method to dynamically alter the ionization states of residues during protein-protein docking simulations to seek improvements in accuracy of prediction of complexes.

How Does a Registry Change in Dynein’s Coiled-Coil Stalk Drive Binding of Dynein to Microtubules?

Dynein is a motor protein that transports cellular cargo along the microtubule (MT) by consuming ATP. Dynein's microtubule-binding domain (MTBD) is separated from the ATP-binding core by a ~15 nm stalk that consists of two α-helices forming an antiparallel coiled coil. It was previously suggested that the coiled-coil stalk creates a registry shift to modulate its binding affinity for MT. A crystal structure of the low-affinity form of MTBD was determined, but that of the high-affinity form with the registry shift is not yet available. In this study, we obtained an all-atom model structure for the high-affinity form of MTBD bound to MT by an anisotropic network model, protein-protein docking, and molecular dynamics simulations. We observe that the magnitude of the coiled-coil helix sliding is dramatically reduced near the two prolines that form the stalk-MTBD boundary and subsequently transformed to cyclic movements of MTBD helices, leading to formation of a new salt bridge with MT at the binding interface. The proposed mechanism explains the roles of highly conserved residues such as the two prolines at the stalk-MTBD boundary, the nonpolar tryptophan and proline residues near the binding interface, and the electropositive residues forming salt bridges with MT.

Structural and Enzymatic Insights into the ATP Binding and Autophosphorylation Mechanism of a Sensor Histidine Kinase

DesK is a sensor histidine kinase (HK) that allows Bacillus subtilis to respond to cold shock, triggering the adaptation of membrane fluidity via transcriptional control of a fatty acid desaturase. It belongs to the HK family HPK7, which includes the nitrogen metabolism regulators NarX/Q and the antibiotic sensor LiaS among other important sensor kinases. Structural information on different HK families is still scarce and several questions remain, particularly concerning the molecular features that determine HK specificity during its catalytic autophosphorylation and subsequent response-regulator phosphotransfer reactions. To analyze the ATP-binding features of HPK7 HKs and dissect their mechanism of autophosphorylation at the molecular level, we have studied DesK in complex with ATP using high resolution structural approaches in combination with biochemical studies. We report the first crystal structure of an HK in complex with its natural nucleotidic substrate. The general fold of the ATP-binding domain of DesK is conserved, compared with well studied members of other families. Yet, DesK displays a far more compact structure at the ATP-binding pocket: the ATP lid loop is much shorter with no secondary structural organization and becomes ordered upon ATP loading. Sequence conservation mapping onto the molecular surface, semi-flexible protein-protein docking simulations, and structure-based point mutagenesis allow us to propose a specific domain-domain geometry during autophosphorylation catalysis. Supporting our hypotheses, we have been able to trap an autophosphorylating intermediate state, by protein engineering at the predicted domain-domain interaction surface.

Molecular mechanism of δ-dendrotoxin-potassium channel recognition explored by docking and molecular dynamic simulations

δ-Dendrotoxin, isolated from mamba snake venom, has 57 residues cross-linked by three disulfide bridges. The protein shares a pharmacological activity with other animal toxins, the potent blockade of potassium channels, but is structurally unrelated to toxins of different species. We employed alanine-scanning mutagenesis to explore the molecular mechanism of δ-dendrotoxin binding to potassium channels, using protein-protein docking and molecular dynamic simulations. In our reasonable model of the δ-dendrotoxin-ShaKv1.1 complex, δ-dendrotoxin interacted mainly with the N-terminal region and the turn of two antiparallel β-sheets of the channel. This binding mode could well explain the functional roles of critical residues in δ-dendrotoxin and the ShaKv1.1 channel. Structural analysis indicated that the critical Lys6 residue of δ-dendrotoxin plugged its side chain into a channel selectivity filter. Another two critical δ-dendrotoxin residues, Lys3 and Arg10, were found to contact channel residues through strong polar and nonpolar interactions, especially salt-bridge interactions. As for the ShaKv1.1 channel, the channel turrets were found in the "half-open state," and two of four Glu423 in the turrets of the channel B and D chains could interact, respectively, with Lys3 and Lys26 of δ-dendrotoxin through electrostatic interactions. The essential Asp431 channel residue was found to associate electrostatically with Arg10 of δ-dendrotoxin, and a critical Tyr449 channel residue was just under the channel-interacting surface of δ-dendrotoxin. Together, these novel data may accelerate the structure-function research of toxins in the dendrotoxin family and be of significant value in revealing the diverse interactions between animal toxins and potassium channels.

An actin-filament-binding interface on the Arp2/3 complex is critical for nucleation and branch stability

The Arp2/3 complex polymerizes new actin filaments from the sides of existing filaments, forming Y-branched networks that are critical for actin-mediated force generation. Binding of the Arp2/3 complex to the sides of actin filaments is therefore central to its actin-nucleating and branching activities. Although a model of the Arp2/3 complex in filament branches has been proposed based on electron microscopy, this model has not been validated using independent approaches, and the functional importance of predicted actin-binding residues has not been extensively tested. Using a combination of molecular dynamics and protein-protein docking simulations, we derived an independent structural model of the interaction between two subunits of the Arp2/3 complex that are key to actin binding, ARPC2 and ARPC4, and the side of an actin filament. This model agreed remarkably well with the previous results from electron microscopy. Complementary mutagenesis experiments revealed numerous residues in ARPC2 and ARPC4 that were required for the biochemical activity of the entire complex. Functionally critical residues clustered together and defined a surface that was predicted by protein-protein docking to be buried in the interaction with actin. Moreover, key residues at this interface were crucial for actin nucleation and Y-branching, high-affinity F-actin binding, and Y-branch stability, demonstrating that the affinity of Arp2/3 complex for F actin independently modulates branch formation and stability. Our results highlight the utility of combining computational and experimental approaches to study protein-protein interactions and provide a basis for further elucidating the role of F-actin binding in Arp2/3 complex activation and function.

Molecular Interaction of Cofilin with F-Actin and Implications for Filament Severing

The regulation of actin polymerization within the cell is critical for many cell functions. Cofilin plays an important part in this process since it binds and severs actin filaments leading to depolymerization as well as the creation of new barbed ends. Although the details of cofilin's interaction with G-actin have been elucidated through a range of experimental studies, the specific interactions with F-actin have remained more elusive. Here we present the results of a detailed computational study involving a combination of protein-protein docking and molecular dynamics simulations. The resulting structural model for the cofilin/F-actin complex matches very well with existing cryoEM and mutagenesis data. Further insight is gained from a sequence and secondary structure alignment of cofilin analogs that strongly supports our binding model. Based on our binding models, we have performed simulations of sparsely and fully decorated filaments in order to gain insight into the mechanism of F-actin twisting and severing. The implications of this binding model in the function and severing action of cofilin are discussed.

Holo‐BtuF stabilizes the open conformation of the vitamin B12 ABC transporter BtuCD

While there is evidence that other ABC transporters can tell between empty and loaded substrate binding protein, reconstitution experiments suggest otherwise for the Escherichia coli vitamin B12 importer BtuCD-F. Here, we address the question of BtuCD-F substrate sensitivity in a combined protein-protein docking and molecular dynamics simulation approach. Starting from the BtuCD and holo-BtuF crystal structures, we model two holo-BtuCD-F docking complexes differing by a 180 degrees orientation of BtuF. One of these is similar to the apo-BtuCD-F crystal structure. Both docking complexes were embedded in a lipid/water environment to investigate their dynamics and BtuCD's conformational response to the presence and absence of BtuF, vitamin B12, and Mg-ATP in a series of 28 independent MD simulations. We find holo-BtuF stabilizing the open conformation of BtuCD, whereas the transporter begins to close again when BtuF or vitamin B12 is removed-suggesting BtuCD-F is capable of substrate sensitivity. We identified BtuC transmembrane helices 3 and 5, the L-loops and the adjacent helices comprised of BtuC residues 170-180 as hotspots of conformational change. We propose the latter to act as substrate sensors. BtuF-Trp44 appears to act as a lid on the vitamin B12 binding cleft in BtuF X-ray structures and protrudes into the BtuCD transport channel in one of our simulations, which might represent an initial step in vitamin B12 uptake. On an average, we observe subunit motions where the nucleotide binding domains approach each other while the transmembrane domains display an opening trend toward the periplasm.

Binding Mode Analysis of Bacillus subtilis Obg with Ribosomal Protein L13 through Computational Docking Study

Introduction: GTPases known as translation factor play a vital role as ribosomal subunit assembly chaperone. The bacterial Obg proteins (<TEX>$Spo{\underline{0B}}$</TEX>-associated <TEX>${\underline{G}}TP$</TEX>-binding protein) belong to the subfamily of P-loop GTPase proteins and now it is considered as one of the new target for antibacterial drug. The majority of bacterial Obgs have been commonly found to be associated with ribosome, implying that these proteins may play a fundamental role in ribosome assembly or maturation. In addition, one of the experimental evidences suggested that Bacillus subtilis Obg (BsObg) protein binds to the L13 ribosomal protein (BsL13) which is known to be one of the early assembly proteins of the 50S ribosomal subunit in Escherichia coli. In order to investigate binding mode between the BsObg and the BsL13, protein-protein docking simulation was carried out after generating 3D structure of the BsL13 structure using homology modeling method. Materials and Methods: Homology model structure of BsL13 was generated using the EcL13 crystal structure as a template. Protein-protein docking of BsObg protein with ribosomal protein BsL13 was performed by DOT, a macro-molecular docking software, in order to predict a reasonable binding mode. The solvated energy minimization calculation of the docked conformation was carried out to refine the structure. Results and Discussion: The possible binding conformation of BsL13 along with activated Obg fold in BsObg was predicted by computational docking study. The final structure is obtained from the solvated energy minimization. From the analysis, three important H-bond interactions between the Obg fold and the L13 were detected: Obg:Tyr27-L13:Glu32, Obg:Asn76-L13:Glu139, and Obg:Ala136-L13:Glu142. The interaction between the BsObg and BsL13 structures were also analyzed by electrostatic potential calculations to examine the interface surfaces. From the results, the key residues for hydrogen bonding and hydrophobic interaction between the two proteins were predicted. Conclusion and Prospects: In this study, we have focused on the binding mode of the BsObg protein with the ribosomal BsL13 protein. The interaction between the activated Obg and target protein was investigated with protein-protein docking calculations. The binding pattern can be further used as a base for structure-based drug design to find a novel antibacterial drug.

Development of a scoring method for predicting protein complex structures

The information about protein-protein interactions increases much more rapidly than the increase of the number of the tertiary structures of those protein complexes. Therefore, precise prediction of protein complex structures by protein-protein docking simulations is required. When the protein complex is rebuilt from its component protomers which derive from experimentally determined complex structure (native structure) by docking, the complex models with rmsd < 10 A from the native structure (near-native model) could be obtained, along with a great number of false positives (decoy). The separation of near-native models from many decoys is therefore needed in the prediction of complex structures by docking. In this study, we developed the method for scoring docking models so that the near-native models were higher in rank than decoys, based on the assumption that the interfaces of near-native models are more complementary in terms of surface properties and shapes compared to those of decoys.