Effect of coarse graining in water models for the study of kinetics and mechanisms of clathrate hydrates nucleation and growth.

Clathrate hydrates are crystalline inclusion compounds wherein a water framework encages small guest atoms/molecules within its cavities. Among the others, methane clathrates are the largest fossil fuel resource still available. They can also be used to safely transport gases and can also form spontaneously under suitable conditions plugging pipelines. Understanding the crystallization mechanism is very important, and given the impossibility of experimentally identifying the atomistic path, simulations played an important role in this field. Given the large computational cost of these simulations, in addition to all-atom force fields, scientists considered coarse-grained water models. Here, we have investigated the effect of coarse-graining, as implemented in the water model mW, on the crystallization characteristics of methane clathrate in comparison with the all-atom TIP4P force field. Our analyses revealed that although the characteristics directly depending on the energetics of the water models are well reproduced, dynamical properties are off by the orders of magnitude. Being crystallization a non-equilibrium process, the altered kinetics of the process results in different characteristics of crystalline nuclei. Both TIP4P and mW water models produce methane clathrate nuclei with some amount of the less stable (in the given thermodynamic conditions) structure II phase and an excess of pentagonal dodecahedral cages over the tetrakaidecahedral ones regarding the ideal ratio in structure I. However, the dependence of this excess on the methane concentration in solution is higher with the former water model, whereas with the latter, the methane concentration in solution dependence is reduced and within the statistical error.

Jarzynski equality on work and free energy: Crystal indentation as a case study.

Mathematical relations concerning particle systems require knowledge of the applicability conditions to become physically relevant and not merely formal. We illustrate this fact through the analysis of the Jarzynski equality (JE), whose derivation for Hamiltonian systems suggests that the equilibrium free-energy variations can be computational or experimentally determined in almost any kind of non-equilibrium processes. This apparent generality is surprising in a mechanical theory. Analytically, we show that the quantity called "work" in the Hamiltonian derivation of the JE is neither a thermodynamic quantity nor mechanical work, except in special circumstances to be singularly assessed. Through molecular dynamics simulations of elastic and plastic deformations induced via nano-indentation of crystalline surfaces that fall within the formal framework of the JE, we illustrate that the JE cannot be verified and that the results of this verification are process dependent.

Adiabatic motion and statistical mechanics via mass-zero constrained dynamics.

In recent work [Coretti et al., J. Chem. Phys., 2018, 149, 191102], a new algorithm to solve numerically the dynamics of the shell model for polarization was presented. The approach, broadly applicable to systems involving adiabatically separated dynamical variables, employs constrained molecular dynamics to strictly enforce the condition that the force on the fast degrees of freedom, modeled as having zero mass, is null at each time step. The algorithm is symplectic and fully time reversible, and results in stable and efficient propagation. In this paper we complete the discussion of the mechanics of mass-zero constrained dynamics by showing how to adapt it to problems where the fast degrees of freedom must satisfy additional conditions. This extension includes, in particular, the important case of first principles molecular dynamics. We then consider the statistical mechanics of the mass-zero constrained dynamical system demonstrating that the marginal probability sampled by the dynamics in the physical phase space recovers the form of the Born-Oppenheimer probability density. The effectiveness of the approach and the favorable scaling of the algorithm with system size are illustrated in test calculations of solid Na via orbital-free density functional dynamics.

Quantum Trajectories for the Dynamics in the Exact Factorization Framework: A Proof-of-Principle Test.

In the framework of the exact factorization of the time-dependent electron-nuclear wave function, we investigate the possibility of solving the nuclear time-dependent Schrödinger equation based on trajectories. The nuclear equation is separated in a Hamilton-Jacobi equation for the phase of the wave function, and a continuity equation for its (squared) modulus. For illustrative adiabatic and nonadiabatic one-dimensional models, we implement a procedure to follow the evolution of the nuclear density along the characteristics of the Hamilton-Jacobi equation. Those characteristics are referred to as quantum trajectories, since they are generated via ordinary differential equations similar to Hamilton's equations, but including the so-called quantum potential, and they can be used to reconstruct exactly the quantum-mechanical nuclear wave function, provided infinite initial conditions are propagated in time.

