Reconciling alternate methods for the determination of charge distributions: A probabilistic approach to high-dimensional least-squares approximations

We propose extensions and improvements of the statistical analysis of distributed multipoles (SADM) algorithm put forth by Chipot et al. in [6] for the derivation of distributed atomic multipoles from the quantum-mechanical electrostatic potential. The method is mathematically extended to general least-squares problems and provides an alternative approximation method in cases where the original least-squares problem is computationally not tractable, either because of its ill-posedness or its high-dimensionality. The solution is approximated employing a Monte Carlo method that takes the average of a random variable defined as the solutions of random small least-squares problems drawn as subsystems of the original problem. The conditions that ensure convergence and consistency of the method are discussed, along with an analysis of the computational cost in specific instances

Subdiffusion in Membrane Permeation of Small Molecules

Citation: Chipot, C. and Comer, J. Subdiffusion in Membrane Permeation of Small Molecules. Sci. Rep. 6, 35913; doi: 10.1038/srep35913 (2016).Within the solubility–diffusion model of passive membrane permeation of small molecules, translocation of the permeant across the biological membrane is traditionally assumed to obey the Smoluchowski diffusion equation, which is germane for classical diffusion on an inhomogeneous free-energy and diffusivity landscape. This equation, however, cannot accommodate subdiffusive regimes, which have long been recognized in lipid bilayer dynamics, notably in the lateral diffusion of individual lipids. Through extensive biased and unbiased molecular dynamics simulations, we show that one-dimensional translocation of methanol across a pure lipid membrane remains subdiffusive on timescales approaching typical permeation times. Analysis of permeant motion within the lipid bilayer reveals that, in the absence of a net force, the mean squared displacement depends on time as t0.7, in stark contrast with the conventional model, which assumes a strictly linear dependence. We further show that an alternate model using a fractional-derivative generalization of the Smoluchowski equation provides a rigorous framework for describing the motion of the permeant molecule on the pico- to nanosecond timescale. The observed subdiffusive behavior appears to emerge from a crossover between small-scale rattling of the permeant around its present position in the membrane and larger-scale displacements precipitated by the formation of transient voids

Mitochondrial ADP/ATP Carrier in Dodecylphosphocholine Binds Cardiolipins with Non-native Affinity.

Biophysical investigation of membrane proteins generally requires their extraction from native sources using detergents, a step that can lead, possibly irreversibly, to protein denaturation. The propensity of dodecylphosphocholine (DPC), a detergent widely utilized in NMR studies of membrane proteins, to distort their structure has been the subject of much controversy. It has been recently proposed that the binding specificity of the yeast mitochondrial ADP/ATP carrier (yAAC3) toward cardiolipins is preserved in DPC, thereby suggesting that DPC is a suitable environment in which to study membrane proteins. In this communication, we used all-atom molecular dynamics simulations to investigate the specific binding of cardiolipins to yAAC3. Our data demonstrate that the interaction interface observed in a native-like environment differs markedly from that inferred from an NMR investigation in DPC, implying that in this detergent, the protein structure is distorted. We further investigated yAAC3 solubilized in DPC and in the milder dodecylmaltoside with thermal-shift assays. The loss of thermal transition observed in DPC confirms that the protein is no longer properly folded in this environment

Structure of HIV-1 gp41 with its membrane anchors targeted by neutralizing antibodies

The HIV-1 gp120/gp41 trimer undergoes a series of conformational changes in order to catalyze gp41-induced fusion of viral and cellular membranes. Here, we present the crystal structure of gp41 locked in a fusion intermediate state by an MPER-specific neutralizing antibody. The structure illustrates the conformational plasticity of the six membrane anchors arranged asymmetrically with the fusion peptides and the transmembrane regions pointing into different directions. Hinge regions located adjacent to the fusion peptide and the transmembrane region facilitate the conformational flexibility that allows high-affinity binding of broadly neutralizing anti-MPER antibodies. Molecular dynamics simulation of the MPER Ab-stabilized gp41 conformation reveals a possible transition pathway into the final post-fusion conformation with the central fusion peptides forming a hydrophobic core with flanking transmembrane regions. This suggests that MPER-specific broadly neutralizing antibodies can block final steps of refolding of the fusion peptide and the transmembrane region, which is required for completing membrane fusion

