Search

Adenylyl cyclases (ACs) play a key role in many signaling cascades. ACs catalyze the production of cyclic AMP from ATP and this function is stimulated or inhibited by the binding of their cognate stimulatory or inhibitory Gα subunits, respectively. Here we used simulation tools to uncover the molecular and subcellular mechanisms of AC function, with a focus on the AC5 isoform, extensively studied experimentally. First, quantum mechanical/molecular mechanical free energy simulations were used to investigate the enzymatic reaction and its changes upon point mutations. Next, molecular dynamics simulations were employed to assess the catalytic state in the presence or absence of Gα subunits. This led to the identification of an inactive state of the enzyme that is present whenever an inhibitory Gα is associated, independent of the presence of a stimulatory Gα. In addition, the use of coevolution‐guided multiscale simulations revealed that the binding of Gα subunits reshapes the free‐energy landscape of the AC5 enzyme by following the classical population‐shift paradigm. Finally, Brownian dynamics simulations provided forward rate constants for the binding of Gα subunits to AC5, consistent with the ability of the protein to perform coincidence detection effectively. Our calculations also pointed to strong similarities between AC5 and other AC isoforms, including AC1 and AC6. Findings from the molecular simulations were used along with experimental data as constraints for systems biology modeling of a specific AC5‐triggered neuronal cascade to investigate how the dynamics of downstream signaling depend on initial receptor activation. This article is categorized under:Structure and Mechanism > Computational Biochemistry and BiophysicsMolecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo MethodsSoftware > Molecular Modeling [ABSTRACT FROM AUTHOR]

Normal mode analysis (NMA) is a fast and inexpensive approach that is largely used to gain insight into functional protein motions, and more recently to create conformations for further computational studies. However, when the protein structure is unknown, the use of computational models is necessary. Here, we analyze the capacity of NMA in internal coordinate space to predict protein motion, its intrinsic flexibility, and atomic displacements, using protein models instead of native structures, and the possibility to use it for model refinement. Our results show that NMA is quite insensitive to modeling errors, but that calculations are strictly reliable only for very accurate models. Our study also suggests that internal NMA is a more suitable tool for the improvement of structural models, and for integrating them with experimental data or in other computational techniques, such as protein docking or more refined molecular dynamics simulations. [ABSTRACT FROM AUTHOR]

Evidence of a special chiral nematic phase is provided using numerical simulation and Onsager theory for systems of hard helical particles. This phase appears at the high density end of the nematic phase, when helices are well aligned, and is characterized by the C2 symmetry axes of the helices spiraling around the nematic director with periodicity equal to the particle pitch. This coupling between translational and rotational degrees of freedom allows a more efficient packing and hence an increase of translational entropy. Suitable order parameters and correlation functions are introduced to identify this screw-like phase, whose main features are then studied as a function of radius and pitch of the helical particles. Our study highlights the physical mechanism underlying a similar ordering observed in colloidal helical flagella [E. Barry, Z. Hensel, Z. Dogic, M. Shribak, and R. Oldenbourg, Phys. Rev. Lett. 96, 018305 (2006)] and raises the question of whether it could be observed in other helical particle systems, such as DNA, at sufficiently high densities. [ABSTRACT FROM AUTHOR]

Using an Onsager-like theory, we have investigated the relationship between the morphology of hard helical particles and the features (pitch and handedness) of the cholesteric phase that they form. We show that right-handed helices can assemble into right- (UGRAPHIC DISPLAY="INLINE" ID="UGT1" SRC="UGT1"/) and left-handed (UGRAPHIC DISPLAY="INLINE" ID="UGT2" SRC="UGT2"/) cholesterics, depending on their curliness, and that the cholesteric pitch is a non-monotonic function of the intrinsic pitch of particles. The theory leads to the definition of a hierarchy of pseudoscalars, which quantify the difference in the average excluded volume between pair configurations of helices having (UGRAPHIC DISPLAY="INLINE" ID="UGT3" SRC="UGT3"/) and (UGRAPHIC DISPLAY="INLINE" ID="UGT4" SRC="UGT4"/)-skewed axes. The predictions of the Onsager-like theory are supported by Monte Carlo simulations of the isotropic phase of hard helices, showing how the cholesteric organization, which develops on scales longer than hundreds of molecular sizes, is encoded in the short-range chiral correlations between the helical axes. [ABSTRACT FROM AUTHOR]