A General Mechanism for Signal Propagation in the Nicotinic Acetylcholine Receptor Family.

Nicotinic acetylcholine receptors (nAChRs) modulate synaptic activity in the central nervous system. The α7 subtype, in particular, has attracted considerable interest in drug discovery as a target for several conditions, including Alzheimer's disease and schizophrenia. Identifying agonist-induced structural changes underlying nAChR activation is fundamentally important for understanding biological function and rational drug design. Here, extensive equilibrium and nonequilibrium molecular dynamics simulations, enabled by cloud-based high-performance computing, reveal the molecular mechanism by which structural changes induced by agonist unbinding are transmitted within the human α7 nAChR. The simulations reveal the sequence of coupled structural changes involved in driving conformational change responsible for biological function. Comparison with simulations of the α4ß2 nAChR subtype identifies features of the dynamical architecture common to both receptors, suggesting a general structural mechanism for signal propagation in this important family of receptors.

The physics of open systems for the simulation of complex molecular environments in soft matter.

Molecular dynamics (MD) has become one of the most powerful tools of investigation in soft matter. Despite such success, simulations of large molecular environments are mostly run using the approximation of closed systems without the possibility of exchange of matter. Due to the molecular complexity of soft matter systems, an optimal simulation strategy would require the application of concurrent multiscale resolution approaches such that each part of a large system can be considered as an open subsystem at a high resolution embedded in a large coarser reservoir of energy and particles. This paper discusses the current capability and the future perspectives of multiscale adaptive resolution MD methods to satisfy the conceptual principles of open systems and to perform simulations of complex molecular environments in soft matter.

Ab initio accelerated molecular dynamics study of the hydride ligands in the ruthenium complex: Ru(H2)2H2(P(C5H9)3)2.

The dihydrogen complex Ru(H2)2H2(P(C5H9)3)2 has been investigated, via ab initio accelerated molecular dynamics, to elucidate the H ligands dynamics and possible reaction paths for H2/H exchange. We have characterized the free energy landscape associated with the H atoms positional exchange around the Ru centre. From the free energy landscape, we have been able to estimate a barrier of 6 kcal mol-1 for the H2/H exchange process. We have also observed a trihydrogen intermediate as a passing state along some of the possible reaction pathways.

On computing the solubility of molecular systems subject to constraints using the extended Einstein crystal method.

A method to compute solubilities for molecular systems using atomistic simulations, based on an extension of the Einstein crystal method, has recently been presented [Li et al., J. Chem. Phys. 146, 214110 (2017)]. This methodology is particularly appealing to compute solubilities in cases of practical importance including, but not limited to, solutions where the solute is sparingly soluble and molecules of importance for the pharmaceutical industry, which are often characterized by strong polar interactions and slow relaxation time scales. The mathematical derivation of this methodology hinges on a factorization of the partition function which is not necessarily applicable in the case of a system subject to holonomic molecular constraints. We show here that, although the mathematical procedure to derive it is slightly different, essentially the same mathematical relation for calculating the solubility can be safely applied for computing the solubility of systems subject to constraints, which are the majority of the systems used for practical molecular simulations.

Solubility of paracetamol in ethanol by molecular dynamics using the extended Einstein crystal method and experiments.

