Molecular dynamics simulation and essential dynamics study of mutated plastocyanin: structural, dynamical and functional effects of a disulfide bridge insertion at the protein surface

Scaffold-based selective SHP2 inhibitors design using core hopping, molecular docking, biological evaluation and molecular simulation

PTPN11 (coding the gene of SHP2), a classic non-receptor protein tyrosine phosphatase, is implicated in multiple cell signaling pathway. Abnormal activation of SHP2 has been shown to contribute to a variety of human diseases, including Juvenile myelomonocytic leukemia (JMML), Noonan syndrome and tumors. Thus, the SHP2 inhibitors have important therapeutic value. Here, based on the compound PubChem CID 8,478,960 (IC₅₀ = 45.01 μM), a series of thiophene [2,3-d] pyrimidine derivatives (IC₅₀ = 0.4-37.87 μM) were discovered as novel and efficient inhibitors of SHP2 through powerful "core hopping" and CDOCKER technology. Furthermore, the SHP2-PTP phosphatase activity assay indicated that Comp#5 (IC₅₀ = 0.4 μM) was the most active SHP2 inhibitor. Subsequently, the effects of Comp#5 on the structure and function of SHP2 were investigated through molecular dynamics (MD) simulation and post-kinetic analysis. The result indicated that Comp#5 enhanced the interaction of residues THR357, ARG362, LYS366, PRO424, CYS459, SER460, ALA461, ILE463, ARG465, THR507 and GLN510 with the surrounding residues, improving the stability of the catalytic active region and the entrance of catalytic active region. In particular, the Comp#5 conjugated with residue ARG362, elevating the efficient and selectivity of SHP2 protein. The study here may pave the way for discovering the novel SHP2 inhibitors for suffering cancer patients.

Effects of Peutz-Jeghers syndrome (PJS) causing missense mutations L67P, L182P, G242V and R297S on the structural dynamics of LKB1 (Liver kinase B1) protein

The liver kinase B1 (LKB1) is encoded by LKB1 gene. Several pathogenic mutations of LKB1 causing Peutz-Jeghers syndrome and also cancers in breast, gastric, pancreas, and colon have been reported. The present study is focused to analyze the effects on the structural dynamics of LKB1 caused by the 4 pathogenic missense mutations (L67P, L182P, G242V, and R297S), which are reported to reduce the catalytic activity. In this study, the structural changes of LKB1 in apo- and in heterotrimeric complex (LKB1-STRADα-MO25α) form with wild and mutated LKB1 are investigated using all atomistic molecular dynamic simulation. The present study reveals that these four mutations initiate local structural distortions and the solvent accessibility of the surrounding regions of ATP-binding pocket such as glycine-rich loop, αB and αC loop, activation and catalytic loops. The mutations of L67P, L182P, and G242 V induce distortions of the secondary structure of β1-β3 sheets, π - π interaction (observed between Phe204 of LKB1 and Phe243 of MO25α), and increase the helical properties (both helical twist and length) of the adjacent αH-helix, respectively. The active kinase features like the conformation of catalytic and activation loops, salt bridge and, finally, the formation of stable R- and C-hydrophobic spines are also found to be perturbed by these mutations. Hence, the observed mutation-induced structural distortions fail to coordinate the essential binding nature of LKB1 with STRADα and MO25α, which eventually affects the native function of LKB1. These observations are in line with the experimentally reported reduced kinase activity of LKB1.

Lid opening and conformational stability of T1 Lipase is mediated by increasing chain length polar solvents

