Search

The secondary transporters of the resistance-nodulation-cell division (RND) superfamily mediate multidrug resistance in Gram-negative bacteria like Pseudomonas aeruginosa. Among these RND transporters, MexB, MexF, and MexY, with partly overlapping specificities, have been implicated in pathogenicity. Only the structure of the former has been resolved experimentally, which together with the lack of data about the functional dynamics of the full set of transporters, limited a systematic investigation of the molecular determinants defining their peculiar and shared features. In a previous work (Ramaswamy et al., Front. Microbiol., 2018, 9, 1144), we compared at an atomistic level the two main putative recognition sites (named access and deep binding pockets) of MexB and MexY. In this work, we expand the comparison by performing extended molecular dynamics (MD) simulations of these transporters and the pathologically relevant transporter MexF. We employed a more realistic model of the inner phospholipid membrane of P. aeruginosa and more accurate force-fields. To elucidate structure/dynamics-activity relationships we performed physico-chemical analyses and mapped the binding propensities of several organic probes on all transporters. Our data revealed the presence, also in MexF, of a few multifunctional sites at locations equivalent to the access and deep binding pockets detected in MexB. Furthermore, we report for the first time about the multidrug binding abilities of two out of five gates of the channels deputed to peripheral (early) recognition of substrates. Overall, our findings help to define a common “recognition topology” characterizing Mex transporters, which can be exploited to optimize transport and inhibition propensities of antimicrobial compounds.

The increasing interest in the molecular mechanism of the binding of different agonists and antagonists to β2-adrenergic receptor (β2AR) inactive and active states has led us to investigate protein–ligand interactions using molecular docking calculations. To perform this study, the 3.2 Å X-ray crystal structure of the active conformation of human β2AR in the complex with the endogenous agonist adrenaline has been used as a template for investigating the binding of two exogenous catecholamines to this adrenergic receptor. Here, we show the derivation of L-DOPA and Droxidopa OPLS all atom (AA) force field (FF) parameters via quantum mechanical (QM) calculations, molecular dynamics (MD) simulations in aqueous solutions of the two catecholamines and the molecular docking of both ligands into rigid and flexible β2AR models. We observe that both ligands share with adrenaline similar experimentally observed binding anchor sites, which are constituted by Asp113/Asn312 and Ser203/Ser204/Ser207 side chains. Moreover, both L-DOPA and Droxidopa molecules exhibit binding affinities comparable to that predicted for adrenaline, which is in good agreement with previous experimental and computational results. L-DOPA and Droxidopa OPLS AA FFs have also been tested by performing MD simulations of these ligands docked into β2AR proteins embedded in lipid membranes. Both hydrogen bonds and hydrophobic interaction networks observed over the 1 μs MD simulation are comparable with those derived from molecular docking calculations and MD simulations performed with the CHARMM FF.

LCAT is activated by apoA-I to form cholesteryl ester. We combined two structures, phospholipase A2 (PLA2) that hydrolyzes the ester bond at the sn-2 position of oxidized (short) acyl chains of phospholipid, and bacteriophage tubulin PhuZ, as C- and N-terminal templates, respectively, to create a novel homology model for human LCAT. The juxtaposition of multiple structural motifs matching experimental data is compelling evidence for the general correctness of many features of the model: i) The N-terminal 10 residues of the model, required for LCAT activity, extend the hydrophobic binding trough for the sn-2 chain 15–20 Å relative to PLA2. ii) The topography of the trough places the ester bond of the sn-2 chain less than 5 Å from the hydroxyl of the catalytic nucleophile, S181. iii) A β-hairpin resembling a lipase lid separates S181 from solvent. iv) S181 interacts with three functionally critical residues: E149, that regulates sn-2 chain specificity, and K128 and R147, whose mutations cause LCAT deficiency. Because the model provides a novel explanation for the complicated thermodynamic problem of the transfer of hydrophobic substrates from HDL to the catalytic triad of LCAT, it is an important step toward understanding the antiatherogenic role of HDL in reverse cholesterol transport.

