Search

We previously molecularly and clinically characterized Mazzanti syndrome, a RASopathy related to Noonan syndrome that is mostly caused by a single recurrent missense variant (c.4A?>?G, p.Ser2Gly) in SHOC2, which encodes a leucine-rich repeat-containing protein facilitating signal flow through the RAS-mitogen-associated protein kinase (MAPK) pathway. We also documented that the pathogenic p.Ser2Gly substitution causes upregulation of MAPK signaling and constitutive targeting of SHOC2 to the plasma membrane due to the introduction of an N-myristoylation recognition motif. The almost invariant occurrence of the pathogenic c.4A?>?G missense change in SHOC2 is mirrored by a relatively homogeneous clinical phenotype of Mazzanti syndrome. Here, we provide new data on the clinical spectrum and molecular diversity of this disorder and functionally characterize new pathogenic variants. The clinical phenotype of six unrelated individuals carrying novel disease-causing SHOC2 variants is delineated, and public and newly collected clinical data are utilized to profile the disorder. In silico, in vitro and in vivo characterization of the newly identified variants provides evidence that the consequences of these missense changes on SHOC2 functional behavior differ from what had been observed for the canonical p.Ser2Gly change but converge toward an enhanced activation of the RAS-MAPK pathway. Our findings expand the molecular spectrum of pathogenic SHOC2 variants, provide a more accurate picture of the phenotypic expression associated with variants in this gene and definitively establish a gain-of-function behavior as the mechanism of disease., European Union (EU); Horizon 2020; European Joint Program on Rare Diseases (EJP-RD); Associazione Italiana per la Ricerca sul Cancro (AIRC); Italian Ministry of Health (Ricerca Corrente 2022); Italian Ministry of Research (FOE 2019); German Federal Ministry of Education and Research; BMBF (German Network for RASopathy Research ""GeNeRARe""); European Reference Network on Rare Congenital Malformations and Rare Intellectual Disability (ERN-ITHACA)

Vesicle biogenesis, trafficking and signaling via the ER-Golgi network support essential processes during development and their disruption can lead to neurodevelopmental disorders and neurodegeneration. We report that de novo missense variants in ARF3, encoding a small GTPase regulating Golgi structure and function, cause a neurodevelopmental disease showing microcephaly and progressive cortical atrophy, with microsomia and rib anomalies in severely affected subjects, suggesting a pleiotropic effect. All microcephaly-associated variants clustered in the guanine nucleotide binding pocket and perturbed the biochemical behavior of the protein by stabilizing it in a GTP-bound state. Functional analysis proved the disruptive consequences of the variants on Golgi integrity, and brain and body plan formation. In-depth analysis in zebrafish embryos expressing ARF3 mutants traced back the developmental alterations to defective gastrulation cell movements as the earliest detectable effect. Our findings document a role of ARF3 in Golgi homeostasis and demonstrate an obligate dependence for early development.

Smyd2 lysine methyltransferase regulates monomethylation of histone and nonhistone lysine residues using S‐adenosylmethionine cofactor as the methyl donor. The nonhistone interactors include several tumorigenic targets, including p53. Understanding this interaction would allow the structural principles that underpin Smyd2‐mediated p53 methylation to be elucidated. Here, we performed μ‐second molecular dynamics (MD) simulations on binary Smyd2‐cofactor and ternary Smyd2‐cofactor‐p53 peptide complexes. We considered both unmethylated and monomethylated p53 peptides (at Lys370 and Lys372). The results indicate that (a) the degree of conformational freedom of the C‐terminal domain of Smyd2 is restricted by the presence of the p53 peptide substrate, (b) the Smyd2 C‐terminal domain shows distinct dynamic properties when interacting with unmethylated and methylated p53 peptides, and (c) Lys372 methylation confines the p53 peptide conformation, with detectable influence on Lys370 accessibility to the cofactor. These MD results are therefore of relevance for studying the biology of p53 in cancer progression.

