A Missense Variant Affecting the N-Terminal Domain of the Laminin-332 β3 Chain Results in a Distinct Form of Junctional Epidermolysis Bullosa With Altered Granulation Tissue Response and No New Blistering: A Second Family Report

Junctional epidermolysis bullosa (JEB) is a rare genodermatosis characterized by fragility of the skin and mucous membranes due to alterations in the dermal epidermal junction. This condition manifests as mechanically induced bullous lesions that heal with hypertrophic granulation tissue and/or atrophic scars. Here, we report two brothers carrying a homozygous LAMB3 missense variant, p.Gly254Asp, which affects the N-terminal end of the laminin-332 (LM332) β3 chain, previously described in another JEB family sharing a common ethnic origin and LAMB3 haplotype with the siblings reported here. Moreover, all affected patients with p.Gly254Asp mutation from both families exhibits a distinct phenotype consisting of a few localized long-standing skin lesions characterized by excessive granulation tissue formation or keloid scars, without new blistering, and associated with amelogenesis imperfecta. Our patients also showed nail dystrophy, expanding the phenotypic spectrum and confirming the peculiar role of the N-terminal end of the β3 chain in regulating the granulation tissue response associated with the wound healing process.

The Putative Role of TIM-3 Variants in Polyendocrine Autoimmunity: Insights from a WES Investigation

Autoimmune polyglandular syndrome (APS) comprises a complex association of autoimmune pathological conditions. APS Type 1 originates from loss-of-function mutations in the autoimmune regulator (AIRE) gene. APS2, APS3 and APS4 are linked to specific HLA alleles within the major histocompatibility complex, with single-nucleotide polymorphisms (SNPs) in non-HLA genes also contributing to disease. In general, variability in the AIRE locus and the presence of heterozygous loss-of-function mutations can impact self-antigen presentation in the thymus. In this study, whole-exome sequencing (WES) was performed on a sixteen-year-old female APS3A/B patient to investigate the genetic basis of her complex phenotype. The analysis identified two variants (p.Arg111Trp and p.Thr101Ile) of the hepatitis A virus cell receptor 2 gene (HAVCR2) encoding for the TIM-3 (T cell immunoglobulin and mucin domain 3) protein. These variants were predicted, through in silico analysis, to impact protein structure and stability, potentially influencing the patient's autoimmune phenotype. While confocal microscopy analysis revealed no alteration in TIM-3 fluorescence intensity between the PBMCs isolated from the patient and those of a healthy donor, RT-qPCR showed reduced TIM-3 expression in the patient's unfractionated PBMCs. A screening conducted on a cohort of thirty APS patients indicated that the p.Thr101Ile and p.Arg111Trp mutations were unique to the proband. This study opens the pathway for the search of TIM-3 variants possibly linked to complex autoimmune phenotypes, highlighting the potential of novel variant discovery in contributing to APS classification and diagnosis.

Functional Study of SNCA p.V15A Variant: Further Linking α-Synuclein and Glucocerebrosidase

BACKGROUND

SNCA p.V15A was reported in five families. In vitro models showed increased aggregation and seeding activity, mitochondrial damage, and apoptosis. Mutant flies had reduced flying ability and survival.

OBJECTIVES

To clinically and functionally evaluate SNCA p.V15A in a large Italian family with Parkinson's disease (PD).

METHODS

Genetic diagnosis was reached through next-generation sequencing. Pathogenicity was assessed by molecular dynamics simulation and biochemical studies on peripheral blood mononuclear cells (PBMCs).

RESULTS

Five siblings carried SNCA p.V15A; three developed bradykinetic-rigid PD in their 50s with rapid motor progression and variable cognitive impairment. A fourth sibling had isolated mood disturbance, whereas the fifth was still unaffected at age 47. The mutant protein showed decreased stability and an unstable folded structure. Proband's PBMCs showed elevated total and phosphorylated α-synuclein (α-syn) levels and significantly reduced glucocerebrosidase activity.

