Role of MDH2 pathogenic variant in pheochromocytoma and paraganglioma patients.

PURPOSE: MDH2 (malate dehydrogenase 2) has recently been proposed as a novel potential pheochromocytoma/paraganglioma (PPGL) susceptibility gene, but its role in the disease has not been addressed. This study aimed to determine the prevalence of MDH2 pathogenic variants among PPGL patients and determine the associated phenotype. METHODS: Eight hundred thirty patients with PPGLs, negative for the main PPGL driver genes, were included in the study. Interpretation of variants of unknown significance (VUS) was performed using an algorithm based on 20 computational predictions, by implementing cell-based enzymatic and immunofluorescence assays, and/or by using a molecular dynamics simulation approach. RESULTS: Five variants with potential involvement in pathogenicity were identified: three missense (p.Arg104Gly, p.Val160Met and p.Ala256Thr), one in-frame deletion (p.Lys314del), and a splice-site variant (c.429+1G>T). All were germline and those with available biochemical data, corresponded to noradrenergic PPGL. CONCLUSION: This study suggests that MDH2 pathogenic variants may play a role in PPGL susceptibility and that they might be responsible for less than 1% of PPGLs in patients without pathogenic variants in other major PPGL driver genes, a prevalence similar to the one recently described for other PPGL genes. However, more epidemiological data are needed to recommend MDH2 testing in patients negative for other major PPGL genes.

Structural basis for interdomain communication in SHIP2 providing high phosphatase activity.

SH2-containing-inositol-5-phosphatases (SHIPs) dephosphorylate the 5-phosphate of phosphatidylinositol-3,4,5-trisphosphate (PI(3,4,5)P3) and play important roles in regulating the PI3K/Akt pathway in physiology and disease. Aiming to uncover interdomain regulatory mechanisms in SHIP2, we determined crystal structures containing the 5-phosphatase and a proximal region adopting a C2 fold. This reveals an extensive interface between the two domains, which results in significant structural changes in the phosphatase domain. Both the phosphatase and C2 domains bind phosphatidylserine lipids, which likely helps to position the active site towards its substrate. Although located distant to the active site, the C2 domain greatly enhances catalytic turnover. Employing molecular dynamics, mutagenesis and cell biology, we identify two distinct allosteric signaling pathways, emanating from hydrophobic or polar interdomain interactions, differentially affecting lipid chain or headgroup moieties of PI(3,4,5)P3. Together, this study reveals details of multilayered C2-mediated effects important for SHIP2 activity and points towards interesting new possibilities for therapeutic interventions.

MD simulations and FRET reveal an environment-sensitive conformational plasticity of importin-ß.

The nuclear pore complex mediates nucleocytoplasmic transport of macromolecules in eukaryotic cells. Transport through the pore is restricted by a hydrophobic selectivity filter comprising disordered phenylalanine-glycine-rich repeats of nuclear pore proteins. Exchange through the pore requires specialized transport receptors, called exportins and importins, that interact with cargo proteins in a RanGTP-dependent manner. These receptors are highly flexible superhelical structures composed of HEAT-repeat motifs that adopt various degrees of extension in crystal structures. Here, we performed molecular-dynamics simulations using crystal structures of Importin-ß in its free form or in complex with nuclear localization signal peptides as the starting conformation. Our simulations predicted that initially compact structures would adopt extended conformations in hydrophilic buffers, while contracted conformations would dominate in more hydrophobic solutions, mimicking the environment of the nuclear pore. We confirmed this experimentally by Förster resonance energy transfer experiments using dual-fluorophore-labeled Importin-ß. These observations explain seemingly contradictory crystal structures and suggest a possible mechanism for cargo protection during passage of the nuclear pore. Such hydrophobic switching may be a general principle for environmental control of protein function.

The SH2 domain regulates c-Abl kinase activation by a cyclin-like mechanism and remodulation of the hinge motion.

Regulation of the c-Abl (ABL1) tyrosine kinase is important because of its role in cellular signaling, and its relevance in the leukemiogenic counterpart (BCR-ABL). Both auto-inhibition and full activation of c-Abl are regulated by the interaction of the catalytic domain with the Src Homology 2 (SH2) domain. The mechanism by which this interaction enhances catalysis is not known. We combined computational simulations with mutagenesis and functional analysis to find that the SH2 domain conveys both local and global effects on the dynamics of the catalytic domain. Locally, it regulates the flexibility of the αC helix in a fashion reminiscent of cyclins in cyclin-dependent kinases, reorienting catalytically important motifs. At a more global level, SH2 binding redirects the hinge motion of the N and C lobes and changes the conformational equilibrium of the activation loop. The complex network of subtle structural shifts that link the SH2 domain with the activation loop and the active site may be partially conserved with other SH2-domain containing kinases and therefore offer additional parameters for the design of conformation-specific inhibitors.