The search for stationary points in the molecular potential energy surfaces (PES) is a problem of increasing relevance in different fields of molecular sciences especially for large, flexible systems characterized by several large-amplitude internal motions leading to shallow minima with comparable energies and separated by small barriers. After structural biology and medicinal chemistry, also high-resolution molecular spectroscopy, which is the focus of our research activity, is nowadays shifting its attention to this kind of molecular systems. In such circumstances, accurate geometrical structures and relative stabilities of all these minima are a mandatory prerequisite for the vis-à-vis comparison between computed and experimental spectra. This task raises, in turn, the problem of the best compromise between accuracy and feasibility. In our opinion, a promising route is offered by a two-level stochastic search in which a relatively inexpensive MM or QM method is used in the initial search, followed by single point energy evaluation at a higher QM level of a relatively large number of low-energy structures in order to select a final short-list of candidates, whose geometries are fully optimized at the higher QM level. Finally, the relative stabilities and properties of the final short-list of energy minima can be computed by a state-of-the-art QM approach. This strategy defines a general two-level search/three-level evaluation approach, which can be finely tuned in terms of the accuracy of the sought results. Setup of the procedure, interface with a general purpose electronic structure code and validation of the most effective low-level methods for some representative molecular systems (three already well characterized and one new) ended up with a general, robust and user-friendly tool, which can be easily used and extended also by non-specialists to aid experimental spectroscopic studies.

Defective protein folding and accumulation of misfolded proteins is associated with neurodegenerative, cardiovascular, secretory, and metabolic disorders. Efforts are being made to identify small-molecule modulators or structural-correctors for conformationally destabilized proteins implicated in various protein aggregation diseases. Using a metastable-reporter-based primary screen, we evaluated pharmacological chaperone activity of a diverse class of natural products. We found that a flavonoid glycoside (C-10, chrysoeriol-7-O-β-D-glucopyranoside) stabilizes metastable proteins, prevents its aggregation, and remodels the oligomers into protease-sensitive species. Data was corroborated with additional secondary screen with disease-specific pathogenic protein. In vitro and cell-based experiments showed that C-10 inhibits α-synuclein aggregation which is implicated in synucleinopathies-related neurodegeneration. C-10 interferes in its structural transition into β-sheeted fibrils and mitigates α-synuclein aggregation-associated cytotoxic effects. Computational modeling suggests that C-10 binds to unique sites in α-synuclein which may interfere in its aggregation amplification. These findings open an avenue for comprehensive SAR development for flavonoid glycosides as pharmacological chaperones for metastable and aggregation-prone proteins implicated in protein conformational diseases.

Multiscale QM/MM approaches have become the most suitable and effective methods for the investigation of spectroscopic properties of medium- or large-size chromophores in condensed phases. On these grounds, we are developing a novel workflow aimed at improving the generality, reliability, and ease of use of the available tools. In the present paper, we report the latest developments of such an approach with specific reference to a general workplan starting with the addition of acetonitrile to the panel of solvents already available in the General Liquid Optimized Boundary (GLOB) model enforcing nonperiodic boundary conditions (NPBC). Next, the solvatochromic shifts induced by acetonitrile on both rigid (uracil and thymine) and flexible (thyrosine) chromophores have been studied introducing in our software a number of new features ranging from rigid-geometry NPBC molecular dynamics based on the quaternion formalism to a full integration of variational (ONIOM) and perturbative (perturbed matrix method (PMM)) approaches for describing different solute-solvent topologies and local fluctuations, respectively. Finally, thymine and uracil have been studied also in methanol to point out the generality of the computational strategy. While further developments are surely needed, the strengths of our integrated approach even in its present version are demonstrated by the accuracy of the results obtained by an unsupervised approach and coupled to a computational cost strongly reduced with respect to that of conventional QM/MM models without any appreciable accuracy deterioration.

Smyd2, the SET and MYND domain containing protein lysine methyltransferase, targets histone and nonhistone substrates. Methylation of nonhistone substrates has direct implications in cancer development and progression. Dynamic regulation of Smyd2 activity and the structural basis of broad substrate specificity still remain elusive. Herein, we report on extensive molecular dynamics simulations on a full length Smyd2 in the presence and absence of AdoMet cofactor (covering together 1.3 μs of sampling), and the accompanying conformational transitions. Additionally, dynamics of the C-terminal domain (CTD) and structural features of substrate crevices of Smyd2 and Smyd3 are compared. The CTD of Smyd2 exhibits conformational flexibility in both states. In the holo form, however, it undergoes larger hinge motions resulting in more opened configurations than the apo form, which is confined around the partially open starting X-ray configuration. AdoMet binding triggers increased elasticity of the CTD leading Smyd2 to adopt fully opened configurations, which completely exposes the substrate binding crevice. These long-range concerted motions highlight Smyd2's ability to target substrates of varying sizes. Substrate crevices of Smyd2 and Smyd3 show distinct features in terms of spatial, hydration, and electrostatic properties that emphasize their characteristic modes of substrates interaction and entry pathways for inhibitor binding. On the whole, our study shows how the elasticity and hinge motion of the CTD regulate its functional role and underpin the basis of broad substrate specificity of Smyd2. We also highlight the specific structural principles that guide substrate and inhibitor binding to Smyd2 and Smyd3.

