NMR refinement and peptide folding using the GROMACS software

Efficient Simulations of Solvent Asymmetry Across Lipid Membranes Using Flat-Bottom Restraints

The routinely employed periodic boundary conditions complicate molecular simulations of physiologically relevant asymmetric lipid membranes together with their distinct solvent environments. Therefore, separating the extracellular fluid from its cytosolic counterpart has often been performed using a costly double-bilayer setup. Here, we demonstrate that the lipid membrane and solvent asymmetry can be efficiently modeled with a single lipid bilayer by applying an inverted flat-bottom potential to ions and other solute molecules, thereby restraining them to only interact with the relevant leaflet. We carefully optimized the parameters of the suggested method so that the results obtained using the flat-bottom and double-bilayer approaches become mutually indistinguishable. Then, we apply the flat-bottom approach to lipid bilayers with various compositions and solvent environments, covering ions and cationic peptides to validate the approach in a realistic use case. We also discuss the possible limitations of the method as well as its computational efficiency and provide a step-by-step guide on how to set up such simulations in a straightforward manner.

Extent of N-Terminus Folding of Semenogelin 1 Cleavage Product Determines Tendency to Amyloid Formation

It is known that four peptide fragments of predominant protein in human semen Semenogelin 1 (SEM1) (SEM1(86-107), SEM1(68-107), SEM1(49-107) and SEM1(45-107)) are involved in fertilization and amyloid formation processes. In this work, the structure and dynamic behavior of SEM1(45-107) and SEM1(49-107) peptides and their N-domains were described. According to ThT fluorescence spectroscopy data, it was shown that the amyloid formation of SEM1(45-107) starts immediately after purification, which is not observed for SEM1(49-107). Seeing that the peptide amino acid sequence of SEM1(45-107) differs from SEM1(49-107) only by the presence of four additional amino acid residues in the N domain, these domains of both peptides were obtained via solid-phase synthesis and the difference in their dynamics and structure was investigated. SEM1(45-67) and SEM1(49-67) showed no principal difference in dynamic behavior in water solution. Furthermore, we obtained mostly disordered structures of SEM1(45-67) and SEM1(49-67). However, SEM1(45-67) contains a helix (E58-K60) and helix-like (S49-Q51) fragments. These helical fragments may rearrange into β-strands during amyloid formation process. Thus, the difference in full-length peptides' (SEM1(45-107) and SEM1(49-107)) amyloid-forming behavior may be explained by the presence of a structured helix at the SEM1(45-107) N-terminus, which contributes to an increased rate of amyloid formation.

Can molecular dynamics be used to simulate biomolecular recognition?

There are many problems in biochemistry that are difficult to study experimentally. Simulation methods are appealing due to direct availability of atomic coordinates as a function of time. However, direct molecular simulations are challenged by the size of systems and the time scales needed to describe relevant motions. In theory, enhanced sampling algorithms can help to overcome some of the limitations of molecular simulations. Here, we discuss a problem in biochemistry that offers a significant challenge for enhanced sampling methods and that could, therefore, serve as a benchmark for comparing approaches that use machine learning to find suitable collective variables. In particular, we study the transitions LacI undergoes upon moving between being non-specifically and specifically bound to DNA. Many degrees of freedom change during this transition and that the transition does not occur reversibly in simulations if only a subset of these degrees of freedom are biased. We also explain why this problem is so important to biologists and the transformative impact that a simulation of it would have on the understanding of DNA regulation.

Application of molecular dynamics simulation in self-assembled cancer nanomedicine

Self-assembled nanomedicine holds great potential in cancer theragnostic. The structures and dynamics of nanomedicine can be affected by a variety of non-covalent interactions, so it is essential to ensure the self-assembly process at atomic level. Molecular dynamics (MD) simulation is a key technology to link microcosm and macroscale. Along with the rapid development of computational power and simulation methods, scientists could simulate the specific process of intermolecular interactions. Thus, some experimental observations could be explained at microscopic level and the nanomedicine synthesis process would have traces to follow. This review not only outlines the concept, basic principle, and the parameter setting of MD simulation, but also highlights the recent progress in MD simulation for self-assembled cancer nanomedicine. In addition, the physicochemical parameters of self-assembly structure and interaction between various assembled molecules under MD simulation are also discussed. Therefore, this review will help advanced and novice researchers to quickly zoom in on fundamental information and gather some thought-provoking ideas to advance this subfield of self-assembled cancer nanomedicine.

Core packing of well-defined X-ray and NMR structures is the same

Numerous studies have investigated the differences and similarities between protein structures determined by solution NMR spectroscopy and those determined by X-ray crystallography. A fundamental question is whether any observed differences are due to differing methodologies or to differences in the behavior of proteins in solution versus in the crystalline state. Here, we compare the properties of the hydrophobic cores of high-resolution protein crystal structures and those in NMR structures, determined using increasing numbers and types of restraints. Prior studies have reported that many NMR structures have denser cores compared with those of high-resolution X-ray crystal structures. Our current work investigates this result in more detail and finds that these NMR structures tend to violate basic features of protein stereochemistry, such as small non-bonded atomic overlaps and few Ramachandran and sidechain dihedral angle outliers. We find that NMR structures solved with more restraints, and which do not significantly violate stereochemistry, have hydrophobic cores that have a similar size and packing fraction as their counterparts determined by X-ray crystallography at high resolution. These results lead us to conclude that, at least regarding the core packing properties, high-quality structures determined by NMR and X-ray crystallography are the same, and the differences reported earlier are most likely a consequence of methodology, rather than fundamental differences between the protein in the two different environments.