Full Enzyme Complex Simulation: Interactions in Human Pyruvate Dehydrogenase Complex.

The pyruvate dehydrogenase complex (PDC) is a large macromolecular machine consisting of dozens of interacting enzymes that are connected and regulated by highly flexible domains, also called swinging arms. The overall structure and function of these domains and how they organize the complex function have not been elucidated in detail to date. This lack of structural and dynamic understanding is frequently observed in multidomain enzymatic complexes. Here we present the first full and dynamic structural model of full human PDC (hPDC), including binding of the linking arms to the surrounding E1 and E3 enzymes via their binding domains with variable stoichiometries. All of the linking domains were modeled at atomistic and coarse-grained levels, and the latter was parametrized to reproduce the same properties of those from the atomistic model. The radii of gyration of the wild-type full complex and functional trimeric subunits were in agreement with available experimental data. Furthermore, the E1 and E3 population effect on the overall structure of the full complex was studied. The results indicated that decreasing the number of E1s increases the flexibility of the now nonoccupied arms. Furthermore, their flexibility depends on the presence of other E1s and E3s in the vicinity, even if they are associated with other arms. As one consequence, the radius of gyration decreases with decreasing number of E1s. This effect also provides an indication of the optimal configuration of E1 and E3 on the basis of the assumption that a certain stability of the enymatic cloud is necessary to avoid free metabolic diffusion of intermediates (metabolic channeling). Our approach and results open a window for future enzyme engineering in a more effective way by evaluating the effect of different linker arm lengths, flexibilities, and combinations of mutations on the activity of other complex enzymes that involve flexible domains, including for example processive enzymes.

Investigation of Core Structure and Stability of Human Pyruvate Dehydrogenase Complex: A Coarse-Grained Approach.

The human pyruvate dehydrogenase complex (hPDC) is a large macromolecular machine, and its unique structural and functional properties make it a versatile target for manipulation aiming for the design of new types of artificial multienzyme cascades. However, model-based and hence systematic understanding of the structure-function relationship of this kind of complexes is yet poor. However, with new structure data, modeling techniques, and increasing computation power available, this shortfall is about to cease. Recently, we have built new atomistic models of E2/E3BP, the two subunits of the massive hPDC core. Here, we present developed coarse-grained models of the same, on the basis of which we built and simulated the full core of hPDC, which is so far the first simulation of the catalytic core of any member in the branched-chain α-keto acid dehydrogenase complex family. We explored the stability of two previously proposed substitutional models of hPDC core: 40E2+20E3BP and 48E2+12E3BP. Our molecular dynamics simulations showed a higher stability and sphericity for the second model. Our core's radius of gyration was found to be in strong agreement with previously published experimental dynamic light scattering (DLS) data. Finally, in the direction of our experimental effort to design active minimized complexes, we simulated C-terminal truncated E2/E3BP cores of different lengths, which clearly showed the instability of the core assembly and symmetry due to subunit separations. This correlated very well with the findings of our experimental investigations by analysis of DLS data for variable truncation lengths. The use of polarizable water and an increased total ion concentration did not lead to a substantially different initial stability of the truncated mutants compared to that of standard MARTINI water; however, a different rearrangement behavior of the fragments was observed. The obtained structure models will serve as a basis for full complex simulations in the future, providing the possibility for the model-assisted targeted manipulation of such a complex enzymatic system.

Reengineering of the human pyruvate dehydrogenase complex: from disintegration to highly active agglomerates.

The pyruvate dehydrogenase complex (PDC) plays a central role in cellular metabolism and regulation. As a metabolite-channeling multi-enzyme complex it acts as a complete nanomachine due to its unique geometry and by coupling a cascade of catalytic reactions using 'swinging arms'. Mammalian and specifically human PDC (hPDC) is assembled from multiple copies of E1 and E3 bound to a large E2/E3BP 60-meric core. A less restrictive and smaller catalytic core, which is still active, is highly desired for both fundamental research on channeling mechanisms and also to create a basis for further modification and engineering of new enzyme cascades. Here, we present the first experimental results of the successful disintegration of the E2/E3BP core while retaining its activity. This was achieved by C-terminal α-helixes double truncations (eight residues from E2 and seven residues from E3BP). Disintegration of the hPDC core via double truncations led to the formation of highly active (approximately 70% of wildtype) apparently unordered clusters or agglomerates and inactive non-agglomerated species (hexamer/trimer). After additional deletion of N-terminal 'swinging arms', the aforementioned C-terminal truncations also caused the formation of agglomerates of minimized E2/E3BP complexes. It is likely that these 'swinging arm' regions are not solely responsible for the formation of the large agglomerates.

Human Pyruvate Dehydrogenase Complex E2 and E3BP Core Subunits: New Models and Insights from Molecular Dynamics Simulations.

