Steered molecular dynamics simulations of a type IV pilus probe initial stages of a force-induced conformational transition

Structure and Dynamics of Type 4a Pili and Type 2 Secretion System Endopili

Within the highly diverse type four filament (TFF or T4F) superfamily, the machineries of type IVa pili (T4aP) and the type 2 secretion system (T2SS) in diderm bacteria exhibit a substantial sequence similarity despite divergent functions and distinct appearances: T4aP can extend micrometers beyond the outer membrane, whereas the endopili in the T2SS are restricted to the periplasm. The determination of the structure of individual components and entire filaments is crucial to understand how their structure enables them to serve different functions. However, the dynamics of these filaments poses a challenge for their high-resolution structure determination. This review presents different approaches that have been used to study the structure and dynamics of T4aP and T2SS endopili by means of integrative structural biology, cryo-electron microscopy (cryo-EM), and molecular dynamics simulations. Their conserved features and differences are presented. The non-helical stretch in the long-conserved N-terminal helix which is characteristic of all members of the TFF and the impact of calcium on structure, function, and dynamics of these filaments are discussed in detail.

Assessing the Stability of Biological Fibrils by Molecular-Scale Simulations

The nanomechanical characterization of several biological fibrils that are the result of protein aggregation via molecular dynamics simulation is nowadays feasible, and together with atomic force microscopy experiments has widened our understanding of the forces in the regime of pN-nN and system sizes of about hundreds of nanometers. Several methodologies have been developed to achieve this target, and they range from the atomistic representation via molecular force fields to coarse-grained strategies that provide comparable results with experiments in a systematic way. In this chapter, we discuss several methodologies for the calculation of mechanical parameters, such as the elastic constants of relevant biological systems. They are presented together with details about parameterization and current limitations. Then, we discuss some of the applications of such methodologies for the description of bacterial filament and β-amyloid systems. Finally, the latest lines of development are discussed.

Computational and biochemical analysis of type IV pilus dynamics and stability

Type IV pili (T4P) are distinctive dynamic filaments at the surface of many bacteria that can rapidly extend and retract and withstand strong forces. T4P are important virulence factors in many human pathogens, including Enterohemorrhagic Escherichia coli (EHEC). The structure of the EHEC T4P has been determined by integrating nuclear magnetic resonance (NMR) and cryo-electron microscopy data. To better understand pilus assembly, stability, and function, we performed a total of 108 ms all-atom molecular dynamics simulations of wild-type and mutant T4P. Extensive characterization of the conformational landscape of T4P in different conditions of temperature, pH, and ionic strength is complemented with targeted mutagenesis and biochemical analyses. Our simulations and NMR experiments reveal a conserved set of residues defining a calcium-binding site at the interface between three pilin subunits. Calcium binding enhances T4P stability ex vivo and in vitro, supporting the role of this binding site as a potential pocket for drug design.

Nanomechanical Stability of Aβ Tetramers and Fibril-like Structures: Molecular Dynamics Simulations

Alzheimer's disease (AD) is a neurodegenerative disorder and one of the main causes of dementia. The disease is associated with amyloid beta (Aβ) peptide aggregation forming initial clusters and then fibril structure and plaques. Other neurodegenerative diseases such as type 2 diabetes, amyotrophic lateral sclerosis, and Parkinson's disease follow a similar mechanism. Therefore, inhibition of Aβ aggregation is considered an effective way to prevent AD. Recent experiments have provided evidence that oligomers are more toxic agents than mature fibrils, prompting researchers to investigate various factors that may influence their properties. One of these factors is nanomechanical stability, which plays an important role in the self-assembly of Aβ and possibly other proteins. This stability is also likely to be related to cell toxicity. In this work, we compare the mechanical stability of Aβ-tetramers and fibrillar structures using a structure-based coarse-grained (CG) approach and all-atom molecular dynamics simulation. Our results support the evidence for an increase in mechanical stability during the Aβ fibrillization process, which is consistent with in vitro AFM characterization of Aβ42 oligomers. Namely, using a CG model, we showed that the Young modulus of tetramers is lower than that of fibrils and, as follows from the experiment, is about 1 GPa. Hydrogen bonds are the dominant contribution to the detachment of one chain from the Aβ fibril fragment. They tend to be more organized along the pulling direction, whereas in the Aβ tetramers no preference is observed.