SHOC2 is a scaffold protein mediating RAS-promoted activation of mitogen-activated protein kinase (MAPK) signaling in response to extracellular stimuli. A recurrent activating mutation in SHOC2 (p.Ser2Gly) causes Mazzanti syndrome, a RASopathy characterized by features resembling Noonan syndrome and distinctive ectodermal abnormalities. A second mutation (p.Met173Ile) supposed to cause loss-of-function was more recently identified in two individuals with milder phenotypes. Here, we report on the third RASopathy-causing SHOC2 mutation (c.807_808delinsTT, p.Gln269_His270delinsHisTyr), which was found associated with prenatal-onset hypertrophic cardiomyopathy. Structural analyses indicated a possible impact of the mutation on the relative orientation of the two SHOC2's leucine-rich repeat domains. Functional studies provided evidence of its activating role, revealing enhanced binding of the mutant protein to MRAS and PPP1CB, and increased signaling through the MAPK cascade. Differing from SHOC2 S2G , SHOC2 Q269_H270delinsHY is not constitutively targeted to the plasma membrane. These data document that diverse mechanisms in SHOC2 functional dysregulation converge toward MAPK signaling upregulation.

Lipid lateral diffusion in membrane bilayers is a fundamental process exploited by cells to enable complex protein structural and dynamic reorganizations. For its importance, lipid mobility in both cellular and model bilayers has been extensively investigated in recent years, especially through the application of time-resolved, fluorescence-based, optical microscopy techniques. However, one caveat of fluorescence techniques is the need to use dye-labeled variants of the lipid of interest, thus potentially perturbing the structural and dynamic properties of the native species. Generally, the effect of the dye/tracer molecule is implicitly assumed to be negligible. Nevertheless, in view of the widespread use of optically modified lipids for studying lipid bilayer dynamics, it is highly desirable to well assess this point. Here, fluorescence correlation spectroscopy (FCS) and molecular dynamics (MD) simulations have been combined together to uncover subtle structural and dynamic effects in DOPC planar membranes enriched with a standard Rhodamine-labeled lipid. Our findings support a non-neutral role of the dye-labeled lipids in diffusion experiments, quantitatively estimating a decrease in lipid mobility of up to 20% with respect to the unlabeled species. Moreover, results highlight the existing interplay between dye concentration, lipid lateral diffusion and membrane permeability, thus suggesting possible implications for future optical microscopy studies of biophysical processes occurring at the membrane level.

Computational simulation of UV/vis spectra in condensed phases can be performed starting from converged molecular dynamics (MD) simulations and then performing quantum mechanical/molecular mechanical (QM/MM) computations for a statistically significant number of snapshots. However, the need of variational solutions (e.g., ONIOM/EE) for a huge number of snapshots makes unpractical the use of state-of-the-art QM Hamiltonians. On the other hand, the effectivity of perturbative approaches (e.g., perturbed matrix method, PMM) comes at the price of poor convergence for configurations strongly different from the reference one. In this paper we introduce an integrated strategy based on a cluster analysis of the MD snapshots. Next, a representative configuration for each cluster is treated at the ONIOM/EE level, whereas local fluctuations within each cluster are described at the PMM level. Some representative systems (uracil in dimethylformamide and in water and tyrosine zwitterion in water) are analyzed to show the effectivity and flexibility of the proposed strategy.