The dynamics and conformational landscape of proteins in organic solvents are events of potential interest in nonaqueous process catalysis. Conformational changes, folding transitions, and stability often correspond to structural rearrangements that alter contacts between solvent molecules and amino acid residues. However, in nonaqueous enzymology, organic solvents limit stability and further application of proteins. In the present study, molecular dynamics (MD) of a thermostable Geobacillus zalihae T1 lipase was performed in different chain length polar organic solvents (methanol, ethanol, propanol, butanol, and pentanol) and water mixture systems to a concentration of 50%. On the basis of the MD results, the structural deviations of the backbone atoms elucidated the dynamic effects of water/organic solvent mixtures on the equilibrium state of the protein simulations in decreasing solvent polarity. The results show that the solvent mixture gives rise to deviations in enzyme structure from the native one simulated in water. The drop in the flexibility in H₂O, MtOH, EtOH and PrOH simulation mixtures shows that greater motions of residues were influenced in BtOH and PtOH simulation mixtures. Comparing the root mean square fluctuations value with the accessible solvent area (SASA) for every residue showed an almost correspondingly high SASA value of residues to high flexibility and low SASA value to low flexibility. The study further revealed that the organic solvents influenced the formation of more hydrogen bonds in MtOH, EtOH and PrOH and thus, it is assumed that increased intraprotein hydrogen bonding is ultimately correlated to the stability of the protein. However, the solvent accessibility analysis showed that in all solvent systems, hydrophobic residues were exposed and polar residues tended to be buried away from the solvent. Distance variation of the tetrahedral intermediate packing of the active pocket was not conserved in organic solvent systems, which could lead to weaknesses in the catalytic H-bond network and most likely a drop in catalytic activity. The conformational variation of the lid domain caused by the solvent molecules influenced its gradual opening. Formation of additional hydrogen bonds and hydrophobic interactions indicates that the contribution of the cooperative network of interactions could retain the stability of the protein in some solvent systems. Time-correlated atomic motions were used to characterize the correlations between the motions of the atoms from atomic coordinates. The resulting cross-correlation map revealed that the organic solvent mixtures performed functional, concerted, correlated motions in regions of residues of the lid domain to other residues. These observations suggest that varying lengths of polar organic solvents play a significant role in introducing dynamic conformational diversity in proteins in a decreasing order of polarity.

A principal component analysis of the dynamics of subdomains and binding sites in human serum albumin

The conformational dynamics of human serum albumin (HSA) was investigated by principal component analysis (PCA) applied to three molecular dynamics trajectories of 200 ns each. The overlap of the essential subspaces spanned by the first 10 principal components (PC) of different trajectories was about 0.3 showing that the PCA based on a trajectory length of 200 ns is not completely convergent for this protein. The contributions of the relative motion of subdomains and of the subdomains (internal) distortion to the first 10 PCs were found to be comparable. Based on the distribution of the first 3 PC, 10 protein conformers are identified showing relative root mean square deviations (RMSD) between 2.3 and 4.6 Å. The main PCs are found to be delocalized over the whole protein structure indicating that the motions of different protein subdomains are coupled. This coupling is considered as being related to the allosteric effects observed upon ligand binding to HSA. On the other hand, the first PC of one of the three trajectories describes a conformational transition of the protein domain I that is close to that experimentally observed upon myristate binding. This is a theoretical support for the older hypothesis stating that changes of the protein onformation favorable to binding can precede the ligand complexation. A detailed all atoms PCA performed on the primary Sites 1 and 2 confirms the multiconformational character of the HSA binding sites as well as the significant coupling of their motions.

The role of electrostatic interactions on klentaq1 insight for domain separation

We investigated the relationship between the thermostability of Klentaq1 and factors stabilizing interdomain interactions. When thermal adaptation of Klentaq1 was analyzed at the atomic level, the protein was stable at 300 and 350 K. It gradually unfolded at 373 K and almost spontaneously unfolded at 400 K. Domain separation was induced by disrupting electrostatic interactions in two salt bridges formed by Lys354-Glu445 and Asp371-Arg435 on the interface domain. The role of these interactions in protein stability was evaluated by comparing free energy solvation (ΔΔG_solv) between wild type and mutants. Substitution of Asp371 by Glu or Asn, and also Glu445 by Asn resulted in a positive value of ΔΔG_solv, suggesting that mutations destabilized the protein structure. Nevertheless, substitution of Glu445 by Asp gave a negative value to ΔΔG_solv reflecting increasing protein stability. Our results demonstrate that interactions at the interface domains of Klentaq1 are essential factors correlated with the Klentaq1 thermostability.

Molecular dynamics of a thermostable multicopper oxidase from Thermus thermophilus HB27: structural differences between the apo and holo forms