CONCLUSION

This study demonstrates accumulation of α-synV15A in PBMCs and strengthens the link between α-syn pathophysiology and glucocerebrosidase dysfunction. © 2024 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Exploiting in silico structural analysis to introduce emerging genotype-phenotype correlations in DHCR24-related sterol biosynthesis disorder: a case study

Desmosterolosis is a rare sterol biosynthesis disorder characterized by multiple congenital anomalies, failure to thrive, severe developmental delay, progressive epileptic encephalopathy, and elevated levels of desmosterol caused by biallelic mutations of DHCR24 encoding 3-β-hydroxysterol Δ-24-reductase. DHCR24 is regarded as the key enzyme of cholesterol synthesis in the metabolism of brain cholesterol as it catalyzes the reduction of the Δ-24 double bond of sterol intermediates during cholesterol biosynthesis. To date, 15 DHCR24 variants, detected in 2 related and 14 unrelated patients, have been associated with the desmosterolosis disorder. Here, we describe a proband harboring the never-described DHCR24 homozygous missense variant NM_014762.4:c.506T>C, NP_055577.1:p.M169T, whose functional validation was confirmed through biochemical assay. By using molecular dynamics simulation techniques, we investigated the impact of this variant on the protein stability and interaction network with the flavin adenine dinucleotide cofactor, thereby providing a preliminary assessment of its mechanistic role in comparison to all known pathogenic variants, the wild-type protein, and a known benign DHCR24 variant. This report expands the clinical and molecular spectra of the DHCR24-related disorder, reports on a novel DHCR24 deleterious variant associated with desmosterolosis, and gives new insights into genotype-phenotype correlations.

Clinical and molecular description of the first Italian cohort of 33 subjects with hypophosphatasia

Introduction

Hypophosphatasia (HPP) is a rare genetic disease caused by inactivating variants of the ALPL gene. Few data are available on the clinical presentation in Italy and/or on Italian HPP surveys.

Methods

There were 30 suspected HPP patients recruited from different Italian tertiary cares. Biological samples and related clinical, biochemical, and anamnestic data were collected and the ALPL gene sequenced. Search for large genomic deletions at the ALPL locus (1p36) was done. Phylogenetic conservation and modeling were applied to infer the effect of the variants on the protein structure.

Results

There were 21 ALPL variants and one large genomic deletion found in 20 out of 30 patients. Unexpectedly, NGS-driven differential diagnosis allowed uncovering three hidden additional HPP cases, for a total of 33 HPP subjects. Eight out of 24 coding variants were novel and classified as "pathogenic", "likely pathogenic", and "variants of uncertain significance". Bioinformatic analysis confirmed that all the variants strongly destabilize the homodimer structure. There were 10 cases with low ALP and high VitB6 that resulted negative to genetic testing, whereas two positive cases have an unexpected normal ALP value. No association was evident with other biochemical/clinical parameters.

Discussion

We present the survey of HPP Italian patients with the highest ALPL mutation rate so far reported and confirm the complexity of a prompt recognition of the syndrome, mostly for HPP in adults. Low ALP and high VitB6 values are mandatory for the genetic screening, this latter remaining the gold standard not only to confirm the clinical diagnosis but also to make differential diagnosis, to identify carriers, to avoid likely dangerous therapy in unrecognized cases.

Protein-Protein Interaction Studies Using Molecular Dynamics Simulation

Protein-protein interaction (PPI) is a crucial event for many biological functions. Studying the molecular details of PPI requires structure determination using X-ray crystallography, nuclear magnetic resistance (NMR), and single particle Cryo-EM. However, sometimes it is not easy to solve the complex structure for various reasons. For example, complex may be unstable, not enough protein expression for structural studies, etc. Further, PPI are intricate processes, and its molecular details cannot be fully explained by experimental observations. Here, we describe a quick and simple method to study the PPI using the combinatorial approach of molecular dynamics simulation and biophysical methods.

