Conformational ensemble of an intrinsically flexible loop in mitochondrial import protein Tim21 studied by modeling and molecular dynamics simulations

BackgroundTim21, a subunit of a highly dynamic translocase of the inner mitochondrial membrane (TIM23) complex, translocates proteins by interacting with subunits in the translocase of the outer membrane (TOM) complex and Tim23 channel in the TIM23 complex. A loop segment in Tim21, which is in close proximity of the binding site of Tim23, has different conformations in X-ray, NMR and new crystal contact-free space (CCFS) structures. MD simulations can provide information on the structure and dynamics of the loop in solution.MethodsThe conformational ensemble of the loop was characterized using loop modeling and molecular dynamics (MD) simulations.ResultsMD simulations confirmed mobility of the loop. Multidimensional scaling and clustering were used to characterize the dynamic conformational ensemble of the loop. Free energy landscape showed that the CCFS crystal structure occupied a low energy region as compared to the conventional X-ray crystal structure. Analysis of crystal packing indicates that the CCFS provides larger conformational space for the motions of the loop.ConclusionsOur work reported the conformational ensemble of the loop in solution, which is in agreement with the structure obtained from CCFS approach. The combination of the experimental techniques and computational methods is beneficial for studying highly flexible regions of proteins.General significanceComputational methods, such as loop modeling and MD simulations, have proved to be useful for studying conformational flexibility of proteins. These methods in integration with experimental techniques such as CCFS has the potential to transform the studies on flexible regions of proteins.     

Comparison of full-atomic and coarse-grained models to examine the molecular fluctuations of c-AMP dependent protein kinase

Molecular fluctuations of the native conformation of c-AMP dependent protein kinase (cAPK) have been investigated with three different approaches. The first approach is the full atomic normal mode analysis (NMA) with empirical force fields. The second and third approaches are based on a coarse-grained model with a single single-parameter- harmonic potential between close residues in the crystal structure of the molecule without any residue specificity. The second method calculates only the magnitude of fluctuations whereas the third method is developed to find the directionality of the fluctuations which are essential to understand the functional importance of biological molecules. The aim, in this study, is to determine whether using such coarse-grained models are appropriate for elucidating the global dynamic characteristics of large proteins which reduces the size of the system at least by a factor of ten. The mean-square fluctuations of C(alpha) atoms and the residue cross-correlations are obtained by three approaches. These results are then compared to test the results of coarse grained models on the overall collective motions. All three of the approaches show that highly flexible regions correspond to the activation and solvent exposed loops, whereas the conserved residues (especially in substrate binding regions) exhibit almost no flexibility, adding stability to the structure. The anti-correlated motions of the two lobes of the catalytic core provide flexibility to the molecule. High similarities among the results of these methods indicate that the slowest modes governing the most global motions are preserved in the coarse grained models for proteins. This finding may suggest that the general shapes of the structures are representative of their dynamic characteristics and the dominant motions of protein structures are robust at coarse-grained levels.     

Mechanics of proteins with a focus on atomic force microscopy

The capacity of proteins to function relies on a balance between molecular stability to maintain their folded state and structural flexibility allowing conformational changes related to biological function. Among many others, four different examples can be chosen. The giant protein titin is stretched and can unfold during muscle contraction providing passive elasticity to muscle tissue; myoglobin adsorbs and releases oxygen molecules thank to conformational changes in its structure; the outer membrane protein G (OmpG) is a bacterial porin with a long and flexible loop that modulates gating; and the proton pump bacteriorhodopsin adapts its cytosolic half to allow proton pumping. All these conformational changes triggered either by chemical or by physical cues, require mechanical flexibility or elasticity of certain protein domains. While the methods to determine protein structure, X-ray crystallography above all, have been dramatically improved over the last decades, the number of tools that directly measure the mechanical flexibility of proteins and protein domains is still limited. In this tutorial, after a brief introduction to protein structure, we present some of the available techniques to estimate protein flexibility, then focusing on atomic force microscopy (AFM). We describe the principles of the technique and its various imaging and force spectroscopy modes of operation that allow probing the elasticity of proteins, protein domains and their surrounding environment.     

Escherichia coli aspartate carbamoyltransferase: the probing of crystal structure analysis via site-specific mutagenesis

