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

Probing Passive Permeation of Tetracycline: Are Simulations Ready for beyond-Rule-of-Five Drug Permeability Calculation?

Passive permeation across lipid membranes is a key determinant of drug bioavailability and efficacy. Accurate computational estimation of drug permeability is essential for rational drug design, yet remains challenging, particularly for ionizable and beyond-Rule-of-Five (bRo5) compounds. In this study, we employ advanced molecular simulations and the inhomogeneous solubility-diffusion model to calculate the effective permeability and elucidate the membrane permeation mechanism of the antibiotic tetracycline (TC), the six hydrogen-bond donors of which violates one of Lipinski’s Rule-of-Five. By integrating the pH-partitioning and Boltzmann-weighted average potential schemes and accounting for both its neutral (TCN) and zwitterionic (TCZ) tautomers, we show that the dominant contribution to the effective permeability of TC arises from TCN, despite its low abundance. This result is attributed to the relatively small permeation barrier of TCN, explaining why the antibiotic exhibits moderate effective permeability even though it predominantly exists in the zwitterionic form at neutral pH. A further systematic investigation reveals that membrane patch size significantly impacts permeability estimates for TC, in contrast to the relative insensitivity observed for three other permeants. This unique sensitivity can be attributed to the hydrogen-bond network formed between TC and its lipid environment, with the smallest 32-POPC patch artificially raising the drug molecule’s permeation barrier and the largest 256-POPC patch exhibiting significant hysteresis that compromises the quality of the one-dimensional free-energy calculation. Overall, our results suggest that while probing the passive permeation of bRo5 drugs by molecular simulations appears increasingly feasible, the protonation and tautomeric states of the permeants, the uncertainty in their microscopic acid dissociation constants, as well as the potential impact of membrane patch-size effects need to be fully considered in permeability predictions for these complex molecules

Probing the formation of a hetero-dimeric membrane transport complex with dual <i>in vitro</i> and <i>in silico</i> mutagenesis

The reversible association of transmembrane helices is a fundamental mechanism in how living cells convey information and respond to physiological events. The cardiac calcium transport regulator phospholamban (PLN) is an example of a single-span transmembrane protein that populates a variety of reversible and competing oligomeric states. PLN primarily forms monomers and pentamers in the membrane, where the PLN pentamer is a storage form and the PLN monomer forms a hetero-dimeric inhibitory complex with SERCA. The binding affinity and free-energy of formation of the SERCA-PLN complex in a membrane have not been determined. As is the case for most transmembrane protein interactions, measuring these quantities experimentally is extremely challenging. In this study, we estimated binding affinities by employing in silico alchemical free-energy calculations for all PLN transmembrane alanine substitutions in a membrane bilayer. The binding affinities were calculated separately for the SERCA-PLN complex, a PLN monomer, and a PLN pentamer and compared to in vitro functional measurements of SERCA regulation by the PLN alanine substitutions. Initially, the changes in SERCA inhibition by PLN alanine substitutions were compared to the changes in free energy for the SERCA-PLN complex formed from the PLN monomer. However, the functional data for the PLN alanine substitutions were better explained by the formation of the SERCA-PLN complex directly from the PLN pentamer. This finding points to an inhibitory mechanism favoring conformational selection of SERCA and the interaction of a PLN pentamer with SERCA for ‘delivery’ of a PLN monomer to the inhibitory site. The implications of these findings suggest that the energetics of helix exchange between homo- and hetero-oligomeric signaling complexes is favored over an intermediate involving a free monomeric helix in the membrane bilayer

Development of Coarse-Grained Lipid Force Fields Based on a Graph Neural NetworkClick to copy article link

Coarse-grained (CG) lipid models enable efficient simulations of large-scale membrane events. However, achieving both speed and atomic-level accuracy remains challenging. Graph neural networks (GNNs) trained on all-atom (AA) simulations can serve as CG force fields, which have demonstrated success in CG simulations of proteins. Herein, we built data sets of AA simulations of DOPC, DOPS, and mixed DOPC/DOPS lipid bilayers and developed the first GNN-based CG lipid models based on the TorchMD-GN architecture. The CG lipid models reproduce the structural correlations of the AA simulations, accelerate the lipid dynamics by 9.4 times, and exhibit some degree of temperature transferability. Moreover, we demonstrate that training CG models on lipid bicelles enhances the performance of models in the lipid self-assembly and vesicle simulations. Our findings indicate that GNN-based CG lipid force fields show promise as a powerful approach for large-scale membrane simulations

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

Elucidation of the noncovalent interactions driving enzyme activity guides branching enzyme engineering for α-glucan modification

Branching enzymes (BEs) confer to α-glucans, the primary energy-storage reservoir in nature, a variety of features, like slow digestion. The full catalytic cycle of BEs can be divided in six steps, namely two covalent catalytic steps involving glycosylation and transglycosylation, and four noncatalytic steps involving substrate binding and transfers (SBTs). Despite the ever-growing wealth of biochemical and structural information on BEs, clear mechanistic insights into SBTs from an industrial-performance perspective are still missing. Here, we report a Rhodothermus profundi BE (RpBE) endowed with twice as much enzymatic activity as the Rhodothermus obamensis BE currently used in industry. Furthermore, we focus on the SBTs for RpBE by means of large-scale computations supported by experiment. Engineering of the crucial positions responsible for the initial substrate-binding step improves enzymatic activity significantly, while offering a possibility to customize product types. In addition, we show that the high-efficiency substrate-transfer steps preceding glycosylation and transglycosylation are the main reason for the remarkable enzymatic activity of RpBE, suggestive of engineering directions for the BE family

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