Molecular Dynamics at Constant Pressure: Allowing the System to Control Volume Fluctuations via a “Shell” Particle

Theory and Monte Carlo simulation of the ideal gas with shell particles in the canonical, isothermal-isobaric, grand canonical, and Gibbs ensembles

Theories of small systems play an important role in the fundamental understanding of finite size effects in statistical mechanics, as well as the validation of molecular simulation results as no computer can simulate fluids in the thermodynamic limit. Previously, a shell particle was included in the isothermal-isobaric ensemble in order to resolve an ambiguity in the resulting partition function. The shell particle removed either redundant volume states or redundant translational degrees of freedom of the system and yielded quantitative differences from traditional simulations in this ensemble. In this work, we investigate the effect of including a shell particle in the canonical, grand canonical, and Gibbs ensembles. For systems comprised of a pure component ideal gas, analytical expressions for various thermodynamic properties are obtained. We also derive the Metropolis Monte Carlo simulation acceptance criteria for these ensembles with shell particles, and the results of the simulations of an ideal gas are in excellent agreement with the theoretical predictions. The system size dependence of various important ensemble averages is also analyzed.

Novel barostat implementation for molecular dynamics

We propose a novel implementation of the extended-dynamics equations for isothermal-isobaric ensemble in molecular dynamics, as the Martyna-Tobias-Klein thermostat and barostat. This method is suitable for systems with constraints and the Verlet-family integrators. Instead of iterations or the Trotter-expansion-based methods, both velocities and box sizes (scaling of bond lengths) are predicted. The algorithm begins with force calculation, requiring neither quarter nor half time steps, and necessitating iterations only inside SHAKE. Several tests demonstrate that the quality is comparable to other implementations. It is found that the formula relating the extended barostat mass to the characteristic time of volume fluctuations is inaccurate for condensed systems, which has consequences for the parameter setup. Emphasis is also put on the verification of the precise isothermal-isobaric ensemble and finite-size effects.

Marine phycocompound screening reveals a potential source of novel senotherapeutics

Cells undergo a controlled and systematic cycle of growth, replication and death. However, the integrity of this process gradually declines, leading to accumulation of senescent cells, a major hallmark of biological ageing. Dietary algae, particularly marine algae, have been long reported to exert anti-ageing benefits as cosmeceuticals and nutraceuticals with limited understanding of the molecular mechanisms underlying their activity. In this study, we have incorporated 1,202 previously reported bioactive small phycocompounds and subjected them to cheminformatic queries to assess these interactions. In-silico ADMET, 2-phase docking, metabolic pathway interaction and molecular dynamics simulations reveal multiple marine phycocompounds to have safe and effective senolytic potentials. We employed a novel deep convolutional neural network driven screening approach to identify (2R*, 3S*, 6R*, 7S*, 10R*, 13R*)-7,13-Dihydroxy-2,6-cyclo-1(9),14-xenicadiene-18,19-dial derived from Dilophus Fasciola, Laurendecumenyne A from Laurencia decumbens and 4-Bromo-3-ethyl-9-[(2E)-2-penten-4-yn-1-yl]-2,8-dioxabicyclo[5.2.1]decan-6-ol from Laurencia sp. to be potent inhibitors of multiple target senescent-cell anti-apoptotic pathway proteins. We simulated the best overall target inhibitors, specific protein inhibitors and molecular pathway regulators with each target protein and found stable interactions with minimum deviations (mean RMSD = 0.17 ± 0.01 nm) and gyrations (mean Rg = 1.64 ± 0.16 nm) of the simulated protein-compound complexes. Finally, molecular mechanics calculation suggests potent (mean ΔG = -69.56 ± 27.19 kCal/mol) and frequent hydrophobic interactions between the top performing marine phycocompounds and target proteins.

In silico and saturation transfer difference NMR approaches to unravel the binding mode of an andrographolide derivative to K-Ras oncoprotein

Background: Andrographolide and its benzylidene derivatives, SRJ09 and SRJ23, potentially bind oncogenic K-Ras to exert anticancer activity. Their molecular interactions with K-Ras oncoproteins that lead to effective biological activity are of major interest. Methods & results: In silico docking and molecular dynamics simulation were performed using Glide and Desmond, respectively; while saturation transfer difference NMR was performed using GDP-bound K-RasG12V. SRJ23 was found to bind strongly and selectively to K-RasG12V, by anchoring to a binding pocket (namely p2) principally via hydrogen bond and hydrophobic interactions. The saturation transfer difference NMR analysis revealed the proximity of protons of functional moieties in SRJ23 to K-RasG12V, suggesting positive binding. Conclusion: SRJ23 binds strongly and interacts stably with K-RasG12V to exhibit its inhibitory activity.

A Novel Approach to Atomistic Molecular Dynamics Simulation of Phenolic Resins Using Symthons

