Structural basis for PHDVC5HCHNSD1-C2HRNizp1 interaction: implications for Sotos syndrome.

Sotos syndrome is an overgrowth syndrome caused by mutations within the functional domains ofNSD1 gene coding for NSD1, a multidomain protein regulating chromatin structure and gene expression. In particular, PHDVC5HCHNSD1 tandem domain, composed by a classical (PHDV) and an atypical (C5HCH) plant homeo-domain (PHD) finger, is target of several pathological missense-mutations. PHDVC5HCHNSD1 is also crucial for NSD1-dependent transcriptional regulation and interacts with the C2HR domain of transcriptional repressor Nizp1 (C2HRNizp1)in vitro To get molecular insights into the mechanisms dictating the patho-physiological relevance of the PHD finger tandem domain, we solved its solution structure and provided a structural rationale for the effects of seven Sotos syndrome point-mutations. To investigate PHDVC5HCHNSD1 role as structural platform for multiple interactions, we characterized its binding to histone H3 peptides and to C2HRNizp1 by ITC and NMR. We observed only very weak electrostatic interactions with histone H3 N-terminal tails, conversely we proved specific binding to C2HRNizp1 We solved C2HRNizp1 solution structure and generated a 3D model of the complex, corroborated by site-directed mutagenesis. We suggest a mechanistic scenario where NSD1 interactions with cofactors such as Nizp1 are impaired by PHDVC5HCHNSD1 pathological mutations, thus impacting on the repression of growth-promoting genes, leading to overgrowth conditions.

Torsion and curvature of FtsZ filaments.

FtsZ filaments participate in bacterial cell division, but it is still not clear how their dynamic polymerization and shape exert force on the underlying membrane. We present a theoretical description of individual filaments that incorporates information from molecular dynamic simulations. The structure of the crystallized Methanococcus jannaschii FtsZ dimer was used to model a FtsZ pentamer that showed a curvature and a twist. The estimated bending and torsion angles between monomers and their fluctuations were included in the theoretical description. The MD data also permitted positioning the curvature with respect to the protein coordinates and allowed us to explore the effect of the relative orientation of the preferred curvature with respect to the surface plane. We find that maximum tension is attained when filaments are firmly attached and oriented with their curvature perpendicular to the surface and that the twist serves as a valve to release or to tighten the tension exerted by the curved filaments on the membrane. The theoretical model also shows that the presence of torsion can explain the shape distribution of short filaments observed by Atomic Force Microscopy in previously published experiments. New experiments with FtsZ covalently attached to lipid membranes show that the filament on-plane curvature depends on lipid head charge, confirming the predicted monomer orientation effects. This new model underlines the fact that the combination of the three elements, filament curvature, twist and the strength and orientation of its surface attachment, can modulate the force exerted on the membrane during cell division.

Simulation of catalytic water activation in mitochondrial F1-ATPase using a hybrid quantum mechanics/molecular mechanics approach: an alternative role for ß-Glu 188.

The use of quantum mechanics/molecular mechanics simulations to study the free energy landscape of the water activation at the catalytic site of mitochondrial F(1)-ATPase affords us insight into the generation of the nucleophile OH(-) prior to ATP hydrolysis. As a result, the ATP molecule was found to be the final proton acceptor. In the simulated pathway, the transfer of a proton to the nucleotide was not direct but occurred via a second water molecule in a manner similar to the Grotthuss mechanism proposed for proton diffusion. Residue ß-Glu 188, previously described as the putative catalytic base, was found to be involved in the stabilization of a transient hydronium ion during water activation. Simulations in the absence of the carboxylate moiety of ß-Glu 188 support this role.

The Role of Gln61 in HRas GTP hydrolysis: a quantum mechanics/molecular mechanics study.

Activation of the water molecule involved in GTP hydrolysis within the HRas·RasGAP system is analyzed using a tailored approach based on hybrid quantum mechanics/molecular mechanics (QM/MM) simulation. A new path emerges: transfer of a proton from the attacking water molecule to a second water molecule, then a different proton is transferred from this second water molecule to the GTP. Gln(61) will stabilize the transient OH(-) and H(3)O(+) molecules thus generated. This newly proposed mechanism was generated by using, for the first time to our knowledge, the entire HRas-RasGAP protein complex in a QM/MM simulation context. It also offers a rational explanation for previous experimental results regarding the decrease of GTPase rate found in the HRas Q61A mutant and the increase exhibited by the HRas Q61E mutant.

Molecular dynamics analysis of conformational change of paramyxovirus F protein during the initial steps of membrane fusion.

The fusion of paramyxovirus to the cell membrane is mediated by fusion protein (F protein) present in the virus envelope, which undergoes a dramatic conformational change during the process. Unlike hemagglutinin in orthomyxovirus, this change is not mediated by an alteration of environmental pH, and its cause remains unknown. Steered molecular dynamics analysis leads us to suggest that the conformational modification is mediated only by stretching mechanical forces once the transmembrane fusion peptide of the protein is anchored to the cell membrane. Such elongating forces will generate major secondary structure rearrangement in the heptad repeat A region of the F protein; from ß-sheet conformation to an elongated coil and then spontaneously to an α-helix. In addition, it is proposed that the heptad repeat A region adopts a final three-helix coiled coil and that this structure appears after the formation of individual helices in each monomer.

Comparing molecular dynamics force fields in the essential subspace.

The continued development and utility of molecular dynamics simulations requires improvements in both the physical models used (force fields) and in our ability to sample the Boltzmann distribution of these models. Recent developments in both areas have made available multi-microsecond simulations of two proteins, ubiquitin and Protein G, using a number of different force fields. Although these force fields mostly share a common mathematical form, they differ in their parameters and in the philosophy by which these were derived, and previous analyses showed varying levels of agreement with experimental NMR data. To complement the comparison to experiments, we have performed a structural analysis of and comparison between these simulations, thereby providing insight into the relationship between force-field parameterization, the resulting ensemble of conformations and the agreement with experiments. In particular, our results show that, at a coarse level, many of the motional properties are preserved across several, though not all, force fields. At a finer level of detail, however, there are distinct differences in both the structure and dynamics of the two proteins, which can, together with comparison with experimental data, help to select force fields for simulations of proteins. A noteworthy observation is that force fields that have been reparameterized and improved to provide a more accurate energetic description of the balance between helical and coil structures are difficult to distinguish from their "unbalanced" counterparts in these simulations. This observation implies that simulations of stable, folded proteins, even those reaching 10 microseconds in length, may provide relatively little information that can be used to modify torsion parameters to achieve an accurate balance between different secondary structural elements.

Molecular dynamics simulation of GTPase activity in polymers of the cell division protein FtsZ.

FtsZ, the prokaryotic ortholog of tubulin, assembles into polymers in the bacterial division ring. The interfaces between monomers contain a GTP molecule, but the relationship between polymerization and GTPase activity is not unequivocally proven. A set of short FtsZ polymers were modelled and the formation of active GTPase structures was monitored using molecular dynamics. Only the interfaces nearest the polymer ends exhibited an adequate geometry for GTP hydrolysis. Simulated conversion of interfaces from close-to-end to internal position and vice versa resulted in their spontaneous rearrangement between active and inactive conformations. This predicted behavior of FtsZ polymer ends was supported by in vitro experiments.