Crystal structures are known for aspartate carbamoyltransferase (ATCase) in the T and R states, with and without the allosteric activator adenosine triphosphate (ATP) or inhibitor cytidine triphosphate (CTP). Visual inspection of X-ray crystal structures does not provide all of the information necessary for the determination of structure--function relationships in protein molecules. This problem is compounded because the crystalline states of the molecule may introduce effects due to crystal packing, restricted flexibility and less than optimum enzymatic conditions. Therefore, alternative techniques are required to test mechanisms conjectured from three-dimensional crystal structures of proteins. The technique of site-specific mutagenesis allows the researcher to test structure--function models based on three-dimensional structures and to obtain further insight into characteristics of the enzyme. Site-specific mutagenesis has been used to probe residues believed to be critical in the structure and function of ATCase. Selection of residues to be mutated has depended extensively on three-dimensional crystal structures of the enzyme. To date, 48 site-specific mutations at 37 different amino acid sites have been published. Although a total of 118 mutants at 58 different sites has been communicated to our laboratory, only published mutants will be considered in this review. In this paper, we compile for the first time, review, and analyze the site-specific mutants of ATCase. Site-specific mutagenesis of proteins has become a powerful technique in modern-day molecular biology, especially in studying a molecule as large as aspartate carbamoyltransferase. In this review, the role of site-specific mutagenesis of ATCase is discussed and improvements in the analysis are suggested.     

Conformational dynamics of human protein kinase CK2α and its effect on function and inhibition

Protein kinase, casein kinase II (CK2), is ubiquitously expressed and highly conserved protein kinase that shows constitutive activity. It phosphorylates a diverse set of proteins and plays crucial role in several cellular processes. The catalytic subunit of this enzyme (CK2α) shows remarkable flexibility as evidenced in numerous crystal structures determined till now. Here, using analysis of multiple crystal structures and long timescale molecular dynamics simulations, we explore the conformational flexibility of CK2α. The enzyme shows considerably higher flexibility in the solution as compared to that observed in crystal structure ensemble. Multiple conformations of hinge region, located near the active site, were observed during the dynamics. We further observed that among these multiple conformations, the most populated conformational state was inadequately represented in the crystal structure ensemble. The catalytic spine, was found to be less dismantled in this state as compared to the “open” hinge/αD state crystal structures. The comparison of dynamics in unbound (Apo) state and inhibitor (CX4945) bound state exhibits inhibitor induced suppression in the overall dynamics of the enzyme. This is especially true for functionally important glycine‐rich loop above the active site. Together, this work gives novel insights into the dynamics of CK2α in solution and relates it to the function. This work also explains the effect of inhibitor on the dynamics of CK2α and paves way for development of better inhibitors.     

Observation of slow dynamic exchange processes in Ras protein crystals by 31P solid state NMR spectroscopy

The folding, structure and biological function of many proteins are inherently dynamic properties of the protein molecule. Often, the respective molecular processes are preserved upon protein crystallization, leading, in X-ray diffraction experiments, to a blurring of the electron density map and reducing the resolution of the derived structure. Nuclear magnetic resonance (NMR) is known to be an alternative method to study molecular structure and dynamics. We designed and built a probe for phosphorus solid state NMR that allows for the first time to study static properties as well as dynamic processes in single-crystals of a protein by NMR spectroscopy. The sensitivity achieved is sufficient to detect the NMR signal from individual phosphorus sites in a 0.3mm(3) size single-crystal of GTPase Ras bound to the nucleotide GppNHp, that is, the signal from approximately 10(15) phosphorus nuclei. The NMR spectra obtained are discussed in terms of the conformational variability of the active center of the Ras-nucleotide complex. We conclude that, in the crystal, the protein complex exists in three different conformations. Magic angle spinning (MAS) NMR spectra of a powder sample of Ras-GppNHp show a splitting of one of the phosphate resonances and thus confirm this conclusion. The MAS spectra provide, furthermore, evidence of a slow, temperature-dependent dynamic exchange process in the Ras protein crystal.     

Comparison of Full-atomic and Coarse-grained Models to Examine the Molecular Fluctuations of c-AMP Dependent Protein Kinase

Abstract