Materials science is beginning to adopt computational simulation to eliminate laboratory trial and error campaigns—much like the pharmaceutical industry of 40 years ago. To further computational materials discovery, new methodology must be developed that enables rapid and accurate testing on accessible computational hardware. To this end, the authors utilise a novel methodology concept of intermediate molecules as a starting point, for which they propose the term ‘symthon’ (The term ‘Symthon’ is being used as a simulation equivalent of the synthon, popularised by Dr Stuart Warren in ‘Organic Synthesis: The Disconnection Approach’, OUP: Oxford, 1983.) rather than conventional monomers. The use of symthons eliminates the initial monomer bonding phase, reducing the number of iterations required in the simulation, thereby reducing the runtime. A novel approach to molecular dynamics, with an NVT (Canonical) ensemble and variable unit cell geometry, was used to generate structures with differing physical and thermal properties. Additional script methods were designed and tested, which enabled a high degree of cure in all sampled structures. This simulation has been trialled on large-scale atomistic models of phenolic resins, based on a range of stoichiometric ratios of formaldehyde and phenol. Density and glass transition temperature values were produced, and found to be in good agreement with empirical data and other simulated values in the literature. The runtime of the simulation was a key consideration in script design; cured models can be produced in under 24 h on modest hardware. The use of symthons has been shown as a viable methodology to reduce simulation runtime whilst generating accurate models.

Unintended specificity of an engineered ligand-binding protein facilitated by unpredicted plasticity of the protein fold

Attempts to create novel ligand-binding proteins often focus on formation of a binding pocket with shape complementarity against the desired ligand (particularly for compounds that lack distinct polar moieties). Although designed proteins often exhibit binding of the desired ligand, in some cases they display unintended recognition behavior. One such designed protein, that was originally intended to bind tetrahydrocannabinol (THC), was found instead to display binding of 25-hydroxy-cholecalciferol (25-D3) and was subjected to biochemical characterization, further selections for enhanced 25-D3 binding affinity and crystallographic analyses. The deviation in specificity is due in part to unexpected altertion of its conformation, corresponding to a significant change of the orientation of an α-helix and an equally large movement of a loop, both of which flank the designed ligand-binding pocket. Those changes led to engineered protein constructs that exhibit significantly more contacts and complementarity towards the 25-D3 ligand than the initial designed protein had been predicted to form towards its intended THC ligand. Molecular dynamics simulations imply that the initial computationally designed mutations may contribute to the movement of the helix. These analyses collectively indicate that accurate prediction and control of backbone dynamics conformation, through a combination of improved conformational sampling and/or de novo structure design, represents a key area of further development for the design and optimization of engineered ligand-binding proteins.

New Method for Constant- NPT Molecular Dynamics

The well-established molecular dynamics simulation methods for constant- NPT ensemble systems such as the Andersen-Nosé-Hoover method and their variants may alter the dynamic properties of the molecules under consideration, because their equations of motion are modified by the coupling with thermostat or barostat. To circumvent this artifact, we propose a new molecular dynamics simulation algorithm, by which only the molecules near the wall of the simulation box are coupled to the thermostat and barostat and the molecules of interest placed in the inner part of the simulation box remain intact. We test the efficiency of our algorithm in attaining the target temperature and pressure and the conformity of the calculated equilibrium and dynamic properties to those of a constant- NPT ensemble system.

An ab initio molecular dynamics method for cocrystal prediction: validation of the approach

Cocrystal formation prediction by ab initio molecular dynamics and validation based on the experimental results of 145 coformers for six drugs.

Frutapin, a lectin from Artocarpus incisa (breadfruit): cloning, expression and molecular insights

Artocarpus incisa (breadfruit) seeds contain three different lectins (Frutalin, Frutapin (FTP) and Frutackin) with distinct carbohydrate specificities. The most abundant lectin is Frutalin, an α-D-galactose-specific carbohydrate-binding glycoprotein with antitumour properties and potential for tumour biomarker discovery as already reported. FTP is the second most abundant, but proved difficult to purify with very low yields and contamination with Frutalin frustrating its characterization. Here, we report for the first time high-level production and isolation of biologically active recombinant FTP in Escherichia coli BL21, optimizing conditions with the best set yielding >40 mg/l culture of soluble active FTP. The minimal concentration for agglutination of red blood cells was 62.5 µg/ml of FTP, a process effectively inhibited by mannose. Apo-FTP, FTP–mannose and FTP–glucose crystals were obtained, and they diffracted X-rays to a resolution of 1.58 (P2₁2₁2₁), 1.70 (P3₁21) and 1.60 (P3₁21) Å respectively. The best solution showed four monomers per asymmetric unit. Molecular dynamics (MD) simulation suggested that FTP displays higher affinity for mannose than glucose. Cell studies revealed that FTP was non-cytotoxic to cultured mouse fibroblast 3T3 cells below 0.5 mg/ml and was also capable of stimulating cell migration at 50 µg/ml. In conclusion, our optimized expression system allowed high amounts of correctly folded soluble FTP to be isolated. This recombinant bioactive lectin will now be tested in future studies for therapeutic potential; for example in wound healing and tissue regeneration.