Li and co-workers [Li et al., J. Chem. Phys. 146, 214110 (2017)] have recently proposed a methodology to compute the solubility of molecular compounds from first principles, using molecular dynamics simulations. We revise and further explore their methodology that was originally applied to naphthalene in water at low concentration. In particular, we compute the solubility of paracetamol in an ethanol solution at ambient conditions. For the simulations, we used a force field that we previously reparameterized to reproduce certain thermodynamic properties of paracetamol but not explicitly its solubility in ethanol. In addition, we have determined the experimental solubility by performing turbidity measurements using a Crystal16 over a range of temperatures. Our work serves a dual purpose: (i) methodologically, we clarify how to compute, with a relatively straightforward procedure, the solubility of molecular compounds and (ii) applying this procedure, we show that the solubility predicted by our force field (0.085 ± 0.014 in mole ratio) is in good agreement with the experimental value obtained from our experiments and those reported in the literature (average 0.0585 ± 0.004), considering typical deviations for predictions from first principle methods. The good agreement between the experimental and the calculated solubility also suggests that the method used to reparameterize the force field can be used as a general strategy to optimize force fields for simulations in solution.

Communication: Constrained molecular dynamics for polarizable models.

A new algorithm to solve numerically the evolution of empirical shell models of polarizable systems is presented. It employs constrained molecular dynamics to satisfy exactly, at each time step, the crucial condition that the gradient of the potential with respect to the shell degrees of freedom is null. The algorithm is efficient, stable, and, contrary to the available alternatives, it is symplectic and time reversible. A proof-of-principle calculation on a polarizable model for NaCl is presented to illustrate its properties in comparison with the current method, which employs a conjugate-gradient procedure to enforce the null gradient condition. The proposed algorithm is applicable to other cases where a minimum condition on a function of an auxiliary set of driven dynamical variables must be satisfied.

Thermal Diffusion in Binary Mixtures: Transient Behavior and Transport Coefficients from Equilibrium and Nonequilibrium Molecular Dynamics.

Equilibrium and nonequilibrium molecular dynamics simulations are combined to compute the full set of coefficients that appear in the phenomenological equations describing thermal transport in a binary mixture subject to a constant thermal gradient. The Dynamical Non-Equilibrium Molecular Dynamics approach (D-NEMD) is employed to obtain the microscopic time evolution of the density and temperature fields, together with that of the mass and energy fluxes. D-NEMD enables one to study not only the steady state, but also the evolution of the fields during the transient that follows the onset of the thermal gradient, up to the establishment of the steady state. This makes it possible to ensure that the system has indeed reached a stationary condition, and to analyze the transient mechanisms and time scales of the mass and energy transport. A local time averaging procedure is applied to each trajectory contributing to the calculation to improve the signal-to-noise ratio in the estimation of the fluxes and to obtain a clear signal with the, relatively limited, statistics available.

On metastability and Markov state models for non-stationary molecular dynamics.

Unlike for systems in equilibrium, a straightforward definition of a metastable set in the non-stationary, non-equilibrium case may only be given case-by-case-and therefore it is not directly useful any more, in particular in cases where the slowest relaxation time scales are comparable to the time scales at which the external field driving the system varies. We generalize the concept of metastability by relying on the theory of coherent sets. A pair of sets A and B is called coherent with respect to the time interval [t1, t2] if (a) most of the trajectories starting in A at t1 end up in B at t2 and (b) most of the trajectories arriving in B at t2 actually started from A at t1. Based on this definition, we can show how to compute coherent sets and then derive finite-time non-stationary Markov state models. We illustrate this concept and its main differences to equilibrium Markov state modeling on simple, one-dimensional examples.

On the properties of a bundle of flexible actin filaments in an optical trap.

