Nanomechanical properties of MscL α helices: A steered molecular dynamics study

Probing membrane deformation energy by KcsA potassium channel gating under varied membrane thickness and tension

This study investigated how membrane thickness and tension modify the gating of KcsA potassium channels when simultaneously varied. The KcsA channel undergoes global conformational changes upon gating: expansion of the cross-sectional area and longitudinal shortening upon opening. Thus, membranes impose differential effects on the open and closed conformations, such as hydrophobic mismatches. Here, the single-channel open probability was recorded in the contact bubble bilayer, by which variable thickness membranes under a defined tension were applied. A fully open channel in thin membranes turned to sporadic openings in thick membranes, where the channel responded moderately to tension increase. Quantitative gating analysis prompted the hypothesis that tension augmented the membrane deformation energy when hydrophobic mismatch was enhanced in thick membranes.

Probing the mechanical properties of ORF3a protein, a transmembrane channel of SARS-CoV-2 virus: Molecular dynamics study

SARS-CoV-2-encoded accessory protein ORF3a was found to be a conserved coronavirus protein that shows crucial roles in apoptosis in cells as well as in virus release and replications. To complete the knowledge and identify the unknown of this protein, further comprehensive research is needed to clarify the leading role of ORF3a in the functioning of the coronavirus. One of the efficient approaches to determining the functionality of this protein is to investigate the mechanical properties and study its structural dynamics in the presence of physical stimuli. Herein, performing all-atom steered molecular dynamics (SMD) simulations, the mechanical properties of the force-bearing components of the ORF3a channel are calculated in different physiological conditions. As variations occurring in ORF3a may lead to alteration in protein structure and function, the G49V mutation was also simulated to clarify the relationship between the mechanical properties and chemical stability of the protein by comparing the behavior of the wild-type and mutant Orf3a. From a physiological conditions point of view, it was observed that in the solvated system, the presence of water molecules reduces Young's modulus of TM1 by ∼30 %. Our results also show that by substitution of Gly49 with valine, Young's modulus of the whole helix increases from 1.61 ± 0.20 to 2.08 ± 0.15 GPa, which is consistent with the calculated difference in free energy of wild-type and mutant helices. In addition to finding a way to fight against Covid-19 disease, understanding the mechanical behavior of these biological nanochannels can lead to the development of the potential applications of the ORF3a protein channel, such as tunable nanovalves in smart drug delivery systems, nanofilters in the new generation of desalination systems, and promising applications in DNA sequencing.

Pocket delipidation induced by membrane tension or modification leads to a structurally analogous mechanosensitive channel state

The mechanosensitive ion channel of large conductance MscL gates in response to membrane tension changes. Lipid removal from transmembrane pockets leads to a concerted structural and functional MscL response, but it remains unknown whether there is a correlation between the tension-mediated state and the state derived by pocket delipidation in the absence of tension. Here, we combined pulsed electron paramagnetic resonance spectroscopy and hydrogen-deuterium exchange mass spectrometry, coupled with molecular dynamics simulations under membrane tension, to investigate the structural changes associated with the distinctively derived states. Whether it is tension- or modification-mediated pocket delipidation, we find that MscL samples a similar expanded subconducting state. This is the final step of the delipidation pathway, but only an intermediate stop on the tension-mediated path, with additional tension triggering further channel opening. Our findings hint at synergistic modes of regulation by lipid molecules in membrane tension-activated mechanosensitive channels.

Principal component analysis of alpha-helix deformations in transmembrane proteins

α-helices are deformable secondary structural components regularly observed in protein folds. The overall flexibility of an α-helix can be resolved into constituent physical deformations such as bending in two orthogonal planes and twisting along the principal axis. We used Principal Component Analysis to identify and quantify the contribution of each of these dominant deformation modes in transmembrane α-helices, extramembrane α-helices, and α-helices in soluble proteins. Using three α-helical samples from Protein Data Bank entries spanning these three cellular contexts, we determined that the relative contributions of these modes towards total deformation are independent of the α-helix's surroundings. This conclusion is supported by the observation that the identities of the top three deformation modes, the scaling behaviours of mode eigenvalues as a function of α-helix length, and the percentage contribution of individual modes on total variance were comparable across all three α-helical samples. These findings highlight that α-helical deformations are independent of cellular location and will prove to be valuable in furthering the development of flexible templates in de novo protein design.

Unfolding mechanism and free energy landscape of single, stable, alpha helices at low pull speeds