Targeted manipulation and exploitation of beneficial properties of multienzyme complexes, especially for the design of novel and efficiently structured enzymatic reaction cascades, require a solid model understanding of mechanistic principles governing the structure and functionality of the complexes. This type of system-level and quantitative knowledge has been very scarce thus far. We utilize the human pyruvate dehydrogenase complex (hPDC) as a versatile template to conduct corresponding studies. Here we present new homology models of the core subunits of the hPDC, namely E2 and E3BP, as the first time effort to elucidate the assembly of hPDC core based on molecular dynamic simulation. New models of E2 and E3BP were generated and validated at atomistic level for different properties of the proteins. The results of the wild type dimer simulations showed a strong hydrophobic interaction between the C-terminal and the hydrophobic pocket which is the main driving force in the intertrimer binding and the core self-assembly. On the contrary, the C-terminal truncated versions exhibited a drastic loss of hydrophobic interaction leading to a dimeric separation. This study represents a significant step toward a model-based understanding of structure and function of large multienzyme systems like PDC for developing highly efficient biocatalyst or bioreaction cascades.

Micellar drug nanocarriers and biomembranes: how do they interact?

Pluronic based formulations are among the most successful nanomedicines and block-copolymer micelles including drugs that are undergoing phase I/II studies as anticancer agents. Using coarse-grained models, molecular dynamics simulations of large-scale systems, modeling Pluronic micelles interacting with DPPC lipid bilayers, on the µs timescale have been performed. Simulations show, in agreement with experiments, the release of Pluronic chains from the micelle to the bilayer. This release changes the size of the micelle. Moreover, the presence of drug molecules inside the core of the micelle has a strong influence on this process. The picture emerging from the simulations is that the micelle stability is a result of an interplay of drug-micelle core and block-copolymer-bilayer interactions. The equilibrium size of the drug vector shows a strong dependency on the hydrophobicity of the drug molecules embedded in the core of the micelle. In particular, the radius of the micelle shows an abrupt increase in a very narrow range of drug molecule hydrophobicity.

Theoretical study of binding and permeation of ether-based polymers through interfaces.

We present a molecular dynamics simulation study on the interactions of poly(ethylene oxide) (PEO), poly(propylene oxide) (PPO), and their ABA-type block copolymer, poloxamers, at water/n-heptane and 1,2-dimyristoyl-sn-glycero-3-phospatidycholine (DMPC) lipid bilayer/water interfaces. The partition coefficients in water/1-octanol of the linear polyethers up to three monomers were calculated. The partition coefficients evidenced a higher hydrophobicity of the PPO in comparison to PEO. At the water/n-heptane interface, the polymers tend to adopt elongated conformations in agreement with similar experimental ellipsometry studies of different poloxamers. In the case of the poloxamers at the n-heptane/water interface, the stronger preference of the PPO block for the hydrophobic phase resulted in bottle-brush-type polymer conformations. At lipid bilayer/water interface, the PEO polymers, as expected from their hydrophilic nature, are weakly adsorbed on the surface of the lipid bilayer and locate in the water phase close to the headgroups. The free energy barriers of permeation calculated for short polymer chains suggest a thermodynamics propensity for the water phase that increase with the chain length. The lower affinity of PEO for the hydrophobic interior of the lipid bilayer resulted in the spontaneous expulsion within the simulation time. On the contrary, PPO chains and poloxamers have a longer residence time inside the bilayer, and they tend to concentrate in the tail region of the bilayer near the polar headgroups. In addition, polymers with PPO unit length comparable to the thickness of the hydrophobic region of the bilayer tend to span across the bilayer.

Understanding the interaction of block copolymers with DMPC lipid bilayer using coarse-grained molecular dynamics simulations.

In this paper, we present a computational model of the adsorption and percolation mechanism of poloxamers (poly(ethylene oxide) (PEO) and poly(propylene oxide) (PPO) triblock copolymers) across a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid bilayer. A coarse-grained model was used to cope with the long time scale of the percolation process. The simulations have provided details of the interaction mechanism of Pluronics with lipid bilayer. In particular, the results have shown that polymer chains containing a PPO block with a length comparable to the DMPC bilayer thickness, such as P85, tends to percolate across the lipid bilayer. On the contrary, Pluronics with a shorter PPO chain, such as L64 and F38, insert partially into the membrane with the PPO block part while the PEO blocks remain in water on one side of the lipid bilayer. The percolation of the polymers into the lipid tail groups reduces the membrane thickness and increases the area per lipid. These effects are more evident for P85 than L64 or F38. Our findings are qualitatively in good agreement with published small-angle X-ray scattering experiments that have evidenced a thinning effect of Pluronics on the lipid bilayer as well as the role of the length of the PPO block on the permeation process of the polymer through the lipid bilayer. Our theoretical results complement the experimental data with a detailed structural and dynamic model of poloxamers at the interface and inside the lipid bilayer.