We establish the statistical mechanics framework for a bundle of Nf living and uncrosslinked actin filaments in a supercritical solution of free monomers pressing against a mobile wall. The filaments are anchored normally to a fixed planar surface at one of their ends and, because of their limited flexibility, they grow almost parallel to each other. Their growing ends hit a moving obstacle, depicted as a second planar wall, parallel to the previous one and subjected to a harmonic compressive force. The force constant is denoted as the trap strength while the distance between the two walls as the trap length to make contact with the experimental optical trap apparatus. For an ideal solution of reactive filaments and free monomers at fixed free monomer chemical potential µ1, we obtain the general expression for the grand potential from which we derive averages and distributions of relevant physical quantities, namely, the obstacle position, the bundle polymerization force, and the number of filaments in direct contact with the wall. The grafted living filaments are modeled as discrete Wormlike chains, with F-actin persistence length ℓp, subject to discrete contour length variations ±d (the monomer size) to model single monomer (de)polymerization steps. Rigid filaments (ℓp = ∞), either isolated or in bundles, all provide average values of the stalling force in agreement with Hill's predictions Fs (H)=NfkBTln(ρ1/ρ1c)/d, independent of the average trap length. Here ρ1 is the density of free monomers in the solution and ρ1c its critical value at which the filament does not grow nor shrink in the absence of external forces. Flexible filaments (ℓp < ∞) instead, for values of the trap strength suitable to prevent their lateral escape, provide an average bundle force and an average trap length slightly larger than the corresponding rigid cases (few percents). Still the stalling force remains nearly independent on the average trap length, but results from the product of two strongly L-dependent contributions: the fraction of touching filaments ∝〈L〉(O.T.) (2) and the single filament buckling force ∝〈L〉(O.T.) (-2).

A semi-flexible model prediction for the polymerization force exerted by a living F-actin filament on a fixed wall.

We consider a single living semi-flexible filament with persistence length ℓp in chemical equilibrium with a solution of free monomers at fixed monomer chemical potential µ1 and fixed temperature T. While one end of the filament is chemically active with single monomer (de)polymerization steps, the other end is grafted normally to a rigid wall to mimic a rigid network from which the filament under consideration emerges. A second rigid wall, parallel to the grafting wall, is fixed at distance L < < ℓp from the filament seed. In supercritical conditions where monomer density ρ1 is higher than the critical density ρ1c, the filament tends to polymerize and impinges onto the second surface which, in suitable conditions (non-escaping filament regime), stops the filament growth. We first establish the grand-potential Ω(µ1, T, L) of this system treated as an ideal reactive mixture, and derive some general properties, in particular the filament size distribution and the force exerted by the living filament on the obstacle wall. We apply this formalism to the semi-flexible, living, discrete Wormlike chain model with step size d and persistence length ℓp, hitting a hard wall. Explicit properties require the computation of the mean force f̄i(L) exerted by the wall at L and associated potential f̄i(L)=-dWi(L)/dL on a filament of fixed size i. By original Monte-Carlo calculations for few filament lengths in a wide range of compression, we justify the use of the weak bending universal expressions of Gholami et al. [Phys. Rev. E 74, 041803 (2006)] over the whole non-escaping filament regime. For a filament of size i with contour length Lc = (i - 1) d, this universal form is rapidly growing from zero (non-compression state) to the buckling value fb(Lc,ℓp)=π(2)kBTℓp4Lc (2) over a compression range much narrower than the size d of a monomer. Employing this universal form for living filaments, we find that the average force exerted by a living filament on a wall at distance L is in practice L independent and very close to the value of the stalling force Fs (H)=(kBT/d)ln(ρˆ1) predicted by Hill, this expression being strictly valid in the rigid filament limit. The average filament force results from the product of the cumulative size fraction x=x(L,ℓp,ρˆ1), where the filament is in contact with the wall, times the buckling force on a filament of size Lc ≈ L, namely, Fs (H)=xfb(L;ℓp). The observed L independence of Fs (H) implies that x ∝ L(-2) for given (ℓp,ρˆ1) and x∝lnρˆ1 for given (ℓp, L). At fixed (L,ρˆ1), one also has x∝ℓp (-1) which indicates that the rigid filament limit ℓp → ∞ is a singular limit in which an infinite force has zero weight. Finally, we derive the physically relevant threshold for filament escaping in the case of actin filaments.

Maximum probability domains for the analysis of the microscopic structure of liquids.

