Dissociation of a Dynamic Protein Complex Studied by All-Atom Molecular Simulations

Abnormal function of EPHA2/p.R957P mutant in congenital cataract

AIM

To identify genetic defects in a Chinese family with congenital posterior polar cataracts and assess the pathogenicity.

METHODS

A four-generation Chinese family affected with autosomal dominant congenital cataract was recruited. Nineteen individuals took part in this study including 5 affected and 14 unaffected individuals. Sanger sequencing targeted hot-spot regions of 27 congenital cataract-causing genes for variant discovery. The pathogenicity of the variant was evaluated by the guidelines of American College of Medical Genetics and InterVar software. Confocal microscopy was applied to detect the subcellular localization of fluorescence-labeled ephrin type-A receptor 2 (EPHA2). Co-immunoprecipitation assay was implemented to estimate the interaction between EphA2 and other lens membrane proteins. The mRNA and protein expression were analyzed by reverse transcription-polymerase chain reaction (qRT-PCR) and Western blotting assay, respectively. The cell migration was analyzed by wound healing assay. Zebrafish model was generated by ectopic expression of human EPHA2/p.R957P mutant to demonstrate whether the mutant could cause lens opacity in vivo.

RESULTS

A novel missense and pathogenic variant c.2870G>C was identified in the sterile alpha motif (SAM) domain of EPHA2. Functional studies demonstrated the variant's impact: reduced EPHA2 protein expression, altered subcellular localization, and disrupted interactions with other lens membrane proteins. This mutant notably enhanced human lens epithelial cell migration, and induced a central cloudy region and roughness in zebrafish lenses with ectopic expression of human EPHA2/p.R957P mutant under differential interference contrast (DIC) optics.

CONCLUSION

Novel pathogenic c.2870G>C variant of EPHA2 in a Chinese congenital cataract family contributes to disease pathogenesis.

Structural basis underlying the synergism of NADase and SLO during group A Streptococcus infection

Group A Streptococcus (GAS) is a strict human pathogen possessing a unique pathogenic trait that utilizes the cooperative activity of NAD+-glycohydrolase (NADase) and Streptolysin O (SLO) to enhance its virulence. How NADase interacts with SLO to synergistically promote GAS cytotoxicity and intracellular survival is a long-standing question. Here, the structure and dynamic nature of the NADase/SLO complex are elucidated by X-ray crystallography and small-angle scattering, illustrating atomic details of the complex interface and functionally relevant conformations. Structure-guided studies reveal a salt-bridge interaction between NADase and SLO is important to cytotoxicity and resistance to phagocytic killing during GAS infection. Furthermore, the biological significance of the NADase/SLO complex in GAS virulence is demonstrated in a murine infection model. Overall, this work delivers the structure-functional relationship of the NADase/SLO complex and pinpoints the key interacting residues that are central to the coordinated actions of NADase and SLO in the pathogenesis of GAS infection.

Integration of immunoinformatics and cheminformatics to design and evaluate a multitope vaccine against Klebsiella pneumoniae and Pseudomonas aeruginosa coinfection

Introduction: Klebsiella pneumoniae (K. pneumoniae) and Pseudomonas aeruginosa (P. aeruginosa) are the most common Gram-negative bacteria associated with pneumonia and coinfecting the same patient. Despite their high virulence, there is no effective vaccine against them. Methods: In the current study, the screening of several proteins from both pathogens highlighted FepA and OmpK35 for K. pneumonia in addition to HasR and OprF from P. aeruginosa as promising candidates for epitope mapping. Those four proteins were linked to form a multitope vaccine, that was formulated with a suitable adjuvant, and PADRE peptides to finalize the multitope vaccine construct. The final vaccine's physicochemical features, antigenicity, toxicity, allergenicity, and solubility were evaluated for use in humans. Results: The output of the computational analysis revealed that the designed multitope construct has passed these assessments with satisfactory scores where, as the last stage, we performed a molecular docking study between the potential vaccine construct and K. pneumonia associated immune receptors, TLR4 and TLR2, showing affinitive to both targets with preferentiality for the TLR4 receptor protein. Validation of the docking studies has proceeded through molecular dynamics simulation, which estimated a strong binding and supported the nomination of the designed vaccine as a putative solution for K. pneumoniae and P. aeruginosa coinfection. Here, we describe the approach for the design and assessment of our potential vaccine.

HBD-2 variants and SARS-CoV-2: New insights into inter-individual susceptibility

Background

A deep understanding of the causes of liability to SARS-CoV-2 is essential to develop new diagnostic tests and therapeutics against this serious virus in order to overcome this pandemic completely. In the light of the discovered role of antimicrobial peptides [such as human b-defensin-2 (hBD-2) and cathelicidin LL-37] in the defense against SARS-CoV-2, it became important to identify the damaging missense mutations in the genes of these molecules and study their role in the pathogenesis of COVID-19.

