New insights into flavivirus biology: the influence of pH over interactions between prM and E proteins.

Diseases caused by flaviviruses, such as dengue and zika, are globally recognized as major threats. During infection, a critical point in their replicative cycle is the maturation step, which occurs throughout the cellular exocytic pathway. This step is a pH-dependent process that involves the modification of the viral envelope by converting prM (pre-membrane) into M (membrane) proteins with the release of a "pr peptide". After this reaction, the pr peptides remain bound to the viral envelope while the virions cross the acidic trans-Golgi network, and are released only at neutral pH after secretion of the virus particles. Despite this current knowledge, the molecular basis of the flavivirus maturation step is largely unknown. Here, based on the crystal structure of the dengue pr-E complex ("pr peptide" bound to virus envelope protein) and using molecular dynamics simulations, we found that the pH shift from acidic to neutral yields considerable structural changes in the system. Dynamic cross correlation maps and root mean square deviation analyses revealed that the pr-E junction is clearly unstable under neutral pH. Secondary structure analysis also revealed that the fusion loop region, present in the E protein, is sensitive to pH and tends to unstructure at a neutral environment. Moreover, we found that five residues present in the E protein, Gly102, His244, Thr70, Thr68 and Asn67 are critical to confer stability to the pr-E complex while inside the Golgi apparatus. This work brings details about the dynamical behavior of the pr-E system, helps to better understand the flavivirus biology and may also be of use in the development of novel antiviral strategies.

The flavivirus capsid protein: Structure, function and perspectives towards drug design.

Flaviviruses, such as dengue and zika viruses, are etiologic agents transmitted to humans mainly by arthropods and are of great epidemiological interest. The flavivirus capsid protein is a structural element required for the viral nucleocapsid assembly that presents the classical function of sheltering the viral genome. After decades of research, many reports have shown its different functionalities and influence over cell normal functioning. The subcellular distribution of this protein, which involves accumulation around lipid droplets and nuclear localization, also corroborates with its multi-functional characteristic. As flavivirus diseases are still in need of global control and in view of the possible key functionalities that the capsid protein promotes over flavivirus biology, novel considerations arise towards anti-flavivirus drug research. This review covers the main aspects concerning structural and functional features of the flavivirus C protein, ultimately, highlighting prospects in drug discovery based on this viral target.

The mechanism by which P250L mutation impairs flavivirus-NS1 dimerization: an investigation based on molecular dynamics simulations.

The flavivirus non-structural protein 1 (NS1) is a conserved glycoprotein with as yet undefined biological function. This protein dimerizes when inside infected cells or associated to cell membranes but also forms lipid-associated hexamers when secreted to the extracellular space. A single amino acid substitution (P250L) is capable of preventing the dimerization of NS1 resulting in lower virulence and slower virus replication. In this work, based on molecular dynamics simulations of the dengue-2 virus NS1 [Formula: see text]-ladder monomer as a core model, we found that this mutation can induce several conformational changes that importantly affect critical monomer-monomer interactions. Based on additional simulations, we suggest a mechanism by which a highly orchestrated sequence of events propagate the local perturbations around the mutation site towards the dimer interface. The elucidation of such a mechanism could potentially support new strategies for rational production of live-attenuated vaccines and highlights a step forward in the development of novel anti-flavivirus measures.

Combining Nuclear Magnetic Resonance Spectroscopy and Density Functional Theory Calculations to Characterize Carvedilol Polymorphs.

The experiments of carvedilol form II, form III, and hydrate by (13)C and (15)N cross-polarization magic-angle spinning (CP MAS) are reported. The GIPAW (gauge-including projector-augmented wave) method from DFT (density functional theory) calculations was used to simulate (13)C and (15)N chemical shifts. A very good agreement was found for the comparison between the global results of experimental and calculated nuclear magnetic resonance (NMR) chemical shifts for carvedilol polymorphs. This work aims a comprehensive understanding of carvedilol crystalline forms employing solution and solid-state NMR as well as DFT calculations.

Association of the anti-tuberculosis drug rifampicin with a PAMAM dendrimer.