As more sequencing data accumulate and novel puzzling genetic regulations are discovered, the need for accurate automated modeling of RNA structure increases. RNA structure modeling from chemical probing experiments has made tremendous progress, however accurately predicting large RNA structures is still challenging for several reasons: RNA are inherently flexible and often adopt many energetically similar structures, which are not reliably distinguished by the available, incomplete thermodynamic model. Moreover, computationally, the problem is aggravated by the relevance of pseudoknots and non-canonical base pairs, which are hardly predicted efficiently. To identify nucleotides involved in pseudoknots and non-canonical interactions, we scrutinized the SHAPE reactivity of each nucleotide of the 188 nt long lariat-capping ribozyme under multiple conditions. Reactivities analyzed in the light of the X-ray structure were shown to report accurately the nucleotide status. Those that seemed paradoxical were rationalized by the nucleotide behavior along molecular dynamic simulations. We show that valuable information on intricate interactions can be deduced from probing with different reagents, and in the presence or absence of Mg2+. Furthermore, probing at increasing temperature was remarkably efficient at pointing to non-canonical interactions and pseudoknot pairings. The possibilities of following such strategies to inform structure modeling software are discussed. [ABSTRACT FROM AUTHOR]

Endoglin (Eng) is an endothelial cell (EC) transmembrane glycoprotein involved in adhesion and angiogenesis. Eng mutations result in vessel abnormalities as observed in hereditary hemorrhagic telangiectasia of type 1. The role of Eng was investigated in endothelial functions and permeability under inflammatory conditions, focusing on the actin dynamic signaling pathway. Endothelial Colony-Forming Cells (ECFC) from human cord blood and mouse lung/aortic EC (MLEC, MAEC) from Eng+/+ and Eng+/− mice were used. ECFC silenced for Eng with Eng-siRNA and ctr-siRNA were used to test tubulogenesis and permeability +/− TNFα and +/− LIM kinase inhibitors (LIMKi). In silico modeling of TNFα–Eng interactions was carried out from PDB IDs 5HZW and 5HZV. Calcium ions (Ca2+) flux was studied by Oregon Green 488 in epifluorescence microscopy. Levels of cofilin phosphorylation and tubulin post-translational modifications were evaluated by Western blot. F-actin and actin–tubulin distribution/co-localization were evaluated in cells by confocal microscopy. Eng silencing in ECFCs resulted in a decrease of cell sprouting by 50 ± 15% (p < 0.05) and an increase in pseudo-tube width (41 ± 4.5%; p < 0.001) compared to control. Upon TNFα stimulation, ECFC Eng–siRNA displayed a significant higher permeability compared to ctr-siRNA (p < 0.01), which is associated to a higher Ca2+ mobilization (p < 0.01). Computational analysis suggested that Eng mitigated TNFα activity. F-actin polymerization was significantly increased in ECFC Eng-siRNA, MAEC+/−, and MLEC+/− compared to controls (p < 0.001, p < 0.01, and p < 0.01, respectively) as well as actin/tubulin distribution (p < 0.01). Furthermore, the inactive form of cofilin (P-cofilin at Ser3) was significantly decreased by 36.7 ± 4.8% in ECFC Eng-siRNA compared to ctr-siRNA (p < 0.001). Interestingly, LIMKi reproduced the absence of Eng on TNFα-induced ECFC-increased permeability. Our data suggest that Eng plays a critical role in the homeostasis regulation of endothelial cells under inflammatory conditions (TNFα), and loss of Eng influences ECFC-related permeability through the LIMK/cofilin/actin rearrangement-signaling pathway. [ABSTRACT FROM AUTHOR]

We investigate the isotropic-to-nematic phase transition in systems of hard helical particles, using Onsager theory and Monte Carlo computer simulations. Motivation of this work resides in the ubiquity of the helical shape motif in many natural and synthetic polymers, as well as in the well known importance that the details of size and shape have in determining the phase behaviour and properties of (soft) condensed matter systems. We discuss the differences with the corresponding spherocylinder phase diagram and find that the helix parameters affect the phase behaviour and the existence of the nematic phase. We find that for high helicity Onsager theory significantly departs from numerical simulations even when a modified form of the Parsons-Lee rescaling is included to account for the non-convexity of particles. [ABSTRACT FROM AUTHOR]

Adenylyl cyclases (ACs) have a crucial role in many signal transduction pathways, in particular in the intricate control of cyclic AMP (cAMP) generation from adenosine triphosphate (ATP). Using homology models developed from existing structural data and docking experiments, we have carried out all-atom, microsecond-scale molecular dynamics simulations on the AC5 isoform of adenylyl cyclase bound to the inhibitory G-protein subunit Gαi in the presence and in the absence of ATP. The results show that Gαi has significant effects on the structure and flexibility of adenylyl cyclase, as observed earlier for the binding of ATP and Gsα. New data on Gαi bound to the C1 domain of AC5 help explain how Gαi inhibits enzyme activity and obtain insight on its regulation. Simulations also suggest a crucial role of ATP in the regulation of the stimulation and inhibition of AC5. [ABSTRACT FROM AUTHOR]