Molecular dynamic (MD) simulations have been performed on Tth-MCO, a hyperthermophilic multicopper oxidase from thermus thermophilus HB27, in the apo as well as the holo form, with the aim of exploring the structural dynamic properties common to the two conformational states. According to structural comparison between this enzyme and other MCOs, the substrate in process to electron transfer in an outer-sphere event seems to transiently occupy a shallow and overall hydrophobic cavity near the Cu type 1 (T1Cu). The linker connecting the β-strands 21 and 24 of the second domain (loop (β21–β24)_D2) has the same conformation in both states, forming a flexible lid at the entrance of the electron-transfer cavity. Loop (β21–β24)_D2 has been tentatively assigned a role occluding the access to the electron-transfer site. The dynamic of the loop (β21–β24)_D2 has been investigated by MD simulation, and results show that the structures of both species have the same secondary and tertiary structure during almost all the MD simulations. In the simulation, loop (β21–β24)_D2 of the holo form undergoes a higher mobility than in the apo form. In fact, loop (β21–β24)_D2 of the holo form experiences a conformational change which enables exposure to the electron-transfer site (open conformation), while in the apo form the opposite effect takes place (closed conformation). To confirm the hypothesis that the open conformation might facilitate the transient electron-donor molecule occupation of the site, the simulation was extended another 40 ns with the electron-donor molecule docked into the protein cavity. Upon electron-donor molecule stabilization, loops near the cavity reduce their mobility. These findings show that coordination between the copper and the protein might play an important role in the general mobility of the enzyme, and that the open conformation seems to be required for the electron transfer process to T1Cu.

How the local geometry of the Cu-binding site determines the thermal stability of blue copper proteins

Identifying the factors that govern the thermal resistance of cupredoxins is essential for understanding their folding and stability, and for improving our ability to design highly stable enzymes with potential biotechnological applications. Here, we show that the thermal unfolding of plastocyanins from two cyanobacteria--the mesophilic Synechocystis and the thermophilic Phormidium--is closely related to the short-range structure around the copper center. Cu K-edge X-ray absorption spectroscopy shows that the bond length between Cu and the S atom from the cysteine ligand is a key structural factor that correlates with the thermal stability of the cupredoxins in both oxidized and reduced states. These findings were confirmed by an additional study of a site-directed mutant of Phormidium plastocyanin showing a reverse effect of the redox state on the thermal stability of the protein.

Molecular dynamics of amicyanin reveals a conserved dynamical core for blue copper proteins

Molecular dynamics simulation has been carried out for the blue copper protein amicyanin from two different sources, Paracoccus denitrificans and Paraccocus versutus, to investigate the structural and dynamical properties common to the two molecules and to identify prominent features shared with proteins of the same family, the monomeric cupredoxins. The two amicyanins have almost identical secondary and tertiary structure. In the simulation, they differ for the number of hydrogen bonds in the main chain and the conformation of some beta-strands. However, they strictly maintain the arrangement of the portions of the beta-barrel that are conserved in the folding architecture of the blue copper proteins. Paracoccus versutus amicyanin equilibrates more rapidly, shows lower atomic deviation values, and is less rigid with respect to Paracoccus denitrificans amicyanin. Principal component analysis reveals that the conformational subspaces corresponding to eigenvectors with the same index for each of the two molecules are not necessarily equivalent. Nevertheless, a core scaffold with constrained dynamics exist for both amicyanins. In addition, two fairly flexible regions that are located on the opposite side with respect to the interaction sites with the partner molecules in the redox process have been evidenced in the protein structure. This description of amicyanin, with a few mobile regions remote from the active site and a rigid scaffold including most of the protein beta-barrel, has a close similarity with that of azurin and plastocyanin, two other cupredoxins previously investigated in simulation. Furthermore, similarities in the distribution of the atomic fluctuations indicate that amicyanin, azurin, and plastocyanin possess common dynamical features, in spite of differences in their structure. On the basis of these findings, we suggest that topological constraints imposed by the folding in correspondence of protein regions that are the most conserved determine the protein dynamics of the cupredoxin family. The dynamical properties of the cupredoxins might be controlled for functional advantages that include the binding mechanism with the biological partners and the collective inner motions of the protein matrix required for the electron transfer, whereas long-range conformational changes in the redox reaction should be excluded.