We introduce the concept of maximum probability domains (MPDs), developed in the context of the analysis of electronic densities, in the study of the microscopic spatial structures of liquids. The idea of locating a particle in a three dimensional region, by determining the domain where the probability of finding that, and only that, particle is maximum, gives an interesting characterization of the local structure of the liquid. The optimization procedure, required for the search of the domain of maximum probability, is carried out by the implementation of the level set method. Results for a couple of case studies are presented, to illustrate the structure of liquid water at ambient conditions and upon increasing pressure from the point of view of MPDs and to compare the information encoded in the solvation shells of sodium in water with, once again, that extracted from the MPDs.

Clathrate structure-type recognition: Application to hydrate nucleation and crystallisation.

For clathrate-hydrate polymorphic structure-type (sI versus sII), geometric recognition criteria have been developed and validated. These are applied to the study of the rich interplay and development of both sI and sII motifs in a variety of hydrate-nucleation events for methane and H2S hydrate studied by direct and enhanced-sampling molecular dynamics (MD) simulations. In the case of nucleation of methane hydrate from enhanced-sampling simulation, we notice that already at the transition state, ∼80% of the enclathrated CH4 molecules are contained in a well-structured (sII) clathrate-like crystallite. For direct MD simulation of nucleation of H2S hydrate, some sI/sII polymorphic diversity was encountered, and it was found that a realistic dissipation of the nucleation energy (in view of non-equilibrium relaxation to either microcanonical (NVE) or isothermal-isobaric (NPT) distributions) is important to determine the relative propensity to form sI versus sII motifs.

Temperature-accelerated molecular dynamics gives insights into globular conformations sampled in the free state of the AC catalytic domain.

The catalytic domain of the adenyl cyclase (AC) toxin from Bordetella pertussis is activated by interaction with calmodulin (CaM), resulting in cAMP overproduction in the infected cell. In the X-ray crystallographic structure of the complex between AC and the C terminal lobe of CaM, the toxin displays a markedly elongated shape. As for the structure of the isolated protein, experimental results support the hypothesis that more globular conformations are sampled, but information at atomic resolution is still lacking. Here, we use temperature-accelerated molecular dynamics (TAMD) simulations to generate putative all-atom models of globular conformations sampled by CaM-free AC. As collective variables, we use centers of mass coordinates of groups of residues selected from the analysis of standard molecular dynamics (MD) simulations. Results show that TAMD allows extended conformational sampling and generates AC conformations that are more globular than in the complexed state. These structures are then refined via energy minimization and further unrestrained MD simulations to optimize inter-domain packing interactions, thus resulting in the identification of a set of hydrogen bonds present in the globular conformations.

Relaxation of a steep density gradient in a simple fluid: comparison between atomistic and continuum modeling.

We compare dynamical nonequilibrium molecular dynamics and continuum simulations of the dynamics of relaxation of a fluid system characterized by a non-uniform density profile. Results match quite well as long as the lengthscale of density nonuniformities are greater than the molecular scale (~10 times the molecular size). In presence of molecular scale features some of the continuum fields (e.g., density and momentum) are in good agreement with atomistic counterparts, but are smoother. On the contrary, other fields, such as the temperature field, present very large difference with respect to reference (atomistic) ones. This is due to the limited accuracy of some of the empirical relations used in continuum models, the equation of state of the fluid in the present example.

Diagnostic and therapeutic utility of transrectal ultrasound in urological office prostatic abscess management: a short report from a single urologic center.