Single alpha helices (SAHs) stable in isolated form are often found in motor proteins where they bridge functional domains. Understanding the mechanical response of SAHs is thus critical to understand their function. The quasi-static force-extension relation of a small number of SAHs is known from single-molecule experiments. Unknown, or still controversial, are the molecular scale details behind those observations. We show that the deformation mechanism of SAHs pulled from the termini at pull speeds approaching the quasi-static limit differs from that of typical helices found in proteins, which are stable only when interacting with other protein domains. Using molecular dynamics simulations with atomistic resolution at low pull speeds previously inaccessible to simulation, we show that SAHs start unfolding from the termini at all pull speeds we investigated. Unfolding proceeds residue-by-residue and hydrogen bond breaking is not the main event determining the barrier to unfolding. We use the molecular simulation data to test the cooperative sticky chain model. This model yields excellent fits of the force-extension curves and quantifies the distance, xE = 0.13 nm, to the transition state, the natural frequency of bond vibration, ν0 = 0.82 ns-1, and the height, V0 = 2.9 kcal mol-1, of the free energy barrier associated with the deformation of single residues. Our results demonstrate that the sticky chain model could advantageously be used to analyze experimental force-extension curves of SAHs and other biopolymers.

Remembering Mechanosensitivity of NMDA Receptors

An increase in post-synaptic Ca²⁺ conductance through activation of the ionotropic N-methyl-D-aspartate receptor (NMDAR) and concomitant structural changes are essential for the initiation of long-term potentiation (LTP) and memory formation. Memories can be initiated by coincident events, as occurs in classical conditioning, where the NMDAR can act as a molecular coincidence detector. Binding of glutamate and glycine, together with depolarization of the postsynaptic cell membrane to remove the Mg²⁺ channel pore block, results in NMDAR opening for Ca²⁺ conductance. Accumulating evidence has implicated both force-from-lipids and protein tethering mechanisms for mechanosensory transduction in NMDAR, which has been demonstrated by both, membrane stretch and application of amphipathic molecules such as arachidonic acid (AA). The contribution of mechanosensitivity to memory formation and consolidation may be to increase activity of the NMDAR leading to facilitated memory formation. In this review we look back at the progress made toward understanding the physiological and pathological role of NMDA receptor channels in mechanobiology of the nervous system and consider these findings in like of their potential functional implications for memory formation. We examine recent studies identifying mechanisms of both NMDAR and other mechanosensitive channels and discuss functional implications including gain control of NMDA opening probability. Mechanobiology is a rapidly growing area of biology with many important implications for understanding form, function and pathology in the nervous system.

Conformational Flexibility of the Protein-Protein Interfaces of the Ebola Virus VP40 Structural Matrix Filament

The Ebola virus (EBOV) is a virulent pathogen that causes severe hemorrhagic fever with a high fatality rate in humans. The EBOV transformer protein VP40 plays crucial roles in viral assembly and budding at the plasma membrane of infected cells. One of VP40's roles is to form the long, flexible, pleomorphic filamentous structural matrix for the virus. Each filament contains three unique interfaces: monomer NTD-NTD to form a dimer, dimer-to-dimer NTD-NTD oligomerization to form a hexamer, and end-to-end hexamer CTD-CTD to build the filament. However, the atomic-level details of conformational flexibility of the VP40 filament are still elusive. In this study, we have performed explicit-solvent, all-atom molecular dynamic simulations to explore the conformational flexibility of the three different interface structures of the filament. Using dynamic network analysis and other calculational methods, we find that the CTD-CTD hexamer interface with weak interdomain amino acid communities is the most flexible, and the NTD-NTD oligomer interface with strong interdomain communities is the least flexible. Our study suggests that the high flexibility of the CTD-CTD interface may be essential for the supple bending of the Ebola filovirus, and such flexibility may present a target for molecular interventions to disrupt the Ebola virus functioning.

Nanomechanical DNA resonators for sensing and structural analysis of DNA-ligand complexes

The effect of direct or indirect binding of intercalant molecules on DNA structure is of fundamental importance in understanding the biological functioning of DNA. Here we report on self-suspended DNA nanobundles as ultrasensitive nanomechanical resonators for structural studies of DNA-ligand complexes. Such vibrating nanostructures represent the smallest mechanical resonator entirely composed of DNA. A correlative analysis between the mechanical and structural properties is exploited to study the intrinsic changes of double strand DNA, when interacting with different intercalant molecules (YOYO-1 and GelRed) and a chemotherapeutic drug (Cisplatin), at different concentrations. Possible implications of our findings are related to the study of interaction mechanism of a wide category of molecules with DNA, and to further applications in medicine, such as optimal titration of chemotherapeutic drugs and environmental studies for the detection of heavy metals in human serum.

Constitutive boost of a K+ channel via inherent bilayer tension and a unique tension-dependent modality