Molecular fluctuations of the native conformation of c-AMP dependent protein kinase (cAPK) have been investigated with three different approaches. The first approach is the full atomic normal mode analysis (NMA) with empirical force fields. The second and third approaches are based on a coarse-grained model with a single single-parameter- harmonic potential between close residues in the crystal structure of the molecule without any residue specificity. The second method calculates only the magnitude of fluctuations whereas the third method is developed to find the directionality of the fluctuations which are essential to understand the functional importance of biological molecules. The aim, in this study, is to determine whether using such coarse-grained models are appropriate for elucidating the global dynamic characteristics of large proteins which reduces the size of the system at least by a factor of ten. The mean-square fluctuations of C^α atoms and the residue cross-correlations are obtained by three approaches. These results are then compared to test the results of coarse grained models on the overall collective motions. AH three of the approaches show that highly flexible regions correspond to the activation and solvent exposed loops, whereas the conserved residues (especially in substrate binding regions) exhibit almost no flexibility, adding stability to the structure. The anti-correlated motions of the two lobes of the catalytic core provide flexibility to the molecule. High similarities among the results of these methods indicate that the slowest modes governing the most global motions are preserved in the coarse grained models for proteins. This finding may suggest that the general shapes of the structures are representative of their dynamic characteristics and the dominant motions of protein structures are robust at coarse-grained levels.     

Stability and rigidity/flexibility—Two sides of the same coin?

Protein molecules require both flexibility and rigidity for functioning. The fast and accurate prediction of protein rigidity/flexibility is one of the important problems in protein science. We have determined flexible regions for four homologous pairs from thermophilic and mesophilic organisms by two methods: the fast FoldUnfold which uses amino acid sequence and the time consuming MDFirst which uses three-dimensional structures. We demonstrate that both methods allow determining flexible regions in protein structure. For three of the four thermophile–mesophile pairs of proteins, FoldUnfold predicts practically the same flexible regions which have been found by the MD/First method. As expected, molecular dynamics simulations show that thermophilic proteins are more rigid in comparison to their mesophilic homologues. Analysis of rigid clusters and their decomposition provides new insights into protein stability. It has been found that the local networks of salt bridges and hydrogen bonds in thermophiles render their structure more stable with respect to fluctuations of individual contacts. Such network includes salt bridge triads Agr-Glu-Lys and Arg-Glu-Arg, or salt bridges (such as Arg-Glu) connected with hydrogen bonds. This ionic network connects alpha helices and rigidifies the structure. Mesophiles can be characterized by stand alone salt bridges and hydrogen bonds or small ionic clusters. Such difference in the network of salt bridges results in different flexibility of homologous proteins. Combining both approaches allows characterizing structural features in atomic detail that determine the rigidity/flexibility of a protein structure. This article is a part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

Effects of environment on the structure of Pyrococcus furiosus rubredoxin: a molecular dynamics study

Molecular dynamics simulations based on a 0.95-A resolution crystal structure of Pyrococcus furiosus have been performed to elucidate the effects of the environment on the structure of rubredoxin, and proteins in general. Three 1-ns simulations are reported here: two crystalline state simulations at 123 and 300 K, and a solution state simulation at 300 K. These simulations show that temperature has a greater impact on the protein structure than the close molecular contacts of the crystal matrix in rubredoxin, although both have an effect on its dynamic properties. These results indicate that differences between NMR solution structures and X-ray crystal structures will be relatively minor if they are done at similar temperatures. In addition, the crystal simulations appears to mimic previous crystallographic experiments on the effects of cryo-temperature on temperature factors, and might provide a useful tool in the structural analysis of protein structures solved at cryo-temperatures.     

Domain-Based Protein Docking with Extremely Large Conformational Changes

Proteins are key components in many processes in living cells, and physical interactions with other proteins and nucleic acids often form key parts of their functions. In many cases, large flexibility of proteins as they interact is key to their function. To understand the mechanisms of these processes, it is necessary to consider the 3D structures of such protein complexes. When such structures are not yet experimentally determined, protein docking has long been present to computationally generate useful structure models. However, protein docking has long had the limitation that the consideration of flexibility is usually limited to very small movements or very small structures. Methods have been developed which handle minor flexibility via normal mode or other structure sampling, but new methods are required to model ordered proteins which undergo large-scale conformational changes to elucidate their function at the molecular level. Here, we present Flex-LZerD, a framework for docking such complexes. Via partial assembly multidomain docking and an iterative normal mode analysis admitting curvilinear motions, we demonstrate the ability to model the assembly of a variety of protein–protein and protein-nucleic acid complexes.     