Background The mechanism of how the hydrophilic threefold channel (C3) of ferritin nanocages facilitates diffusion of diverse metal ions into the internal cavity remains poorly explored. Methods Computational modeling and free energy estimations were carried out on R. catesbeiana H´ ferritin. Transit features and associated energetics for Fe2+, Mg2+, Zn2+ ions through the C3 channel have been examined. Results We highlight that iron conduction requires the involvement of two Fe2+ ions in the channel. In such doubly occupied configuration, as observed in X-ray structures, Fe2+ is displaced from the internal site (stabilized by D127) at lower energetic cost. Moreover, comparison of Fe2+, Mg2+ and Zn2+ transit features shows that E130 geometric constriction provides not only an electrostatic anchor to the incoming ions but also differentially influence their diffusion kinetics. Conclusions Overall, the study provides insights into Fe2+ entry mechanism and characteristic features of metal-protein interactions that influence the metal ions passage. The dynamics data suggest that E130 may act as a metal selectivity gate. This implicates an ion-specific entry mechanism through the channel with the distinct diffusion kinetics being the discriminating factor. General Significance Ferritin nanocages not only act as biological iron reservoirs but also have gained importance in material science as template scaffolds for synthesizing metal nanoparticles. This study provides mechanistic understanding on the conduction of different metal ions through the channel.

The fruitful interplay of high-resolution spectroscopy and quantum chemistry has a long history, especially in the field of small, semi-rigid molecules. However, in recent years, the targets of spectroscopic studies are shifting toward flexible molecules, characterized by a large number of closely spaced energy minima, all contributing to the overall spectrum. Here, artificial intelligence comes into play since it is at the basis of powerful unsupervised techniques for the exploration of soft degrees of freedom. Integration of such algorithms with a two-stage QM/QM′ (Quantum Mechanical) exploration/refinement strategy driven by a user-friendly graphical interface is the topic of the present paper. We will address in particular: (i) the performances of different semi-empirical methods for the exploration step and (ii) the comparison between stochastic and meta-heuristic algorithms in achieving a cheap yet complete exploration of the conformational space for medium sized chromophores. As test cases, we choose three amino acids of increasing complexity, whose full conformer enumeration has been reached only very recently. Next, we show that systems in condensed phases can be treated at the same level and with the same efficiency when employing a polarizable continuum description of the solvent. Finally, the challenging issue represented by the vibrational circular dichroism spectra of some rhodium complexes with flexible ligands has been addressed, showing that our fully unsupervised approach leads to remarkable agreement with the experiment. The fruitful interplay of high-resolution spectroscopy and quantum chemistry has a long history, especially in the field of small, semi-rigid molecules. However, in recent years, the targets of spectroscopic studies are shifting toward flexible molecules, characterized by a large number of closely spaced energy minima, all contributing to the overall spectrum. Here, artificial intelligence comes into play since it is at the basis of powerful unsupervised techniques for the exploration of soft degrees of freedom. Integration of such algorithms with a two-stage QM/QM′ (Quantum Mechanical) exploration/refinement strategy driven by a user-friendly graphical interface is the topic of the present paper. We will address in particular: (i) the performances of different semi-empirical methods for the exploration step and (ii) the comparison between stochastic and meta-heuristic algorithms in achieving a cheap yet complete exploration of the conformational space for medium sized chromophores. As test cases, we choose three amino acids of increasing complexity, whose full conformer enumeration has been reached only very recently. Next, we show that systems in condensed phases can be treated at the same level and with the same efficiency when employing a polarizable continuum description of the solvent. Finally, the challenging issue represented by the vibrational circular dichroism spectra of some rhodium complexes with flexible ligands has been addressed, showing that our fully unsupervised approach leads to remarkable agreement with the experiment.

Computational modeling involving Quantum Mechanics (QM) and Molecular Mechanics (MM) calculations are widely utilized to unveil the atomic-molecular properties that underpin their inherent characteristic features. The choice over the either of the QM and MM methods or a multiscale composite approach is driven by the target property of interest, and of course, the molecular size. Often, tailor-made schemes need to be devised to match the specific study purpose. Herein, we provide a perspective of these approaches addressing their effectiveness in terms of the delicate balance between the accuracy and computational feasibility. We focus on representative examples to highlight how different approaches can be fruitfully exploited for modeling the conformational landscape, and possibly, the spectroscopic behavior of biochemical molecules, especially amino acids building blocks.

