A Practical Guide to All-Atom and Coarse-Grained Molecular Dynamics Simulations Using Amber and Gromacs: A Case Study of Disulfide-Bond Impact on the Intrinsically Disordered Amyloid Beta.

Intrinsically disordered proteins (IDPs) pose challenges to conventional experimental techniques due to their large-scale conformational fluctuations and transient structural elements. This work presents computational methods for studying IDPs at various resolutions using the Amber and Gromacs packages with both all-atom (Amber ff19SB with the OPC water model) and coarse-grained (Martini 3 and SIRAH) approaches. The effectiveness of these methodologies is demonstrated by examining the monomeric form of amyloid-ß (Aß42), an IDP, with and without disulfide bonds at different resolutions. Our results clearly show that the addition of a disulfide bond decreases the ß-content of Aß42; however, it increases the tendency of the monomeric Aß42 to form fibril-like conformations, explaining the various aggregation rates observed in experiments. Moreover, analysis of the monomeric Aß42 compactness, secondary structure content, and comparison between calculated and experimental chemical shifts demonstrates that all three methods provide a reasonable choice to study IDPs; however, coarse-grained approaches may lack some atomistic details, such as secondary structure recognition, due to the simplifications used. In general, this study not only explains the role of disulfide bonds in Aß42 but also provides a step-by-step protocol for setting up, conducting, and analyzing molecular dynamics (MD) simulations, which is adaptable for studying other biomacromolecules, including folded and disordered proteins and peptides.

Scrutinising the Conformational Ensemble of the Intrinsically Mixed-Folded Protein Galectin-3.

Galectin-3 is a protein involved in many intra- and extra-cellular processes. It has been identified as a diagnostic or prognostic biomarker for certain types of heart disease, kidney disease and cancer. Galectin-3 comprises a carbohydrate recognition domain (CRD) and an N-terminal domain (NTD), which is unstructured and contains eight collagen-like Pro-Gly-rich tandem repeats. While the structure of the CRD has been solved using protein crystallography, current knowledge about conformations of full-length galectin-3 is limited. To fill in this knowledge gap, we performed molecular dynamics (MD) simulations of full-length galectin-3. We systematically re-scaled the solute-solvent interactions in the Martini 3 force field to obtain the best possible agreement between available data from SAXS experiments and the ensemble of conformations generated in the MD simulations. The simulation conformations were found to be very diverse, as reflected, e.g., by (i) large fluctuations in the radius of gyration, ranging from about 2 to 5 nm, and (ii) multiple transient contacts made by amino acid residues in the NTD. Consistent with evidence from NMR experiments, contacts between the CRD and NTD were observed to not involve the carbohydrate-binding site on the CRD surface. Contacts within the NTD were found to be made most frequently by aromatic residues. Formation of fuzzy complexes with unspecific stoichiometry was observed to be mediated mostly by the NTD. Taken together, we offer a detailed picture of the conformational ensemble of full-length galectin-3, which will be important for explaining the biological functions of this protein at the molecular level.

Membrane curvature sensing by model biomolecular condensates.

Biomolecular condensates (BCs) are fluid droplets that form in biological cells by liquid-liquid phase separation. Their major components are intrinsically disordered proteins. Vast attention has been given in recent years to BCs inside the cytosol and nucleus. BCs at the cell membrane have not been studied to the same extent so far. However, recent studies provide increasingly more examples of interfaces between BCs and membranes which function as platforms for diverse biomolecular processes. Galectin-3, for example, is known to mediate clathrin-independent endocytosis and has been recently shown to undergo liquid-liquid phase separation, but the function of BCs of galectin-3 in endocytic pit formation is unknown. Here, we use dissipative particle dynamics simulations to study a generic coarse-grained model for BCs interacting with lipid membranes. In analogy to galectin-3, we consider polymers comprising two segments - one of them mediates multivalent attractive interactions between the polymers, and the other one has affinity for association with specific lipid head groups. When these polymers are brought into contact with a multi-component membrane, they spontaneously assemble into droplets and, simultaneously, induce lateral separation of lipids within the membrane. Interestingly, we find that if the membrane is bent, the polymer droplets localize at membrane regions curved inward. Although the polymers have no particular shape or intrinsic curvature, they appear to sense membrane curvature when clustered at the membrane. Our results indicate toward a generic mechanism of membrane curvature sensing by BCs involved in such processes as endocytosis.

Crystal Structure of the ORP8 Lipid Transport ORD Domain: Model of Lipid Transport.

ORPs are lipid-transport proteins belonging to the oxysterol-binding protein family. They facilitate the transfer of lipids between different intracellular membranes, such as the ER and plasma membrane. We have solved the crystal structure of the ORP8 lipid transport domain (ORD8). The ORD8 exhibited a ß-barrel fold composed of anti-parallel ß-strands, with three α-helices replacing ß-strands on one side. This mixed alpha-beta structure was consistent with previously solved structures of ORP2 and ORP3. A large cavity (≈1860 Å3) within the barrel was identified as the lipid-binding site. Although we were not able to obtain a lipid-bound structure, we used computer simulations based on our crystal structure to dock PS and PI4P molecules into the putative lipid-binding site of the ORD8. Comparative experiments between the short ORD8ΔLid (used for crystallography) and the full-length ORD8 (lid containing) revealed the lid's importance for stable lipid binding. Fluorescence assays revealed different transport efficiencies for PS and PI4P, with the lid slowing down transport and stabilizing cargo. Coarse-grained simulations highlighted surface-exposed regions and hydrophobic interactions facilitating lipid bilayer insertion. These findings enhance our comprehension of ORD8, its structure, and lipid transport mechanisms, as well as provide a structural basis for the design of potential inhibitors.