Tropomyosin crystal structure and muscle regulation

The crystal structure of tropomyosin filaments has been solved to 15 A resolution by refinement of models against the diffraction data and heavy atom labeling of cysteine residues. These results confirm and extend earlier findings. The improved maps reveal the pitch of the coiled coil, the location of the cysteine residues, and the location and features of the overlapping molecular ends in the filaments. A correlation can now be made between regions of the amino acid sequence and key features of the molecule, such as contact sites in the lattice and departures from regularity along the coiled coil. The crystal shows remarkable dynamic features and the relative flexibility of different parts of the molecule as well as its anisotropic character have been determined. The structure and motions of tropomyosin in the crystal provide information on the structure of tropomyosin in muscle and its possible role in regulation. An atomic model of the molecule has been constructed, based on the low resolution X-ray results, together with the stereochemistry of alpha-helical coiled coils. In contrast to previous views, the molecule appears to display but one set of seven alpha-sites that permit weak linkages of the flexible tropomyosin filament to the actin helix. Correspondingly, we picture that in the "off" state of ATPase activity, the alpha-sites are not occupied; in the "on" state, they are only partly occupied; and in the "potentiated" state, they are more completely saturated. Control of contraction is therefore seen as a statistical mechanism requiring at least three distinct average conformations for the tropomyosin molecule on the actin helix.     

Conformational Transitions upon Ligand Binding: Holo-Structure Prediction from Apo Conformations

Biological function of proteins is frequently associated with the formation of complexes with small-molecule ligands. Experimental structure determination of such complexes at atomic resolution, however, can be time-consuming and costly. Computational methods for structure prediction of protein/ligand complexes, particularly docking, are as yet restricted by their limited consideration of receptor flexibility, rendering them not applicable for predicting protein/ligand complexes if large conformational changes of the receptor upon ligand binding are involved. Accurate receptor models in the ligand-bound state (holo structures), however, are a prerequisite for successful structure-based drug design. Hence, if only an unbound (apo) structure is available distinct from the ligand-bound conformation, structure-based drug design is severely limited. We present a method to predict the structure of protein/ligand complexes based solely on the apo structure, the ligand and the radius of gyration of the holo structure. The method is applied to ten cases in which proteins undergo structural rearrangements of up to 7.1 Å backbone RMSD upon ligand binding. In all cases, receptor models within 1.6 Å backbone RMSD to the target were predicted and close-to-native ligand binding poses were obtained for 8 of 10 cases in the top-ranked complex models. A protocol is presented that is expected to enable structure modeling of protein/ligand complexes and structure-based drug design for cases where crystal structures of ligand-bound conformations are not available.     

Recent developments in structural proteomics for protein structure determination

The major challenges in structural proteomics include identifying all the proteins on the genome-wide scale, determining their structure-function relationships, and outlining the precise three-dimensional structures of the proteins. Protein structures are typically determined by experimental approaches such as X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. However, the knowledge of three-dimensional space by these techniques is still limited. Thus, computational methods such as comparative and de novo approaches and molecular dynamic simulations are intensively used as alternative tools to predict the three-dimensional structures and dynamic behavior of proteins. This review summarizes recent developments in structural proteomics for protein structure determination; including instrumental methods such as X-ray crystallography and NMR spectroscopy, and computational methods such as comparative and de novo structure prediction and molecular dynamics simulations.     

The Differential Response of Proteins to Macromolecular Crowding

The habitat in which proteins exert their function contains up to 400 g/L of macromolecules, most of which are proteins. The repercussions of this dense environment on protein behavior are often overlooked or addressed using synthetic agents such as poly(ethylene glycol), whose ability to mimic protein crowders has not been demonstrated. Here we performed a comprehensive atomistic molecular dynamic analysis of the effect of protein crowders on the structure and dynamics of three proteins, namely an intrinsically disordered protein (ACTR), a molten globule conformation (NCBD), and a one-fold structure (IRF-3) protein. We found that crowding does not stabilize the native compact structure, and, in fact, often prevents structural collapse. Poly(ethylene glycol) PEG500 failed to reproduce many aspects of the physiologically-relevant protein crowders, thus indicating its unsuitability to mimic the cell interior. Instead, the impact of protein crowding on the structure and dynamics of a protein depends on its degree of disorder and results from two competing effects: the excluded volume, which favors compact states, and quinary interactions, which favor extended conformers. Such a viscous environment slows down protein flexibility and restricts the conformational landscape, often biasing it towards bioactive conformations but hindering biologically relevant protein-protein contacts. Overall, the protein crowders used here act as unspecific chaperons that modulate the protein conformational space, thus having relevant consequences for disordered proteins.     

