Martini 3 Coarse-Grained Model for the Cofactors Involved in Photosynthesis.

As a critical step in advancing the simulation of photosynthetic complexes, we present the Martini 3 coarse-grained (CG) models of key cofactors associated with light harvesting (LHCII) proteins and the photosystem II (PSII) core complex. Our work focuses on the parametrization of beta-carotene, plastoquinone/quinol, violaxanthin, lutein, neoxanthin, chlorophyll A, chlorophyll B, and heme. We derived the CG parameters to match the all-atom reference simulations, while structural and thermodynamic properties of the cofactors were compared to experimental values when available. To further assess the reliability of the parameterization, we tested the behavior of these cofactors within their physiological environments, specifically in a lipid bilayer and bound to photosynthetic complexes. The results demonstrate that our CG models maintain the essential features required for realistic simulations. This work lays the groundwork for detailed simulations of the PSII-LHCII super-complex, providing a robust parameter set for future studies.

Molecular basis of PIP2-dependent conformational switching of phosphorylated CD44 in binding FERM.

Association of the cellular adhesive protein CD44 and the N-terminal (FERM) domain of cytoskeleton adaptors is critical for cell proliferation, migration, and signaling. Phosphorylation of the cytoplasmic domain (CTD) of CD44 acts as an important regulator of the protein association, but the structural transformation and dynamics mechanism remain enigmatic. In this study, extensive coarse-grained simulations were employed to explore the molecular details in the formation of CD44-FERM complex under S291 and S325 phosphorylation, a modification path known to exert reciprocal effects on the protein association. We find that phosphorylation of S291 inhibits complexation by causing the CTD of CD44 to adopt a more closed structure. In contrast, S325 phosphorylation liberates the CD44-CTD from the membrane surface and promotes the linkage with FERM. The phosphorylation-driven transformation is found to occur in a PIP2-dependent manner, with PIP2 effecting the relative stability of the closed and open conformation, and a replacement of PIP2 by POPS greatly abrogates this effect. The revealed interdependent regulation mechanism by phosphorylation and PIP2 in the association of CD44 and FERM further strengthens our understanding of the molecular basis of cellular signaling and migration.

Adaption to glucose limitation is modulated by the pleotropic regulator CcpA, independent of selection pressure strength.

BACKGROUND: A central theme in (micro)biology is understanding the molecular basis of fitness i.e. which strategies are successful under which conditions; how do organisms implement such strategies at the molecular level; and which constraints shape the trade-offs between alternative strategies. Highly standardized microbial laboratory evolution experiments are ideally suited to approach these questions. For example, prolonged chemostats provide a constant environment in which the growth rate can be set, and the adaptive process of the organism to such environment can be subsequently characterized. RESULTS: We performed parallel laboratory evolution of Lactococcus lactis in chemostats varying the quantitative value of the selective pressure by imposing two different growth rates. A mutation in one specific amino acid residue of the global transcriptional regulator of carbon metabolism, CcpA, was selected in all of the evolution experiments performed. We subsequently showed that this mutation confers predictable fitness improvements at other glucose-limited growth rates as well. In silico protein structural analysis of wild type and evolved CcpA, as well as biochemical and phenotypic assays, provided the underpinning molecular mechanisms that resulted in the specific reprogramming favored in constant environments. CONCLUSION: This study provides a comprehensive understanding of a case of microbial evolution and hints at the wide dynamic range that a single fitness-enhancing mutation may display. It demonstrates how the modulation of a pleiotropic regulator can be used by cells to improve one trait while simultaneously work around other limiting constraints, by fine-tuning the expression of a wide range of cellular processes.

Altered secondary structure of Dynorphin A associates with loss of opioid signalling and NMDA-mediated excitotoxicity in SCA23.