Since spheroidal HDL particles (sHDL) are highly dynamic, molecular dynamics (MD) simulations are useful for obtaining structural models. Here we use MD to simulate sHDL with stoichiometries of reconstituted and circulating particles. The hydrophobic effect during simulations rapidly remodels discoidal HDL containing mixed lipids to sHDL containing a cholesteryl ester/triglyceride (CE/TG) core. We compare the results of simulations of previously characterized reconstituted sHDL particles containing two or three apoA-I created in the absence of phospholipid transfer protein (PLTP) with simulations of circulating human HDL containing two or three apoA-I without apoA-II. We find that circulating sHDL compared with reconstituted sHDL with the same number of apoA-I per particle contain approximately equal volumes of core lipid but significantly less surface lipid monolayers. We conclude that in vitro reconstituted sHDL particles contain kinetically trapped excess phospholipid and are less than ideal models for circulating sHDL particles. In the circulation, phospholipid transfer via PLTP decreases the ratio of phospholipid to apolipoprotein for all sHDL particles. Further, sHDL containing two or three apoA-I adapt to changes in surface area by condensation of common conformational motifs. These results represent an important step toward resolving the complicated issue of the protein and lipid stoichiometry of circulating HDL.

HDL is a population of apoA-I-containing particles inversely correlated with heart disease. Because HDL is a soft form of matter deformable by thermal fluctuations, structure determination has been difficult. Here, we compare the recently published crystal structure of lipid-free (Δ185-243)apoA-I with apoA-I structure from models and molecular dynamics (MD) simulations of discoidal HDL. These analyses validate four of our previous structural findings for apoA-I: i) a baseline double belt diameter of 105 Å ii) central α helixes with an 11/3 pitch; iii) a “presentation tunnel” gap between pairwise helix 5 repeats hypothesized to move acyl chains and unesterified cholesterol from the lipid bilayer to the active sites of LCAT; and iv) interchain salt bridges hypothesized to stabilize the LL5/5 chain registry. These analyses are also consistent with our finding that multiple salt bridge-forming residues in the N-terminus of apoA-I render that conserved domain “sticky.” Additionally, our crystal MD comparisons led to two new hypotheses: i) the interchain leucine-zippers previously reported between the pair-wise helix 5 repeats drive lipid-free apoA-I registration; ii) lipidation induces rotations of helix 5 to allow formation of interchain salt bridges, creating the LCAT presentation tunnel and “zip-locking” apoA-I into its full LL5/5 registration.

We study the structure and dynamics of spherical high density lipoprotein (HDL) particles through coarse-grained multi-microsecond molecular dynamics simulations. We simulate both a lipid droplet without the apolipoprotein A-I (apoA-I) and the full HDL particle including two apoA-I molecules surrounding the lipid compartment. The present models are the first ones among computational studies where the size and lipid composition of HDL are realistic, corresponding to human serum HDL. We focus on the role of lipids in HDL structure and dynamics. Particular attention is paid to the assembly of lipids and the influence of lipid-protein interactions on HDL properties. We find that the properties of lipids depend significantly on their location in the particle (core, intermediate region, surface). Unlike the hydrophobic core, the intermediate and surface regions are characterized by prominent conformational lipid order. Yet, not only the conformations but also the dynamics of lipids are found to be distinctly different in the different regions of HDL, highlighting the importance of dynamics in considering the functionalization of HDL. The structure of the lipid droplet close to the HDL-water interface is altered by the presence of apoA-Is, with most prominent changes being observed for cholesterol and polar lipids. For cholesterol, slow trafficking between the surface layer and the regimes underneath is observed. The lipid-protein interactions are strongest for cholesterol, in particular its interaction with hydrophobic residues of apoA-I. Our results reveal that not only hydrophobicity but also conformational entropy of the molecules are the driving forces in the formation of HDL structure. The results provide the first detailed structural model for HDL and its dynamics with and without apoA-I, and indicate how the interplay and competition between entropy and detailed interactions may be used in nanoparticle and drug design through self-assembly.