Biased coarse-grained molecular dynamics simulation approach for flexible fitting of X-ray structure into cryo electron microscopy maps

Several approaches have been introduced to interpret, in terms of high-resolution structure, low-resolution structural data as obtained from cryo-EM. As conformational changes are often observed in biological molecules, these techniques need to take into account the flexibility of proteins. Flexibility has been described in terms of movement between rigid domains and between rigid secondary structure elements, which present some limitations for studying dynamical properties. Normal mode analysis has also been used, but is limited to medium resolution data. All-atom molecular dynamics fitting techniques are more appropriate to fit structures into higher-resolution data as full protein flexibility is considered, but are cumbersome in terms of computational time. Here, we introduce a coarse-grained approach; a Go-model was used to represent biological molecules, combined with biased molecular dynamics to reproduce accurately conformational transitions. Illustrative examples on simulated data are shown. Accurate fittings can be obtained for resolution ranging from 5 to 20 Å. The approach was also tested on experimental data of Elongation Factor G and Escherichia coli RNA polymerase, where its validity is compared to previous models obtained from different techniques. This comparison demonstrates that quantitative flexible techniques, as opposed to manual docking, need to be considered to interpret low-resolution data.     

Determination of protein oligomeric structure from small‐angle X‐ray scattering

Small‐angle X‐ray scattering (SAXS) is useful for determining the oligomeric states and quaternary structures of proteins in solution. The average molecular mass in solution can be calculated directly from a single SAXS curve collected on an arbitrary scale from a sample of unknown protein concentration without the need for beamline calibration or protein standards. The quaternary structure in solution can be deduced by comparing the experimental SAXS curve to theoretical curves calculated from proposed models of the oligomer. This approach is especially robust when the crystal structure of the target protein is known, and the candidate oligomer models are derived from the crystal lattice. When SAXS data are obtained at multiple protein concentrations, this analysis can provide insight into dynamic self‐association equilibria. Herein, we summarize the computational methods that are used to determine protein molecular mass and quaternary structure from SAXS data. These methods are organized into a workflow and demonstrated with four case studies using experimental SAXS data from the published literature.     

Towards a complete description of the structural and dynamic properties of the denatured state of barnase and the role of residual structure in folding

The detailed characterization of denatured proteins remains elusive due to their mobility and conformational heterogeneity. NMR studies are beginning to provide clues regarding residual structure in the denatured state but the resulting data are too sparse to be transformed into molecular models using conventional techniques. Molecular dynamics simulations can complement NMR by providing detailed structural information for components of the denatured ensemble. Here, we describe three independent 4 ns high-temperature molecular dynamics simulations of barnase in water. The simulated denatured state was conformationally heterogeneous with respect to the conformations populated both within a single simulation and between simulations. Nonetheless, there were some persistent interactions that occurred to varying degrees in all simulations and primarily involved the formation of fluid hydrophobic clusters with participating residues changing over time. The region of the beta(3-4) hairpin contained a particularly high degree of such side-chain interactions but it lacked beta-structure in two of the three denatured ensembles: beta(3-4) was the only portion of the beta-structure to contain significant residual structure in the denatured state. The two principal alpha-helices (alpha1 and alpha2) adopted dynamic helical structure. In addition, there were persistent contacts that pinched off core 2 from the body of the protein. The rest of the protein was unstructured, aside from transient and mostly local side-chain interactions. Overall, the simulated denatured state contains residual structure in the form of dynamic, fluctuating secondary structure in alpha1 and alpha2, as well as fluctuating tertiary contacts in the beta(3-4) region, and between alpha1 and beta(3-4), in agreement with previous NMR studies. Here, we also show that these regions containing residual structure display impaired mobility by both molecular dynamics and NMR relaxation experiments. The residual structure was important in decreasing the conformational states available to the chain and in repairing disrupted regions. For example, tertiary contacts between beta(3-4) and alpha1 assisted in the refolding of alpha1. This contact-assisted helix formation was confirmed in fragment simulations of beta(3-4) and alpha1 alone and complexed, and, as such, alpha1 and beta(3-4) appear to be folding initiation sites. The role of these sites in folding was investigated by working backwards and considering the simulation in reverse, noting that earlier time-points from the simulations provide models of the major intermediate and transition states in quantitative agreement with data from both unfolding and refolding experiments. Both beta(3-4) and alpha1 are dynamic in the denatured state but when they collide and make enough contacts, they provide a loose structural scaffold onto which further beta-strands pack. The beta-structure condenses about beta(3-4), while alpha1 aids in stabilizing beta(3-4) and maintaining its orientation. The resulting beta-structure is relatively planar and loose in the major intermediate. Further packing ensues, and as a result the beta-sheet twists, leading to the major transition state. The structure is still expanded and loops are not well formed at this point. Fine-tuning of the packing interactions and the final condensation of the structure then occurs to yield the native state.     