Molecular dynamics simulation study of solvent effects on conformation and dynamics of polyethylene oxide and polypropylene oxide chains in water and in common organic solvents.

In this paper, the conformation and dynamics properties of polyethylene oxide (PEO) and polypropylene oxide (PPO) polymer chains at 298 K have been studied in the melt and at infinite dilution condition in water, methanol, chloroform, carbon tetrachloride, and n-heptane using molecular dynamics simulations. The calculated density of PEO melt with chain lengths of n = 2, 3, 4, 5 and, for PPO, n = 7 are in good agreement with the available experimental data. The conformational properties of PEO and PPO show an increasing gauche preference for the O-C-C-O dihedral in the following order water>methanol>chloroform>carbon tetrachloride = n-heptane. On the contrary, the preference for trans conformation has a maximum in carbon tetrachloride and n-heptane followed in the order by chloroform, methanol, and water. The PEO conformational preferences are in qualitative agreement with results of NMR studies. PEO chains formed different types of hydrogen bonds with polar solvent molecules. In particular, the occurrence of bifurcated hydrogen bonding in chloroform was also observed. Radii of gyration of PEO chains of length larger than n = 9 monomers showed a good agreement with light scattering data in water and in methanol. For the shorter chains the observed deviations are probably due to the enhanced hydrophobic effects caused by the terminal methyl groups. For PEO the fitting of end-to-end distance distributions with the semi-flexible chain model at 298 K provided persistence lengths of 0.375 and 0.387 nm in water and methanol, respectively. Finally, the radius of gyration of Pluronic P85 turned out to be 2.25 ± 0.4 nm at 293 K in water in agreement with experimental data.

Diffusion of 1,2-dimethoxyethane and 1,2-dimethoxypropane through phosphatidycholine bilayers: a molecular dynamics study.

In this paper, a theoretical study of 1,2-dimethoxyethane (DME) and 1,2-dimethoxypropane (DMP) at water/n-heptane and 1,2-dimyristoyl-sn-glycero-3-phospatidycholine (DMPC) lipid bilayer/water interfaces using the umbrella sampling method is reported. Recently proposed GROMOS96/OPLS compatible models for DME and DMP have been used for the simulation studies. The percolation free energy barrier of one DME and DMP molecule from water to n-heptane phase calculated using the umbrella sampling method turned out to be equal to ~18.5 kJ/mol and ~6 kJ/mol, respectively. In the case of the DMPC lipid bilayer, overall free energy barriers of ~20 kJ/mol and ∼12 kJ/mol were obtained for DME and DMP, respectively. The spontaneous diffusion of DME and DMP in the lipid bilayer has also been investigated using unconstrained molecular dynamics simulations at the water/DMPC interface and inside the lipid bilayer. As expected from the estimated percolation barriers, simulation results show that DME, contrary to DMP, spontaneously diffuse into the aqueous solution from the lipid interior. In addition, simulations with multiple DME or DMP molecules at the interface show spontaneous diffusion within 50 ns inside the DMPC layer only for DMP.

Structure and dynamics of 1,2-dimethoxyethane and 1,2-dimethoxypropane in aqueous and non-aqueous solutions: a molecular dynamics study.

Herein, we report a comparative modelling study of 1,2-dimethoxyethane (DME) and 1,2-dimethoxypropane (DMP) at 298 K and 318 K in the liquid state, water mixtures, and at infinite dilution condition in water, methanol, carbon tetrachloride, and n-heptane. Both DME and DMP are united-atom models compatible with GROMOS∕OPLS force fields. Calculated thermodynamic and structural properties of the pure DME and DMP liquids resulted in excellent agreement with the experimental data. In aqueous solutions, densities, diffusion coefficients, and concentration dependent conformers of DME, were in agreement with experimental data. The calculated free energy of solvation (ΔG(hyd)) at 298 K is equal to -22.1 ± 0.8 kJ mol(-1) in good agreement with the experimental value of 20.2 kJ mol(-1). In addition, the free energy of solvation of DME in non-aqueous solvents follows the trend methanol ≈ water < carbon tetrachloride < n-heptane, consistently with the dielectric constant of the solvents. On contrary, the presence of an extra methyl group on chiral carbon makes DMP less soluble than DME in water (ΔG(hyd) = -16.0 ± 1.1 kJ mol(-1)) but more soluble in non-polar solvents as n-heptane. Finally, for the DMP the chiral discrimination of the two enantiomers was calculated as solvation free energy difference of one DMP isomer in the solution of the other. The obtained value of ΔΔG(RS) = -3.7 ± 1.4 kJ mol(-1) indicates a net chiral discrimination of the two enantiomers.