Luca Maragliano

String method in collective variables: Minimum free energy paths and isocommittor surfaces

A computational technique is proposed which combines the string method with a sampling technique to determine minimum free energy paths. The technique only requires to compute the mean force and another conditional expectation locally along the string, and therefore can be applied even if the number of collective variables kept in the free energy calculation is large. This is in contrast with other free energy sampling techniques which aim at mapping the full free energy landscape and whose cost increases exponentially with the number of collective variables kept in the free energy. Provided that the number of collective variables is large enough, the new technique captures the mechanism of transition in that it allows to determine the committor function for the reaction and, in particular, the transition state region. The new technique is illustrated on the example of alanine dipeptide, in which we compute the minimum free energy path for the isomerization transition using either two or four dihedral angles as collective variables. It is shown that the mechanism of transition can be captured using the four dihedral angles, but it cannot be captured using only two of them.

Mapping the Network of Pathways of CO Diffusion in Myoglobin

The pathways of diffusion of a CO molecule inside a myoglobin protein and toward the solvent are investigated. Specifically, the three-dimensional potential of mean force (PMF or free energy) of the CO molecule position inside the protein is calculated by using the single-sweep method in concert with fully resolved atomistic simulations in explicit solvent. The results are interpreted under the assumption that the diffusion of the ligand can be modeled as a navigation on the PMF in which the ligand hops between the PMF local minima following the minimum free energy paths (MFEPs) with rates set by the free energy barriers that need to be crossed. Here, all the local minima of the PMF, the MFEPs, and the barriers along them are calculated. The positions of the local minima are in good agreement with all the known binding cavities inside the protein, which indicates that these cavities may indeed serve as dynamical traps inside the protein and thereby influence the binding process. In addition, the MFEPs connecting the local PMF minima show a complicated network of possible pathways of exit of the dissociated CO starting from the primary docking site, in which the histidine gate is the closest exit from the binding site for the ligand but it is not the only possible one.

Single-sweep methods for free energy calculations

A simple, efficient, and accurate method is proposed to map multidimensional free energy landscapes. The method combines the temperature-accelerated molecular dynamics (TAMD) proposed in [L. Maragliano and E. Vanden-Eijnden, Chem. Phys. Lett. 426, 168 (2006)] with a variational reconstruction method using radial-basis functions for the representation of the free energy. TAMD is used to rapidly sweep through the important regions of the free energy landscape and to compute the gradient of the free energy locally at points in these regions. The variational method is then used to reconstruct the free energy globally from the mean force at these points. The algorithmic aspects of the single-sweep method are explained in detail, and the method is tested on simple examples and used to compute the free energy of the solvated alanine dipeptide in two and four dihedral angles.

Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer.

An extremely scalable computational strategy is described for calculations of the potential of mean force (PMF) in multidimensions on massively distributed supercomputers. The approach involves coupling thousands of umbrella sampling (US) simulation windows distributed to cover the space of order parameters with a Hamiltonian molecular dynamics replica-exchange (H-REMD) algorithm to enhance the sampling of each simulation. In the present application, US/H-REMD is carried out in a two-dimensional (2D) space and exchanges are attempted alternatively along the two axes corresponding to the two order parameters. The US/H-REMD strategy is implemented on the basis of parallel/parallel multiple copy protocol at the MPI level, and therefore can fully exploit computing power of large-scale supercomputers. Here the novel technique is illustrated using the leadership supercomputer IBM Blue Gene/P with an application to a typical biomolecular calculation of general interest, namely the binding of calcium ions to the small protein Calbindin D9k. The free energy landscape associated with two order parameters, the distance between the ion and its binding pocket and the root-mean-square deviation (rmsd) of the binding pocket relative the crystal structure, was calculated using the US/H-REMD method. The results are then used to estimate the absolute binding free energy of calcium ion to Calbindin D9k. The tests demonstrate that the 2D US/H-REMD scheme greatly accelerates the configurational sampling of the binding pocket, thereby improving the convergence of the potential of mean force calculation.

Direct imaging of DNA fibers: the visage of double helix.

Direct imaging becomes important when the knowledge at few/single molecule level is requested and where the diffraction does not allow to get structural and functional information. Here we report on the direct imaging of double stranded (ds) λ-DNA in the A conformation, obtained by combining a novel sample preparation method based on super hydrophobic DNA molecules self-aggregation process with transmission electron microscopy (TEM). The experimental breakthrough is the production of robust and highly ordered paired DNA nanofibers that allowed its direct TEM imaging and the double helix structure revealing.

Intermediate state trapping of a voltage sensor.

Voltage sensor domains (VSDs) regulate ion channels and enzymes by undergoing conformational changes depending on membrane electrical signals. The molecular mechanisms underlying the VSD transitions are not fully understood. Here, we show that some mutations of I241 in the S1 segment of the Shaker Kv channel positively shift the voltage dependence of the VSD movement and alter the functional coupling between VSD and pore domains. Among the I241 mutants, I241W immobilized the VSD movement during activation and deactivation, approximately halfway between the resting and active states, and drastically shifted the voltage activation of the ionic conductance. This phenotype, which is consistent with a stabilization of an intermediate VSD conformation by the I241W mutation, was diminished by the charge-conserving R2K mutation but not by the charge-neutralizing R2Q mutation. Interestingly, most of these effects were reproduced by the F244W mutation located one helical turn above I241. Electrophysiology recordings using nonnatural indole derivatives ruled out the involvement of cation-Π interactions for the effects of the Trp inserted at positions I241 and F244 on the channel’s conductance, but showed that the indole nitrogen was important for the I241W phenotype. Insight into the molecular mechanisms responsible for the stabilization of the intermediate state were investigated by creating in silico the mutations I241W, I241W/R2K, and F244W in intermediate conformations obtained from a computational VSD transition pathway determined using the string method. The experimental results and computational analysis suggest that the phenotype of I241W may originate in the formation of a hydrogen bond between the indole nitrogen atom and the backbone carbonyl of R2. This work provides new information on intermediate states in voltage-gated ion channels with an approach that produces minimum chemical perturbation.