Plasmodium falciparum encounters frequent environmental challenges during its life cycle which makes productive protein folding immensely challenging for its metastable proteome. To identify the important components of protein folding machinery involved in maintaining P. falciparum proteome, we performed a proteome-wide phylogenetic profiling across various species. We found that except HSP110, the parasite lost all other cytosolic nucleotide exchange factors essential for regulating HSP70 which is the centrum of the protein folding network. Evolutionary and structural analysis shows that besides its canonical interaction with HSP70, PfHSP110 has acquired sequence insertions for additional dynamic interactions. Molecular co-evolution profile depicts that the co-evolving proteins of PfHSP110 belong to distinct pathways like genetic variation, DNA repair, fatty acid biosynthesis, protein modification/trafficking, molecular motions, and apoptosis. These proteins exhibit unique physiochemical properties like large size, high iso-electric point, low solubility, and antigenicity, hence require PfHSP110 chaperoning to attain functional state. Co-evolving protein interaction network suggests that PfHSP110 serves as an important hub to coordinate protein quality control, survival, and immune evasion pathways in the parasite. Overall, our findings highlight potential accessory roles of PfHSP110 that may provide survival advantage to the parasite during its lifecycle and febrile conditions. The data also open avenues for experimental validation of auxiliary functions of PfHSP110 and their exploration for design of better antimalarial strategies. Proteins 2015; 83:1513–1525. © 2015 Wiley Periodicals, Inc.

We extend the library of solvents that can be treated using the GLOB (general liquid optimized boundary) method, that allows to perform MD simulations under non-periodic boundary conditions (NPBC) optimizing effective potentials between explicit molecules and the boundary for four organic solvents: CHCl3, CCl4, CH3OH and CH3CN. We show that GLOB allows reducing the number of explicit solvent shells to be included, while yielding results comparable with PBC and significant advantages over simulations without explicit boundaries. Finally, we provide polynomial fittings for all available GLOB effective potentials (including SPC water) to simplify their implementation in NPBC MD simulations.

Water solvation has a central role in several biochemical processes ranging from protein folding to biomolecular recognition and enzyme catalysis. Because of its importance, the structure and dynamics of hydration layers around biological macromolecules have been the targets of a great number of experimental and computational studies. In the present contribution, we have investigated the effects of periodic boundary conditions (PBCs), as used in conjunction with molecular dynamics (MD) simulations, on the dynamic and electric properties of water layers. In particular, we have systematically performed MD simulations of neat water and biomolecules in aqueous solutions by imposing a different external dielectric constant, a generally overlooked parameter in PBC simulations. The effect of the system size has also been addressed. Overall, our results consistently indicate that the dipole moment properties of water layers, and specifically the dipole moment fluctuations and the reorientational correlation functions, can be sensitive to the choice of the external boundary conditions, whereas other molecular properties, such as the self-diffusion coefficient and the reorientational relaxation times, are not affected. We think that our investigation may help to assess appropriate simulation conditions for modeling the aqueous environment of relevant biochemical systems and processes.

Dysfunctional organellar protein quality control machinery leads to protein misfolding associated cardiovascular, neurodegenerative, metabolic and secretory disorders. To understand organellar homeostasis, suitable tools are required which can sense changes in their respective protein folding capacity upon exposure to environmental and pharmacological perturbations. Herein, we have assessed protein folding capacity of cellular organelles using a metastable sensor selectively targeted to cytosol, nucleus, mitochondria, endoplasmic reticulum, golgi and peroxisomes. Microscopy and biochemical data revealed that these sensors report both acute and organelle-specific cellular insults. It also provided insights into contrasting refolding capacities of cellular organelles to recover from proteotoxic challenges. Further, we used these metastable sensors to evaluate pharmacological modulation of organellar protein folding capacity by small molecules. We observed pyrazole based scaffolds increased organellar protein folding capacity through upregulation of chaperones, mainly HSP90 and its co-chaperone HOP which coordinate refolding of misfolded/aggregated species. Overall, our data highlights the potential use of organelle-specific metastable sensors to understand protein folding capacity of sub-cellular compartments and assess pharmacological correction of their proteostasis imbalance. This study also provides additional avenue for use of these organelle-specific metastable sensors in drug discovery programs for identification of novel pharmacophores and drug repositioning of promising scaffolds for protein conformational diseases associated with different cellular organelles.