Impact of an alpha helix and a cysteine-cysteine disulfide bond on the resistance of bacterial adhesion pili to stress

Escherichia coli express adhesion pili that mediate attachment to host cell surfaces and are exposed to body fluids in the urinary and gastrointestinal tracts. Pilin subunits are organized into helical polymers, with a tip adhesin for specific host binding. Pili can elastically unwind when exposed to fluid flow forces, reducing the adhesin load, thereby facilitating sustained attachment. Here we investigate biophysical and structural differences of pili commonly expressed on bacteria that inhabit the urinary and intestinal tracts. Optical tweezers measurements reveal that class 1a pili of uropathogenic E. coli (UPEC), as well as class 1b of enterotoxigenic E. coli (ETEC), undergo an additional conformational change beyond pilus unwinding, providing significantly more elasticity to their structure than ETEC class 5 pili. Examining structural and steered molecular dynamics simulation data, we find that this difference in class 1 pili subunit behavior originates from an α-helical motif that can unfold when exposed to force. A disulfide bond cross-linking β-strands in class 1 pili stabilizes subunits, allowing them to tolerate higher forces than class 5 pili that lack this covalent bond. We suggest that these extra contributions to pilus resiliency are relevant for the UPEC niche, since resident bacteria are exposed to stronger, more transient drag forces compared to those experienced by ETEC bacteria in the mucosa of the intestinal tract. Interestingly, class 1b ETEC pili include the same structural features seen in UPEC pili, while requiring lower unwinding forces that are more similar to those of class 5 ETEC pili.

Unravelling the molecular basis of AM-6494 high potency at BACE1 in Alzheimer's disease: an integrated dynamic interaction investigation

β-amyloid precursor protein cleaving enzyme1 (BACE1) has prominently been an important drug design target implicated in Alzheimer's disease pathway. The failure rate of most of the already tested drugs at different clinical phases remains a major concern. Recently, AM-6494 was reported as a novel potent, highly selective, and orally effective inhibitor against BACE1. AM-6494 displayed no alteration of skin/fur colour in animal studies, an adverse effect common to previous BACE1 inhibitors. However, the atomistic molecular mechanism of BACE1 inhibition by AM-6494 remains unclear. To elucidate the binding mechanism of AM-6494 relative to umibecestat (CNP-520) as well as the structural changes when bound to BACE1, advanced computational techniques such as accelerated MD simulation and principal component analysis have been utilised. The results demonstrated higher binding affinity of AM-6494 at BACE1 with van der Waals as dominant energy contributor compared to umibecestat. Conformational monitoring of the β-hairpin flap covering the active site revealed an effective flap closure when bound with AM-6494 compared to CNP-520, which predominantly alternates between semi-open and closed conformations. The observed effective flap closure of AM-6494 explains its higher inhibitory power towards BACE1. Besides the catalytic Asp32/228 dyad, Tyr14, Leu30, Tyr71 and Gly230 represent critical residues in the potency of these inhibitors at BACE1 binding interface. The findings highlighted in this research provide a basis to explain AM-6494 high inhibitory potency and might assist in the design of new inhibitors with improved selectivity and potency for BACE1.

Perturb-Scan-Pull: A Novel Method Facilitating Conformational Transitions in Proteins