Comparison between mean forces and swarms-of-trajectories string methods

The original formulation of the string method in collective variable space is compared with a recent variant called string method with swarms-of-trajectories. The assumptions made in the original method are revisited and the significance of the minimum free energy path (MFEP) is discussed in the context of reactive events. These assumptions are compared to those made in the string method with swarms-of-trajectories, and shown to be equivalent in a certain regime: in particular an expression for the path identified by the swarms-of-trajectories method is given and shown to be closely related to the MFEP. Finally, the algorithmic aspects of both methods are compared.

Regulation of neural gene transcription by optogenetic inhibition of the RE1-silencing transcription factor.

Significance Repressor element 1-silencing transcription factor (REST) is a transcriptional repressor that regulates nervous system development. Normally expressed at low levels by mature neurons, REST is up-regulated in various brain pathologies. Using a light-sensitive domain from the oat plant (Avena sativa), we engineered novel optogenetic proteins that inhibit REST activity when illuminated by blue light, thus obtaining the spatial and temporal control of the transcription of REST target genes. This approach may have an impact in the development of new therapies for all the diseases in which REST is dysregulated, such as epilepsy, ischemia, and cancers of various origin, as such therapies may counteract the long-term changes in gene expression that take place in the context of the pathological brain. Optogenetics provides new ways to activate gene transcription; however, no attempts have been made as yet to modulate mammalian transcription factors. We report the light-mediated regulation of the repressor element 1 (RE1)-silencing transcription factor (REST), a master regulator of neural genes. To tune REST activity, we selected two protein domains that impair REST-DNA binding or recruitment of the cofactor mSin3a. Computational modeling guided the fusion of the inhibitory domains to the light-sensitive Avena sativa light–oxygen–voltage-sensing (LOV) 2-phototrophin 1 (AsLOV2). By expressing AsLOV2 chimeras in Neuro2a cells, we achieved light-dependent modulation of REST target genes that was associated with an improved neural differentiation. In primary neurons, light-mediated REST inhibition increased Na+-channel 1.2 and brain-derived neurotrophic factor transcription and boosted Na+ currents and neuronal firing. This optogenetic approach allows the coordinated expression of a cluster of genes impinging on neuronal activity, providing a tool for studying neuronal physiology and correcting gene expression changes taking place in brain diseases.

A Novel Topology of Proline-rich Transmembrane Protein 2 (PRRT2) HINTS FOR AN INTRACELLULAR FUNCTION AT THE SYNAPSE

Proline-rich transmembrane protein 2 (PRRT2) has been identified as the single causative gene for a group of paroxysmal syndromes of infancy, including epilepsy, paroxysmal movement disorders, and migraine. On the basis of topology predictions, PRRT2 has been assigned to the recently characterized family of Dispanins, whose members share the two-transmembrane domain topology with a large N terminus and short C terminus oriented toward the outside of the cell. Because PRRT2 plays a role at the synapse, it is important to confirm the exact orientation of its N and C termini with respect to the plasma membrane to get clues regarding its possible function. Using a combination of different experimental approaches, including live immunolabeling, immunogold electron microscopy, surface biotinylation and computational modeling, we demonstrate a novel topology for this protein. PRRT2 is a type II transmembrane protein in which only the second hydrophobic segment spans the plasma membrane, whereas the first one is associated with the internal surface of the membrane and forms a helix-loop-helix structure without crossing it. Most importantly, the large proline-rich N-terminal domain is not exposed to the extracellular space but is localized intracellularly, and only the short C terminus is extracellular (Ncyt/Cexo topology). Accordingly, we show that PRRT2 interacts with the Src homology 3 domain-bearing protein Intersectin 1, an intracellular protein involved in synaptic vesicle cycling. These findings will contribute to the clarification of the role of PRRT2 at the synapse and the understanding of pathogenic mechanisms on the basis of PRRT2-related neurological disorders.

Identification and Expression of Acetylcholinesterase in Octopus vulgaris Arm Development and Regeneration: a Conserved Role for ACHE?

Acetylcholinesterase (ACHE) is a glycoprotein with a key role in terminating synaptic transmission in cholinergic neurons of both vertebrates and invertebrates. ACHE is also involved in the regulation of cell growth and morphogenesis during embryogenesis and regeneration acting through its non-cholinergic sites. The mollusk Octopus vulgaris provides a powerful model for investigating the mechanisms underlying tissue morphogenesis due to its high regenerative power. Here, we performed a comparative investigation of arm morphogenesis during adult arm regeneration and embryonic arm development which may provide insights on the conserved ACHE pathways. In this study, we cloned and characterized O. vulgaris ACHE, finding a single highly conserved ACHE hydrophobic variant, characterized by prototypical catalytic sites and a putative consensus region for a glycosylphosphatidylinositol (GPI)-anchor attachment at the COOH-terminus. We then show that its expression level is correlated to the stage of morphogenesis in both adult and embryonic arm. In particular, ACHE is localized in typical neuronal sites when adult-like arm morphology is established and in differentiating cell locations during the early stages of arm morphogenesis. This possibility is also supported by the presence in the ACHE sequence and model structure of both cholinergic and non-cholinergic sites. This study provides insights into ACHE conserved roles during processes of arm morphogenesis. In addition, our modeling study offers a solid basis for predicting the interaction of the ACHE domains with pharmacological blockers for in vivo investigations. We therefore suggest ACHE as a target for the regulation of tissue morphogenesis.