Methods

We conducted a comprehensive analysis with multiple in silico approaches to identify the damaging missense SNPs for hBD-2 and LL-37; moreover, we applied docking methods and molecular dynamics analysis to study the impact of the filtered mutations.

Results

The comprehensive analysis reveals the presence of three damaging SNPs in hBD-2; these SNPs were predicted to decrease the stability of hBD-2 with a damaging impact on hBD-2 structure as well. G51D and C53G mutations were located in highly conserved positions and were associated with differences in the secondary structures of hBD-2. Docking-coupled molecular dynamics simulation analysis revealed compromised binding affinity for hBD-2 SNPs towards the SARS-CoV-2 spike domain. Different protein-protein binding profiles for hBD-2 SNPs, in relation to their native form, were guided through residue-wise levels and differential adopted conformation/orientation.

Conclusions

The presented model paves the way for identifying patients prone to COVID-19 in a way that would guide the personalization of both the diagnostic and management protocols for this serious disease.

Hunting for Novel Routes in Anticancer Drug Discovery: Peptides against Sam-Sam Interactions

Among the diverse protein binding modules, Sam (Sterile alpha motif) domains attract attention due to their versatility. They are present in different organisms and play many functions in physiological and pathological processes by binding multiple partners. The EphA2 receptor contains a Sam domain at the C-terminus (EphA2-Sam) that is able to engage protein regulators of receptor stability (including the lipid phosphatase Ship2 and the adaptor Odin). Ship2 and Odin are recruited by EphA2-Sam through heterotypic Sam-Sam interactions. Ship2 decreases EphA2 endocytosis and consequent degradation, producing chiefly pro-oncogenic outcomes in a cellular milieu. Odin, through its Sam domains, contributes to receptor stability by possibly interfering with ubiquitination. As EphA2 is upregulated in many types of tumors, peptide inhibitors of Sam-Sam interactions by hindering receptor stability could function as anticancer therapeutics. This review describes EphA2-Sam and its interactome from a structural and functional perspective. The diverse design strategies that have thus far been employed to obtain peptides targeting EphA2-mediated Sam-Sam interactions are summarized as well. The generated peptides represent good initial lead compounds, but surely many efforts need to be devoted in the close future to improve interaction affinities towards Sam domains and consequently validate their anticancer properties.

Dynamic stability of salt stable cowpea chlorotic mottle virus capsid protein dimers and pentamers of dimers

Intermediates of the self-assembly process of the salt stable cowpea chlorotic mottle virus (ss-CCMV) capsid can be modelled atomistically on realistic computational timescales either by studying oligomers in equilibrium or by focusing on their dissociation instead of their association. Our previous studies showed that among the three possible dimer interfaces in the icosahedral capsid, two are thermodynamically relevant for capsid formation. The aim of the current study is to evaluate the relative structural stabilities of the three different ss-CCMV dimers and to find and understand the conditions that lead to their dissociation. Long timescale molecular dynamics simulations at 300 K of the various dimers and of the pentamer of dimers underscore the importance of large contact surfaces on stabilizing the capsid subunits within an oligomer. Simulations in implicit solvent show that at higher temperature (350 K), the N-terminal tails of the protein units act as tethers, delaying dissociation for all but the most stable interface. The pentamer of dimers is also found to be stable on long timescales at 300 K, with an inherent flexibility of the outer protein chains.

Molecular Dynamics Simulation Investigation of the Binding and Interaction of the EphA6-Odin Protein Complex

Protein-protein interaction plays an important role in the development of almost all cells. Elucidating the dynamic binding and affinity of a protein-protein complex is essential for understanding the biological functions of proteins. EphA6 and Odin proteins are members of the Eph (erythropoietin-producing hepatocyte) receptor family and the Anks (ankyrin repeat and sterile α motif domain-containing) family, respectively. Odin significantly functions in regulating endocytosis, degradation, and stability of EphA receptors. In this work, the key residues of the interaction interface were determined through a hydrogen bond, root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), and dynamic correlation analysis of the conventional molecular dynamics (MD) simulations. The calculated standard binding free energy, -7.92 kcal/mol, between EphA6 and Odin is quite consistent with the experimental measurement value, -8.73 kcal/mol. By the combination of several MD simulation techniques, our investigation of the binding process reveals the detailed representative characteristics of the entire binding pathway at the molecular level. Based on the obtained potential of the mean force (PMF) curve, the analysis of the simulation trajectories shows that the residue Arg1013 in the receptor EphA6 is responsible for capturing Asp739 and Asp740 in the ligand Odin during the initial stage of binding. In the later stage of binding, the hydrogen bonds and salt bridges between a series of residues Lys973, Leu1007, Gly1009, His1010, and Arg1012 in the receptor and residues Leu735, Asn736, Asp739, Asp740, and Asp753 in the ligand mainly contribute to the stability of the protein complex. In addition, the specific change process of the receptor-ligand-binding mode is also clarified during the binding process. Our present simulation will promote a deep understanding of the protein-protein interaction, and the identified key interresidue interaction will be theoretical guidance for the design of protein drugs.