Conformational transitions in proteins facilitate precise physiological functions. Therefore, it is crucial to understand the mechanisms underlying these processes to modulate protein function. Yet, studying structural and dynamical properties of proteins is notoriously challenging due to the complexity of the underlying potential energy surfaces (PES). We have previously developed the perturbation-response scanning (PRS) method to identify key residues that participate in the communication network responsible for specific conformational transitions. PRS is based on a residue-by-residue scan of the protein to determine the subset of residues/forces which provide the closest conformational change leading to a target conformational state, inasmuch as linear response theory applies to these motions. Here, we develop a novel method to further evaluate if conformational transitions may be triggered on the PES. We aim to study functionally relevant conformational transitions in proteins by using results obtained from PRS and feeding them as inputs to steered molecular dynamics simulations. The success and the transferability of the method are evaluated on three protein systems having different complexities of motion on the PES: calmodulin, adenylate kinase, and bacterial ferric binding protein. We find that the method captures the target conformation, while providing key residues and the optimum paths with relatively low free energy profiles.

The Role of Glycans in Bacterial Adhesion to Mucosal Surfaces: How Can Single-Molecule Techniques Advance Our Understanding?

Bacterial adhesion is currently the subject of increased interest from the research community, leading to fast progress in our understanding of this complex phenomenon. Resent research within this field has documented the important roles played by glycans for bacterial surface adhesion, either through interaction with lectins or with other glycans. In parallel with this increased interest for and understanding of bacterial adhesion, there has been a growth in the sophistication and use of sensitive force probes for single-molecule and single cell studies. In this review, we highlight how the sensitive force probes atomic force microscopy (AFM) and optical tweezers (OT) have contributed to clarifying the mechanisms underlying bacterial adhesion to glycosylated surfaces in general and mucosal surfaces in particular. We also describe research areas where these techniques have not yet been applied, but where their capabilities appear appropriate to advance our understanding.

Strength of Neisseria meningitidis binding to endothelial cells requires highly-ordered CD147/β2-adrenoceptor clusters assembled by alpha-actinin-4

Neisseria meningitidis (meningococcus) is an invasive bacterial pathogen that colonizes human vessels, causing thrombotic lesions and meningitis. Establishment of tight interactions with endothelial cells is crucial for meningococci to resist haemodynamic forces. Two endothelial receptors, CD147 and the β2-adrenergic receptor (β₂AR), are sequentially engaged by meningococci to adhere and promote signalling events leading to vascular colonization, but their spatiotemporal coordination is unknown. Here we report that CD147 and β₂AR form constitutive hetero-oligomeric complexes. The scaffolding protein α-actinin-4 directly binds to the cytosolic tail of CD147 and governs the assembly of CD147–β₂AR complexes in highly ordered clusters at bacterial adhesion sites. This multimolecular assembly process increases the binding strength of meningococci to endothelial cells under shear stress, and creates molecular platforms for the elongation of membrane protrusions surrounding adherent bacteria. Thus, the specific organization of cellular receptors has major impacts on host–pathogen interaction.

Neisseria meningitidis bacteria bind to host proteins CD147 and β2-adrenergic receptor on the surface of endothelial cells. Here, Maïssa et al. show that the two proteins interact with each other forming clusters that increase the binding strength of the bacteria to endothelial cells.

Mechanical response and deformation mechanics of Type IV pili investigated using steered coarse-grained molecular dynamics simulation

Type IV pili are long filamentous structures on the surface of bacteria, which can be rapidly assembled or disassembled with pilin subunits by molecular motors. They can generate force during retraction and are involved in many bacterial functions. Steered molecular dynamics simulations with coarse-grained MARTINI models are carried out to investigate the mechanical behaviors of pili under tension. Our study is the first to report a Young's modulus of 0.80±0.07GPa and a spring constant of 1294.6±116.5kJmol^-1nm^-2 for pilus. Our results show the mechanical responses of pili are different from those described by the worm-like chain model and the van der Waal's interactions play a critical role in the mechanical responses. Moreover, the effects of pulling rates and virtual spring constants of pilus on Young's modulus are studied and two distinct morphological stages with the conformational changes appear during the extension of pilus are observed. This work provide insight into the mechanics and the deformation mechanism of pilus assembly.

Emergence of complex behavior in pili-based motility in early stages of P. aeruginosa surface adaptation