Background The bacterial proteins IscS, IscU and CyaY, the bacterial orthologue of frataxin, play an essential role in the biological machine that assembles the prosthetic Fe—S cluster groups on proteins. They form functionally binary and ternary complexes both in vivo and in vitro. Yet, the mechanism by which they work remains unclear. Methods We carried out extensive molecular dynamics simulations to understand the nature of their interactions and the role of dynamics starting from the crystal structure of a IscS-IscU complex and the experimentally-based model of a ternary IscS-IscU-CyaY complex and used nuclear magnetic resonance to experimentally test the interface. Results We show that, while being firmly anchored to IscS, IscU has a pivotal motion around the interface. Our results also describe how the catalytic loop of IscS can flip conformation to allow Fe—S cluster assembly. This motion is hampered in the ternary complex explaining its inhibitory properties in cluster formation. Conclusions We conclude that the observed ‘fluid’ IscS-IscU interface provides the binary complex with a functional adaptability exploited in partner recognition and unravels the molecular determinants of the reported inhibitory action of CyaY in the IscS-IscU-CyaY complex explained in terms of the hampering effect on specific IscU-IscS movements. General significance Our study provides the first mechanistic basis to explain how the IscS-IscU complex selects its binding partners and supports the inhibitory role of CyaY in the ternary complex., Graphical abstract Image 1, Highlights • IscU protomers show a pivotal motion around the interface with IscS, in complex. • Both IscS catalytic loop and IscU motion seem correlated with Fe—S cluster assembly. • Shuttling of the IscS catalytic loop is accompanied by a structural transition. • CyaY inhibits cluster formation by hampering the dynamics of the IscS-IscU complex.

It is widely accepted that endocytosis mediates the uptake of cationic cell penetrating peptides (CPPs) at relatively low concentrations (i.e. nano- to micromolar), while direct transduction across the plasma membrane comes into play at higher concentrations (i.e. micro- to millimolar). This latter process appears to depend on peptide-driven cellular processes, which in turn may induce local perturbations of plasma-membrane composition and/or integrity, and to be favored by peptide aggregation, especially into dimers. Besides, in most studies CPPs are tethered to fluorescent dyes in order to track peptide transduction events under the microscope, although often overlooking the possible role played by the dyes in assisting translocation. In an effort to provide some insights into the transduction process, here we report on a molecular dynamics (MD) simulation study of a prototype of the CPP family, namely the Tat11 arginine-rich motif. To be specific, the translocation of Tat11 across a purposely-created membrane pore, either or not covalently-linked to the tetramethylrhodamine-5-maleimide (TAMRA) dye and in both its monomeric and dimeric form, is analyzed in some detail. Results from several unconstrained and steered MD simulations, as well as energy decomposition analysis, nicely support the latest experimental evidence and help to shed light on key factors enabling peptide transduction. In particular, our study highlights the much slower translocation kinetics of Tat11 dimer in comparison to the single peptide, and therefore its enhanced capability to stabilize membrane pores. Notably, it also shows how TAMRA has overall negligible kinetic and energetic effects on peptide transduction, yet it promotes this process indirectly by favoring peptide aggregation.

Fluorescent molecular rotors (FMRs) belong to an important class of environment-sensitive dyes capable of acting as nanoprobes in the measurement of viscosity and polarity of their micro-environment. FMRs have found widespread applications in various research fields, ranging from analytical to biochemical sciences, for example in intracellular imaging studies or in volatile organic compound detection. Here, a computational investigation of a recently proposed FMR, namely 4-(diphenylamino)phthalonitrile (DPAP), in various chemical environments is presented. A purposely developed molecular mechanics force field is proposed and then applied to simulate the rotor in a high- and low-polar solvent (i.e., acetonitrile, tetrahydrofuran, o-xylene and cyclohexane), a polymer matrix and a lipid membrane. Subtle effects of the molecular interactions with the embedding medium, the structural fluctuations of the rotor and its rotational dynamics are analyzed in some detail. The results correlate with a previous work, thus supporting the reliability of the model, and provide further insights into the environment-specific properties of the dye. In particular, it is shown how molecular diffusion and rotational correlation times of the FMR are affected by the surrounding medium and how the molecular orientation of the dye becomes anisotropic once immersed in the lipid bilayer. Moreover, a qualitative correlation between the FMR rotational dynamics and the fluorescence lifetime is detected, a result in line with the observed viscosity dependence of its emission. Finally, optical absorption spectra are computed and successfully compared with their experimental counterparts.