Connecting the dots: A practical evaluation of web-tools for describing protein dynamics as networks

Protein Structure Networks (PSNs) are a well-known mathematical model for estimation and analysis of the three-dimensional protein structure. Investigating the topological architecture of PSNs may help identify the crucial amino acid residues for protein stability and protein-protein interactions, as well as deduce any possible mutational effects. But because proteins go through conformational changes to give rise to essential biological functions, this has to be done dynamically over time. The most effective method to describe protein dynamics is molecular dynamics simulation, with the most popular software programs for manipulating simulations to infer interaction networks being RING, MD-TASK, and NAPS. Here, we compare the computational approaches used by these three tools-all of which are accessible as web servers-to understand the pathogenicity of missense mutations and talk about their potential applications as well as their advantages and disadvantages.

Design, synthesis and biological evaluation of N-(4-alkoxy-3-(1H-tetrazol-1-yl)phenyl) heterocyclic aromatic amide derivatives as xanthine oxidase inhibitors

Xanthine oxidase (XO) is a flavoprotein that exists in various organisms and can catalyze the uric acid formation in the human body. Based on the amide framework of N-(4-((3-cyanobenzyl)oxy)-3-(1H-tetrazol-1-yl)phenyl)isonicotinamide (compound 1) reported in our previous work, a series of N-(4-alkoxy-3-(1H-tetrazol-1-yl)phenyl) heterocyclic aromatic amide derivatives were designed, synthesized and evaluated as novel amide-based XO inhibitors. Structure-activity relationship campaign identified the most promising compound g25 (IC₅₀ = 0.022 μM), which possesses a special 1H-imidazole-5-carboxamide scaffold and presented comparable XO inhibitory potency to topiroxostat (IC₅₀ = 0.017 μM). Enzyme kinetic studies revealed that compound g25 acted as a mixed-type XO inhibitor. Molecular docking and molecular dynamics indicated that imidazole NH of g25 formed two stable hydrogen bonds with Glu1261 residue of XO that provided a vital contribution for the binding affinity. In addition, in vivo activity evaluation demonstrated that compound g25 exhibited obviously hypouricemic effect on a potassium oxonate induced hyperuricemic rat model.

Gain of Function of Malate Dehydrogenase 2 and Familial Hyperglycemia

CONTEXT

Genes causing familial forms of diabetes mellitus are only partially known.

OBJECTIVE

We set out to identify the genetic cause of hyperglycemia in multigenerational families with an apparent autosomal dominant form of adult-onset diabetes not due to mutations in known monogenic diabetes genes.

METHODS

Existing whole-exome sequencing (WES) data were used to identify exonic variants segregating with diabetes in 60 families from the United States and Italy. Functional studies were carried out in vitro (transduced MIN6-K8 cells) and in vivo (Caenorhabditis elegans) to assess the diabetogenic potential of 2 variants in the malate dehydrogenase 2 (MDH2) gene linked with hyperglycemia in 2 of the families.

RESULTS

A very rare mutation (p.Arg52Cys) in MDH2 strongly segregated with hyperglycemia in 1 family from the United States. An infrequent MDH2 missense variant (p.Val160Met) also showed disease cosegregation in a family from Italy, although with reduced penetrance. In silico, both Arg52Cys and Val160Met were shown to affect MDH2 protein structure and function. In transfected HepG2 cells, both variants significantly increased MDH2 enzymatic activity, thereby decreasing the NAD+/NADH ratio-a change known to affect insulin signaling and secretion. Stable expression of human wild-type MDH2 in MIN6-K8 cell lines enhanced glucose- and GLP-1-stimulated insulin secretion. This effect was blunted by the Cys52 or Met160 substitutions. Nematodes carrying equivalent changes at the orthologous positions of the mdh-2 gene showed impaired glucose-stimulated insulin secretion.

CONCLUSION