Pseudomonas aeruginosa move across surfaces by using multiple Type IV Pili (TFP), motorized appendages capable of force generation via linear extension/retraction cycles, to generate surface motions collectively known as twitching motility. Pseudomonas cells arrive at a surface with low levels of piliation and TFP activity, which both progressively increase as the cells sense the presence of a surface. At present, it is not clear how twitching motility emerges from these initial minimal conditions. Here, we build a simple model for TFP-driven surface motility without complications from viscous and solid friction on surfaces. We discover the unanticipated structural requirement that TFP motors need to have a minimal amount of effective angular rigidity in order for cells to perform the various classes of experimentally-observed motions. Moreover, a surprisingly small number of TFP are needed to recapitulate movement signatures associated with twitching: Two TFP can already produce movements reminiscent of recently observed slingshot type motion. Interestingly, jerky slingshot motions characteristic of twitching motility comprise the transition region between different types of observed crawling behavior in the dynamical phase diagram, such as self-trapped localized motion, 2-D diffusive exploration, and super-diffusive persistent motion.

Revealing Atomic-Level Mechanisms of Protein Allostery with Molecular Dynamics Simulations

Molecular dynamics (MD) simulations have become a powerful and popular method for the study of protein allostery, the widespread phenomenon in which a stimulus at one site on a protein influences the properties of another site on the protein. By capturing the motions of a protein's constituent atoms, simulations can enable the discovery of allosteric binding sites and the determination of the mechanistic basis for allostery. These results can provide a foundation for applications including rational drug design and protein engineering. Here, we provide an introduction to the investigation of protein allostery using molecular dynamics simulation. We emphasize the importance of designing simulations that include appropriate perturbations to the molecular system, such as the addition or removal of ligands or the application of mechanical force. We also demonstrate how the bidirectional nature of allostery-the fact that the two sites involved influence one another in a symmetrical manner-can facilitate such investigations. Through a series of case studies, we illustrate how these concepts have been used to reveal the structural basis for allostery in several proteins and protein complexes of biological and pharmaceutical interest.

How Bacteria Use Type IV Pili Machinery on Surfaces

The bacterial type IV pilus (T4P) is a versatile molecular machine with a broad range of functions. Recent advances revealed that the molecular components and the biophysical properties of the machine are well conserved among phylogenetically distant bacterial species. However, its functions are diverse, and include adhesion, motility, and horizontal gene transfer. This review focusses on the role of T4P in surface motility and bacterial interactions. Different species have evolved distinct mechanisms for intracellular coordination of multiple pili and of pili with other motility machines, ranging from physical coordination to biochemical clocks. Coordinated behavior between multiple bacteria on a surface is achieved by active manipulation of surfaces and modulation of pilus-pilus interactions. An emerging picture is that the T4P actively senses and responds to environmental conditions.

Structural and functional insights into the conductive pili of Geobacter sulfurreducens revealed in molecular dynamics simulations

Geobacter sulfurreducens (GS) electronically connects with extracellular electron acceptors using conductive protein filaments or pili. To gain insights into their role as biological nanowires, we investigated the structural dynamics of the GS pilus in solution via molecular dynamics simulations. In the model, all of the pilin's aromatics clustered as a right-handed helical band along the pilus, maintaining inter-aromatic distances and dimer configurations optimal for multistep hopping. The aromatics were interspersed within the regions of highest negative potential, which influenced the type and configuration of the aromatic contacts and the rates of electron transfer. Small foci of positive potential were also present but were neutralized within uncharged regions, thus minimizing charge trapping. Consistent with the model predictions, mutant strains with reduced aromatic contacts or negative potentials had defects in pili functions such as the reduction of Fe(III) oxides and electrodes. The results therefore support the notion of a pilus fiber evolved to function as an electronic conduit between the cell and extracellular electron acceptors.

Accelerated molecular dynamics and protein conformational change: a theoretical and practical guide using a membrane embedded model neurotransmitter transporter