A study of collective atomic fluctuations and cooperativity in the U1A-RNA complex based on molecular dynamics simulations

Cooperative interactions play an important role in recognition and binding in macromolecular systems. In this study, we find that cross-correlated atomic fluctuations can be used to identify cooperative networks in a protein-RNA system. The dynamics of the RRM-containing protein U1A-stem loop 2 RNA complex have been calculated theoretically from a 10 ns molecular dynamics (MD) simulation. The simulation was analyzed by calculating the covariance matrix of all atomic fluctuations. These matrix elements are then presented in the form of a two-dimensional grid, which displays fluctuations on a per residue basis. The results indicate the presence of strong, selective cross-correlated fluctuations throughout the RRM in U1A-RNA. The atomic fluctuations correspond well with previous biophysical studies in which a multiplicity of cooperative networks have been reported and indicate that the various networks identified in separate individual experiments are fluctuationally correlated into a hyper-network encompassing most of the RRM. The calculated results also correspond well with independent results from a statistical covariance analysis of 330 aligned RRM sequences. This method has significant implications as a predictive tool regarding cooperativity in the protein-nucleic acid recognition process.

Role of disulfide bonds in modulating internal motions of proteins to tune their function: molecular dynamics simulation of scorpion toxin Lqh III

A series of 1-ns MD simulations were performed on the scorpion toxin Lqh III in native and disulfide bond broken states. The removal of disulfide bonds has caused hydrogen bond network alteration in the five-residue turn, the long loop, the alpha-helix, the loop connecting strands II and III, and the C-terminal region. In addition and more importantly, it has influenced the amplitude of the fluctuations of five-residue turn, loops, and C-terminal region with a minor effect on the fluctuations of the cysteines in the broken bond sites. These findings suggest that disulfide bonds are not the most important factors in rigidifying their own locations, while they have more important effects at a global scale. Furthermore, our results reveal that disulfide bonds have considerable influence on the functionally important essential modes of motions and the correlations between the motions of the binding site residues. Therefore, we can conclude that disulfide bonds have a crucial role in modulating the function via adjusting the dynamics of scorpion toxin molecules. Although this conclusion cannot be generalized to all peptides and proteins, it demonstrates the importance of more investigations on this aspect of disulfide bond efficacy.

Molecular Dynamics of Apo-Adenylate Kinase: A Principal Component Analysis

Adenylate kinase from E. coli (AKE) is studied with molecular dynamics. AKE undergoes large-scale motions of its Lid and AMP-binding domains when its open form closes over its substrates, AMP and Mg2+-ATP. The third domain, the Core, is relatively stable during closing. The resulting trajectory is analyzed with a principal component analysis method that decomposes the atom motions into modes ordered by their decreasing contributions to the total protein fluctuation. Simulations at 303 K (normal T) and 500 K (high T) reveal that at both temperatures the first three modes account for 70% of the total fluctuation. The residues that contribute the most to these three modes are concentrated in the Lid and AMP-binding domains. Analysis of the normal T modes indicates that the Lid and AMP-binding domains sample a broad distribution of conformations indicating that AKE is designed to provide its substrates with a large set of conformations. The high T results show that the Lid initially closes toward the Core. Subsequently, the Lid rotates to a new stable conformation that is different from what is observed in the substrate-bound AKE. These results are discussed in the context of experimental data that indicate that adenylate kinases do sample more than one conformational state in solution and that each of these conformational states undergoes substantial fluctuations. A pair of residues is suggested for labeling that would be useful for monitoring distance fluctuations by energy transfer experiments.

Conservation and specialization in PAS domain dynamics

