The Rift Valley fever virus (RVFV) is transmitted by infected mosquitoes, causing severe disease in humans and livestock across Africa. We determined the x-ray structure of the RVFV class II fusion protein Gc in its postfusion form and in complex with a glycerophospholipid (GPL) bound in a conserved cavity next to the fusion loop. Site-directed mutagenesis and molecular dynamics simulations further revealed a built-in motif allowing en bloc insertion of the fusion loop into membranes, making few nonpolar side-chain interactions with the aliphatic moiety and multiple polar interactions with lipid head groups upon membrane restructuring. The GPL head-group recognition pocket is conserved in the fusion proteins of other arthropod-borne viruses, such as Zika and chikungunya viruses, which have recently caused major epidemics worldwide.

Molecular dynamics simulations are capable of predicting the permeability of lipid membranes for drug-like solutes, but the calculations have remained prohibitively expensive for high-throughput studies. Here, we analyze simple measures for accelerating potential of mean force (PMF) calculations of membrane permeation, namely, (i) using smaller simulation systems, (ii) simulating multiple solutes per system, and (iii) using shorter cutoffs for the Lennard-Jones interactions. We find that PMFs for membrane permeation are remarkably robust against alterations of such parameters, suggesting that accurate PMF calculations are possible at strongly reduced computational cost. In addition, we evaluated the influence of the definition of the membrane center of mass (COM), used to define the transmembrane reaction coordinate. Membrane-COM definitions based on all lipid atoms lead to artifacts due to undulations and, consequently, to PMFs dependent on membrane size. In contrast, COM definitions based on a cylinder around the solute lead to size-independent PMFs, down to systems of only 16 lipids per monolayer. In summary, compared to popular setups that simulate a single solute in a membrane of 128 lipids with a Lennard-Jones cutoff of 1.2 nm, the measures applied here yield a speedup in sampling by factor of ∼40, without reducing the accuracy of the calculated PMF.

Formate-nitrite transporters (FNTs) facilitate the translocation of monovalent polyatomic anions, such as formate and nitrite, across biological membranes. FNTs are widely distributed among pathogenic bacteria and eukaryotic parasites, but they lack human homologues, making them attractive drug targets. The mechanisms and energetics involved in anion permeation across the FNTs have remained largely unclear. Both, channel and transporter mode of function have been proposed, with strong indication of proton coupling to the permeation process. We combine molecular dynamics simulations, quantum mechanical calculations, and pK a calculations, to compute the energetics of the complete permeation cycle of an FNT. We find that anions as such, are not able to traverse the FNT pore. Instead, anion binding into the pore is energetically coupled to protonation of a centrally located histidine. In turn, the histidine can protonate the permeating anion, thereby enabling its release. Such mechanism can accommodate the functional diversity among the FNTs, as it may facilitate both, export and import of substrates, with or without proton co-transport. The mechanism excludes proton leakage via the Grotthuss mechanism, and it rationalises the selectivity for weak acids.

Steroids have numerous physiological functions associated with cellular signaling or modulation of the lipid membrane structure and dynamics, and as such, they have found broad pharmacological applications. Steroid–membrane interactions are relevant to multiple steps of steroid biosynthesis and action, as steroids are known to interact with neurotransmitter or membrane steroid receptors, and steroids must cross lipid membranes to exert their physiological functions. Therefore, rationalizing steroid function requires understanding of steroid–membrane interactions. We combined molecular dynamics simulations and isothermal titration calorimetry to characterize the conformations and the energetics of partitioning, in addition to the kinetics of flip–flop transitions and membrane exit, of 26 representative steroid compounds in a model lipid membrane. The steroid classes covered in this study include birth control and anabolic drugs, sex and corticosteroid hormones, neuroactive steroids, as well as steroids modulating the lipid membrane structure. We found that the conformational ensembles adopted by different steroids vary greatly, as quantified by their distributions of tilt angles and insertion depths into the membrane, ranging from well-defined steroid conformations with orientations either parallel or normal to the membrane, to wide conformational distributions. Surprisingly, despite their chemical diversity, the membrane/water partition coefficient is similar among most steroids, except for structural steroids such as cholesterol, leading to similar rates for exiting the membrane. By contrast, the rates of steroid flip–flop vary by at least 9 orders of magnitude, revealing that flip–flop is the rate-limiting step during cellular uptake of polar steroids. This study lays the ground for a quantitative understanding of steroid–membrane interactions, and it will hence be of use for studies of steroid biosynthesis and function as well as for the development and usage of steroids in a pharmacological context.

Molecular docking has been extensively applied in virtual screening of small molecule libraries for lead identification and optimization. A necessary prerequisite for successful differentiation between active and non-active ligands is the accurate prediction of their binding affinities in the complex by use of docking scoring functions. However, many studies have shown rather poor correlations between docking scores and experimental binding affinities. Our work aimed to improve this correlation by implementing a multipose binding concept in the docking scoring scheme. Multipose binding, i.e., the property of certain protein-ligand complexes to exhibit different ligand binding modes, has been shown to occur in nature for a variety of molecules. We conducted a high-throughput docking study and implemented multipose binding in the scoring procedure by considering multiple docking solutions in binding affinity prediction. In general, improvement of the agreement between docking scores and experimental data was observed, and this was most pronounced in complexes with large and flexible ligands and high binding affinities. Further developments of the selection criteria for docking solutions for each individual complex are still necessary for a general utilization of the multipose binding concept for accurate binding affinity prediction by molecular docking.