Structural determinants and mechanism of mammalian CRM1 allostery.

Proteins carrying nuclear export signals cooperatively assemble with the export factor CRM1 and the effector protein RanGTP. In lower eukaryotes, this cooperativity is coupled to CRM1 conformational changes; however, it is unknown if mammalian CRM1 maintains its compact conformation or shows similar structural flexibility. Here, combinations of small-angle X-ray solution scattering and electron microscopy experiments with molecular dynamics simulations reveal pronounced conformational flexibility in mammalian CRM1 and demonstrate that RanGTP binding induces association of its N- and C-terminal regions to form a toroid structure. The CRM1 toroid is stabilized mainly by local interactions between the terminal regions, rather than by global strain. The CRM1 acidic loop is key in transmitting the effect of this RanGTP-induced global conformational change to the NES-binding cleft by shifting its population to the open state, which displays enhanced cargo affinity. Cooperative CRM1 export complex assembly thus constitutes a highly dynamic process, encompassing an intricate interplay of global and local structural changes.

Universal relaxation governs the nonequilibrium elasticity of biomolecules.

Experimental and computational dynamic force spectroscopy is widely used to determine the mechanical properties of single biomolecules. Whereas so far the focus has mainly been on rupture or unfolding forces, recent force-probe molecular dynamics simulations have revealed a strong loading rate dependence of biomolecular elasticities, which cannot be explained by the established one-dimensional transition-state treatments. We show that this nonequilibrium behavior can be explained by a theory that includes relaxation effects. For three structurally and mechanically quite diverse systems, a single relaxation mode suffices to quantitatively describe their loading-rate-dependent elastic behavior. Atomistic simulations of these systems revealed the microscopic nature of the respective relaxation modes. This result suggests a new type of "elasticity spectroscopy" experiment, which should render nonequilibrium properties of structured macromolecules accessible to single-molecule force spectroscopy.

The different flexibility of c-Src and c-Abl kinases regulates the accessibility of a druggable inactive conformation.

c-Src and c-Abl are two closely related protein kinases that constitute important anticancer targets. Despite their high sequence identity, they show different sensitivities to the anticancer drug imatinib, which binds specifically to a particular inactive conformation in which the Asp of the conserved DFG motif points outward (DFG-out). We have analyzed the DFG conformational transition of the two kinases using massive molecular dynamics simulations, free energy calculations, and isothermal titration calorimetry. On the basis of the reconstruction of the free energy surfaces for the DFG-in to DFG-out conformational changes of c-Src and c-Abl, we propose that the different flexibility of the two kinases results in a different stability of the DFG-out conformation and might be the main determinant of imatinib selectivity.

An unusual hydrophobic core confers extreme flexibility to HEAT repeat proteins.

Alpha-solenoid proteins are suggested to constitute highly flexible macromolecules, whose structural variability and large surface area is instrumental in many important protein-protein binding processes. By equilibrium and nonequilibrium molecular dynamics simulations, we show that importin-beta, an archetypical alpha-solenoid, displays unprecedentedly large and fully reversible elasticity. Our stretching molecular dynamics simulations reveal full elasticity over up to twofold end-to-end extensions compared to its bound state. Despite the absence of any long-range intramolecular contacts, the protein can return to its equilibrium structure to within 3 A backbone RMSD after the release of mechanical stress. We find that this extreme degree of flexibility is based on an unusually flexible hydrophobic core that differs substantially from that of structurally similar but more rigid globular proteins. In that respect, the core of importin-beta resembles molten globules. The elastic behavior is dominated by nonpolar interactions between HEAT repeats, combined with conformational entropic effects. Our results suggest that alpha-solenoid structures such as importin-beta may bridge the molecular gap between completely structured and intrinsically disordered proteins.

Hydrophilic linkers and polar contacts affect aggregation of FG repeat peptides.