HBD-2 binds SARS-CoV-2 RBD and blocks viral entry: Strategy to combat COVID-19

New approaches to complement vaccination are needed to combat the spread of SARS-CoV-2 and stop COVID-19-related deaths and medical complications. Human beta defensin 2 (hBD-2) is a naturally occurring epithelial cell-derived host defense peptide that has anti-viral properties. Our comprehensive in-silico studies demonstrate that hBD-2 binds the site on the CoV-2-RBD that docks with the ACE2 receptor. Biophysical measurements confirm that hBD-2 indeed binds to the CoV-2-receptor-binding domain (RBD) (K_D ∼ 2μM by surface plasmon resonance), preventing it from binding to ACE2-expressing cells. Importantly, hBD-2 shows specificity by blocking CoV-2/spike pseudoviral infection, but not VSVG-mediated infection, of ACE2-expressing human cells with an IC₅₀ of 2.8 ± 0.4 μM. These promising findings offer opportunities to develop hBD-2 and/or its derivatives and mimetics to safely and effectively use as agents to prevent SARS-CoV-2 infection.

•

HBD-2 binds spike-RBD at the ACE2 interaction site in silico

•

Biophysical and biological assays confirm hBD-2 binding to spike-RBD

•

HBD-2 blocks spike-RBD:ACE2 binding

•

HBD-2 prevents CoV-2/spike pseudovirions from infecting ACE2-expressing human cells

Protein-Protein Interaction-Gaussian Accelerated Molecular Dynamics (PPI-GaMD): Characterization of Protein Binding Thermodynamics and Kinetics

Protein-protein interactions (PPIs) play key roles in many fundamental biological processes such as cellular signaling and immune responses. However, it has proven challenging to simulate repetitive protein association and dissociation in order to calculate binding free energies and kinetics of PPIs due to long biological timescales and complex protein dynamics. To address this challenge, we have developed a new computational approach to all-atom simulations of PPIs based on a robust Gaussian accelerated molecular dynamics (GaMD) technique. The method, termed "PPI-GaMD", selectively boosts interaction potential energy between protein partners to facilitate their slow dissociation. Meanwhile, another boost potential is applied to the remaining potential energy of the entire system to effectively model the protein's flexibility and rebinding. PPI-GaMD has been demonstrated on a model system of the ribonuclease barnase interactions with its inhibitor barstar. Six independent 2 μs PPI-GaMD simulations have captured repetitive barstar dissociation and rebinding events, which enable calculations of the protein binding thermodynamics and kinetics simultaneously. The calculated binding free energies and kinetic rate constants agree well with the experimental data. Furthermore, PPI-GaMD simulations have provided mechanistic insights into barstar binding to barnase, which involves long-range electrostatic interactions and multiple binding pathways, being consistent with previous experimental and computational findings of this model system. In summary, PPI-GaMD provides a highly efficient and easy-to-use approach for binding free energy and kinetics calculations of PPIs.

Letting go: Deep computational modeling insights into pH-dependent calcium affinity

Calcium and other cofactors can feature as key additions to a molecular interface, to the extent that the cofactor is completely buried in the bound state. How can such an interaction be regulated then? The answer: By facilitating a switch through an allosteric network. Although a number of unbinding mechanisms are being characterized, an extensive computational study by Joswig et al. reveals a detailed model for the pattern recognition receptor langerin.

Reaction Pathway Sampling and Free-Energy Analyses for Multimeric Protein Complex Disassembly by Employing Hybrid Configuration Bias Monte Carlo/Molecular Dynamics Simulation

Physicochemical characterization of multimeric biomacromolecule assembly and disassembly processes is a milestone to understand the mechanisms for biological phenomena at the molecular level. Mass spectroscopy (MS) and structural bioinformatics (SB) approaches have become feasible to identify subcomplexes involved in assembly and disassembly, while they cannot provide atomic information sufficient for free-energy calculation to characterize transition mechanism between two different sets of subcomplexes. To combine observations derived from MS and SB approaches with conventional free-energy calculation protocols, we here designed a new reaction pathway sampling method by employing hybrid configuration bias Monte Carlo/molecular dynamics (hcbMC/MD) scheme and applied it to simulate the disassembly process of serum amyloid P component (SAP) pentamer. The results we obtained are consistent with those of the earlier MS and SB studies with respect to SAP subcomplex species and the initial stage of SAP disassembly processes. Furthermore, we observed a novel dissociation event, ring-opening reaction of SAP pentamer. Employing free-energy calculation combined with the hcbMC/MD reaction pathway trajectories, we moreover obtained experimentally testable observations on (1) reaction time of the ring-opening reaction and (2) importance of Asp42 and Lys117 for stable formation of SAP oligomer.