Molecular dynamics simulations of bovine rhodopsin: influence of protonation states and different membrane-mimicking environments

G-protein coupled receptors (GPCRs) are a protein family of outstanding pharmaceutical interest. GPCR homology models, based on the crystal structure of bovine rhodopsin, have been shown to be valuable tools in the drug-design process. The initial model is often refined by molecular dynamics (MD) simulations, a procedure that has been recently discussed controversially. We therefore analyzed MD simulations of bovine rhodopsin in order to identify contacts that could serve as constraints in the simulation of homology models. Additionally, the effect of an N-terminal truncation, the nature of the membrane mimic, the influence of varying protonation states of buried residues and the importance of internal water molecules was analyzed. All simulations were carried out using the program-package GROMACS. While N-terminal truncation negatively influenced the overall protein stability, a stable simulation was possible in both solvent environments. As regards the protonation state of titratable sites, the experimental data could be reproduced by the program UHBD (University of Houston Brownian Dynamics), suggesting its application for studying homology models of GPCRs. A high flexibility was observed for internal water molecules at some sites. Finally, interhelical hydrogen-bonding interactions could be derived, which can now serve as constraints in the simulations of GPCR homology models.     

Molecular dynamics study of a calmodulin-like protein with an IQ peptide: spontaneous refolding of the protein around the peptide

The Calmodulin (CaM) is a small (16.7 kDa), highly acidic protein that is crucial to all eukaryotes by serving as a prototypical calcium sensor. In the present study, we investigated, through molecular dynamics simulations, the dynamics of a complex between the Mlc1p protein, which is a CaM-like protein, and the IQ4 peptide. This protein-peptide interaction is of high importance because IQ motifs are widely distributed among different kinds of CaM-binding proteins. The Mlc1p-IQ4 complex, which had been resolved by crystallography to 2.1 A, confers to a Ca(+2)-independent stable structure. During the simulations, the complex undergoes a complicated modulation process, which involves bending of the angles between the alpha-helices of the protein, breaking of the alpha-helical structure of the IQ4 peptide into two sections, and formation of new contact points between the protein and the peptide. The dynamics of the process consist of fast sub picosecond events and much slower ones that take a few nanoseconds to completion. Our study expands the information embedded in the crystal structure of the Mlc1p-IQ4 complex by describing its dynamic behavior as it evolves from the crystal structure to a form stable in solution. The article shows that careful application of molecular dynamics simulations can be used for extending the structural information presented by the crystal structure, thereby revealing the dynamic configuration of the protein in its physiological environment.     

Membrane simulations of OpcA: gating in the loops?

Mobility of extracellular loops may play an important role in the function of outer membrane proteins from Gram-negative bacteria. Molecular dynamics simulations of OpcA from Neisseria meningitidis, embedded in a lipid bilayer, have been used to explore the relationship between the crystal structure and dynamic function of this protein. The results suggest that the crystal environment may constrain the membrane protein structure in a nonphysiological state. The presence of lipids and physiological salt concentrations result in changes in the conformation of the extracellular loops of OpcA, leading to opening of a pore, and to modulation of the molecular surface implicated in recognition of proteoglycan. These changes may be related to the role of OpcA in pathogenesis via modulation of the conformation of a possible sialic acid binding site.