The PAS (Per-ARNT-Sim) superfamily is presented as a well-suited study case to demonstrate how comparison of functional motions among distant homologous proteins with conserved fold characteristics may give insight into their functional specialization. Based on the importance of structural flexibility of the receptive structures in anticipating the signal-induced conformational changes of these sensory systems, the dynamics of these structures were analysed. Molecular dynamics was proved to be an effective method to obtain a reliable picture of the dynamics of the crystal structures of HERG, phy3, PYP and FixL, provided that an extensive conformational space sampling is performed. Other reliable sources of dynamic information were the ensembles of NMR structures of hPASK, HIF-2alpha and PYP. Essential dynamics analysis was successfully employed to extract the relevant information from the sampled conformational spaces. Comparison of motion patterns in the essential subspaces, based on the structural alignment, allowed identification of the specialized region in each domain. This appears to be evolved in the superfamily by following a specific trend, that also suggests the presence of a limited number of general solutions adopted by the PAS domains to sense external signals. These findings may give insight into unknown mechanisms of PAS domains and guide further experimental studies.

Thermal stability of wild type and disulfide bridge containing mutant of poplar plastocyanin

A comparative study of the thermal stability of wild type poplar plastocyanin and of a mutant form containing a disulfide bridge between residues 21 and 25 was performed using differential scanning calorimetry and optical spectroscopic techniques. For wild type plastocyanin the transition temperature, determined from the calorimetric profiles, is 62.7 degrees C at the scan rate of 60 degrees C/h, whereas for the mutant it is reduced to 58.0 degrees C. In both cases, the endothermic peak is followed by an exothermic one at higher temperatures. The unfolding process monitored by optical absorption at 596 nm also reveals a reduced thermal stability of the mutated plastocyanin compared to the wild type protein, with transition temperatures of 54.8 and 58.0 degrees C, respectively. For both proteins, the denaturation process was found to be irreversible and dependent on the scan rate preventing the thermodynamic analysis of the unfolding process. In parallel, small conformational changes between wild type and mutant plastocyanin emerge from fluorescence spectroscopy measurements. Here, a difference in the interaction of the two proteins between the microenvironment surrounding the fluorophores and the solvent was proposed. The destabilization observed in the disulfide containing mutant of plastocyanin suggests that the double mutation, Ile21Cys and Glu25Cys, introduces strain into the protein which offsets the stabilizing effect expected from the formation of a covalent crosslink.

In silico evidence for DNA polymerase-beta's substrate-induced conformational change

Structural information for mammalian DNA pol-beta combined with molecular and essential dynamics studies have provided atomistically detailed views of functionally important conformational rearrangements that occur during DNA repair and replication. This conformational closing before the chemical reaction is explored in this work as a function of the bound substrate. Anchors for our study are available in crystallographic structures of the DNA pol-beta in "open" (polymerase bound to gapped DNA) and "closed" (polymerase bound to gapped DNA and substrate, dCTP) forms; these different states have long been used to deduce that a large-scale conformational change may help the polymerase choose the correct nucleotide, and hence monitor DNA synthesis fidelity, through an "induced-fit" mechanism. However, the existence of open states with bound substrate and closed states without substrates suggest that substrate-induced conformational closing may be more subtle. Our dynamics simulations of two pol-beta/DNA systems (with/without substrates at the active site) reveal the large-scale closing motions of the thumb and 8-kDa subdomains in the presence of the correct substrate--leading to nearly perfect rearrangement of residues in the active site for the subsequent chemical step of nucleotidyl transfer--in contrast to an opening trend when the substrate is absent, leading to complete disassembly of the active site residues. These studies thus provide in silico evidence for the substrate-induced conformational rearrangements, as widely assumed based on a variety of crystallographic open and closed complexes. Further details gleaned from essential dynamics analyses clarify functionally relevant global motions of the polymerase-beta/DNA complex as required to prepare the system for the chemical reaction of nucleotide extension.

Exploring the electronic and mechanical properties of protein using conducting atomic force microscopy