Our findings suggest a central role of MDH2 in human glucose homeostasis and indicate that gain of function variants in this gene may be involved in the etiology of familial forms of diabetes.

KDM6A missense variants hamper H3 histone demethylation in lung squamous cell carcinoma

KDM6A is the disease causative gene of type 2 Kabuki Syndrome, a rare multisystem disease; it is also a known cancer driver gene, with multiple somatic mutations found in a few cancer types. In this study, we looked at eleven missense variants in lung squamous cell carcinoma, one of the most common lung cancer subtypes, to see how they affect the KDM6A catalytic mechanisms. We found that they influence the interaction with histone H3 and the exposure of the trimethylated Lys27, which is critical for wild-type physiological function to varying degrees, by altering the conformational transition.

Inference of dynamic spatial GRN models with multi-GPU evolutionary computation

Reverse engineering mechanistic gene regulatory network (GRN) models with a specific dynamic spatial behavior is an inverse problem without analytical solutions in general. Instead, heuristic machine learning algorithms have been proposed to infer the structure and parameters of a system of equations able to recapitulate a given gene expression pattern. However, these algorithms are computationally intensive as they need to simulate millions of candidate models, which limits their applicability and requires high computational resources. Graphics processing unit (GPU) computing is an affordable alternative for accelerating large-scale scientific computation, yet no method is currently available to exploit GPU technology for the reverse engineering of mechanistic GRNs from spatial phenotypes. Here we present an efficient methodology to parallelize evolutionary algorithms using GPU computing for the inference of mechanistic GRNs that can develop a given gene expression pattern in a multicellular tissue area or cell culture. The proposed approach is based on multi-CPU threads running the lightweight crossover, mutation and selection operators and launching GPU kernels asynchronously. Kernels can run in parallel in a single or multiple GPUs and each kernel simulates and scores the error of a model using the thread parallelism of the GPU. We tested this methodology for the inference of spatiotemporal mechanistic gene regulatory networks (GRNs)-including topology and parameters-that can develop a given 2D gene expression pattern. The results show a 700-fold speedup with respect to a single CPU implementation. This approach can streamline the extraction of knowledge from biological and medical datasets and accelerate the automatic design of GRNs for synthetic biology applications.

MitImpact 3: modeling the residue interaction network of the Respiratory Chain subunits

Numerous lines of evidence have shown that the interaction between the nuclear and mitochondrial genomes ensures the efficient functioning of the OXPHOS complexes, with substantial implications in bioenergetics, adaptation, and disease. Their interaction is a fascinating and complex trait of the eukaryotic cell that MitImpact explores with its third major release. MitImpact expands its collection of genomic, clinical, and functional annotations of all non-synonymous substitutions of the human mitochondrial genome with new information on putative Compensated Pathogenic Deviations and co-varying amino acid sites of the Respiratory Chain subunits. It further provides evidence of energetic and structural residue compensation by techniques of molecular dynamics simulation. MitImpact is freely accessible at http://mitimpact.css-mendel.it.

Characterization of the adipogenic protein E4orf1 from adenovirus 36 through an in silico approach

Adenovirus 36 (Ad-36) is related to human obesity due to its adipogenic activity mediated by the early 4 open reading frame 1 (E4orf1) protein. Mechanisms underlying the adipogenic effect of E4orf1 are not completely understood; however, the proliferation and differentiation of fat cells are increased through the activation of the phosphatidyl inositol 3 kinase pathway by binding proteins containing PDZ domain. This study characterized E4orf1 tridimensional structure and analyzed its interactions with PDZ domain-containing proteins in order to provide new information about the behavior of this viral protein and its targets, which could provide an interesting druggable target for obesity-related cardiometabolic alterations. In silico strategies such as homology modeling, docking, and molecular dynamics (MD) were used to study the interaction of E4orf1 with five PDZ domains of disk large homolog 1 (PDZ-1 and PDZ-2), membrane-associated guanylate kinase 1 (PDZ-3), and multi-PDZ domain protein 1 (PDZ-7 and PDZ-10). Mutagenesis analysis of selected residues was performed to evaluate their effects on the stabilization of E4orf1:PDZ complexes. MD simulations showed that the E4orf1:PDZ10 complex was more stable than the others ones. The highly hydrophobic residues at the C-terminal region (114-125) of the E4orf1 are essential in the initial phase stabilization of the complexes. Moreover, the residues 80-85 in the core region contribute to longer stabilization of the E4orf1:PDZ10 complex, a result that was confirmed by in silico mutagenesis. In conclusion, E4orf1 forms a stable complex with PDZ10 domain, and the residues 80-85 are of particular importance. The characterization of E4orf1 interactions with PDZ domains provides an initial approach to discover druggable targets for Ad-36-induced obesity.