Spinocerebellar ataxia type 23 (SCA23) is caused by missense mutations in prodynorphin, encoding the precursor protein for the opioid neuropeptides α-neoendorphin, Dynorphin (Dyn) A and Dyn B, leading to neurotoxic elevated mutant Dyn A levels. Dyn A acts on opioid receptors to reduce pain in the spinal cord, but its cerebellar function remains largely unknown. Increased concentration of or prolonged exposure to Dyn A is neurotoxic and these deleterious effects are very likely caused by an N-methyl-d-aspartate-mediated non-opioid mechanism as Dyn A peptides were shown to bind NMDA receptors and potentiate their glutamate-evoked currents. In the present study, we investigated the cellular mechanisms underlying SCA23-mutant Dyn A neurotoxicity. We show that SCA23 mutations in the Dyn A-coding region disrupted peptide secondary structure leading to a loss of the N-terminal α-helix associated with decreased κ-opioid receptor affinity. Additionally, the altered secondary structure led to increased peptide stability of R6W and R9C Dyn A, as these peptides showed marked degradation resistance, which coincided with decreased peptide solubility. Notably, L5S Dyn A displayed increased degradation and no aggregation. R6W and wt Dyn A peptides were most toxic to primary cerebellar neurons. For R6W Dyn A, this is likely because of a switch from opioid to NMDA- receptor signalling, while for wt Dyn A, this switch was not observed. We propose that the pathology of SCA23 results from converging mechanisms of loss of opioid-mediated neuroprotection and NMDA-mediated excitotoxicity.

Role of Charge and Hydrophobicity in Liprotide Formation: A Molecular Dynamics Study with Experimental Constraints.

Bovine α-lactalbumin (aLA) and oleate (OA) form a complex that has been intensively studied for its tumoricidal activity. Small-angle X-ray scattering (SAXS) has revealed that this complex consists of a lipid core surrounded by partially unfolded protein. We call this type of complex a liprotide. Little is known of the molecular interactions between OA and aLA, and no technique has so far provided any high-resolution structure of a liprotide. Here we have used coarse-grained (CG) molecular dynamics (MD) simulations, isothermal titration calorimetry (ITC) and SAXS to investigate the interactions between aLA and OA during the process of liprotide formation. With ITC we found that the strongest enthalpic interactions occurred at a molar ratio of 12.0±1.4:1 OA/aLA. Liprotides formed between OA and aLA at several OA/aLA ratios in silico were stable both in CG and in all-atom simulations. From the simulated structures we calculated SAXS spectra that show good agreement with experimentally measured patterns of matching liprotides. The simulations showed that aLA assumes a molten globular (MG) state, exposing several hydrophobic patches involved in interactions with OA. Initial binding of aLA to OA occurs in an area of aLA in which a high amount of positive charge is located, and only later do hydrophobic interactions become important. The results reveal how unfolding of aLA to expose hydrophobic residues is important for complex formation between aLA and OA. Our findings suggest a general mechanism for liprotide formation and might explain the ability of a large number of proteins to form liprotides with OA.

The structural basis for peptide selection by the transport receptor OppA.

Oligopeptide-binding protein A (OppA) from Lactococcus lactis binds peptides of an exceptionally wide range of lengths (4-35 residues), with no apparent sequence preference. Here, we present the crystal structures of OppA in the open- and closed-liganded conformations. The structures directly explain the protein's phenomenal promiscuity. A huge cavity allows binding of very long peptides, and a lack of constraints for the position of the N and C termini of the ligand is compatible with binding of peptides with varying lengths. Unexpectedly, the peptide's amino-acid composition (but not the exact sequence) appears to have a function in selection, with a preference for proline-rich peptides containing at least one isoleucine. These properties can be related to the physiology of the organism: L. lactis is auxotrophic for branched chain amino acids and favours proline-rich caseins as a source of amino acids. We propose a new mechanism for peptide selection based on amino-acid composition rather than sequence.

Rhodopsin forms a dimer with cytoplasmic helix 8 contacts in native membranes.

G protein-coupled receptors form dimers and higher-order oligomers in membranes, but the precise mode of receptor-receptor interaction remains unknown. To probe the intradimeric proximity of helix 8 (H8), we conducted chemical cross-linking of endogenous cysteines in rhodopsin in disk membranes. We identified a Cys316-Cys316 cross-link using partial proteolysis and liquid chromatography with mass spectrometry. These results show that a symmetric dimer interface mediated by H1 and H8 contacts is present in native membranes.

Immobilization of the plug domain inside the SecY channel allows unrestricted protein translocation.