Transport of large proteins into the nucleus involves two events, binding of the cargo protein to a transport receptor in the cytoplasm and passage of the cargo-transporter complex through the selective permeability barrier of the nuclear pore complex. The permeability barrier is formed by largely disordered polypeptides, each containing a number of conserved hydrophobic phenylalanine-glycine (FG) sequence motifs, connected by hydrophilic linkers of varying sequence (FG nups). How the motifs interact to form the permeability barrier, however, is not yet known. We have, therefore, carried out molecular dynamics simulations on various model FG repeat peptides to study the aggregation propensity of FG nups and the specific roles of the hydrophobic FG motifs and the hydrophilic linkers. Our simulations show spontaneous aggregation of the model nups into hydrated aggregates, which exhibit structural features assumed to be part of the permeability barrier. Our simulations suggest that short beta-sheets are an important structural feature of the aggregates and that Phe residues are sufficiently exposed to allow rapid binding of transport receptors. A surprisingly large influence of the amino acid composition of the hydrophilic linkers on aggregation is seen, as well as a major contribution of hydrogen-bonding patterns.

Influence of the g- conformation of Ser and Thr on the structure of transmembrane helices.

In order to study the influence of Ser and Thr on the structure of transmembrane helices we have analyzed a database of helix stretches extracted from crystal structures of membrane proteins and an ensemble of model helices generated by molecular dynamics simulations. Both complementary analyses show that Ser and Thr in the g- conformation induce and/or stabilize a structural distortion in the helix backbone. Using quantum mechanical calculations, we have attributed this effect to the electrostatic repulsion between the side chain Ogamma atom of Ser and Thr and the backbone carbonyl oxygen at position i-3. In order to minimize the repulsive force between these negatively charged oxygens, there is a modest increase of the helix bend angle as well as a local opening of the helix turn preceding Ser/Thr. This small distortion can be amplified through the helix, resulting in a significant displacement of the residues located at the other side of the helix. The crystal structures of aquaporin Z and the beta(2)-adrenergic receptor are used to illustrate these effects. Ser/Thr-induced structural distortions can be implicated in processes as diverse as ligand recognition, protein function and protein folding.

Structural models of class a G protein-coupled receptors as a tool for drug design: insights on transmembrane bundle plasticity.

G protein-coupled receptors (GPCRs) interact with an extraordinary diversity of ligands by means of their extracellular domains and/or the extracellular part of the transmembrane (TM) segments. Each receptor subfamily has developed specific sequence motifs to adjust the structural characteristics of its cognate ligands to a common set of conformational rearrangements of the TM segments near the G protein binding domains during the activation process. Thus, GPCRs have fulfilled this adaptation during their evolution by customizing a preserved 7TM scaffold through conformational plasticity. We use this term to describe the structural differences near the binding site crevices among different receptor subfamilies, responsible for the selective recognition of diverse ligands among different receptor subfamilies. By comparing the sequence of rhodopsin at specific key regions of the TM bundle with the sequences of other GPCRs we have found that the extracellular region of TMs 2 and 3 provides a remarkable example of conformational plasticity within Class A GPCRs. Thus, rhodopsin-based molecular models need to include the plasticity of the binding sites among GPCR families, since the "quality" of these homology models is intimately linked with the success in the processes of rational drug-design or virtual screening of chemical databases.

Computational study on the difference between the Co-C bond dissociation energy in methylcobalamin and adenosylcobalamin.

The bond dissociation energies of the Co-C bonds in the cobalamin cofactors methylcobalamin and adenosylcobalamin were calculated using the hybrid quantum mechanics/molecular mechanics method IMOMM (integrated molecular orbital and molecular mechanics). Calculations were performed on models of differing complexities as well as on the full systems. We investigated the origin of the different experimental values for the Co-C bond dissociation energies in methylcobalamin and adenosylcobalamin, and have provided an explanation for the difficulties encountered when we attempt to reproduce this difference in quantum chemistry. Additional calculations have been performed using the Miertus-Scrocco-Tomasi method in order to estimate the influence of solvent effects on the homolytic Co-C bond cleavage. Introduction of these solvation effects is shown to be necessary for the correct reproduction of experimental trends in bond dissociation energies in solution, which consequently have no direct correlation with dissociation processes in the enzyme.

Some critical issues in the application of quantum mechanics/molecular mechanics methods to the study of transition metal complexes.

The application of quantum mechanics/molecular mechanics (QM/MM) methods in transition metal chemistry is growing steadily. It becomes therefore appropriate to assess the importance of a number of technical issues associated to their implementation. This work presents the discussion of several of these issues, including the eventual need for conformational searches, the choice of the MM force field and the possibility of its tuning. The examples presented here prove that a proper handling of these technical aspects can lead to an improvement in the efficiency and quality of QM/MM calculations.