Structural and energetic study of cation-p-cation interactions in proteins

Cation-pi interactions of aromatic rings and positively charged groups are among the most important interactions in structural biology. The role and energetic characteristics of these interactions are well established. However, the occurrence of cation-pi-cation interactions is an unexpected motif, which raises intriguing questions about its functional role in proteins. We present a statistical analysis of the occurrence, composition and geometrical preferences of cation-pi-cation interactions identified in a set of non-redundant protein structures taken from the Protein Data Bank. Our results demonstrate that this structural motif is observed at a small, albeit non-negligible frequency in proteins, and suggest a preference to establish cation-pi-cation motifs with Trp, followed by Tyr and Phe. Furthermore, we have found that cation-pi-cation interactions tend to be highly conserved, which supports their structural or functional role. Finally, we have performed an energetic analysis of a representative subset of cation-pi-cation complexes combining quantum-chemical and continuum solvation calculations. Our results point out that the protein environment can strongly screen the cation-cation repulsion, leading to an attractive interaction in 64% of the complexes analyzed. Together with the high degree of conservation observed, these results suggest a potential stabilizing role in the protein fold, as demonstrated recently for a miniature protein (Craven et al., J. Am. Chem. Soc. 2016, 138, 1543). From a computational point of view, the significant contribution of non-additive three-body terms challenges the suitability of standard additive force fields for describing cation-p-cation motifs in molecular simulations. Keywords: Cation-π−cation complexes; noncovalent interactions; cooperativity; protein structure

NMR structure and ion channel activity of the p7 protein from hepatitis C virus

The small membrane protein p7 of hepatitis C virus forms oligomers and exhibits ion channel activity essential for virus infectivity. These viroporin features render p7 an attractive target for antiviral drug development. In this study, p7 from strain HCV-J (genotype 1b) was chemically synthesized and purified for ion channel activity measurements and structure analyses. p7 forms cation-selective ion channels in planar lipid bilayers and at the single-channel level by the patch clamp technique. Ion channel activity was shown to be inhibited by hexamethylene amiloride but not by amantadine. Circular dichroism analyses revealed that the structure of p7 is mainly α-helical, irrespective of the membrane mimetic medium (e.g. lysolipids, detergents, or organic solvent/water mixtures). The secondary structure elements of the monomeric form of p7 were determined by 1H and 13C NMR in trifluoroethanol/water mixtures. Molecular dynamics simulations in a model membrane were combined synergistically with structural data obtained from NMR experiments. This approach allowed us to determine the secondary structure elements of p7, which significantly differ from predictions, and to propose a three-dimensional model of the monomeric form of p7 associated with the phospholipid bilayer. These studies revealed the presence of a turn connecting an unexpected N-terminal α-helix to the first transmembrane helix, TM1, and a long cytosolic loop bearing the dibasic motif and connecting TM1 to TM2. These results provide the first detailed experimental structural framework for a better understanding of p7 processing, oligomerization, and ion channel gating mechanism.Instituto de Física La Plat

How Detergent Impacts Membrane Proteins: Atomic-Level Views of Mitochondrial Carriers in Dodecylphosphocholine.

Characterizing the structure of membrane proteins (MPs) generally requires extraction from their native environment, most commonly with detergents. Yet, the physicochemical properties of detergent micelles and lipid bilayers differ markedly and could alter the structural organization of MPs, albeit without general rules. Dodecylphosphocholine (DPC) is the most widely used detergent for MP structure determination by NMR, but the physiological relevance of several prominent structures has been questioned, though indirectly, by other biophysical techniques, e.g., functional/thermostability assay (TSA) and molecular dynamics (MD) simulations. Here, we resolve unambiguously this controversy by probing the functional relevance of three different mitochondrial carriers (MCs) in DPC at the atomic level, using an exhaustive set of solution-NMR experiments, complemented by functional/TSA and MD data. Our results provide atomic-level insight into the structure, substrate interaction and dynamics of the detergent-membrane protein complexes and demonstrates cogently that, while high-resolution NMR signals can be obtained for MCs in DPC, they systematically correspond to nonfunctional states