The SecYEG complex forms a protein-conducting channel in the inner membrane of Escherichia coli to support the translocation of secretory proteins in their unfolded state. The SecY channel is closed at the periplasmic face of the membrane by a small re-entrance loop that connects transmembrane segment 1 with 2b. This helical domain 2a is termed the plug domain. By the introduction of pairs of cysteines and crosslinkers, the plug domain was immobilized inside the channel and connected to transmembrane segment 10. Translocation was inhibited to various degrees depending on the position and crosslinker spacer length. With one of the crosslinked mutants translocation occurred unrestricted. Biochemical characterization of this mutant as well as molecular dynamics simulations suggest that only a limited movement of the plug domain suffices for translocation.

Curvature-dependent elastic properties of liquid-ordered domains result in inverted domain sorting on uniaxially compressed vesicles.

Using a coarse-grained molecular model we study the spatial distribution of lipid domains on a 20-nm-sized vesicle. The lipid mixture laterally phase separates into a raftlike, liquid-ordered (l(o)) phase and a liquid-disordered phase. As we uniaxially compress the mixed vesicle keeping the enclosed volume constant, we impart tension onto the membrane. The vesicle adopts a barrel shape, which is composed of two flat contact zones and a curved edge. The l(o) domain, which exhibits a higher bending rigidity, segregates to the highly curved edge. This inverted domain sorting switches to normal domain sorting, where the l(o) domain prefers the flat contact zone, when we release the contents of the vesicle. We rationalize this domain sorting by a pronounced reduction of the bending rigidity and area compressibility of the l(o) phase upon bending.

Sphere-to-rod transitions of nonionic surfactant micelles in aqueous solution modeled by molecular dynamics simulations.

Control of the size and agglomeration of micellar systems is important for pharmaceutical applications such as drug delivery. Although shape-related transitions in surfactant solutions are studied experimentally, their molecular mechanisms are still not well understood. In this study, we use coarse-grained molecular dynamics simulations to describe micellar assemblies of pentaethylene glycol monododecyl ether (C(12)E(5)) in aqueous solution at different concentrations. The obtained size and aggregation numbers of the aggregates formed are in very good agreement with the available experimental data. Importantly, increase of the concentration leads to a second critical micelle concentration where a transition to rod-like aggregates is observed. This transition is quantified in terms of shape anisotropy, together with a detailed structural analysis of the micelles as a function of aggregation number.

Role of lipids in spheroidal high density lipoproteins.

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.

Thermostable D-amino acid decarboxylases derived from Thermotoga maritima diaminopimelate decarboxylase.

Diaminopimelate decarboxylases (DAPDCs) are highly selective enzymes that catalyze the common final step in different lysine biosynthetic pathways, i.e. the conversion of meso-diaminopimelate (DAP) to L-lysine. We examined the modification of the substrate specificity of the thermostable decarboxylase from Thermotoga maritima with the aim to introduce activity with 2-aminopimelic acid (2-APA) since its decarboxylation leads to 6-aminocaproic acid (6-ACA), a building block for the synthesis of nylon-6. Structure-based mutagenesis of the distal carboxylate binding site resulted in a set of enzyme variants with new activities toward different D-amino acids. One of the mutants (E315T) had lost most of its activity toward DAP and primarily acted as a 2-APA decarboxylase. We next used computational modeling to explain the observed shift in catalytic activities of the mutants. The results suggest that predictive computational protocols can support the redesign of the catalytic properties of this class of decarboxylating PLP-dependent enzymes.

Bottom-up fabrication of a proteasome-nanopore that unravels and processes single proteins.

The precise assembly and engineering of molecular machines capable of handling biomolecules play crucial roles in most single-molecule methods. In this work we use components from all three domains of life to fabricate an integrated multiprotein complex that controls the unfolding and threading of individual proteins across a nanopore. This 900 kDa multicomponent device was made in two steps. First, we designed a stable and low-noise ß-barrel nanopore sensor by linking the transmembrane region of bacterial protective antigen to a mammalian proteasome activator. An archaeal 20S proteasome was then built into the artificial nanopore to control the unfolding and linearized transport of proteins across the nanopore. This multicomponent molecular machine opens the door to two approaches in single-molecule protein analysis, in which selected substrate proteins are unfolded, fed to into the proteasomal chamber and then addressed either as fragmented peptides or intact polypeptides.

Solvent-exposed tails as prestalk transition states for membrane fusion at low hydration.