Molecular dynamics simulations and functional studies reveal that hBD-2 binds SARS-CoV-2 spike RBD and blocks viral entry into ACE2 expressing cells

New approaches to complement vaccination are needed to combat the spread of SARS-CoV-2 and stop COVID-19 related deaths and long-term medical complications. Human beta defensin 2 (hBD-2) is a naturally occurring epithelial cell derived host defense peptide that has antiviral properties. Our comprehensive in-silico studies demonstrate that hBD-2 binds the site on the CoV-2-RBD that docks with the ACE2 receptor. Biophysical and biochemical assays confirm that hBD-2 indeed binds to the CoV-2-receptor binding domain (RBD) (KD ~ 300 nM), preventing it from binding to ACE2 expressing cells. Importantly, hBD-2 shows specificity by blocking CoV-2/spike pseudoviral infection, but not VSV-G mediated infection, of ACE2 expressing human cells with an IC50 of 2.4± 0.1 μM. These promising findings offer opportunities to develop hBD-2 and/or its derivatives and mimetics to safely and effectively use as novel agents to prevent SARS-CoV-2 infection.

Conformational Entropy from Mobile Bond Vectors in Proteins: A Viewpoint that Unifies NMR Relaxation Theory and Molecular Dynamics Simulation Approaches

A new method for determining conformational entropy in proteins is reported. Proteins prevail as conformational ensembles, p ∝ exp(-u). By selecting a bond vector (e.g., N-H) as a conformation representative, molecular dynamics simulations can provide (relative to a reference structure) p as approximate Boltzmann probability density and u as N-H potential of mean force (POMF). The latter is as accurate as implied by the force field but statistical in character; this limits the insights it can provide and its utilization. Conformational entropy is given exclusively by u. Deriving it from POMFs renders it accurate but statistical in character. Previously, we devised explicit (i.e., analytical but not exact) potentials made of Wigner functions, D_0K^L, with L ≤ 4, which closely resemble the corresponding POMFs in form; hence, they also approach the latter in accuracy. Such potentials can be beneficially characterized/compared in terms of composition, symmetry, and associated order parameters. In this study, we develop a method for deriving conformational entropy from them, which also features the benefits specified above. The method developed is applied to the dimerization of the Rho GTPase-binding domain of plexin-B1. Insights into local ordering, entropy compensation, and features of allostery are gained. In previous work, we developed the slowly relaxing local structure (SRLS) approach for the analysis of NMR relaxation from restricted bond vector motion in proteins. SRLS comprises explicit (restricting) potentials of the kind developed here. It also comprises diffusion tensors describing the local motion and related features of local geometry. The complete model fits experimental data. In future work, the explicit potentials developed here will be inserted unchanged in SRLS-based data fitting, thereby improving the picture of structural dynamics. Given that SRLS is unique in featuring potentials that can closely approach the corresponding POMFs in accuracy, the present study is an important step toward generally improving protein dynamics by NMR relaxation.

Disulfide Bonds Affect the Binding Sites of Human β Defensin Type 3 on Negatively Charged Lipid Membranes

Human β defensin type 3 (hBD-3) is a small natural antimicrobiotic. It is strongly cationic and has six cysteine residues which can form three pairs of intramolecular disulfide bonds under oxidized condition. Those disulfide bonds can break under reducing condition. However, the antibacterial activities of hBD-3 in its wild-type and analog forms are similar. In this project, the structure and dynamics of hBD-3 were investigated by running simulations on hBD-3 in its wild-type and analog forms in solvent, binding to negatively charged lipid bilayers, and self-assembly with POPG lipids. It was found that the RMSFs of hBD-3 in both its wild-type and analog forms are similar in solvent, while they are very diverse depending on the binding sites of hBD-3 with negatively charged bilayers. Calculating both the distance map and insertion depths for 18 hBD-3 molecules binding on the POPG bilayer, hBD-3 in its analog form binds stably with the POPG bilayer through the head and loop regions, while hBD-3 wild-type binds with the POPG bilayer on the two loop regions stably. hBD-3 analog caused membrane thinning and disrupted the POPG lipids more significantly than the wildtype. Based on the self-assembly simulations, hBD-3 monomer can bind with and embed inside the negatively charged POPG lipid membrane and have more contacts with the POPG lipid heads than with tails. The current work emphasized the structural diversity of hBD-3 interacting with negatively charged lipid membrane affected by the disulfide bonding states.

A Systematic Test of Receptor Binding Kinetics for Ligands in Tumor Necrosis Factor Superfamily by Computational Simulations