Mechanisms of pathogenesis of missense mutations on the KDM6A-H3 interaction in type 2 Kabuki Syndrome

Mutations in genes encoding for histone methylation proteins are associated with several developmental disorders. Among them, KDM6A is the disease causative gene of type 2 Kabuki Syndrome, a rare multisystem disease. While nonsense mutations and short insertions/deletions are known to trigger pathogenic mechanisms, the functional effects of missense mutations are still uncharacterized. In this study, we demonstrate that a selected set of missense mutations significantly hamper the interaction between KDM6A and the histone H3, by modifying the dynamics of the linker domain, and then causing a loss of function effect.

Insights into the molecular pathogenesis of cardiospondylocarpofacial syndrome: MAP3K7 c.737-7A > G variant alters the TGFβ-mediated α-SMA cytoskeleton assembly and autophagy

Transforming growth factor beta-activated kinase 1 (TAK1) is a highly conserved kinase protein encoded by MAP3K7, and activated by multiple extracellular stimuli, growth factors and cytokines. Heterozygous variants in MAP3K7 cause the cardiospondylocarpofacial syndrome (CSCFS) which is characterized by short stature, dysmorphic facial features, cardiac septal defects with valve dysplasia, and skeletal anomalies. CSCFS has been described in seven patients to date and its molecular pathogenesis is only partially understood. Here, the functional effects of the MAP3K7 c.737-7A > G variant, previously identified in a girl with CSCFS and additional soft connective tissue features, were explored. This splice variant generates an in-frame insertion of 2 amino acid residues in the kinase domain of TAK1. Computational analysis revealed that this in-frame insertion alters protein dynamics in the kinase activation loop responsible for TAK1 autophosphorylation after binding with its interactor TAB1. Co-immunoprecipitation studies demonstrate that the ectopic expression of TAK1-mutated protein impairs its ability to physically bind TAB1. In patient's fibroblasts, MAP3K7 c.737-7A > G variant results in reduced TAK1 autophosphorylation and dysregulation of the downstream TAK1-dependent signaling pathway. TAK1 loss-of-function is associated with an impaired TGFβ-mediated α-SMA cytoskeleton assembly and cell migration, and defective autophagy process. These findings contribute to our understanding of the molecular pathogenesis of CSCFS and might offer the rationale for the design of novel therapeutic targets.

Molecular modelling and dynamic simulations of sequestosome 1 (SQSTM1) missense mutations linked to Paget disease of bone