Membrane fusion is a key step in intracellular trafficking and viral infection. The underlying molecular mechanism is poorly understood. We have used molecular dynamics simulations in conjunction with a coarse grained model to study early metastable and transition states during the fusion of two planar palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers separated by five waters per lipid in the cis leaflets at zero tension. This system mimics the contact area between two vesicles with large diameters compared to the membrane thickness at conditions where fusion may start in the core of the contact area. At elevated temperatures, the two proximal leaflets become connected via multiple lipid molecules and form a stalklike structure. At room temperature, this structure has a free energy of 3k(B)T and is separated from the unconnected state by a significant free energy barrier of 20k(B)T. Stalk formation is initiated by the establishment of a localized hydrophobic contact between the bilayers. This contact is either formed by two partially splayed lipids or a single fully splayed one leading to the formation of a (metastable) splayed lipid bond intermediate. These findings indicate that, for low hydration, early membrane fusion kinetics is not determined by the stalk energy but by the energy of prestalk transition states involving solvent-exposed lipid tails.

Simulations of the c-subunit of ATP-synthase reveal helix rearrangements.

The c-subunit of the enzyme, ATP synthase couples the proton movement through the a-subunit with its own rotation and subsequent rotation of the F1 ring to drive ATP synthesis. Here, we perform mus time-scale coarse-grained molecular dynamics simulations of the c-subunit to characterize its structure and dynamics. Two different helix-helix interfaces, albeit with similar interfacial characteristics, are sampled in the simulations. The helix-2 of the c-subunit monomer rotates around the axis of helix-1 bringing about a change in the interface. Previous models have also proposed such a change in the helix interface but postulated that helix-2 swivels around its own axis. Such large-scale changes in helix packing motifs have not been observed before. The helix-swirling persists even in the c-subunit ring but the dynamics is much slower. The cooperative behavior in the ring appears to stabilize a conformation less-populated in the monomer. Analyzing the stability of the c-subunit ring, it was found that six lipid molecules are necessary to fill the central cavity of the ring. These lipid molecules were not aligned with the surrounding bilayer but protruded towards the periplasmic side. The characterization of the monomer and ring presented in this work sheds light into the structural dynamics of the c-subunit and its functional relevance.

Alternative mechanisms for the interaction of the cell-penetrating peptides penetratin and the TAT peptide with lipid bilayers.

Cell-penetrating peptides (CPPs) have recently attracted much interest due to their apparent ability to penetrate cell membranes in an energy-independent manner. Here molecular-dynamics simulation techniques were used to study the interaction of two CPPs: penetratin and the TAT peptide with 1,2-Dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) phospolipid bilayers shed light on alternative mechanisms by which these peptides might cross biological membranes. In contrast to previous simulation studies of charged peptides interacting with lipid bilayers, no spontaneous formation of transmembrane pores was observed. Instead, the simulations suggest that the peptides may enter the cell by micropinocytosis, whereby the peptides induce curvature in the membrane, ultimately leading to the formation of small vesicles within the cell that encapsulate the peptides. Specifically, multiple peptides were observed to induce large deformations in the lipid bilayer that persisted throughout the timescale of the simulations (hundreds of nanoseconds). Pore formation could be induced in simulations in which an external potential was used to pull a single penetratin or TAT peptide into the membrane. With the use of umbrella-sampling techniques, the free energy of inserting a single penetratin peptide into a DPPC bilayer was estimated to be approximately 75 kJmol(-1), which suggests that the spontaneous penetration of single peptides would require a timescale of at least seconds to minutes. This work also illustrates the extent to which the results of such simulations can depend on the initial conditions, the extent of equilibration, the size of the system, and the conditions under which the simulations are performed. The implications of this with respect to the current systems and to simulations of membrane-peptide interactions in general are discussed.

Areas of monounsaturated diacylphosphatidylcholines.

We have studied the structural properties of monounsaturated diacylphosphatidylcholine lipid bilayers (i.e., diCn:1PC, where n = 14, 16, 18, 20, 22, and 24 is the number of acyl chain carbons). High-resolution x-ray scattering data were analyzed in conjunction with contrast-varied neutron scattering data using a technique we recently developed. Analyses of the data show that the manner by which bilayer thickness increases with increasing n in monounsaturated diacylphosphatidylcholines is dependent on the double bond's position. For commonly available monounsaturated diacylphosphatidylcholines, this results in the nonlinear behavior of both bilayer thickness and lipid area, whereas for diC18:1PC bilayers, lipid area assumes a maximum value. It is worthwhile to note that compared to previous data, our results indicate that lipid areas are smaller by approximately 10%. This observation highlights the need to revisit lipid areas, as they are often used in comparisons with molecular dynamics simulations. Moreover, simulators are encouraged to compare their results not only to x-ray scattering data, but to neutron data as well.