The association of the anti-tuberculosis drug rifampicin (RIF) with a 4th-generation poly(amidoamine) (G4-PAMAM) dendrimer was investigated by means of molecular dynamics simulations. The RIF load capacity was estimated to be around 20 RIF per G4-PAMAM at neutral pH. The complex formed by 20 RIF molecules and the dendrimer (RIF20-PAMAM) was subjected to 100 ns molecular dynamics (MD) simulations at two different pH conditions (neutral and acidic). The complex was found to be significantly more stable in the simulation at neutral pH compared to the simulation at low pH in which the RIF molecules were rapidly and almost simultaneously expelled to the solvent bulk. The high stability of the RIF-PAMAM complex under physiological pH and the rapid release of RIF molecules under acidic medium provide an interesting switch for drug targeting since the Mycobacterium resides within acidic domains of the macrophage. Altogether, these results suggest that, at least in terms of stability and pH-dependent release, PAMAM-like dendrimers may be considered suitable drug delivery systems for RIF and derivatives.

Design of inhibitors of thymidylate kinase from Variola virus as new selective drugs against smallpox.

Recently we constructed a homology model of the enzyme thymidylate kinase from Variola virus (VarTMPK) and proposed it as a new target to the drug design against smallpox. In the present work, we used the antivirals cidofovir and acyclovir as reference compounds to choose eleven compounds as leads to the drug design of inhibitors for VarTMPK. Docking and molecular dynamics (MD) studies of the interactions of these compounds inside VarTMPK and human TMPK (HssTMPK) suggest that they compete for the binding region of the substrate and were used to propose the structures of ten new inhibitors for VarTMPK. Further docking and MD simulations of these compounds, inside VarTMPK and HssTMPK, suggest that nine among ten are potential selective inhibitors of VarTMPK.

Mixed Monte Carlo/Molecular Dynamics simulations of the prion protein.

In this paper we present the results of mixed Monte Carlo/Molecular Dynamics (MC/MD) simulations of the D178N mutant of the human prion protein. We have used the MC moves for polypeptide sampling known as Concerted Rotations with Angles (CRA) to selectively sample the region of the prion protein comprising the ß-sheet and one of the α-helices. The results indicate that the MC/MD simulations sample the phase space substantially faster than regular Molecular Dynamics simulations starting with the same initial conditions. This work further indicates the MC/MD technique as a potentially powerful simulation tool, allowing the selective sampling of a region of a physical system that is deemed important.

CoMFA/CoMSIA 3D-QSAR of pyrimidine inhibitors of Pneumocystis carinii dihydrofolate reductase.

Pneumocystis carinii is typically a non-pathogenic fungus found in the respiratory tract of healthy humans. However, it may cause P. carinii pneumonia (PCP) in people with immune deficiency, affecting mainly premature babies, cancer patients and transplant recipients, and people with acquired immunodeficiency syndrome (AIDS). In the latter group, PCP occurs in approximately 80% of patients, a major cause of death. Currently, there are many available therapies to treat PCP patients, including P. carinii dihydrofolate reductase (PcDHFR) inhibitors, such as trimetrexate (TMX), piritrexim (PTX), trimethoprim (TMP), and pyrimethamine (PMT). Nevertheless, the high percentage of adverse side effects and the limited therapeutic success of the current drug therapy justify the search for new drugs rationally planned against PCP. This work focuses on the study of pyrimidine inhibitors of PcDHFR, using both CoMFA and CoMSIA 3D-QSAR methods.

Mixed Monte Carlo/molecular dynamics simulations in explicit solvent.

A mixed Monte Carlo/Molecular Dynamics method using the trial moves for peptide backbone sampling known as Concerted Rotations with Angles was implemented. The algorithm was used to study polyalanine systems. Equivalent results to conventional Molecular Dynamics were obtained for simulations of Ala(6) in implicit solvent. To test the efficiency of the implemented method, several 150 ns simulations of Ala(12) in explicit water were performed. The results show that the present method yields significantly faster formation of secondary structure than the conventional Molecular Dynamics simulations. This opens the possibility to selectively sample alanine-rich regions of larger peptides or proteins. It remains to be established whether hydrophilic amino acid residues can be successfully treated with the present methodology.

Engineered monomeric human histidine triad nucleotide-binding protein 1 hydrolyzes fluorogenic acyl-adenylate and lysyl-tRNA synthetase-generated lysyl-adenylate.