Molecular mechanisms underlying channel-membrane interplay have been extensively studied. Cholesterol, as a major component of the cell membrane, participates either in specific binding to channels or via modification of membrane physical features. Here, we examined the action of various sterols (cholesterol, epicholesterol, etc.) on a prototypical potassium channel (KcsA). Single-channel current recordings of the KcsA channel were performed in a water-in-oil droplet bilayer (contact bubble bilayer) with a mixed phospholipid composition (azolectin). Upon membrane perfusion of sterols, the activated gate at acidic pH closed immediately, irrespective of the sterol species. During perfusion, we found that the contacting bubbles changed their shapes, indicating alterations in membrane physical features. Absolute bilayer tension was measured according to the principle of surface chemistry, and inherent bilayer tension was ∼5 mN/m. All tested sterols decreased the tension, and the nonspecific sterol action to the channel was likely mediated by the bilayer tension. Purely mechanical manipulation that reduced bilayer tension also closed the gate, whereas the resting channel at neutral pH never activated upon increased tension. Thus, rather than conventional stretch activation, the channel, once ready to activate by acidic pH, changes the open probability through the action of bilayer tension. This constitutes a channel regulating modality by two successive stimuli. In the contact bubble bilayer, inherent bilayer tension was high, and the channel remained boosted. In the cell membrane, resting tension is low, and it is anticipated that the ready-to-activate channel remains closed until bilayer tension reaches a few millinewton/meter during physiological and pathological cellular activities.

Structural Dynamics of the MscL C-terminal Domain

The large conductance mechanosensitive channel (MscL), acts as an osmoprotective emergency valve in bacteria by opening a large, water-filled pore in response to changes in membrane tension. In its closed configuration, the last 36 residues at the C-terminus form a bundle of five α-helices co-linear with the five-fold axis of symmetry. Here, we examined the structural dynamics of the C-terminus of EcMscL using site-directed spin labelling electron paramagnetic resonance (SDSL EPR) spectroscopy. These experiments were complemented with computational modelling including molecular dynamics (MD) simulations and finite element (FE) modelling. Our results show that under physiological conditions, the C-terminus is indeed an α-helical bundle, located near the five-fold symmetry axis of the molecule. Both experiments and computational modelling demonstrate that only the top part of the C-terminal domain (from the residue A110 to E118) dissociates during the channel gating, while the rest of the C-terminus stays assembled. This result is consistent with the view that the C-terminus functions as a molecular sieve and stabilizer of the oligomeric MscL structure as previously suggested.

Pulling MscL open via N-terminal and TM1 helices: A computational study towards engineering an MscL nanovalve

There are great opportunities in the manipulation of bacterial mechanosensitive (MS) ion channels for specific and targeted drug delivery purposes. Recent research has shown that these ion channels have the potential to be converted into nanovalves through clever use of magnetic nanoparticles and magnetic fields. Using a combination of molecular dynamics (MD) simulations and the finite element (FE) modelling, this study investigates the theoretical feasibility of opening the MscL channel (MS channel of large conductance of E. coli) by applying mechanical force directly to its N-terminus. This region has already been reported to function as a major mechanosensor in this channel. The stress-strain behaviour of each MscL helix was obtained using all atom MD simulations. Using the same method, we simulated two models, the wild-type (WT) MscL and the G22N mutant MscL, both embedded in a POPE lipid bilayer. In addition to indicating the main interacting residues at the hydrophobic pore, their pairwise interaction energies were monitored during the channel gating. We implemented these inputs into our FE model of MscL using curve-fitting codes and continuum mechanics equations. In the FE model, the channel could be fully opened via pulling directly on the N-terminus and bottom of TM1 by mutating dominant van der Waals interactions in the channel pore; otherwise the stress generated on the channel protein can irreversibly unravel the N-secondary structure. This is a significant finding suggesting that applying force in this manner is sufficient to open an MscL nanovalve delivering various drugs used, for example, in cancer chemotherapy. More importantly, the FE model indicates that to fully operate an MscL nanovalve by pulling directly on the N-terminus and bottom of TM1, gain-of-function (GOF) mutants (e.g., G22N MscL) would have to be employed rather than the WT MscL channel.

Nanomechanical properties of MscL α helices: A steered molecular dynamics study

Gating of mechanosensitive (MS) channels is driven by a hierarchical cascade of movements and deformations of transmembrane helices in response to bilayer tension. Determining the intrinsic mechanical properties of the individual transmembrane helices is therefore central to understanding the intricacies of the gating mechanism of MS channels. We used a constant-force steered molecular dynamics (SMD) approach to perform unidirectional pulling tests on all the helices of MscL in M. tuberculosis and E. coli homologs. Using this method, we could overcome the issues encountered with the commonly used constant-velocity SMD simulations, such as low mechanical stability of the helix during stretching and high dependency of the elastic properties on the pulling rate. We estimated Young's moduli of the α-helices of MscL to vary between 0.2 and 12.5 GPa with TM2 helix being the stiffest. We also studied the effect of water on the properties of the pore-lining TM1 helix. In the absence of water, this helix exhibited a much stiffer response. By monitoring the number of hydrogen bonds, it appears that water acts like a 'lubricant' (softener) during TM1 helix elongation. These data shed light on another physical aspect underlying hydrophobic gating of MS channels, in particular MscL.