Entry of HIV virus into the host cell is initiated by the interaction of its surface exposed gp120 protein with the cell surface CD4 receptor and a co-receptor that can be either CCR5 or CXCR4. The third variable region (V3 loop) of gp120 has an important role in co-receptor selection by gp120 and forms an epitope for neutralizing antibodies. In this work the dynamical behavior of the V3 loop in a trimeric environment has been investigated by generating an atomistic trimer model of gp120-CD4 complex and has been compared with the result of a monomeric gp120-CD4 complex. The main results coming from this work are that the three V3 loops belonging to the three subunits of the trimer display a different dynamical behavior in terms of its flexibility, spatial orientation, motion along the principal modes, conformations, solvent exposure and electrostatic potential distribution. We propose that the ability of the V3 loop to present, in the trimeric environment, simultaneous multiple alternative conformations that increase its capability of co-receptor recognition, is at least in part due to the effect of electrostatic potential generated by two subunits over the third one.

Non-coded alpha-amino acids, originally exploited by nature, have been successfully reproduced by recent synthetic strategies to confer special structural and functional properties to small peptides. The most known and well-studied atypical residue is alpha-aminoisobutyric acid (Aib), which is contained in a fairly large number of peptides with known antibiotic effects. Here, we report on a molecular dynamics (MD) study of a series of homooligopeptides based on alpha-aminoisobutyric acid (Aib) with increasing length (Ac-(Aib)(n)-NMe, n = 5, 6, 7 and 10) and at various temperatures, employing a recent extension of the AMBER force field tailored for the Aib residue. Solvent effects have been analyzed by comparative MD simulations of a heptapeptide in water and dimethylsulfoxide at different temperatures. Our results show that the preference for the 3(10)- and/or alpha-helix structures, which typically characterize Aib based peptides, is finely tuned by several factors including the chain length, temperature and solvent nature. While the transitions between intra-molecular i -> i + 3 and i -> i + 4 hydrogen bonds characterizing 3(10) and alpha-helices, respectively, are rather fast in small peptides (in the picosecond timescale), our analysis shows that the above physical and chemical factors modulate the relative equilibrium populations of the two helical structures. The obtained results nicely agree with available experimental data and support the use of the new force field for modeling Aib containing peptides.

Ferritin molecular cages are marvelous 24-mer supramolecular architectures that enable massive iron storage (>2000 iron atoms) within their inner cavity. This cavity is connected to the outer environment by two channels at C3 and C4 symmetry axes of the assembly. Ferritins can also be exploited as carriers for in vivo imaging and therapeutic applications, owing to their capability to effectively protect synthetic non-endogenous agents within the cage cavity and deliver them to targeted tissue cells without stimulating adverse immune responses. Recently, X-ray crystal structures of Fe(2+)-loaded ferritins provided important information on the pathways followed by iron ions toward the ferritin cavity and the catalytic centers within the protein. However, the specific mechanisms enabling Fe(2+) uptake through wild-type and mutant ferritin channels is largely unknown. To shed light on this question, we report extensive molecular dynamics simulations, site-directed mutagenesis, and kinetic measurements that characterize the transport properties and translocation mechanism of Fe(2+) through the two ferritin channels, using the wild-type bullfrog Rana catesbeiana H' protein and some of its variants as case studies. We describe the structural features that determine Fe(2+) translocation with atomistic detail, and we propose a putative mechanism for Fe(2+) transport through the channel at the C3 symmetry axis, which is the only iron-permeable channel in vertebrate ferritins. Our findings have important implications for understanding how ion permeation occurs, and further how it may be controlled via purposely engineered channels for novel biomedical applications based on ferritin.

Background The human lysine methyltransferase Smyd3, a member of the SET and MYND domain containing protein family, harbors methylation activity on both histone and non-histone targets in a tightly regulated manner. The mechanism of how Smyd3 dynamically regulates substrate recognition is still not fully unveiled. Methods Here, we employed molecular dynamics simulations on full length human Smyd3, performed to a total of 1.2 μ-second, in the presence (holo) and absence (apo) of the S-Adenosyl methionine (AdoMet) cofactor. The dynamical features of Smyd3 in apo and holo states have been examined and compared via examining geometrical and electrostatic properties. Results The results show a distinct dynamics of the C-terminal domain (CTD) in the two states. In the apo state, the CTD undergoes a large hinge like motion and samples more opened configurations, thus acting like a loosened clamp and resulting in expanded substrate binding crevice. In the holo state, the CTD exhibits a restricted motion while the overall structure remains compact, mimicking a closed clamp. This leads to a localized increase in the negative potential at the substrate binding cleft. Further, solvent accessibility of critical residues at the target lysine access channel, important for methylation activity, is increased. Conclusions We postulate that AdoMet cofactor acts like a key and locks Smyd3 in a closed conformation. In effect, the cofactor binding restricts the elasticity of the CTD, presenting a compact substrate binding cleft with high negative potential, which may have implications on substrate recruitment via long range electrostatics. General significance The deletion of the CTD from Smyd3 has been shown to abolish the basal histone methylation activity. Our study highlights the importance of the CTD elasticity in shaping the substrate binding site for recognition and supports the previously proposed role of the CTD in stabilizing the active site for methylation activity.