The Paget disease (PDB; OMIM is 167250) is a chronic bone disease caused by pathogenic mutations in Sequestome1/p62 (SQSTM1) gene. This study has aimed to interpret the relationship of PDB linked SQSTM1 mutations with protein structure and its molecular dynamic features. The disease causative missense mutations were initially collected, and then analyzed for their, exonic and domain distribution, impact on secondary and tertiary structures, and their ability on protein-ligand interactions, using a combination of systems biology approaches. Our results show that most PDB linked SQSTM1 missense mutations affect amino acid residues clustered within or near the UBA domain (aa 389-434), which participates in the ubiquitination of substrates. We also report that the majority mutations occurred in α-helices over β-strands but their effects on the secondary structure were mostly neutral. Global tertiary structure deviations were minimal; however, at amino acid residue level minor structural changes were evident. The molecular dynamics simulation analysis showed that both PB1 and UBA domains were under constant structural fluctuations resulting in closed form conformation of SQSMT1 protein structure, when it is bound to PRKCI ligand. We also found salt bridge conformation changes in the UBA domain of SQSTM1 mutants when they bound to the PRKCI interactor protein. This finding suggests the possibility that mutations in SQSTM1 could impair its ability to ubiquitinate the substrates, eventually affecting autophagy and apoptosis, especially in mature osteoclasts. This study presents the additional insight into structure and function relationship between SQSTM1 mutations and PDB pathogenesis. Communicated by Ramaswamy H. Sarma.

Structural and dynamic analysis of G558R mutation in chicken TSHR gene shows altered signal transduction and corroborates its role as a domestication gene

The thyroid-stimulating hormone receptor (TSHR) has been indicated as a putative domestication gene in chicken. Comparison of WGS identified a variant in residue 558 of the transmembrane domain (TM) of TSHR, where the domestic chicken (GGD) presents an arginine, whereas the red jungle fowl (RJF) shares a conserved glycine with other vertebrates. This variant has been demonstrated to be associated with phenotypes that are important for domestication and related to thyroid regulation, such as less fearful behavior, reduced aggressive behavior and reduced dependence on seasonal reproduction in GGD as compared with RJF. By means of molecular dynamics simulations, we highlighted the structural and dynamic differences of variant Gly558Arg in the TSHR TM domain. Alterations in TM helix flexibility, structure and protein overall motion are described. The so-called 'arginine snorkeling' of residue 568 in GGD is observed and we hypothesize it as the originating force that produces the observed whole-protein perturbation in the helix bundle dynamics, capable of altering the TSHR signal transduction. The results are discussed in the context of their implications for a better understanding of biological mechanisms in chicken under control of the thyroid, such as body metabolism, as well as for their usefulness in biomedical research.

Agenesis of the putamen and globus pallidus caused by recessive mutations in the homeobox gene GSX2

Basal ganglia are subcortical grey nuclei that play essential roles in controlling voluntary movements, cognition and emotion. While basal ganglia dysfunction is observed in many neurodegenerative or metabolic disorders, congenital malformations are rare. In particular, dysplastic basal ganglia are part of the malformative spectrum of tubulinopathies and X-linked lissencephaly with abnormal genitalia, but neurodevelopmental syndromes characterized by basal ganglia agenesis are not known to date. We ascertained two unrelated children (both female) presenting with spastic tetraparesis, severe generalized dystonia and intellectual impairment, sharing a unique brain malformation characterized by agenesis of putamina and globi pallidi, dysgenesis of the caudate nuclei, olfactory bulbs hypoplasia, and anomaly of the diencephalic-mesencephalic junction with abnormal corticospinal tract course. Whole-exome sequencing identified two novel homozygous variants, c.26C>A; p.(S9*) and c.752A>G; p.(Q251R) in the GSX2 gene, a member of the family of homeobox transcription factors, which are key regulators of embryonic development. GSX2 is highly expressed in neural progenitors of the lateral and median ganglionic eminences, two protrusions of the ventral telencephalon from which the basal ganglia and olfactory tubercles originate, where it promotes neurogenesis while negatively regulating oligodendrogenesis. The truncating variant resulted in complete loss of protein expression, while the missense variant affected a highly conserved residue of the homeobox domain, was consistently predicted as pathogenic by bioinformatic tools, resulted in reduced protein expression and caused impaired structural stability of the homeobox domain and weaker interaction with DNA according to molecular dynamic simulations. Moreover, the nuclear localization of the mutant protein in transfected cells was significantly reduced compared to the wild-type protein. Expression studies on both patients' fibroblasts demonstrated reduced expression of GSX2 itself, likely due to altered transcriptional self-regulation, as well as significant expression changes of related genes such as ASCL1 and PAX6. Whole transcriptome analysis revealed a global deregulation in genes implicated in apoptosis and immunity, two broad pathways known to be involved in brain development. This is the first report of the clinical phenotype and molecular basis associated to basal ganglia agenesis in humans.