In interfacing man-made electronic components with specifically folded biomacromolecules, the perturbative effects of junction structure on any signal generated should be considered. We report herein on the electron-transfer characteristics of the blue copper metalloprotein, azurin, as characterized at a refined level by conducting atomic force microscopy (C-AFM). Specifically, the modulation of current-voltage (I-V) behavior with compressional force has been examined. In the absence of assignable resonant electron tunneling within the confined bias region, from -1 to 1 V, the I-V behavior was analyzed with a modified Simmons formula. To interpret the variation of tunneling barrier height and barrier length obtained by fitting with the modified Simmons formula, an atom packing density model associated with protein mechanical deformation was proposed and simulated by molecular dynamics. The barrier heights determined at the minimum forces necessary for stable electrical contact correlate reasonably well with those estimated from bulk biophysical (electroanalytical and photochemical) experiments previously reported. At higher forces, the tunnel barrier decreases to fall within the range observed with saturated organic systems. Molecular dynamics simulations revealed changes in secondary structure and atomic density of the protein with respect to compression. At low compression, where transport measurements are made, secondary structure is retained, and atomic packing density is observed to increase linearly with force. These predictions, and those made at higher compression, are consistent with both experimentally observed modulations of tunneling barrier height with applied force and the applicability of the atom packing density model of electron tunneling in proteins to molecular-level analyses.

Excited state charge-transfer dynamics study of poplar plastocyanin by ultrafast pump-probe spectroscopy and molecular dynamics simulation

We have applied ultrafast pump-probe spectroscopy to investigate the excited state dynamics of the blue copper protein poplar plastocyanin, by exciting in the blue side of its 600-nm absorption band. The decay of the charge-transfer excited state occurs exponentially with a time constant of approximately 280 fs and is modulated by well visible oscillations. The Fourier transform of the oscillatory component, besides providing most of the vibrational modes found by conventional resonance Raman, presents additional bands in the low frequency region modes, which are reminiscent of collective motions of biological relevance. Notably, a high frequency mode at approximately 508 cm(-1), whose dynamics are consistent with that of the excited state and already observed for other blue copper proteins, is shown to be present also in poplar plastocyanin. This vibrational mode is reproduced by a molecular dynamics simulation involving the excited state of the copper site.

MD simulation of a plastocyanin mutant adsorbed onto a gold surface

MD simulation of plastocyanin, an electron transfer protein, adsorbed onto a gold surface, has been performed for 10 ns. Starting from the crystallographic structure of a poplar plastocyanin mutant engineered with the insertion of a disulfide bridge, the protein has been anchored to a gold substrate modeled by a cluster of three layers in the Au<111> configuration. A number of significant structural and dynamical properties of the protein molecule, covalently bound through either one or two sulfur atoms to the gold surface, has been extracted and compared with those of the free protein. Attention has been paid to investigate the dynamical aspects putatively related to the electron transfer process. In particular, the cross-correlation function between specific active site vibrations and all the other protein atom motions and the principal component analysis have been calculated in order to put into evidence dynamical correlation of some functional relevance. The results are discussed also in connection with related experiments.

Mutagenesis of prochlorothrix plastocyanin reveals additional features in photosystem I interactions

Three surface residues of plastocyanin from Prochlorothrix hollandica have been modified by site-directed mutagenesis. Changes have been made in methionine 33, located in the hydrophobic patch of the copper protein, and in arginine 86 and proline 53, both located in the eastern hydrophilic area. The reactivity toward photosystem I of single mutants M33N, P53A, P53E, R86Q, R86E, and the double mutant M33N/P14L has been studied by laser flash absorption spectroscopy. All the mutations yield increased reactivity of plastocyanin toward photosystem I as compared with wild type plastocyanin, thus indicating that in Prochlorothrix electron donation to photosystem I is not optimized. The most drastic increases in the intracomplex electron transfer rate are obtained with mutants in methionine 33, whereas replacing arginine 86 only modestly affects the plastocyanin-photosystem I equilibrium constant for complex formation. Mutations at position 53 also promote major changes in the association of plastocyanin with photosystem I, yielding a change from a mechanism involving complex formation to a simpler collisional interaction. Molecular dynamics calculations indicate that mutations at position 33 promote changes in the H-bond network around the copper center. The comparative kinetic analysis of the reactivity of Prochlorothrix plastocyanin mutants toward photosystem I from other cyanobacteria reveals that mutations M33N, P53A, and P53E result in enhanced general reactivity.