Ligands in the tumor necrosis factor (TNF) superfamily are one major class of cytokines that bind to their corresponding receptors in the tumor necrosis factor receptor (TNFR) superfamily and initiate multiple intracellular signaling pathways during inflammation, tissue homeostasis, and cell differentiation. Mutations in the genes that encode TNF ligands or TNFR receptors result in a large variety of diseases. The development of therapeutic treatment for these diseases can be greatly benefitted from the knowledge on binding properties of these ligand–receptor interactions. In order to complement the limitations in the current experimental methods that measure the binding constants of TNF/TNFR interactions, we developed a new simulation strategy to computationally estimate the association and dissociation between a ligand and its receptor. We systematically tested this strategy to a comprehensive dataset that contained structures of diverse complexes between TNF ligands and their corresponding receptors in the TNFR superfamily. We demonstrated that the binding stabilities inferred from our simulation results were compatible with existing experimental data. We further compared the binding kinetics of different TNF/TNFR systems, and explored their potential functional implication. We suggest that the transient binding between ligands and cell surface receptors leads into a dynamic nature of cross-membrane signal transduction, whereas the slow but strong binding of these ligands to the soluble decoy receptors is naturally designed to fulfill their functions as inhibitors of signal activation. Therefore, our computational approach serves as a useful addition to current experimental techniques for the quantitatively comparison of interactions across different members in the TNF and TNFR superfamily. It also provides a mechanistic understanding to the functions of TNF-associated cell signaling pathways.

Sam Domains in Multiple Diseases

BACKGROUND

The sterile alpha motif (Sam) domain is a small helical protein module, able to undergo homo- and hetero-oligomerization, as well as polymerization, thus forming different types of protein architectures. A few Sam domains are involved in pathological processes and consequently, they represent valuable targets for the development of new potential therapeutic routes. This study intends to collect state-of-the-art knowledge on the different modes by which Sam domains can favor disease onset and progression.

METHODS

This review was build up by searching throughout the literature, for: a) the structural properties of Sam domains, b) interactions mediated by a Sam module, c) presence of a Sam domain in proteins relevant for a specific disease.

RESULTS

Sam domains appear crucial in many diseases including cancer, renal disorders, cataracts. Often pathologies are linked to mutations directly positioned in the Sam domains that alter their stability and/or affect interactions that are crucial for proper protein functions. In only a few diseases, the Sam motif plays a kind of "side role" and cooperates to the pathological event by enhancing the action of a different protein domain.

CONCLUSION

Considering the many roles of the Sam domain into a significant variety of diseases, more efforts and novel drug discovery campaigns need to be engaged to find out small molecules and/or peptides targeting Sam domains. Such compounds may represent the pillars on which to build novel therapeutic strategies to cure different pathologies.

ATP Converts Aβ42 Oligomer into Off-Pathway Species by Making Contact with Its Backbone Atoms Using Hydrophobic Adenosine

Adenosine triphosphate (ATP) is newly expected to be involved in the clearance of amyloid β 1-42 (Aβ₄₂) fibril and its precursors, Aβ₄₂ oligomer. Meanwhile, the microscopic mechanism of the role in dissolving the protein aggregate still remains elusive. Aiming to elucidate the mechanism, we examined effects of ATP on the conformational change and thermodynamic stability of the protomer dimer of Aβ₄₂ pentamer and tetramer, Aβ₄₂(9), by employing all-atom molecular dynamics simulations. We observed interprotomer twisting and intraprotomer peeling of Aβ₄₂(9). These conformational changes remarkably accelerate dissociation of the protomer dimer. However, the presence of ATP itself has no positive effect on dissociation processes of the protomer dimer and a monomer from the dimer, indicating its irrelevance to decomposition of the Aβ₄₂ oligomer. Rather, it could be supposed that ATP prevents additional binding and rebinding of Aβ₄₂ monomers to the Aβ₄₂ oligomer and it then converts Aβ₄₂ oligomer into an off-pathway species which is excluded from Aβ₄₂ fibril growth processes. Interestingly, hydrophobic adenosine in ATP makes contact with Aβ₄₂(9) on its backbone atoms, with respect to both Aβ₄₂ monomers on the edge of Aβ₄₂(9) and dissociated Aβ₄₂ monomers in Aβ₄₂(9). These roles of ATP would be applied without regard to the structural polymorphism of the Aβ₄₂ fibril.

AMPK Signaling Regulates the Age-Related Decline of Hippocampal Neurogenesis