Association of a homozygous GCK missense mutation with mild diabetes

Background

Homozygous inactivating GCK mutations have been repeatedly reported to cause severe hyperglycemia, presenting as permanent neonatal diabetes mellitus (PNDM). Conversely, only two cases of GCK homozygous mutations causing mild hyperglycemia have been so far described. We here report a novel GCK mutation (c.1116G>C, p.E372D), in a family with one homozygous member showing mild hyperglycemia.

Methods

GCK mutational screening was carried out by Sanger sequencing. Computational analyses to investigate pathogenicity and molecular dynamics (MD) were performed for GCK‐E372D and for previously described homozygous mutations associated with mild (n = 2) or severe (n = 1) hyperglycemia, used as references.

Results

Of four mildly hyperglycemic family‐members, three were heterozygous and one, diagnosed in the adulthood, was homozygous for GCK‐E372D. Two nondiabetic family members carried no mutations. Fasting glucose (p = 0.016) and HbA1c (p = 0.035) correlated with the number of mutated alleles (0–2).

In‐silico predicted pathogenicity was not correlated with the four mutations’ severity. At MD, GCK‐E372D conferred protein structure flexibility intermediate between mild and severe GCK mutations.

Conclusions

We present the third case of homozygous GCK mutations associated with mild hyperglycemia, rather than PNDM. Our in‐silico analyses support previous evidences suggesting that protein stability plays a role in determining clinical severity of GCK mutations.

Are Gaming-Enabled Graphic Processing Unit Cards Convenient for Molecular Dynamics Simulation?

In several fields of research, molecular dynamics simulation techniques are exploited to evaluate the temporal motion of particles constituting water, ions, small molecules, macromolecules, or more complex systems over time. These techniques are considered difficult to setup, computationally demanding and require high specialization and scientific skills. Moreover, they need specialized computing infrastructures to run faster and make the simulation of big systems feasible. Here, we have simulated 3 systems of increasing sizes on scientific- and gaming-enabled graphic processing unit (GPU) cards with Amber, GROMACS, and NAMD and measured their performance accounting also for the market prices of the GPU cards where they were run on.

Integration of network models and evolutionary analysis into high-throughput modeling of protein dynamics and allosteric regulation: theory, tools and applications

Proteins are dynamical entities that undergo a plethora of conformational changes, accomplishing their biological functions. Molecular dynamics simulation and normal mode analysis methods have become the gold standard for studying protein dynamics, analyzing molecular mechanism and allosteric regulation of biological systems. The enormous amount of the ensemble-based experimental and computational data on protein structure and dynamics has presented a major challenge for the high-throughput modeling of protein regulation and molecular mechanisms. In parallel, bioinformatics and systems biology approaches including genomic analysis, coevolution and network-based modeling have provided an array of powerful tools that complemented and enriched biophysical insights by enabling high-throughput analysis of biological data and dissection of global molecular signatures underlying mechanisms of protein function and interactions in the cellular environment. These developments have provided a powerful interdisciplinary framework for quantifying the relationships between protein dynamics and allosteric regulation, allowing for high-throughput modeling and engineering of molecular mechanisms. Here, we review fundamental advances in protein dynamics, network theory and coevolutionary analysis that have provided foundation for rapidly growing computational tools for modeling of allosteric regulation. We discuss recent developments in these interdisciplinary areas bridging computational biophysics and network biology, focusing on promising applications in allosteric regulations, including the investigation of allosteric communication pathways, protein–DNA/RNA interactions and disease mutations in genomic medicine. We conclude by formulating and discussing future directions and potential challenges facing quantitative computational investigations of allosteric regulatory mechanisms in protein systems.