OBJECTIVES: Prostatic abscess (PA) is an infrequent condition in the modern antibiotic era. The everyday use of transrectal ultrasound (TRUS) during diagnostic work-up and the widespread recurrence to prostatic biopsies may lead to an increase of PA diagnosis. In this short report we analyze the patients characteristics and the management of seven recent cases of PA diagnosed in our institution. MATERIALS AND METHODS: The records of 7 patients admitted to our Center for LUTS associated to septic fever or acute urinary retention, was prospectively collected. Suspect of PA was done on digital rectal examination (DRE) and confirmed by TRUS performed after urinary system ultrasound (UUS) evaluation. Patients were admitted to hospital only in case of septic signs. A sovrapubic (SPC) or urethral catheter (UC) was placed depending on symptoms. A TRUS-guided aspiration of PA was performed with patient in lithotomic position, using a 18 gauge two-part needle, side/end fire needle access. Patient was discharged with antibiotic therapy and followed up until complete resolution of the PA and symptoms. RESULTS: Mean age was 62 years (range 24-82). Two patients were diabetics and one was affected by the immunodeficiency acquired syndrome (HIV). In one case, PA was detected after a persistent fever post TRUS guided prostate biopsy. Average prostate volume was 69 ml (range 19-118 ml). DRE was able to diagnose PA only in 2 cases (29%), UUS evaluation in 1 case (14%). All cases were confirmed by TRUS as hypo-anechoic areas with or without internal echoes in all patients. Mean PA dimension was 3.64 cm (range 1.5-8). SPC was placed in 3 cases (43%), UC in 3 patients (43%). Only 1 patient refused catheterization. Side fire needle aspiration was performed in all cases and in combination with end fire access in case of particular location of abscess cavities. Second look was needed in 2 cases (29%). Antibiotics were administered in all cases. The aspirated pus showed a positive culture for Escherichia coli (43%), Klebsiella pneumoniae (29%), Pseudomonas aeruginosa (14%) and Enterococcus faecalis (14%). PA resolution time mean was 9 days (range 3-24). CONCLUSIONS: TRUS evaluation in case of persistent LUTS associated with fever or acute urinary retention is determinant in the diagnosis of PA. Office or institutional management with TRUS needle aspiration is a good option in these cases.

Diagnostic ultrasound-guided excisional testicular biopsy for small (< 1 cm) incidental nodules. A single institution experience.

INTRODUCTION AND OBJECTIVES: The widespread use of scrotal ultrasound (SUS) has led to a marked increase in the number of incidentally detected testicular lesions. A small incidental nodule (STN) has defined as a non palpable (< 10 mm), asymptomatic solid lesion with normal levels of oncological testicular markers. Nowadays the lack of agreement on the topic causes managing problems to andrologists. We present our experience consisting in 8 cases of STN discovered by SUS performed for different clinical indications. MATHERIAL AND METHODS: We retrieved from our ultrasonographic files the clinical information about 717 patients evaluated for andrological problems. Patients with STN underwent to a complete clinical history and physical examination as well as oncological testicular markers measurement and ormonal assessment and then received a diagnostic ultrasound guided excisional biopsy (DEB). Surgical approach was performed through an inguinal incision. Using the coordinates previously obtained from preoperative SUS, STN was localized by intraoperative SUS. The lesion was enucleated and sent to the Pathology department for frozen section examination (FSE). Biopsies of affected testis (TB) were also performed. Post-excision ultrasound has been used to confirm the complete removal of the nodule. Whether pathological findings were benign, testis sparing surgery (TSS) was performed. Immediate radical orchidectomy (IRO) was performed if FSE and TB findings suggested a malignant lesion. RESULTS: STNs were discovered in 8 patients (1,1%). Very small lesions (< 5 mm) were detected in 50% of cases. We performed four IRO and four DEB with consequent TSS. In one case we performed a delayed radical orchidectomy (DRO). At FSE pathologist reported 3 Leydig cell tumor and 3 seminoma and an inflammatory regressive lesion in one case. FSE on TB reported intratesticular neoplasia (TIN) in three cases. In one case nodule wasn't sent to FSE. We observed a concordance between FSE and definitive pathologic report in six cases (75%). CONCLUSIONS: The management of STN is still a challenge for the surgical andrologist. A correct diagnosis has a crucial role in making the best treatment and patients outcome. Ultrasound guided excisional biopsy and the close collaboration with a dedicated pathologist are very useful in reducing errors.