Recognition of methylated DNA through methyl-CpG binding domain proteins

DNA methylation is a key regulatory control route in epigenetics, involving gene silencing and chromosome inactivation. It has been recognized that methyl-CpG binding domain (MBD) proteins play an important role in interpreting the genetic information encoded by methylated DNA (mDNA). Although the function of MBD proteins has attracted considerable attention and is well characterized, the mechanism underlying mDNA recognition by MBD proteins is still poorly understood. In this article, we demonstrate that the methyl-CpG dinucleotides are recognized at the MBD–mDNA interface by two MBD arginines through an interplay of hydrogen bonding and cation-π interaction. Through molecular dynamics and quantum-chemistry calculations we investigate the methyl-cytosine recognition process and demonstrate that methylation enhances MBD–mDNA binding by increasing the hydrophobic interfacial area and by strengthening the interaction between mDNA and MBD proteins. Free-energy perturbation calculations also show that methylation yields favorable contribution to the binding free energy for MBD–mDNA complex

NMR structure and ion channel activity of the p7 protein from hepatitis C virus

The small membrane protein p7 of hepatitis C virus forms oligomers and exhibits ion channel activity essential for virus infectivity. These viroporin features render p7 an attractive target for antiviral drug development. In this study, p7 from strain HCV-J (genotype 1b) was chemically synthesized and purified for ion channel activity measurements and structure analyses. p7 forms cation-selective ion channels in planar lipid bilayers and at the single-channel level by the patch clamp technique. Ion channel activity was shown to be inhibited by hexamethylene amiloride but not by amantadine. Circular dichroism analyses revealed that the structure of p7 is mainly α-helical, irrespective of the membrane mimetic medium (e.g. lysolipids, detergents, or organic solvent/water mixtures). The secondary structure elements of the monomeric form of p7 were determined by 1H and 13C NMR in trifluoroethanol/water mixtures. Molecular dynamics simulations in a model membrane were combined synergistically with structural data obtained from NMR experiments. This approach allowed us to determine the secondary structure elements of p7, which significantly differ from predictions, and to propose a three-dimensional model of the monomeric form of p7 associated with the phospholipid bilayer. These studies revealed the presence of a turn connecting an unexpected N-terminal α-helix to the first transmembrane helix, TM1, and a long cytosolic loop bearing the dibasic motif and connecting TM1 to TM2. These results provide the first detailed experimental structural framework for a better understanding of p7 processing, oligomerization, and ion channel gating mechanism.Instituto de Física La Plat

Gating of a pH-Sensitive K 2P Potassium Channel by an Electrostatic Effect of Basic Sensor Residues on the Selectivity Filter

Abstract K + channels share common selectivity characteristics but exhibit a wide diversity in how they are gated open. Leak K 2P K + channels TASK-2, TALK-1 and TALK-2 are gated open by extracellular alkalinization. The mechanism for this alkalinizationdependent gating has been proposed to be the neutralization of the side chain of a single arginine (lysine in TALK-2) residue near the pore of TASK-2, which occurs with the unusual pK a of 8.0. We now corroborate this hypothesis by transplanting the TASK-2 extracellular pH (pH o ) sensor in the background of a pH o -insensitive TASK-3 channel, which leads to the restitution of pH o -gating. Using a concatenated channel approach, we also demonstrate that for TASK-2 to open, pH o sensors must be neutralized in each of the two subunits forming these dimeric channels with no apparent cross-talk between the sensors. These results are consistent with adaptive biasing force analysis of K + permeation using a model selectivity filter in wild-type and mutated channels. The underlying free-energy profiles confirm that either a doubly or a singly charged pH o sensor is sufficient to abolish ion flow. Atomic detail of the associated mechanism reveals that, rather than a collapse of the pore, as proposed for other K 2P channels gated at the selectivity filter, an increased height of the energetic barriers for ion translocation accounts for channel blockade at acid pH o . Our data, therefore, strongly suggest that a cycle of protonation/deprotonation of pH o -sensing arginine 224 side chain gates the TASK-2 channel by electrostatically tuning the conformational stability of its selectivity filter