Hint1 is a homodimeric protein and member of the ubiquitous HIT superfamily. Hint1 catalyzes the hydrolysis of purine phosphoramidates and lysyl-adenylate generated by lysyl-tRNA synthetase (LysRS). To determine the importance of homodimerization on the biological and catalytic activity of Hint1, the dimer interface of human Hint1 (hHint1) was destabilized by replacement of Val(97) of hHint1 with Asp, Glu, or Arg. The mutants were shown to exist as monomers in solution by a combination of size exclusion chromatograph, static light scattering, and chemically induced dimerization studies. Circular dichroism studies revealed little difference between the stability of the V97D, V97E, and wild-type hHint1. Relative to wild-type and the V97E mutant, however, significant perturbation of the V97D mutant structure was observed. hHint1 was shown to prefer 3-indolepropionic acyl-adenylate (AIPA) over tryptamine adenosine phosphoramidate monoester (TpAd). Wild-type hHint1 was found to be 277- and 1000-fold more efficient (k(cat)/K(m) values) than the V97E and V97D mutants, respectively. Adenylation of wild-type, V97D, and V97E hHint1 by human LysRS was shown to correlate with the mutant k(cat)/K(m) values using 3-indolepropionic acyl-adenylate as a substrate, but not tryptamine adenosine phosphoramidate monoester. Significant perturbations of the active site residues were not detected by molecular dynamics simulations of the hHint1s. Taken together, these results demonstrate that for hHint1; 1) the efficiency (k(cat)/K(m)) of acylated AMP hydrolysis, but not maximal catalytic turnover (k(cat)), is dependent on homodimerization and 2) the hydrolysis of lysyl-AMP generated by LysRS is not dependent on homodimerization if the monomer structure is similar to the wild-type structure.

Dynamical behavior of the vascular endothelial growth factor: biological implications.

The vascular endothelial growth factor (VEGF) seems to be the most important regulator of physiological and pathological angiogenesis, being, for this reason, a favorite target for therapies against angiogenesis-related diseases. VEGF is a homodimer in which the monomers are formed by beta-strands interconnected on the poles by three loops. A recent work showed that an intimate relationship between loops-1 and -3 is required for high affinity binding to the receptors (Kiba et al., J Biol Chem 2003;278:13453-13461). In this work, we report the results of a 10-ns molecular dynamics simulation of VEGF. We analyzed the dynamical behavior of the protein (using a dynamical cross-correlation map) and found that it is governed by a high degree of correlation between the motions of the loops. We also performed a principal component analysis and found an overall motion in which the opposite poles are projected against each other, just like the movement of the wings of a butterfly. From the biological point of view, it is likely that this motion would facilitate receptor binding since VEGF must enter a restricted cavity formed by the two subunits of the receptor.

The Na+ binding channel of human coagulation proteases: novel insights on the structure and allosteric modulation revealed by molecular surface analysis.

Thrombovascular diseases result from imbalanced haemostasis and comprise important health problems in the aging population worldwide. The activity of enzymes pertaining to the coagulation cascade of mammalians exhibit several control mechanisms in order to maintain a proper balance between bleeding and thrombosis. For instance, human coagulation serine proteases carrying a F225 or Y225 are allosteric modulated by the binding of Na+ in a water-filled channel connected to the primary specificity pocket (S1 subsite) of these enzymes. We have characterized the structure, topography and lipophilicity of this channel in the ligand-free fast (sodium-bound) and slow (sodium-free) forms of thrombin, in the sole available structure of activated protein C and in several structures of the coagulation factors VIIa, IXa and Xa, differing in the nature of the bound inhibitor and in the occupancy of exosite-I as well as the Ca2+ and Na+ binding sites. Opposite to thrombin, the aqueous channels in all other coagulation enzymes sheltering a Na+ binding site do not have an aperture on the enzyme surface opposite to the S1 subsite entrance. In these enzymes, the lack of the three-residue insertion in loop 1 (183-189) as found in thrombin allied to compensatory mutations in the positions 187-185 and 222 effects a constriction in the water-filled channel that ends up by segregating the ion binding site from the S1 subsite. We also disclosed major topographical changes on the thrombin's surface upon sodium release and transition to the slow form that culminate in the narrowing of the S1 subsite entrance and, strikingly, in the loss of communication between the primary specificity pocket and the exosite-I. Such observation is in accordance with existing experimental data demonstrating thermodynamic linkage between these distant regions on the thrombin surface. Conformational changes in F34, L40, R73 and T74 were the main responsible for this effect. A path by which these changes in the vicinity of exosite-I could be transmitted to the S1 subsite and, consequently, to the sodium binding site is proposed.