Protein-protein interactions are at the basis of many protein functions, and the knowledge of 3D structures of protein-protein complexes provides structural, mechanical and dynamical pieces of information essential to understand these functions. Protein-protein interfaces can be seen as stable, organized regions where residues from different partners form non-covalent interactions that are responsible for interaction specificity and strength. They are commonly described as a peripheral region, whose role is to protect the core region that concentrates the most contributing interactions, from the solvent. To get insights into the dynamics of protein-protein complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on eight different protein-protein complexes of different functional class and interface size by taking into account the bound and unbound forms. On the one hand, we characterized structural changes upon binding of the proteins, and on the other hand we extensively analyzed the interfaces and the structural waters involved in the binding. Based on our analysis, in 6 cases out of 8, the interfaces rearranged during the simulation time, in stable and long-lived substates with alternative residue-residue contacts. These rearrangements are not restricted to side-chain fluctuations in the periphery but also affect the core interface. Finally, the analysis of the waters at the interface and involved in the binding pointed out the importance to take into account their role in the estimation of the interaction strength., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Martin and Frezza.)

RNA molecules can easily adopt alternative structures in response to different environmental conditions. As a result, a molecule's energy landscape is rough and can exhibit a multitude of deep basins. In the absence of a high-resolution structure, small-angle X-ray scattering data (SAXS) can narrow down the conformational space available to the molecule and be used in conjunction with physical modeling to obtain high-resolution putative structures to be further tested by experiments. Because of the low resolution of these data, it is natural to implement the integration of SAXS data into simulations using a coarse-grained representation of the molecule, allowing for much wider searches and faster evaluation of SAXS theoretical intensity curves than with atomistic models. We present here the theoretical framework and the implementation of a simulation approach based on our coarse-grained model HiRE-RNA combined with SAXS evaluations "on-the-fly" leading the simulation toward conformations agreeing with the scattering data, starting from partially folded structures as the ones that can easily be obtained from secondary structure prediction-based tools. We show on three benchmark systems how our approach can successfully achieve high-resolution structures with remarkable similarity with the native structure recovering not only the overall shape, as imposed by SAXS data, but also the details of initially missing base pairs.

Unveiling the subtle rules that control the buildup of macroscopic chirality starting from chiral molecular elements is a challenge for theory and computations. In this context, a remarkable phenomenon is the formation of helically twisted nematic (cholesteric) phases, with pitch in the micrometer range, driven by self-assembly of relatively small chiral species into supramolecular semiflexible polymers. We have developed a theoretical framework to connect the cholesteric organization to the shape and chirality of the constituents, described with molecular detail, in this kind of system. The theory has been tested against new accurate measurements for solutions of short DNA duplexes. We show that the cholesteric organization is determined by steric repulsion between duplexes, and we identify distinctive features of linear self-assembly in the temperature and concentration dependence of the pitch.

We analyze the capacity of normal modes to predict observed protein conformational changes, and, notably, those induced by the formation of protein-protein complexes. We show that normal modes calculated in internal coordinate space (ICS) provide better predictions. For a large test set, using the ICS approach describes the conformational changes more completely, and with fewer low-frequency modes than the equivalent Cartesian coordinate modes, despite the fact that the internal coordinate calculations were restricted to torsional angles. This can be attributed to the fact that the use of ICS extends the range over which movements along the corresponding eigenvectors remain close to the true conformational energy hypersurface. We also show that the PaLaCe coarse-grain protein model performs better than a simple elastic network model. We apply ICS normal-mode analysis to protein complexes and, by extending the approach of Sunada and Go̅, [Sunada, S.; Go̅, N. J. Comput. Chem. 1995, 16, 328-336], we show that we can couple an accurate view of the Cartesian coordinate movements induced by ICS modes with the detection of the key residues responsible for the movements.

The chirality transfer in liquid crystals induced by two helical perylenequinones (namely, the natural compounds cercosporin and phleichrome) was investigated by integrating measurements of helical twisting power with a conformational analysis by DFT calculations and with the prediction of their twisting ability by the surface-chirality method. The two quasi-enantiomeric derivatives induce oppositely handed cholesteric phases when introduced as dopants in nematic solvents. We evaluated the role of the different conformations of the chiral hydroxyalkyl side chains in determining the helical twisting power: They were found to affect the strength of the chirality transfer, although the handedness of the induced cholesteric phase is essentially determined by the axial chirality (helicity) of the core of the perylenequinones.