Classical Molecular Dynamics (MD) simulations can provide insights at the nanoscopic scale into protein dynamics. Currently, simulations of large proteins and complexes can be routinely carried out in the ns-μs time regime. Clustering of MD trajectories is often performed to identify selective conformations and to compare simulation and experimental data coming from different sources on closely related systems. However, clustering techniques are usually applied without a careful validation of results and benchmark studies involving the application of different algorithms to MD data often deal with relatively small peptides instead of average or large proteins; finally clustering is often applied as a means to analyze refined data and also as a way to simplify further analysis of trajectories. Herein, we propose a strategy to classify MD data while carefully benchmarking the performance of clustering algorithms and internal validation criteria for such methods. We demonstrate the method on two showcase systems with different features, and compare the classification of trajectories in real and PCA space. We posit that the prototype procedure adopted here could be highly fruitful in clustering large trajectories of multiple systems or that resulting especially from enhanced sampling techniques like replica exchange simulations.

Item does not contain fulltext Microtubules are dynamic cytoskeletal elements coordinating and supporting a variety of neuronal processes, including cell division, migration, polarity, intracellular trafficking, and signal transduction. Mutations in genes encoding tubulins and microtubule-associated proteins are known to cause neurodevelopmental and neurodegenerative disorders. Growing evidence suggests that altered microtubule dynamics may also underlie or contribute to neurodevelopmental disorders and neurodegeneration. We report that biallelic mutations in TBCD, encoding one of the five co-chaperones required for assembly and disassembly of the alphabeta-tubulin heterodimer, the structural unit of microtubules, cause a disease with neurodevelopmental and neurodegenerative features characterized by early-onset cortical atrophy, secondary hypomyelination, microcephaly, thin corpus callosum, developmental delay, intellectual disability, seizures, optic atrophy, and spastic quadriplegia. Molecular dynamics simulations predicted long-range and/or local structural perturbations associated with the disease-causing mutations. Biochemical analyses documented variably reduced levels of TBCD, indicating relative instability of mutant proteins, and defective beta-tubulin binding in a subset of the tested mutants. Reduced or defective TBCD function resulted in decreased soluble alpha/beta-tubulin levels and accelerated microtubule polymerization in fibroblasts from affected subjects, demonstrating an overall shift toward a more rapidly growing and stable microtubule population. These cells displayed an aberrant mitotic spindle with disorganized, tangle-shaped microtubules and reduced aster formation, which however did not alter appreciably the rate of cell proliferation. Our findings establish that defective TBCD function underlies a recognizable encephalopathy and drives accelerated microtubule polymerization and enhanced microtubule stability, underscoring an additional cause of altered microtubule dynamics with impact on neuronal function and survival in the developing brain.

In recent years, engineered biological pores responsive to external stimuli have been fruitfully used for various biotechnological applications. Moreover, the strategy of tethering photo-switchable moieties into biomolecules has provided an unprecedented temporal control of purposely designed nanodevices, as demonstrated, for example, by the light-mediated regulation of the activity of enzymes and biochannels. Inspired by these advancements, we propose here a de novo designed nanodevice featuring the α-hemolysin (αHL) membrane channel purposely functionalized by an artificial "on/off" molecular switch. The switch, which is based on the photo-isomerization of the azobenzene moiety, introduces a smart nano-valve into the natural non-gated pore to confer tunable transport properties. We validated through molecular dynamics simulations and free energy calculations the effective inter-conversion of the engineered αHL pore between two configurations corresponding to an "open" and a "closed" form. The reported switchable translocation of a single-stranded DNA fragment under applied voltage supports the promising capabilities of this nanopore prototype in view of molecular sensing, detection and delivery applications at single-molecule level.