The global incidence of age-associated neurological diseases is expected to rise with increasingly greying societies. In the aged brain, there is a dramatic decrease in the number of stem cells, which is a main cause for the decrease in brain function. Intrinsic factors, such as cell metabolism, have been studied but its role in neurogenesis is still unknown. Therefore, this study sought to establish whether AMP-activated protein kinase (AMPK) signaling does indeed regulate hippocampal neurogenesis in the aged brain. We found that i) AMPKα2 was the predominant catalytic subunit in the subgranular and subventricular zones; ii) AMPK activation was at a significantly higher level in the aged vs. young hippocampus; iii) short term (7 days) treatment with selective AMPK signaling inhibitor Compound C (10 mg/kg/day, i.p.) significantly increased the numbers of newborn (BrdU+), Type 2 (MCM2+), and Type 3 (DCX+) neural stem cells, but not Type 1 (GFAP+/Sox2+) cells, in the aged hippocampus. Taken together, our results demonstrate that AMPK signaling plays a critical role in the age-related decline of hippocampal neurogenesis.

Interaction specificity of clustered protocadherins inferred from sequence covariation and structural analysis

Clustered protocadherins, a large family of paralogous proteins that play important roles in neuronal development, provide an important case study of interaction specificity in a large eukaryotic protein family. A mammalian genome has more than 50 clustered protocadherin isoforms, which have remarkable homophilic specificity for interactions between cellular surfaces. A large antiparallel dimer interface formed by the first 4 extracellular cadherin (EC) domains controls this interaction. To understand how specificity is achieved between the numerous paralogs, we used a combination of structural and computational approaches. Molecular dynamics simulations revealed that individual EC interactions are weak and undergo binding and unbinding events, but together they form a stable complex through polyvalency. Strongly evolutionarily coupled residue pairs interacted more frequently in our simulations, suggesting that sequence coevolution can inform the frequency of interaction and biochemical nature of a residue interaction. With these simulations and sequence coevolution, we generated a statistical model of interaction energy for the clustered protocadherin family that measures the contributions of all amino acid pairs at the interface. Our interaction energy model assesses specificity for all possible pairs of isoforms, recapitulating known pairings and predicting the effects of experimental changes in isoform specificity that are consistent with literature results. Our results show that sequence coevolution can be used to understand specificity determinants in a protein family and prioritize interface amino acid substitutions to reprogram specific protein-protein interactions.

Modified Potential Functions Result in Enhanced Predictions of a Protein Complex by All-Atom Molecular Dynamics Simulations, Confirming a Stepwise Association Process for Native Protein-Protein Interactions

The relative prevalence of native protein-protein interactions (PPIs) are the cornerstone for understanding the structure, dynamics and mechanisms of function of protein complexes. In this study, we develop a scheme for scaling the protein-water interaction in the CHARMM36 force field, in order to better fit the solvation free energy of amino acids side-chain analogues. We find that the molecular dynamics simulation with the scaled force field, CHARMM36s, as well as a recently released version, CHARMM36m, effectively improve on the overly sticky association of proteins, such as ubiquitin. We investigate the formation of a heterodimer protein complex between the SAM domains of the EphA2 receptor and the SHIP2 enzyme by performing a combined total of 48 μs simulations with the different potential functions. While the native SAM heterodimer is only predicted at a low rate of 6.7% with the original CHARMM36 force field, the yield is increased to 16.7% with CHARMM36s, and to 18.3% with CHARMM36m. By analyzing the 25 native SAM complexes formed in the simulations, we find that their formation involves a preorientation guided by Coulomb interactions, consistent with an electrostatic steering mechanism. In 12 cases, the complex could directly transform to the native protein interaction surfaces with only small adjustments in domain orientation. In the other 13 cases, orientational and/or translational adjustments are needed to reach the native complex. Although the tendency for non-native complexes to dissociate has nearly doubled with the modified potential functions, a dissociation followed by a reassociation to the correct complex structure is still rare. Instead, the remaining non-native complexes undergo configurational changes/surface searching, which, however, rarely leads to native structures on a time scale of 250 ns. These observations provide a rich picture of the mechanisms of protein-protein complex formation and suggest that computational predictions of native complex PPIs could be improved further.

Protein-protein binding pathways and calculations of rate constants using fully-continuous, explicit-solvent simulations

A grand challenge in the field of biophysics has been the complete characterization of protein-protein binding processes at atomic resolution. This characterization requires the direct simulation of binding pathways starting from the initial, unbound state and proceeding through states that are too transient to be captured by experiment. Here, we applied the weighted ensemble path sampling strategy to orchestrate atomistic simulation of protein-protein binding pathways. Our simulation generated 203 fully-continuous and independent pathways along with rate constants for the binding process involving the barnase and barstar proteins. Results reveal multiple binding pathways along a "funnel-like" free energy landscape in which the formation of the "encounter complex" intermediate is rate-limiting followed by a relatively rapid rearrangement of the encounter complex to the bound state. Among all diffusional collisions, only ∼11% were productive. In the most probable binding pathways, the proteins rotated to a large extent (likely via electrostatic steering) in order to collide productively followed by "rolling" of the proteins along each other's binding interfaces to reach the bound state. Consistent with experiment, R59 was identified as the most kinetically important barnase residue for the binding process. Furthermore, protein desolvation occurs late in the binding process during the rearrangement of the encounter complex to the bound state. Notably, the positions of crystallographic water molecules that bridge hydrogen bonds between barnase and barstar are occupied in the bound-state ensemble. Our simulation was completed in a month using 1600 CPU cores at a time, demonstrating that it is now practical to carry out atomistic simulations of protein-protein binding.