Molecular dynamics simulation provides a powerful and accurate method to model protein conformational change, yet timescale limitations often prevent direct assessment of the kinetic properties of interest. A large number of molecular dynamic steps are necessary for rare events to occur, which allow a system to overcome energy barriers and conformationally transition from one potential energy minimum to another. For many proteins, the energy landscape is further complicated by a multitude of potential energy wells, each separated by high free-energy barriers and each potentially representative of a functionally important protein conformation. To overcome these obstacles, accelerated molecular dynamics utilizes a robust bias potential function to simulate the transition between different potential energy minima. This straightforward approach more efficiently samples conformational space in comparison to classical molecular dynamics simulation, does not require advanced knowledge of the potential energy landscape and converges to the proper canonical distribution. Here, we review the theory behind accelerated molecular dynamics and discuss the approach in the context of modeling protein conformational change. As a practical example, we provide a detailed, step-by-step explanation of how to perform an accelerated molecular dynamics simulation using a model neurotransmitter transporter embedded in a lipid cell membrane. Changes in protein conformation of relevance to the substrate transport cycle are then examined using principle component analysis.

Exceptionally widespread nanomachines composed of type IV pilins: the prokaryotic Swiss Army knives

Prokaryotes have engineered sophisticated surface nanomachines that have allowed them to colonize Earth and thrive even in extreme environments. Filamentous machineries composed of type IV pilins, which are associated with an amazing array of properties ranging from motility to electric conductance, are arguably the most widespread since distinctive proteins dedicated to their biogenesis are found in most known species of prokaryotes. Several decades of investigations, starting with type IV pili and then a variety of related systems both in bacteria and archaea, have outlined common molecular and structural bases for these nanomachines. Using type IV pili as a paradigm, we will highlight in this review common aspects and key biological differences of this group of filamentous structures.

Using type IV pili as a paradigm, we review common genetic, structural and mechanistic features (many) as well as differences (few) of the exceptionally widespread and functionally versatile prokaryotic nano-machines composed of type IV pilins.

The hypervariable region of meningococcal major pilin PilE controls the host cell response via antigenic variation

Type IV pili (Tfp) are expressed by many Gram-negative bacteria to promote aggregation, adhesion, internalization, twitching motility, or natural transformation. Tfp of Neisseria meningitidis, the causative agent of cerebrospinal meningitis, are involved in the colonization of human nasopharynx. After invasion of the bloodstream, Tfp allow adhesion of N. meningitidis to human endothelial cells, which leads to the opening of the blood-brain barrier and meningitis. To achieve firm adhesion, N. meningitidis induces a host cell response that results in elongation of microvilli surrounding the meningococcal colony. Here we study the role of the major pilin subunit PilE during host cell response using human dermal microvascular endothelial cells and the pharynx carcinoma-derived FaDu epithelial cell line. We first show that some PilE variants are unable to induce a host cell response. By engineering PilE mutants, we observed that the PilE C-terminus domain, which contains a disulfide bonded region (D-region), is critical for the host cell response and that hypervariable regions confer different host cell specificities. Moreover, the study of point mutants of the pilin D-region combined with structural modeling of PilE revealed that the D-region contains two independent regions involved in signaling to human dermal microvascular endothelial cells (HDMECs) or FaDu cells. Our results indicate that the diversity of the PilE D-region sequence allows the induction of the host cell response via several receptors. This suggests that Neisseria meningitidis has evolved a powerful tool to adapt easily to many niches by modifying its ability to interact with host cells.

Type IV pili (Tfp) are long appendages expressed by many Gram-negative bacteria, including Neisseria meningitidis, the causative agent of cerebrospinal meningitis. These pili are involved in many aspects of pathogenesis: natural competence, aggregation, adhesion, and twitching motility. More specifically, Neisseria meningitidis, which is devoid of a secretion system to manipulate its host, has evolved its Tfp to signal to brain endothelial cells and open the blood-brain barrier. In this report, we investigate, at the molecular level, the involvement of the major pilin subunit PilE in host cell response. Our results indicate that the PilE C-terminal domain, which contains a disulfide bonded region (D-region), is critical for the host cell response and contains two independent regions involved in host cell signaling.