Toroidal pores formed by antimicrobial peptides show significant disorder.

A large variety of antimicrobial peptides have been shown to act, at least in vitro, by poration of the lipid membrane. The nanometre size of these pores, however, complicates their structural characterization by experimental techniques. Here we use molecular dynamics simulations, to study the interaction of a specific class of antimicrobial peptides, melittin, with a dipalmitoylphosphatidylcholine bilayer in atomic detail. We show that transmembrane pores spontaneously form above a critical peptide to lipid ratio. The lipid molecules bend inwards to form a toroidally shaped pore but with only one or two peptides lining the pore. This is in strong contrast to the traditional models of toroidal pores in which the peptides are assumed to adopt a transmembrane orientation. We find that peptide aggregation, either prior or after binding to the membrane surface, is a prerequisite to pore formation. The presence of a stable helical secondary structure of the peptide, however is not. Furthermore, results obtained with modified peptides point to the importance of electrostatic interactions in the poration process. Removing the charges of the basic amino-acid residues of melittin prevents pore formation. It was also found that in the absence of counter ions pores not only form more rapidly but lead to membrane rupture. The rupture process occurs via a novel recursive poration pathway, which we coin the Droste mechanism.

A coarse-grained model for polyethylene oxide and polyethylene glycol: conformation and hydrodynamics.

A coarse-grained (CG) model for polyethylene oxide (PEO) and polyethylene glycol (PEG) developed within the framework of the MARTINI CG force field (FF) using the distributions of bonds, angles, and dihedrals from the CHARMM all-atom FF is presented. Densities of neat low molecular weight PEO agree with experiment, and the radius of gyration R(g) = 19.1 A +/- 0.7 for 76-mers of PEO (M(w) approximately 3400), in excellent agreement with neutron scattering results for an equal sized PEG. Simulations of 9, 18, 27, 36, 44, 67, 76, 90, 112, 135, and 158-mers of the CG PEO (442 < M(w) < 6998) at low concentration in water show the experimentally observed transition from ideal chain to real chain behavior at 1600 < M(w) < 2000, in excellent agreement with the dependence of experimentally observed hydrodynamic radii of PEG. Hydrodynamic radii of PEO calculated from diffusion coefficients of the higher M(w) PEO also agree well with experiment. R(g) calculated from both all-atom and CG simulations of PEO76 at 21 and 148 mg/cm(3) are found to be nearly equal. This lack of concentration dependence implies that apparent R(g) from scattering experiments at high concentration should not be taken to be the chain dimension. Simulations of PEO grafted to a nonadsorbing surface yield a mushroom to brush transition that is well described by the Alexander-de Gennes formalism.

In silico study of full-length amyloid beta 1-42 tri- and penta-oligomers in solution.

Amyloid oligomers are considered to play causal roles in the pathogenesis of amyloid-related degenerative diseases including Alzheimer's disease. Using MD simulation techniques, we explored the contributions of the different structural elements of trimeric and pentameric full-length Abeta1-42 aggregates in solution to their stability and conformational dynamics. We found that our models are stable at a temperature of 310 K, and converge toward an interdigitated side-chain packing for intermolecular contacts within the two beta-sheet regions of the aggregates: beta1 (residues 18-26) and beta2 (residues 31-42). MD simulations reveal that the beta-strand twist is a characteristic element of Abeta-aggregates, permitting a compact, interdigitated packing of side chains from neighboring beta-sheets. The beta2 portion formed a tightly organized beta-helix, whereas the beta1 portion did not show such a firm structural organization, although it maintained its beta-sheet conformation. Our simulations indicate that the hydrophobic core comprising the beta2 portion of the aggregate is a crucial stabilizing element in the Abeta aggregation process. On the basis of these structure-stability findings, the beta2 portion emerges as an optimal target for further antiamyloid drug design.