Computational studies of protein-protein dissociation by statistical potential and coarse-grained simulations: a case study on interactions between colicin E9 endonuclease and immunity proteins

Proteins carry out their diverse functions in cells by forming interactions with each other. The dynamics of these interactions are quantified by the measurement of association and dissociation rate constants. Relative to the efforts made to model the association of biomolecules, little has been studied to understand the principles of protein complex dissociation. Using the interaction between colicin E9 endonucleases and immunity proteins as a test system, here we develop a coarse-grained simulation method to explore the dissociation mechanisms of protein complexes. The interactions between proteins in the complex are described by the knowledge-based potential that was constructed by the statistics from available protein complexes in the structural database. Our study provides the supportive evidences to the dual recognition mechanism for the specificity of binding between E9 DNase and immunity proteins, in which the conserved residues of helix III of Im2 and Im9 proteins act as the anchor for binding, while the sequence variations in helix II make positive or negative contributions to specificity. Beyond that, we further suggest that this binding specificity is rooted in the process of complex dissociation instead of association. While we increased the flexibility of protein complexes, we further found that they are less prone to dissociation, suggesting that conformational fluctuations of protein complexes play important functional roles in regulating their binding and dissociation. Our studies therefore bring new insights to the molecule mechanisms of protein-protein interactions, while the method can serve as a new addition to a suite of existing computational tools for the simulations of protein complexes.

Phosphorylation promotes binding affinity of Rap-Raf complex by allosteric modulation of switch loop dynamics

The effects of phosphorylation of a serine residue on the structural and dynamic properties of Ras-like protein, Rap, and its interactions with effector protein Ras binding domain (RBD) of Raf kinase, in the presence of GTP, are investigated via molecular dynamics simulations. The simulations show that phosphorylation significantly effects the dynamics of functional loops of Rap which participate in the stability of the complex with effector proteins. The effects of phosphorylation on Rap are significant and detailed conformational analysis suggest that the Rap protein, when phosphorylated and with GTP ligand, samples different conformational space as compared to non-phosphorylated protein. In addition, phosphorylation of SER11 opens up a new cavity in the Rap protein which can be further explored for possible drug interactions. Residue network analysis shows that the phosphorylation of Rap results in a community spanning both Rap and RBD and strongly suggests transmission of allosteric effects of local alterations in Rap to distal regions of RBD, potentially affecting the downstream signalling. Binding free energy calculations suggest that phosphorylation of SER11 residue increases the binding between Rap and Raf corroborating the network analysis results. The increased binding of the Rap-Raf complex can have cascading effects along the signalling pathways where availability of Raf can influence the oncogenic effects of Ras proteins. These simulations underscore the importance of post translational modifications like phosphorylation on the functional dynamics in proteins and can be an alternative to drug-targeting, especially in notoriously undruggable oncoproteins belonging to Ras-like GTPase family.

The Sam-Sam interaction between Ship2 and the EphA2 receptor: design and analysis of peptide inhibitors

The lipid phosphatase Ship2 represents a drug discovery target for the treatment of different diseases, including cancer. Its C-terminal sterile alpha motif domain (Ship2-Sam) associates with the Sam domain from the EphA2 receptor (EphA2-Sam). This interaction is expected to mainly induce pro-oncogenic effects in cells therefore, inhibition of the Ship2-Sam/EphA2-Sam complex may represent an innovative route to discover anti-cancer therapeutics. In the present work, we designed and analyzed several peptide sequences encompassing the interaction interface of EphA2-Sam for Ship2-Sam. Peptide conformational analyses and interaction assays with Ship2-Sam conducted through diverse techniques (CD, NMR, SPR and MST), identified a positively charged penta-amino acid native motif in EphA2-Sam, that once repeated three times in tandem, binds Ship2-Sam. NMR experiments show that the peptide targets the negatively charged binding site of Ship2-Sam for EphA2-Sam. Preliminary in vitro cell-based assays indicate that -at 50 µM concentration- it induces necrosis of PC-3 prostate cancer cells with more cytotoxic effect on cancer cells than on normal dermal fibroblasts. This work represents a pioneering study that opens further opportunities for the development of inhibitors of the Ship2-Sam/EphA2-Sam complex for therapeutic applications.

A role of the SAM domain in EphA2 receptor activation

Among the 20 subfamilies of protein receptor tyrosine kinases (RTKs), Eph receptors are unique in possessing a sterile alpha motif (SAM domain) at their C-terminal ends. However, the functions of SAM domains in Eph receptors remain elusive. Here we report on a combined cell biology and quantitative fluorescence study to investigate the role of the SAM domain in EphA2 function. We observed elevated tyrosine autophosphorylation levels upon deletion of the EphA2 SAM domain (EphA2ΔS) in DU145 and PC3 prostate cancer cells and a skin tumor cell line derived from EphA1/A2 knockout mice. These results suggest that SAM domain deletion induced constitutive activation of EphA2 kinase activity. In order to explain these effects, we applied fluorescence correlation spectroscopy to investigate the lateral molecular organization of EphA2. Our results indicate that SAM domain deletion (EphA2ΔS-GFP) increases oligomerization compared to the full length receptor (EphA2FL-GFP). Stimulation with ephrinA1, a ligand for EphA2, induced further oligomerization and activation of EphA2FL-GFP. The SAM domain deletion mutant, EphA2ΔS-GFP, also underwent further oligomerization upon ephrinA1 stimulation, but the oligomers were larger than those observed for EphA2FL-GFP. Based on these results, we conclude that the EphA2 SAM domain inhibits kinase activity by reducing receptor oligomerization.

Different dynamics and pathway of disulfide bonds reduction of two human defensins, a molecular dynamics simulation study

Human defensins are a class of antimicrobial peptides that are crucial components of the innate immune system. Both human α defensin type 5 (HD5) and human β defensin type 3 (hBD-3) have 6 cysteine residues which form 3 pairs of disulfide bonds in oxidizing condition. Disulfide bond linking is important to the protein structure stabilization, and the disulfide bond linking and breaking order have been shown to influence protein function. In this project, microsecond long molecular dynamics simulations were performed to study the structure and dynamics of HD5 and hBD-3 wildtype and analogs which have all 3 disulfide bonds released in reducing condition. The structure of hBD-3 was found to be more dynamic and flexible than HD5, based on RMSD, RMSF, and radius of gyration calculations. The disulfide bridge breaking order of HD5 and hBD-3 in reducing condition was predicted by two kinds of methods, which gave consistent results. It was found that the disulfide bonds breaking pathways for HD5 and hBD-3 are very different. The breaking of disulfide bonds can influence the dimer interface by making the dimer structure less stable for both kinds of defensin. In order to understand the difference in dynamics and disulfide bond breaking pathway, hydrophilic and hydrophobic accessible surface areas (ASA), buried surface area between cysteine pairs, entropy of cysteine pairs, and internal energy were calculated. Comparing to the wildtype, hBD-3 analog is more hydrophobic, while HD5 is more hydrophilic. For hBD-3, the disulfide breaking is mainly entropy driven, while other factors such as the solvation effects may take the major role in controlling HD5 disulfide breaking pathway. Proteins 2017; 85:665-681. © 2016 Wiley Periodicals, Inc.

Molecular Dynamics Simulations Reveal Isoform Specific Contact Dynamics between the Plexin Rho GTPase Binding Domain (RBD) and Small Rho GTPases Rac1 and Rnd1

The Plexin family of transmembrane receptors are unique in that their intracellular region interacts directly with small GTPases of the Rho family. The Rho GTPase binding domain of plexin (RBD)-which is responsible for these interactions-can bind with Rac1 as well as Rnd1 GTPases. GTPase complexes have been crystallized with the RBDs of plexinA1, -A2, and -B1. The protein association is thought to elicit different functional responses in a GTPase and plexin isoform specific manner, but the origin of this is unknown. In this project, we investigated complexes between several RBD and Rac1/Rnd1 GTPases using multimicrosecond length all atom molecular dynamics simulations, also with reference to the free forms of the RBDs and GTPases. In accord with the crystallographic data, the RBDs experience more structural changes than Rho-GTPases upon complex formation. Changes in protein dynamics and networks of correlated motions are revealed by analyzing dihedral angle fluctuations in the proteins. The extent of these changes differs between the different RBDs and also between the Rac1 and Rnd1 GTPases. While the RBDs in the free and bound states have similar-if not decreased-correlations, correlations within the GTPases are increased upon binding. Mapping highly correlated residues to the structures, it is found that the plexinA1, -B1, and -A2 RBDs all have similar communication pathways within the ubiquitin fold, but that different residues are involved. Dynamic network analyses indicate that plexinA1 and -B1 RBDs interact with small GTPases in a similar manner, whereas complexes with the plexinA2 RBD display different features. Importantly complexes with Rnd1 have a considerable number of dynamic correlations and network connections between the proteins, whereas such features are missing in the RBD-Rac1 complexes. Overall, the simulations suggest mechanisms that are consistent with the experimental data on plexinB1 and indicate RBD and GTPase isoform specific changes in protein dynamics upon complex formation.